Cách tính sai số cho phụ kiện đa thức (theo độ dốc và chặn)

Question

Xin chào Tôi muốn tính toán sai số về độ dốc và đánh chặn được tính bằng hàm scipy.polyfit. Tôi có (+/-) sự không chắc chắn đối với ydata vì vậy làm thế nào tôi có thể bao gồm nó để tính toán độ không chắc chắn vào độ dốc và chặn? Mã của tôi là,

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

Bạn có thể sử dụng scipy.optimize.curve_fit thay vì polyfit. Nó có thông số sigma vì lỗi của ydata. Nếu bạn có lỗi của bạn cho mỗi giá trị y trong một chuỗi yerror (để yerror có chiều dài tương tự như chuỗi y_d bạn), bạn có thể làm:

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

Đối với một sự lựa chọn xem đoạn Lắp một cường quốc-pháp luật để dữ liệu có lỗi trong số Scipy Cookbook.