[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