# NOT RUN {
##### data = seizure
seiz1 <- cold(y ~ lage + lbase + v4 + trt + trt:lbase, data = seizure,
start = NULL, dependence = "AR1")
seiz2 <-cold(y ~ lage + lbase + v4 + trt, data = seizure, start = NULL,
dependence = "AR1")
anova(seiz1, seiz2)
##### data = datacold
mod0 <- cold(z ~ Time * Treatment, data = datacold, time = "Time",
id = "Subject", dependence = "ind")
mod0R <- cold(z ~ Time * Treatment, random = ~ 1, data = datacold,
time = "Time", id = "Subject", dependence = "indR")
summary(mod0R)
anova(mod0, mod0R)
# }
Run the code above in your browser using DataLab