[R-sig-ME] estimating variance components for arbitrarily definedvar/covar matrices

David Duffy David.Duffy at qimr.edu.au
Thu Feb 26 01:52:54 CET 2015

Hi Matt.

It depends if you want to do GLMMs or not. Look at the AnimalINLA, kinship 
(now kinship2) with nlme, and regress packages for starters, since you can 
feed them an arbitrary matrix out of the box.  I have used glmmPQL too. 
WRT Ben's comment about pedigreemm, you would have to look at their code 
but it should be _relatively_ straightforward.

CHeers, David.

| David Duffy (MBBS PhD)
| email: David.Duffy at qimrberghofer.edu.au  ph: INT+61+7+3362-0217 fax: -0101
| Genetic Epidemiology, QIMR Berghofer Institute of Medical Research
| 300 Herston Rd, Brisbane, Queensland 4006, Australia  GPG 4D0B994A

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