# NOT RUN {
data(imps)
fit <- wgee(Y ~ Drug+Sex+Time, data=imps, id=imps$ID, family="binomial",
corstr="exchangeable",scale=NULL,mismodel= R ~ Drug+Time)
WRsquare.gee(fit, weight_mean=TRUE)
# }
Run the code above in your browser using DataLab