Skip to content
Prev 687 / 20628 Next

use of pedigree structures

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).