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

[David Duffy]

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