1, 数据
这次使用一个PPT里面的数据, 用R语言演示一下如何做BLUP值计算....ID Sire Dam
1 NA NA
2 NA NA
3 1 NA
4 1 2
5 3 2
如果是计算逆矩阵的矩阵形式, 可以使用makeAinv(pped)$Ainv
Ainv = makeAinv...-1.0 -1.0000000 . 2
如果是计算逆矩阵的行列形式, 可以使用makeAinv(pped)$listAinv
makeAinv(pped)$listAinv
row..., 两者一样
3, 构建模型
$$ y = Xb + Zu + e $$
构建固定因子矩阵
这里使用函数model.matrix构建矩阵, 比较方便
for(i in 1:4) dat[,i] 值
solve(LHS)%*%RHS
7 x 1 Matrix of class "dgeMatrix"
[,1]
[1,] 142.842105
[2,] 151.118421