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