如何在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.删除缺失值,将因子扩展为虚拟变量等因素也可能导致差异.