H.BAR {lmDiallel}

Average heterosis effect

Description

H.BAR effect to fit GE2 and GE3 models with lm function

Usage

H.BAR(P1, P2, type = "fix", data)

Arguments

P1

a variable for the first parent

P2

a variable for the second parent

type

a variable for model selection. May be "fix" (fixed model) or "random" (random model).

data

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

Value

A design matrix for the H.BAR effect

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("zhang05")
dMod <- lm(Yield ~ Env/Block + H.BAR(Par1, Par2) + VEi(Par1, Par2) +
                   Hi(Par1, Par2) + SCA(Par1, Par2) +
                   H.BAR(Par1, Par2):Env + VEi(Par1, Par2):Env +
                   Hi(Par1, Par2):Env + SCA(Par1, Par2):Env, data = zhang05)
anova(dMod)