model.matrixDiallel {lmDiallel}

Incidence matrices for Diallel model parametrisation

Description

model.matrixDiallel is useful to build design matrices, according to the user-defined (or default) parameterisation for lm function. It shares the same syntax of the lm.diallel function.

Arguments

formula

an object of class "formula" (or one that can be coerced to that class): a symbolic description of the model to be fitted.'formula' uses the regular R syntax to specify the response variable and the two variables for parentals (e.g., Yield ~ Par1 + Par2)

Block

used to specify an optional variable coding for blocks

Env

used to specify an optional variable coding for environments

data

a 'data.frame' where to look for explanatory variables

fct

a string variable coding for the selected model. 6 main diallel models: Hayman's model 1 (="HAYMAN1"), Hayman's model 2 (="HAYMAN2"), Griffing's model 1 (="GRIFFING1"), Griffing's model 2 (="GRIFFING2"), Gardner-Eberhart model 2 (="GE2") and Gardner-Eberhart model 3 (="GE3"). The strings "GE2r" and "GE3r" can be used to specify the 'enhanced' GE2 and GE3 models, including the effect of reciprocals (REC).

Details

model.matrixDiallel creates a design matrix for a diallel model, as specified in the 'fct' argument.

Value

The design matrix for a diallel model as specified in the 'fct' argument.

Author(s)

Andrea Onofri, Niccolo' Terzaroli, Luigi Russi

References

Onofri, A., Terzaroli, N. & Russi, L. Linear models for diallel crosses: a review with R functions. Theor Appl Genet (2020). https://doi.org/10.1007/s00122-020-03716-8

Examples

data("diallelMET")
ModMat <- model.matrixDiallel(Yield ~ Par1 + Par2,
                   Env, Block, fct= "GE3",
                   data = diallelMET)