Hi R-Help.
I am working on a data set with a 3-level nested structure. I have
individuals nested in households and multiple observations on each
individual. I assume that the individuals inside a given household are
correlated and that the individuals are correlated with themselves
over time. The data is not balanced.
I have computed a GLMM with logit link function and two random normal
distributed intercepts attributed to households and individuals. Now I
want to compare it with a population average (marginal) model using
GEE.
I have looked at the geepack library but I'm unsure if it is possible
to use it in my case where I have the additional household clustering.
I am willing to assume an exchangable working correlation stucture
inside households, similarly for invidividuals over waves and
otherwise independence.
Can such a thing be done with the packages available?
Thanks in advance,
Andreas Jensen
