Entiendo que una de las ventajas de GEE es que usted no haga suposiciones acerca de la distribución de probabilidad conjunta, confiando en cambio en la media, la varianza y las asociaciones (corr). Entonces, ¿por qué hemos de especificar una familia en el comando gee
?
model <- geeglm(y~x, data=data, id=id, family=gaussian('identity'), corstr='ar1')