【发布时间】：2014-03-02 12:08:26
【问题描述】：
我想用lm 来预测y by a/b（a 乘 b 的商）当我做lm(y~a/b) 我得到a 和b 交互时，我该怎么做？
【问题讨论】：
【发布时间】：2014-03-02 12:08:26
【问题描述】：
lm(y~I(a/b))
?formula 对此进行了介绍。
【讨论】：
-
I(something)告诉 R 将该部分评估为数学运算，而不是公式命令，例如I(log(x))获取x的对数，然后将系数拟合到x的对数。 -
@statquant a:b 与 seq(a,b,1) 完全不同
-
@lebatsnok 不...这是
a:b的标准解释，但在公式内部（如lm），它指定a和b的交互。a/b，当不在I()内部时，指定a及其与b的交互，但没有单独表示b的效果。 -
每当我对这些东西感到困惑时，我都会设置一些随机数据来测试它。试试
y <- rnorm(100); a <- rnorm(100); b <- rnorm(100); summary(lm(y ~ a +b)); summary(lm(y ~ a * b)); summary(lm(y ~ a/b)); summary(lm(y ~ I(a/b)))之类的东西 - 比较摘要会让您快速了解正在发生的事情。 -
比较
a和/或b的值与y的模型预测值也很有帮助。