Multi-trait GBLUP via vec-stacked mixed-model equations

Usage

morie_multi_trait_gblup(x, y, markers, Sigma_g = NULL, Sigma_e = NULL)

Arguments

x: Fixed-effect design (vector or matrix).
y: Multi-trait response (n x t).
markers: Genotype matrix (n x m).
Sigma_g: Optional t x t genetic covariance.
Sigma_e: Optional t x t residual covariance.

Value

list(estimate, G_hat, B_hat, Sigma_g, Sigma_e, n, t, method).

References

Montesinos Lopez Ch 10.

Examples

morie_multi_trait_gblup(
  x = rnorm(50), y = rnorm(50),
  markers = matrix(sample(0:2, 200, TRUE), 50, 4)
)
#> $estimate
#> [1] 3.441244e-17
#> 
#> $G_hat
#>                [,1]
#>  [1,]  0.0016614242
#>  [2,]  0.0327492645
#>  [3,]  0.0362280324
#>  [4,] -0.0322057049
#>  [5,]  0.0361337488
#>  [6,]  0.0334452445
#>  [7,] -0.0350193491
#>  [8,]  0.0052700024
#>  [9,]  0.0670960028
#> [10,] -0.0353705225
#> [11,]  0.0644075078
#> [12,] -0.0260391010
#> [13,]  0.0052676642
#> [14,]  0.0733888679
#> [15,]  0.0393928663
#> [16,]  0.0984056148
#> [17,]  0.0357870963
#> [18,]  0.0016642302
#> [19,] -0.1004535850
#> [20,] -0.0292984292
#> [21,]  0.0046988798
#> [22,] -0.0263898036
#> [23,] -0.0322056191
#> [24,]  0.0023612591
#> [25,]  0.0701314941
#> [26,] -0.0663308914
#> [27,]  0.0392667048
#> [28,]  0.0046993501
#> [29,]  0.0390439904
#> [30,] -0.0664562200
#> [31,]  0.0733885159
#> [32,]  0.0360091529
#> [33,] -0.0256921449
#> [34,]  0.0300630289
#> [35,] -0.0319842682
#> [36,] -0.0350226828
#> [37,] -0.0664568125
#> [38,]  0.0670950732
#> [39,] -0.0632930618
#> [40,] -0.0321129343
#> [41,] -0.0006762566
#> [42,] -0.0603831807
#> [43,]  0.0392656971
#> [44,] -0.0350205887
#> [45,]  0.0304105640
#> [46,] -0.0008988025
#> [47,] -0.0039340843
#> [48,] -0.0634221366
#> [49,] -0.0693689096
#> [50,] -0.0292961881
#> 
#> $B_hat
#>             [,1]
#> [1,] -0.03625982
#> [2,]  0.07579978
#> 
#> $Sigma_g
#>           [,1]
#> [1,] 0.5541933
#> 
#> $Sigma_e
#>           [,1]
#> [1,] 0.5541933
#> 
#> $n
#> [1] 50
#> 
#> $t
#> [1] 1
#> 
#> $method
#> [1] "Multi-trait GBLUP (vec-stacked MME)"
#>