Re: [R] multiple nesting levels in GEE

About this list Date view Thread view Subject view Author view Attachment view

From: Thomas Lumley (tlumley@u.washington.edu)
Date: Sun 30 May 2004 - 02:12:03 EST


Message-id: <Pine.A41.4.58.0405290905390.148958@homer07.u.washington.edu>

On Sat, 29 May 2004 Aurelie.Cohas@univ-lyon1.fr wrote:

> Hello,
>
> I'm actually trying to fit a gee model with 2 nesting levels since I expect a
> correlation between all members of a litter at a first level and between all
> individuals sharing a mother at a second superior level with an exchangeable
> matrix. I order my dataframe by both mother and litter

The gee() function does not compute this correlation structure (because
no-one has implemented it). I don't think other software implements it
either, although people keep pointing out that it would be useful.

  This need not matter -- the point of GEE is that the correlation
structure is not important for validity of inference (although it may
affect efficiency). If you are fitting a linear model I would recommend
lme(), otherwise you should just be able to use an exchangeable model at
the largest clustering size.

        -thomas

> I try several syntaxes:
> id= mother*litter which give the same correlation matrix as id= litter*mother
> and id=litter with a matrix size of 7*7 which correspond to the maximum number of
> young per litter
>
> id=litter|mother which give the same correlation matrix as id=mother|litter with
> a correlation matrix size of 235*235 which correspond to the number of
> observations
>
> id=litter/mother and id=mother/litter give different result with a strange matrix
> of 11*11 with no 1 on the diagonal
>
> If anybody know if it is possible to construct nesting levels in R and what is
> the good syntax for this kind of model?
>
> Thanks in advance
>
>
> Aurélie Cohas
>
> ______________________________________________
> R-help@stat.math.ethz.ch mailing list
> https://www.stat.math.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
>

Thomas Lumley Assoc. Professor, Biostatistics
tlumley@u.washington.edu University of Washington, Seattle

______________________________________________
R-help@stat.math.ethz.ch mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


About this list Date view Thread view Subject view Author view Attachment view

This archive was generated by hypermail 2.1.3 : Mon 31 May 2004 - 23:05:14 EST