R中多元线性回归斜率的比较

我正在调查尺寸随时间的变化。在一个多元线性模型中，我有五个大小变量。我在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)

names（iris）斜率（和截距）包含在
mod.iris$系数中

要比较SL和SW的斜率，请执行以下操作：

mod.iris$coefficients[1,'SL'] / mod.iris$coefficients[1,'S']

要比较SL和SW的斜率，请执行以下操作：

mod.iris$coefficients[1,'SL'] / iris$SW

谢谢我想我不清楚。我想测试这些坡度是否存在显著差异，对每个组合进行测试，并（可选）调整p值。我正在寻求多元回归中斜率的事后检验。我现在根据这个编辑了这个问题。