Fatal编程技术网
  • r/
  • 我如何以与R'相同的方式计算效果；s lm.fit（）？

我如何以与R'相同的方式计算效果；s lm.fit（）？

r

我如何以与R'相同的方式计算效果；s lm.fit（）？,r,lm,R,Lm,这里是非统计学家，正在将一些R代码移植到Python 我需要Python中的Rlm.fit 移植R代码很简单，但由于lm.fit调用C函数C_Cdqrls我遇到了麻烦 我已经开始使用numpy.linalg.lstsq复制lm.fit： lstsq_sol, residuals, rank, singulars = lstsq(x, y, 1e-07) 这给了我最小二乘解，残差和，秩和奇点 但是，我想计算fit$effects——我该怎么做 代码如下： （在lm.fit代码中引用的z$pivo

这里是非统计学家，正在将一些R代码移植到Python

我需要Python中的R
lm.fit

移植R代码很简单，但由于
lm.fit
调用C函数
C_Cdqrls
我遇到了麻烦

我已经开始使用
numpy.linalg.lstsq
复制
lm.fit
：

lstsq_sol, residuals, rank, singulars = lstsq(x, y, 1e-07)
这给了我最小二乘解，残差和，秩和奇点

但是，我想计算
fit$effects
——我该怎么做

代码如下：

（在
lm.fit
代码中引用的
z$pivots
和
z$pivoted
是什么？我需要担心这些吗

还有一个问题看起来很相似，但问题和答案似乎不包括如何计算影响（至少我不知道如何：）

链接问题中提到的Q可能来自这里：你找到方法了吗？因为我想在scala中做同样的事情，但我不明白它们是如何计算残差的。对不起，没有，但如果你找到方法，请告诉我：）该函数是用Fortran语言编写的：我认为自从您引用
trunk
以来，
lm.R
行号已经过时了。。。请将链接更新到正确的当前行，然后通过添加commit ref使链接永久化，好吗？