本文主要是介绍R语言求解混合线性方程组(有系谱),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!
用矩阵混合线性方程组计算:方差组分已知
数据
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 200
2 5 3 2 170
3 6 3 4 180
A-1
> Ainv <- makeAinv(pped)$Ainv;Ainv
6 x 6 sparse Matrix of class "dgCMatrix"1 1.5 . 0.5
这篇关于R语言求解混合线性方程组(有系谱)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!