在统计信息编程语言R中,以下公式(例如lm()
或glm()
)
z ~ (x+y)^2
相当于
z ~ x + y + x:y
假设我只有连续的预测变量,是否有一种简明的方法来获得
z ~ I(x^2) + I(y^2) + I(x) + I(y) + I(x*y)
一个对因子预测变量做正确的事情的公式是加号。
一种可能的解决方案是
z ~ (poly(x,2) + poly(y,2))^2
我正在寻找更优雅的东西。
最佳答案
我不知道它是否更优雅,但是poly函数可以采用多个向量:
z ~ poly(x, y, degree=2)
这将创建您要求的所有组合,而无需其他组合。请注意,以这种方式进行操作时,您需要指定度数= 2,而不仅仅是2。
关于r - 所有一阶和二阶预测变量的公式,包括R中的相互作用,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/3661302/