r - 如何在 lm 上使用 lm.fit

Question

我想使用 lm.fit 来提高速度，但第二个版本提供了 NA

 sum <-   summary(lm(y~x))
 slope <-  sum$coefficients[2]

或者

 sum <-   lm.fit(as.matrix(x,ncol=1),y)
 slope <-  sum$coefficients[2]

编辑 1

我现在看到 sum$coefficients 只有 1 个值。为什么会这样，它是什么？帮助没有解释这一点

Answer 1

提供可重现的示例会更容易。

但是，我的猜测是您在第二种情况下错过了拦截。尝试lm.fit(cbind(1,x),y)看看这是否为您提供了您正在寻找的比较。

其他差异可能是由于lm在调用之前进行的其他预处理lm.fit，但我们不知道那是什么x。诸如删除缺失值、将因子扩展为虚拟变量之类的事情也可能导致差异。