[R-sig-ME] Mixed-model-binary logistic model with dependence between individual repeated measures

[Ben Bolker]

-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1

On 11-01-07 06:59 AM, Anna Ekman wrote:
> Hi, I am a novice R user and do not know how to properly mail to this
> list. I apologies if I do it in the wrong way.
> 
> I want to analyze my data using a random intercept (later extended
> also to random slope) logistic model for a binary outcome (later
> extended to a ordinal outcome). This I have been able to do in SAS if
> assuming the repeated measurements within an individual to be
> independent, but I want to be able to choose different covariance
> structures for the individual measurements. This I cannot do directly
> in either SAS or STATA, and therefore now turn to R. How can I do
> this in R?
> 
> Anna

  I'm surprised that you can't do this in SAS (PROC MIXED, NLMIXED, or
GLIMMIX?) or Stata <http://www.gllamm.org/>, but: if you want to do it
in R, your choices are glmmPQL in the MASS package or possibly one of
the generalized estimating equation packages (geese, geepack?)  I would
recommend the following references for getting started:

  Zuur et al Mixed models (Springer)
  Pinheiro and Bates 2000 (Springer), especially the material on
temporal autocorrelation models

  Extending to a ordinal outcome with temporal autocorrelation could be
tricky ...




-----BEGIN PGP SIGNATURE-----
Version: GnuPG v1.4.10 (GNU/Linux)
Comment: Using GnuPG with Mozilla - http://enigmail.mozdev.org/

iEYEARECAAYFAk0nNJsACgkQc5UpGjwzenMymACfbqx+gtOjyhoX9hHpPyO/vbVg
NXUAniLvM+8voyzKA6axKLyJLclqeBhY
=wBwK
-----END PGP SIGNATURE-----