r - lm.fit() 的困难：未找到“C_Cdqrls”

Question

我试图重现 lm.fit 函数正在做什么。但我有一个问题，它调用'C_Cdqrls'例程：

n<-100
p<-5
x<-cbind(1,matrix(rnorm(n*(p-1)),n,p-1))
res<-rnorm(n)
y<-x%*%bet+res
tol<-1e-07
z<-.Call(C_Cdqrls, x, y,tol)

我得到：

Error: object 'C_Cdqrls' not found

使用 'C_Cdqrls' 代替产生：

the Error in .Call("C_Cdqrls", x, y, tol) : 
  C symbol name "C_Cdqrls" not in load table

同样，summary.lm() 在 lm-object 上调用一个名为 qr.lm 的函数，但执行以下操作：

qr.lm(lm(y~x-1))

产量：

Error: could not find function "qr.lm"

我的问题是：如何让 .Call(C_Cdqrls, x, y,tol) 和 qr.lm() 分别给出与 lm.fit() 和 summary.lm() 调用时相同的结果？

Answer 1

C_Cdqrls是stats包命名空间中未导出的对象。要查看它，请使用stats:::C_Cdqrls.