我的R公式是否与我想到的统计模型相当？

Question

问题：

使用formula构建统计模型是R语言强大而优雅的功能。 我没有尽可能多地使用formula的一个原因是语法有点令人困惑（例如x*y并不仅仅意味着“x和{y的产品{1}}“）。

问题：

我正在寻找一种方法来确保我正确使用formula语法并且我输入的formula确实实现了我想到的统计模型。理想情况下，我希望在实际拟合模型之前得到此确认。

示例：

说，我想通过线性回归找到模型a的参数b和y = a + b*(x1*x2)。天真地，我在R

中输入

df <- data.frame(y=seq(5), x1=runif(5), x2=runif(5)) # toy data
lm(y~x1*x2, data=df)    # this is wrong

我可以从lm的输出中看出，由于x1和x2的额外系数，这不是我想要的。但是在调用拟合函数之前应该可以调试公式。 （适合此模型的正确方法是lm(y~x1:x2, data=df)）

Answer 1

在运行模型之前调试公式的一种方法是使用formula和update：

f <- formula( y ~ x1*x2)
update( f , terms( f ) )
# y ~ x1 + x2 + x1:x2

f <- formula( y ~ x1:x2)
update( f , terms( f ) )
# y ~ x1:x2

巧合的是，您还可以通过包含1（1 * a = a）来指定模型中的截距项（即a的系数），这是等效的：< / p>

f <- formula( y ~ 1 + x1:x2)
update( f , terms( f ) )
# y ~ x1:x2