fcm: fixed effect constraint indication matrix

In sommer: Solving Mixed Model Equations in R

fixed effect constraint indication matrix

Description

fcm creates a matrix with the correct number of columns to specify a constraint in the fixed effects using the Gtc argument of the vsr function.

Usage

  fcm(x, reps=NULL)

Arguments

x	vector of 1's and 0's corresponding to the traits for which this fixed effect should be fitted. For example, for a trivariate model if the fixed effect "x" wants to be fitted only for trait 1 and 2 but not for the 3rd trait then you would use fcm(c(1,1,0)) in the Gtc argument of the vsr() function.
reps	integer specifying the number of times the matrix should be repeated in a list format to provide easily the constraints in complex models that use the ds(), us() or cs() structures.

Value

$res

a matrix or a list of matrices with the constraints to be provided in the Gtc argument of the vsr function.

Author(s)

Giovanny Covarrubias-Pazaran

References

Covarrubias-Pazaran G (2016) Genome assisted prediction of quantitative traits using the R package sommer. PLoS ONE 11(6): doi:10.1371/journal.pone.0156744

See Also

The function vsr to know how to use fcm in the mmer solver.

Examples

fcm(c(1,1,0))
fcm(c(0,1,1))
fcm(c(1,1,1))

fcm(c(1,1,1),2)

# ## model with Env estimated for both traits
# data(DT_example)
# DT <- DT_example
# A <- A_example
# ans4 <- mmer(cbind(Yield, Weight) ~ Env,
#               random= ~ vsr(Name) + vsr(Env:Name),
#               rcov= ~ vsr(units),
#               data=DT)
# summary(ans4)$betas
# ## model with Env only estimated for Yield
# ans4b <- mmer(cbind(Yield, Weight) ~ vsr(Env, Gtc=fcm(c(1,0))),
#              random= ~ vsr(Name) + vsr(Env:Name),
#              rcov= ~ vsr(units),
#              data=DT)
# summary(ans4b)$betas

sommer documentation built on Sept. 11, 2024, 6:22 p.m.