R中多元線性回歸的斜率比較

Question

我正在調查隨時間變化的尺寸。 在針對年份的多元線性模型中，我有五個大小變量。 我在汽車包裝中使用了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

Answer 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