बहुपद फिटिंग (ढलान और अवरोध में) के लिए त्रुटि की गणना कैसे करें

Question

हाय मैं ढलान और अवरोध में त्रुटियों की गणना करना चाहता हूं जिन्हें 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() 

धन्यवाद एक बहुत

Answer 1

आप scipy.optimize.curve_fit बजाय polyfit इस्तेमाल कर सकते हैं। Ydata की त्रुटियों के लिए इसमें पैरामीटर sigma है। आप एक दृश्य yerror में हर y मूल्य के लिए अपने त्रुटि है, तो (ताकि yerror अपने y_d अनुक्रम के रूप में एक ही लंबाई है) आप कर सकते हैं:

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

एक विकल्प के लिए एक शक्ति जी के लिए फिटिंग पैरा देखना Scipy Cookbook में त्रुटियों वाले डेटा।