Сравнение наклонов многомерной линейной регрессии в 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 и т.д., значительно отличается или нет, сделать это для каждой комбинации, и (при необходимости) регулировать р- значения. Я ищу некоторый пост-hoc-тест для склонов в многомерной регрессии.

Благодаря,

Габора

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