# NOT RUN {
# data (for demonstration only)
library(MASS)
data(birthwt)
# impose an artificial (not fully appropriate) grid structure
x<-rep(1:14,14)
y<-as.integer(gl(14,14))
coords<-cbind(x[-(190:196)],y[-(190:196)])
# }
# NOT RUN {
formula<-formula(low ~ race + smoke + bwt)
mgee<-GEE(formula, family = "gaussian", data = birthwt,
coord=coords, corstr="fixed",scale.fix=TRUE)
mmi<-mmiGEE(mgee,birthwt)
# }
Run the code above in your browser using DataLab