非统计学家,将一些R代码移植到Python。
我需要在Python中使用R {#1}}。
移植R代码非常简单,但由于lm.fit
调用了C函数lm.fit
,我遇到了一个问题。
我已开始使用C_Cdqrls
复制lm.fit
:
numpy.linalg.lstsq
这给了我最小二乘解,残差之和,等级和单数。
但是,我想计算lstsq_sol, residuals, rank, singulars = lstsq(x, y, 1e-07)
- 我该怎么做?
代码在这里:https://github.com/wch/r-source/blob/trunk/src/library/stats/R/lm.R#L114
(fit$effects
代码中提到的z$pivots
和z$pivoted
是什么 - 我需要担心这些吗?
聚苯乙烯。有一个问题似乎非常相似,但问题和答案似乎并没有包括如何计算效果(至少我不能弄清楚如何:https://stats.stackexchange.com/questions/130704/how-the-effectsqy-in-lm-fit-calculated。)