如何以与R&m; lm.fit()相同的方式计算效果?

时间:2015-09-22 08:43:39

标签: r lm

非统计学家,将一些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$pivotsz$pivoted是什么 - 我需要担心这些吗?

聚苯乙烯。有一个问题似乎非常相似,但问题和答案似乎并没有包括如何计算效果(至少我不能弄清楚如何:https://stats.stackexchange.com/questions/130704/how-the-effectsqy-in-lm-fit-calculated。)

0 个答案:

没有答案
相关问题