我正在调查尺寸随时间的变化。我在针对年份的多元线性模型中有五个大小变量。我使用 car 包中的 Anova() 函数来测试每个尺寸测量的斜率是否相等。由于它们被证明是显着不同的，我想对大小测量的每个组合进行成对测试（记住，这些是响应变量），但我没有设法弄清楚如何去做。到目前为止，我的代码与下面的代码几乎相同，但不是分类变量，而是数字作为解释变量。
names(iris) <- c("SL", "SW", "PL", "PW", "SPP")
mod.iris <- lm(cbind(SL, SW, PL, PW) ~ SPP, data=iris)
summary(mod.iris)
manova.iris <- Anova(mod.iris)
summary(manova.iris)
因此，使用此示例，我想测试 SL 到 SW、SL 到 PL、SL 到 PW 等的斜率是否显着不同，对每个组合都这样做，并（可选地）调整 p 值。我正在为多元回归中的斜率寻找一些事后检验。
谢谢，
嘉宝