9
如果我在R中有一個線性模型的彙總表，我怎樣才能得到與相互作用估計值相關聯的p值，或者僅僅是羣組截取值等，而不必計算行號？在R線性模型中，僅爲相互作用係數獲得p值
例如，與這樣的模型如lm(y ~ x + group)與x如連續和group作爲分類，對於lm對象彙總表具有用於估算值：
- 截距
- X，所有的斜率組
- 5組內差異從整體截距
- 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)
謝謝，這是一個很好的方法來做到這一點。請注意，如果我使用與默認值不同的對比度，則行名稱爲group1，group2，group3等，而不是groupA，groupB，groupC等。如果它們是不依賴於其他方法的附加方法瞭解組級別名稱和正在使用的對比。 – Jdub
我不確定我是否正確理解你。如果你想讓它工作而不考慮因素級別的名稱，你可以嘗試類似於％paste0（'group（％）'）的％summary（myMod1）$ coefficients [names（summary（myMod1）$ coefficients [，4]）％ ' '，levels（myData $ group）），4]' – RoyalTS
如果這回答你的問題，你會好好接受它（點擊答案旁邊的綠色複選標記）？ – RoyalTS