我的首要问题是，在 WLS 案例中，R 如何计算 R^2？它不仅对观察结果进行加权，然后计算 R^2。为了弄清楚这一点，我一直在查看源代码，直到在lm.wfit 代码中遇到这个问题：
z <- .Call(C_Cdqrls, x *wts, y*wts, tol)
这里正在做什么？有谁知道我如何访问此代码以获取详细信息？即，返回给z 的内容是什么？ C_Cdqrls、x*wts、y*wts、tol是如何使用的？
到目前为止我的理解（我不确定它是否正确）是.Call 表示 R 正在 C 中执行此代码。但是，我想看看这是如何在 C 中完成的，如果可能的。
谢谢！
【问题讨论】：
R squared的值实际上是在调用summary.lm时计算出来的，你可以查看任何函数的源代码，无论是在实际的svn存储库中（https://svn.r-project.org/R/），还是在github上的this read only mirror。
在https://github.com/wch/r-source/blob/trunk/src/library/stats/R/lm.R 寻找
summary.lm
我们看到以下对权重的计算 (w)
r <- z$residuals
f <- z$fitted.values
w <- z$weights
if (is.null(w)) {
mss <- if (attr(z$terms, "intercept"))
sum((f - mean(f))^2) else sum(f^2)
rss <- sum(r^2)
} else {
mss <- if (attr(z$terms, "intercept")) {
m <- sum(w * f /sum(w))
sum(w * (f - m)^2)
} else sum(w * f^2)
rss <- sum(w * r^2)
r <- sqrt(w) * r
}
# ..... some other code
# ... then this definition
ans$r.squared <- mss/(mss + rss)
【讨论】：
【讨论】：