Wie berechnet man Fehler für polynomische Anpassung (in Steigung und Schnittpunkt)

Question

Hallo Ich möchte Fehler in Steigung und Achsenabschnitt berechnen, die von scipy.polyfit Funktion berechnet werden. Ich habe (+/-) Ungewissheit für Ydata, also wie kann ich es für die Berechnung der Unsicherheit in Steigung und Schnittpunkt einbeziehen? Mein Code ist,

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

Vielen Dank

Answer 1

Sie scipy.optimize.curve_fit statt polyfit nutzen könnten. Es hat einen Parameter sigma für Fehler von YData. Wenn Sie Ihren Fehler für jeden y-Wert in einer Folge haben yerror (so dass yerror hat die gleiche Länge wie Ihre y_d Sequenz) können Sie tun:

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

Für eine Alternative siehe Absatz Montage ein Potenzgesetz zu Daten mit Fehlern in der Scipy Cookbook.