methods: Methods for diallel model fitting

In lmDiallel: Linear Fixed/Mixed Effects Models for Diallel Crosses

Methods for diallel model fitting

Description

The object returned by the 'lm.diallel()' function is of classes 'lm' and 'diallel'. Specific methods were devised to explore the 'diallel' object.

Usage

## S3 method for class 'diallel'
summary(object, MSE, dfr, ...)
## S3 method for class 'diallel'
vcov(object, MSE, ...)
## S3 method for class 'diallel'
anova(object, MSE, dfr, type = 1, ...)
## S3 method for class 'diallel'
model.matrix(object, ...)

Arguments

object	an object of class diallel.
MSE	Mean Square Error, when it cannot be derived from model fit
dfr	Residual degrees of freedom, when they cannot be derived from model fit
type	It is used to select between Type I (sequential) or Type III (marginal) F tests in ANOVA
...	Other optional arguments

Details

To be defined

Value

vcov.diallel: a variance-covariance matrix summary.diallel: a data.frame of estimated parameters with standard errors anova.diallel: an ANOVA table predict.diallel: a vector of predictions from a diallel model model.matrix.diallel: a design matrix for the fitted diallel model

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

library(lmDiallel)
library(multcomp)
# Data with replicates
data("hayman54")

fit <- lm.diallel(Ftime ~ Par1 + Par2, data = hayman54,
                  fct = "HAYMAN1")
summary(fit)
anova(fit)
gh <- glht(linfct = diallel.eff(fit), adjust = "none")

lmDiallel documentation built on April 19, 2023, 5:10 p.m.