# NOT RUN {
##### data = seizure
### AR1
seiz1M <- cold(y ~ lage + lbase + v4 + trt + trt:lbase, data = seizure,
start = NULL, dependence = "AR1")
getAIC(seiz1M)
### indR
seiz0R <- cold(y ~ lage + lbase + trt + trt:lbase + v4, random = ~ 1,
data = seizure, dependence = "indR")
getAIC(seiz0R)
# }
Run the code above in your browser using DataLab