data(shoulder)
if (FALSE) {
fit <- nparLD(resp ~ group1 * group2 * time,
data = shoulder,
subject = "subject",
hypothesis = "H0F")
fit
}
if (FALSE) {
fit <- nparLD(resp ~ group1 * group2 * time,
data = shoulder,
subject = "subject",
hypothesis = "H0p",
contrast = list("group1:time"),
Factor.Information=TRUE)
fit
plot(fit)
plot(fit$MCTP)
}
Run the code above in your browser using DataLab