5

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

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()

どうもありがとう

4

1 に答える 1

4

scipy.optimize.curve_fitの代わりに使用できますpolyfit。sigmaydata のエラー用のパラメータがあります。シーケンス内のすべての y 値にエラーがある場合yerror(したがって、シーケンスyerrorと同じ長さy_d)、次のことができます。

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

別の方法については、Scipy クックブックのパラグラフFitting a power-law to data with errors を参照してください。

于 2012-10-03T10:56:19.827 に答える