如何在matplotlib图的底部显示残差

Question

我想重现这个情节.错误显示在图的底部.你能分享一下它的成功吗？

我在这里找到了一个关于stackoverflow的例子,但是它在R. 如何创建一个图表来显示R中的预测模型,数据和残差

Answer 1

您只能使用add_axes在Matplotlib中创建此类绘图.这是一个例子.

from scipy.optimize import curve_fit
#Data
x = arange(1,10,0.2)
ynoise = x*numpy.random.rand(len(x)) 
#Noise; noise is scaled by x, in order to it be noticable on a x-squared function
ydata = x**2 + ynoise #Noisy data

#Model
Fofx = lambda x,a,b,c: a*x**2+b*x+c
#Best fit parameters
p, cov = curve_fit(Fofx,x,ydata)

#PLOT
fig1 = figure(1)
#Plot Data-model
frame1=fig1.add_axes((.1,.3,.8,.6))
#xstart, ystart, xend, yend [units are fraction of the image frame, from bottom left corner]
plot(x,ydata,'.b') #Noisy data
plot(x,Fofx(x,*p),'-r') #Best fit model
frame1.set_xticklabels([]) #Remove x-tic labels for the first frame
grid()

#Residual plot
difference = Fofx(x,*p) - ydata
frame2=fig1.add_axes((.1,.1,.8,.2))        
plot(x,difference,'or')
grid()

Answer 2

我认为您正在寻找像这样的pylab_examples 示例代码的错误栏：errorbar_demo.py

您可以添加额外的子图并用误差线绘制点。

编辑：图之间没有边界：

from pylab import *
subplots_adjust(hspace=0.,wspace=0.)
subplot(211)
imshow(rand(100,100), cmap=cm.BuPu_r)
subplot(212)
imshow(rand(100,100), cmap=cm.BuPu_r)
show()