如何在r中为多类自变量和响应变量拟合GEE模型？超过3类

我正试图在

R

中拟合GEE模型，即：

fit1 <-geeglm (crime~month+date+time,data=data,family =gaussian,=subject,corstr="ind")
Fit2 <-update(fit1,corstr ="exch")
Fit3<-update(fit1,corstr="ar1")
QIC(fit1,Fit2,Fit3)

---

fit1对不起，我是否应该用in在r中做了什么来回答这个问题？最好有一个可复制的示例，因为其他人可以测试它并检查潜在的问题。