如何在lm上使用lm.fit

Question

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

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

或

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

编辑1

我现在看到总和$系数只有1个值。为什么这是什么呢？帮助没有解释这个

Answer 1

通过可重复的示例提供帮助会更容易。

然而，我的猜测是你错过了第二种情况下的拦截。试试lm.fit(cbind(1,x),y)，看看是否能为您提供所需的比较结果。

其他差异可能是由于lm在调用lm.fit之前执行的其他预处理，但我们不知道在没有看到x的情况下是什么。删除缺失值，将因子扩展为虚拟变量等事情也可能导致差异。