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。