В линейной модели R получите значения p -только для коэффициентов взаимодействия

] Если у меня есть сводная таблица для линейной модели в R, как я могу получить значения p -, связанные только с оценками взаимодействия или только с групповыми пересечениями и т. д., без подсчета номеров строк?

Например, в такой модели, как lm(y ~ x + group)с xкак непрерывной и groupкак категориальной, сводная таблица для объекта lmимеет оценки для:

  1. перехват
  2. x, наклон по всем группам
  3. 5 внутригрупповых отличий от общего перехвата
  4. 5 внутригрупповых отличий от общего наклона.

Я хотел бы найти способ получить каждое из них как группу значений 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)
10
задан Sam Firke 1 June 2015 в 21:46
поделиться