多項式近似（勾配と切片）の誤差を計算する方法

Question

こんにちは、scipy.polyfit関数で計算される勾配と切片の誤差を計算したいと思います。私はydataについて（+/-）不確実性を持っているので、それを勾配と迎撃に不確実性を計算するためにどのように含めることができますか？私のコードは、あなたがscipy.optimize.curve_fit代わりのpolyfitを使用することができ、

from scipy import polyfit 
import pylab as plt 
from numpy import * 

data = loadtxt("data.txt") 
xdata,ydata = data[:,0],data[:,1] 


x_d,y_d = log10(xdata),log10(ydata) 
polycoef = polyfit(x_d, y_d, 1) 
yfit = 10**(polycoef[0]*x_d+polycoef[1]) 


plt.subplot(111) 
plt.loglog(xdata,ydata,'.k',xdata,yfit,'-r') 
plt.show() 

どうもありがとう

Answer 1

です。 ydataのエラーにはパラメータsigmaがあります。あなたが順序yerror内のすべてのy値のためのあなたのエラーを持っている場合（yerrorあなたy_dシーケンスと同じ長さを有するように）あなたが行うことができます。代替のために

polycoef, _ = scipy.optimize.curve_fit(lambda x, a, b: a*x+b, x_d, y_d, sigma=yerror) 

はべき乗則にフィッティング段落を見ますエラーのデータはScipy Cookbookにあります。