Sorry, I sent it quickly and forgot to thank in advance
Marcio
Marcio Resende wrote:
>
> hello guys, I need to do a BLUP in the simplest model
> y = Xm + Zg + e
> however I have missing data in the analysis which I can´t consider as
> 0(zero). So I need to generate the matrix X'Z, Z'X and Z'Z
hello guys, I need to do a BLUP in the simplest model
y = Xm + Zg + e
however I have missing data in the analysis which I can´t consider as
0(zero). So I need to generate the matrix X'Z, Z'X and Z'Z step by step; I
can´t use
crossprod(x) #neither
X'X <- t(x)%*%x
because I should skip the element
2 matches
Mail list logo