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