#-------------------------------------------------------------
## Example 1.3.2 p-16
#-------------------------------------------------------------
# PROC GLM DATA=ex124;
# CLASS herd dose drug;
# MODEL PCVdif=drug herd(drug) dose dose*drug;
# RANDOM herd(drug);
# RUN;
str(ex124)
summary(ex124)
ex124$herd1 <- factor(ex124$herd)
ex124$drug1 <- factor(ex124$drug)
ex124$dose1 <- factor(ex124$dose)
fm1.1 <-
aov(
formula = PCVdif ~ drug1 + Error(herd1:drug1) + dose1 + dose1:drug1
, data = ex124
, projections = FALSE
, qr = TRUE
, contrasts = NULL
# , ...
)
if (requireNamespace("report", quietly = TRUE)) {
fm1.1 |>
report::report()
}
summary(fm1.1)
Run the code above in your browser using DataLab