Skip to contents

Ordinal predictor with environment and genetic effects

Source: R/msm098.R
Msm098.Rd

Formula: L = X_E beta_E + Z_L g + eps (eq. 7.5), environment effects and genetic effects without the line-by-environment interaction. With L_g the Cholesky factor of G, the genetic block enters as Z_L L_g, the design used in Table 7.6.

Usage

Msm098(n, X_E, Z_L, L_g = NULL)

Arguments

n

Number of rows to take from each block.

X_E

Environment design block.

Z_L

Line incidence matrix.

L_g

Cholesky factor of G, so that G = L_g L_g'.

Value

List with estimate (the total column count), design, widths, method.

References

Montesinos Lopez, Montesinos Lopez & Crossa (2022), Multivariate Statistical Machine Learning Methods for Genomic Prediction, Springer, eq. (7.5) p.220 and Table 7.6 p.233. DOI 10.1007/978-3-030-89010-0.

Examples

set.seed(17)
G <- crossprod(matrix(rnorm(16), 4)) / 4 + diag(0.1, 4)
XE <- model.matrix(~ factor(c(1, 1, 2, 2)) - 1)
Msm098(n = 4L, X_E = XE, Z_L = diag(4), L_g = t(chol(G)))$estimate
#> [1] 6