我正在研究尺寸随时间的变化。我在多变量线性模型中有五个大小变量与年份。我在汽车包装中使用了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值。我正在寻求多变量回归中斜率的一些事后检验。
谢谢，
的Gabor
答案 0 :(得分：1)
斜率（和截距）包含在mod.iris$coefficients中。
为了将SL的斜率与SW的斜率进行比较，请执行：
mod.iris$coefficients[1,'SL'] / mod.iris$coefficients[1,'S']
并且比较SL与SW的斜率，执行：
mod.iris$coefficients[1,'SL'] / iris$SW