Learn R Programming

lmDiallel (version 1.0.1)

zhang05: Data for diallel analysis from Zhang (2005)

Description

Data collected in XXX with 5 parents, 2 reps and 2 environments

Usage

data("zhang05")

Arguments

Format

A data.frame with 60 observations on the following 6 variables.

Par1

male parent, a factor with 5 levels

Par2

female parent, a factor with 5 levels

Block

block, a factor with 2 levels

Combination

combination between environment and block, an integer vector

Env

environment, a factor with 2 levels

Yield

yield, a numeric vector

Author

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

Run this code
data("zhang05")
contrasts(zhang05$Block) <- c("contr.sum")
contrasts(zhang05$Env) <- c("contr.sum")
dMod <- lm(Yield ~ Env/Block + H.BAR(Par1, Par2) + VEi(Par1, Par2) +
                   Hi(Par1, Par2) + SCA(Par1, Par2) +
                   H.BAR(Par1, Par2):Env + GCA(Par1, Par2):Env +
                   Hi(Par1, Par2):Env + SCA(Par1, Par2):Env, data = zhang05)
anova(dMod)
#or
dMod2 <- lm(Yield ~ Env/Block + H.BAR(Par1, Par2) + VEi(Par1, Par2) +
                   Hi(Par1, Par2) + SCA(Par1, Par2), data = zhang05)
summary(dMod2)$coefficients

Run the code above in your browser using DataLab