在R中指示交互的问题

Question

leaps4 <- regsubsets(Y~X1+X2+X3+X4+X5+X6+X7 + X1:X2 + X1:X3 + X1:X4 + X1:X5 + 
                       X1:X6 + X1:X7 + X2:X3 + X2:X4 + X2:X5 + X2:X6 + X2:X7 + 
                       X3:X4 + X3:X5 + X3:X6 + X3:X7 + X4:X5 + X4:X6 + X4:X7 + 
                       X5:X6 + X5:X7 + X6:X7, 
                     data=prostateData, nbest=1)

我们有97个观察值和7个变量：

p = rowSums(summary(leaps4)$which)
print (p)
1 2 3 4 5 6 7 8
2 3 4 5 6 7 8 9

我们希望看到不同变量之间的所有交互作用，但我们写的p值只表示所有变量之间的一个交互作用。那么，我们如何编写代码，以便R程序可以指示所有的交互。(我们所写的只适用于summary()和anova()，而不计算PRESSp BICp AICp和R^2p。)

Answer 1

看起来，如果您想要显示比“最佳”子集更多的内容，您应该增加nbest参数。如果希望公式包含所有单向和双向交互，只需使用"^2“

library(leaps)
b<-regsubsets(Y  ~  (X1+X2+X3+X4+X5+X6+X7 )^2, 
                 data=prostateData, nbest=10)

Answer 2

这个公式适用于lm和glm模型，我尝试将其成功地与regsubsets中的内置数据集swiss一起使用。

leaps4 <- regsubsets(Y~.*., data=prostateData, nbest=1)

您可能还想看看this post，它用glm描述它。

This site还对R中的公式调用进行了很好的细分。

原则证明：

b<-regsubsets(肥力~.*，data=swiss，nbest=2)