data(impsdata)
id<-impsdata$id
r<-impsdata$r
data.real <- data.frame(id=id,y=impsdata$y, impsdata$x)
head(data.real,n=10)
# each participant has three records
time<-4
# the outcome is binary
family=binomial()
models <- list(y~Time, y~Drug,y~Time+Drug,
y~Time*Drug,y~Time+Sex+Drug,
y~Time+Sex+Drug+Time:Sex+Sex:Drug+Drug:Time)
model_mis<-r~Drug+Time+Sex
candidate.cor.sets<-c("exchangeable","independence","ar1")
#not run
#output_ELCIC<-ELCICwgee(models, candidate.cor.sets,data=data.real,model_mis,family,r,id,time)
#output_MLIC<-MLICwgee(models, candidate.cor.sets,data=data.real,model_mis,family,r,id,time)
#output_QICW<-QICWwgee(models, candidate.cor.sets,data=data.real,model_mis,family,r,id,time)
Run the code above in your browser using DataLab