[R-sig-ME] use of pedigree structures
ian white
i.m.s.white at ed.ac.uk
Thu Mar 27 13:40:48 CET 2008
The lmer function includes random terms as
(expression for a model matrix|grouping factor)
and I wonder how general the expression can be. For example, a pedigree
object can be converted to an A matrix (in the language of animal
breeding), but can the A matrix be incorporated into an lmer analysis?
I thought that
(chol(pedigree)|ones)
might be OK, where 'ones' is a constant vector of value 1. With the only
example I have tried, the syntax is accepted but the analysis ends with
a forced exit from R (core dump).
--
I White
University of Edinburgh
Ashworth Laboratories, West Mains Road
Edinburgh EH9 3JT (Tel 0131 650 5490)
More information about the R-sig-mixed-models
mailing list