# NOT RUN {
# First generate a dataset with 200 patients, rho =0.25 and tau = 0.5 and
# then estimate the parameters using estimcov.
set.seed(2015)
dataset <- r.gee.1subgroup(n=200, reg=list(c(0,0,0,0.1),c(0,0,0,0.1)), sigma=c(3,2.5),
tau=0.5, rho=0.25, theta=1, k=1.5, Time=c(0:5), OD=0)
estimations <- estimcov(data=dataset,Time=c(0:5))
estimations[[1]]
estimations[[2]]
# }
Run the code above in your browser using DataLab