GBLUP – full mixed-model implementation — morie_gblup

Solves Henderson's MME with VanRaden G.

Usage

morie_gblup_full(x, y, markers, lambda_gblup = NULL)

Arguments

x: Fixed-effect design (vector or matrix).
y: Numeric response.
markers: Genotype matrix (n x m).
lambda_gblup: Optional ratio sigma_e^2 / sigma_g^2 (default h^2=0.5).

Value

Named list (estimate, g_hat, beta, se, lambda_gblup, n, method).

References

Montesinos Lopez Ch 3.

Examples

morie_gblup_full(x = rnorm(50), y = rnorm(50), markers = matrix(sample(0:2, 200, TRUE), 50, 4))
#> $estimate
#> [1] -2.456694e-17
#> 
#> $g_hat
#>  [1] -0.17858809 -0.14278764 -0.03389322  0.66054445  0.19185823  0.15726347
#>  [7]  0.66054546 -0.44284104  0.44275347  0.21309953 -0.03389191 -0.65793495
#> [13]  0.05106459  0.15996227 -0.26624370 -0.44284050 -0.59023538  0.66054525
#> [19] -0.10819244 -0.22774460 -0.11089011  0.40695264 -0.73493045 -0.80532766
#> [25] -0.35788268  0.32199717  0.14270105  0.21309907 -0.76682884  0.26615997
#> [31]  0.25159921  0.12146100 -0.30204493  0.73094306  0.23036014 -0.66063038
#> [37] -0.48133909  0.25159678 -0.22774363 -0.11089240  0.44275591  0.10690133
#> [43]  0.44275259 -0.44283897  0.25159842  0.44545142  0.07230328 -0.15734776
#> [49]  0.23035777  0.15726282
#> 
#> $beta
#> [1] 0.1770434 0.2167859
#> 
#> $se
#> [1] 1.081721
#> 
#> $y_hat
#>  [1] -0.0227388137 -0.6673470137  0.3506280504  0.9588435787  0.2610071162
#>  [6]  0.2922546543  0.8563050274 -0.0452163914  0.4095323115  0.6213130606
#> [11]  0.3069351904 -0.6333467466  0.1659225982  0.7361316307 -0.1231846181
#> [16] -0.5670867215 -0.7325401611  0.8225179057  0.1207151020  0.0035862770
#> [21]  0.0088298815  0.1557571903 -0.6554138440 -0.4274396047 -0.6837762995
#> [26]  0.6321326783  0.0003884412  0.3426556734 -0.8940263063  0.6211040063
#> [31]  0.3001478902  0.3899812574 -0.2984632390  0.6514851284  0.0801528012
#> [36] -0.7341791531 -0.3812232875  0.7289864206  0.0727939673 -0.0722657438
#> [41]  0.6717855092  0.0669909333  0.4864105068 -0.2316019159  1.0105082929
#> [46]  0.6722400976  0.2710804499 -0.4794566042  0.6623185580  0.6049351164
#> 
#> $lambda_gblup
#> [1] 1
#> 
#> $n
#> [1] 50
#> 
#> $method
#> [1] "GBLUP with VanRaden G"
#>