擅长:python、mysql、java
<p>Numpy的<code>polyfit</code>至少在1.3版之后支持<code>full</code>参数。如果设置为<code>True</code>,<code>polyfit</code>将返回更多的值,包括残差的平方。它是卡方(不受自由度的规范化)。</p>
<p>一个简单的例子是</p>
<pre><code>p, residuals, _, _, _ = numpy.polyfit(x, y, 2, full=True)
chisq_dof = residuals / (len(x) - 3)
</code></pre>
<p>我自己并没有尝试使用权重,但是我假设<code>polyfit</code>在这里做的是正确的(因为numpy 1.7,<code>polyfit</code>接受一个参数<code>w</code>来提供适合的权重)。</p>