超级会员免费看
这篇博客介绍了如何在R语言中解决混合线性方程组，特别是当方差组分已知的情况。通过矩阵运算求解方程，并使用asreml库来计算最佳线性无偏预测(BLUP)值，同时详细说明了方差组分的定义和计算过程，以及查看拟合模型的结果和BLUP值的标准误计算方法。
用矩阵混合线性方程组计算：方差组分已知
数据
Y = Xb + Za + e
> dat <- data.frame(id=c(4,5,6),sire = c(1,3,3),dam=c(2,2,4),y=c(200,170,180))
> dat
id sire dam y
1 4 1 2 200
2 5 3 2 170
3 6 3 4 180
> for( i in 1:3) dat[,i] <- as.factor(dat[,i])
> str(dat)
'data.frame': 3 obs. of 4 variables:
$ id : Factor w/ 3 levels "4","5","6": 1 2 3
$ sire: Factor w/ 2 levels "1","3": 1 2 2
$ dam : Factor w/ 2 levels "2","4": 1 1 2
$ y : num 200 170 180
> dat
id sire dam y
1 4 1 2 20
订阅专栏 解锁全文
扫一扫
657
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
