r - lm の代わりに lm.fit を使用する方法

Question

速度のために lm.fit を使用したいのですが、2 番目のバージョンでは 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

再現可能な例を支援する方が簡単です。

ただし、2 番目のケースではインターセプトが欠落していると思います。lm.fit(cbind(1,x),y)探している比較が得られるかどうか試してみてください。

他の違いは、 をlm呼び出す前に行う他の前処理が原因である可能性がありますlm.fitが、 を見なければそれが何であるかはわかりませんx。欠損値を削除したり、因子をダミー変数に展開したりすることでも、違いが生じる可能性があります。