# NOT RUN {
data(epilepsy)
fit <- gcmr(counts ~ offset(log(time)) + visit + trt + visit:trt, data = epilepsy,
subset = (id != 49), marginal = negbin.marg, cormat = cluster.cormat(id, "ar1"),
options=gcmr.options(seed=123, nrep=c(25,100) ))
summary(fit)
# }
Run the code above in your browser using DataCamp Workspace