data(plasma)
if (FALSE) {
fit <- nparLD(resp ~ group * time,
data = plasma,
subject = "subject",
hypothesis = "H0F", effect="weighted")
fit
}
if (FALSE) {
fit <- nparLD(resp ~ group * time,
data = plasma,
subject = "subject",
hypothesis = "H0p",
contrast = list("group", "Tukey"))
fit
plot(fit)
plot(fit$MCTP)
}
Run the code above in your browser using DataLab