Skip to contents

Formula: GY_ij = beta-bar + (beta_i - beta-bar) + e_ij (eq. 1.4), the fixed-effects fit rewritten around the average level.

Usage

Msm005(groups)

Arguments

groups

List of equal-length observation vectors, one per level.

Value

List with estimate, beta_bar, deviations, deviations_sum, method.

References

Montesinos Lopez, Montesinos Lopez & Crossa (2022), Multivariate Statistical Machine Learning Methods for Genomic Prediction, Springer, eq. (1.4) p.17. DOI 10.1007/978-3-030-89010-0.

Examples

D <- data.frame(x = c(1, 2, 3, 4), y = c(2, 4, 5, 9))
Msm005(D)
#> $estimate
#> [1] 3.75
#> 
#> $beta_bar
#> [1] 3.75
#> 
#> $deviations
#>     x     y 
#> -1.25  1.25 
#> 
#> $deviations_sum
#> [1] 0
#> 
#> $method
#> [1] "reparameterized one-way model (MVSML 2022 eq. 1.4)"
#>