muscatine$cage <- muscatine$age - 12
muscatine$cage2 <- muscatine$cage^2
f1 <- numobese ~ gender
f2 <- numobese ~ gender + cage + cage2 +
gender:cage + gender:cage2
gee1 <- geeglm(formula = f1, id = id,
waves = occasion, data = muscatine, family = binomial(),
corstr = "independence")
gee2 <- geeglm(formula = f2, id = id,
waves = occasion, data = muscatine, family = binomial(),
corstr = "independence")
tidy(gee1)
tidy(gee2)
QIC(gee1)
QIC(gee2)
Run the code above in your browser using DataLab