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

Jarrod Hadfield j.hadfield at ed.ac.uk
Fri Jan 7 19:46:51 CET 2011


Hi All,

I think one issue is going to be constraining the within-individual  
residual structure to a correlation matrix, since the diagonal  
elements of the covariance matrix cannot be estimated from binary  
data. A residual structure of the form rcov=~cor(time):units  does  
this in MCMCglmm, but I would not advocate it if there are many time  
points because there will be N*(N-1)/2 parameters, rather than 1  
parameter with an appropriate AR structure. N is the number of time  
points.

Cheers,

Jarrod



Quoting Gavin Simpson <gavin.simpson at ucl.ac.uk>:

> On Fri, 2011-01-07 at 11:49 -0500, Ben Bolker wrote:
>> -----BEGIN PGP SIGNED MESSAGE-----
>> Hash: SHA1
> <snip />
>>  I do not
>> > want to assume that. In addition I would like to be able to chose
>> > other distributions than the normal for my random effect, which is
>> > not possible in SAS (proc NLMIXED).
>>
>>   It's not possible in R either as far as I know.
>
> I was reading the article in the latest issue of the R Journal on the
> hglm package, and although I was only giving it a cursory scan over
> lunch it looks like it might be able to fit the sort of model implied
> here; random effects distributed as a member of the exponential family.
>
> G
>
>>
>> The generalized estimating
>> > equation packages are probably not an option as I do not whant
>> > marginal models. I will look at the references you suggested. Thank
>> > you. /Anna
>> >
>>
>>   If you want a non-marginal model with non-normal random effects and
>> within-individual correlation structures other than compound symmetry
>> (i.e. simple block structures), you are probably going to have to
>> construct your own solution with WinBUGS or AD Model Builder or ... ? If
>> you're lucky, MCMCglmm may be able to do what you want -- I would check
>> it out. (Molenbergh and Verbeke's book on longitudinal models describes
>> approaches for non-normal random effects, but in the context of LMMs
>> (i.e. normally distributed errors) -- they may have done something to
>> extend this stuff to GLMMs more recently.  It's possible that someone
>> out there has done what you want and encapsulated it into a canned
>> package, but I doubt it.
>>
>>    cheers
>>     Ben Bolker
>> -----BEGIN PGP SIGNATURE-----
>> Version: GnuPG v1.4.10 (GNU/Linux)
>> Comment: Using GnuPG with Mozilla - http://enigmail.mozdev.org/
>>
>> iEYEARECAAYFAk0nRBsACgkQc5UpGjwzenPr6wCfb7wixKYRGfglmge7Ejg0+i26
>> k5QAoJpuc+QiWC8iaMXC6aDnW75GfrYE
>> =292C
>> -----END PGP SIGNATURE-----
>>
>> _______________________________________________
>> R-sig-mixed-models at r-project.org mailing list
>> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>
> --
> %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%
>  Dr. Gavin Simpson             [t] +44 (0)20 7679 0522
>  ECRC, UCL Geography,          [f] +44 (0)20 7679 0565
>  Pearson Building,             [e] gavin.simpsonATNOSPAMucl.ac.uk
>  Gower Street, London          [w] http://www.ucl.ac.uk/~ucfagls/
>  UK. WC1E 6BT.                 [w] http://www.freshwaters.org.uk
> %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%
>
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>
>



-- 
The University of Edinburgh is a charitable body, registered in
Scotland, with registration number SC005336.




More information about the R-sig-mixed-models mailing list