Polinom uydurma hatası (eğim ve kesişme noktası) nasıl hesaplanır

Question

Merhaba Scipy.polyfit işlevi tarafından hesaplanan eğim ve kesişme noktalarındaki hataları hesaplamak istiyorum. Ydata için (+/-) belirsizliğim var, bu yüzden belirsizliği eğime ve kesişme noktasına hesaplamak için nasıl ekleyebilirim? Kodum,

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

Teşekkür Sen scipy.optimize.curve_fit yerine polyfit kullanabilirsiniz

Answer 1

çok. Ydata hataları için sigma parametresine sahiptir. Eğer bir dizinin yerror her y değerinin için hata varsa (yerror sizin y_d dizisi ile aynı uzunluğa sahiptir, böylece) Yapabileceğiniz: bir alternatifi için

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

güç kanunu göre donanım parçası paragraf bkz Scipy Cookbook numaralı telefondan numaralı hatalarla ilgili veriler.