Jak obliczyć błąd dla wielomianu (nachylenie i przecięcie)

Question

Witam Chcę obliczyć błędy w nachyleniu i przecięciu, które są obliczane przez funkcję scipy.polyfit. Mam (+/-) niepewność dla ydata, więc jak mogę włączyć ją do obliczenia niepewności na nachylenie i przechwycenie? Mój kod jest

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

Thanks a lot

Answer 1

Można użyć scipy.optimize.curve_fit zamiast polyfit. Ma parametr sigma dla błędów ydata. Jeśli masz swój błąd dla każdej wartości y w sekwencji yerror (tak, że yerror ma taką samą długość jak swojej sekwencji y_d) można zrobić:

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

alternatywy patrz paragraf Zamontowanie potęgowego do dane z błędami w Scipy Cookbook.