Я использую пакет effects
, чтобы найти влияние переменных в моей линейной модели.
library(effects)
data(iris)
lm1 <- lm(Sepal.Length ~ Sepal.Width + Petal.Length*Petal.Width,data=iris)
Для простого термина в модели я могу получить эффекты для каждой точки данных, используя
effect("Sepal.Width", lm1, xlevels=iris['Sepal.Width'])
Как получить аналогичный одномерный вектор значений для моегосрок взаимодействия в каждой точке?Имеет ли это смысл?Все, что я пробовал, это возвращение двумерной матрицы, например
effect("Petal.Length:Petal.Width", lm1 ,xlevels=iris['Petal.Length']*iris['Petal.Width'])
Я не уверен, что в этом случае следует использовать для аргумента xlevels
, чтобы дать мне больше, чем просто значение по умолчанию5 одинаково расположенных точек.