在统计信息编程语言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/

10-12 23:33