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")
Run the code above in your browser using DataLab