] Если у меня есть сводная таблица для линейной модели в R, как я могу получить значения p -, связанные только с оценками взаимодействия или только с групповыми пересечениями и т. д., без подсчета номеров строк?
Например, в такой модели, как lm(y ~ x + group)
с x
как непрерывной и group
как категориальной, сводная таблица для объекта lm
имеет оценки для:
Я хотел бы найти способ получить каждое из них как группу значений p -, даже если количество групп или формула модели изменятся. Может быть, сводная таблица каким-то образом использует информацию для группировки строк?
Ниже приведен пример набора данных с двумя разными моделями. В первой модели есть четыре разных набора значений p-, которые я мог бы получить по отдельности, тогда как во второй модели есть только два набора значений p-.
x <- 1:100
groupA <-.5*x + 10 + rnorm(length(x), 0, 1)
groupB <-.5*x + 20 + rnorm(length(x), 0, 1)
groupC <-.5*x + 30 + rnorm(length(x), 0, 1)
groupD <-.5*x + 40 + rnorm(length(x), 0, 1)
groupE <-.5*x + 50 + rnorm(length(x), 0, 1)
groupF <-.5*x + 60 + rnorm(length(x), 0, 1)
myData <- data.frame(x = x,
y = c(groupA, groupB, groupC, groupD, groupE, groupF),
group = rep(c("A","B","C","D","E","F"), each = length(x))
)
myMod1 <- lm(y ~ x + group + x:group, data = myData)
myMod2 <- lm(y ~ group + x:group - 1, data = myData)
summary(myMod1)
summary(myMod2)