# NOT RUN {
data(datacoldM)
str(datacoldM)
mod0.M<- cold(z~Time*Treatment, data=datacoldM, time="Time",
id="Subject", dependence="ind")
summary (mod0.M)
mod1.M<- cold(z~Time*Treatment, data=datacoldM, time="Time",
id="Subject", dependence="AR1")
summary (mod1.M)
anova(mod0.M,mod1.M)
# }
Run the code above in your browser using DataLab