MeteoInfoLab脚本示例：误差棒 - errorbar

MeteoInfo

更新了MeteoInfoLab的功能，可以使用errorbar()函数绘制误差棒（error bar）图。

固定x, y误差：

1. x = arange(0.1, 4, 0.5)
   
2. y = exp(-x)
   
3. errorbar(x, y, fmt='b', xerr=0.2, yerr=0.4)
   
4. title('Fixed error values example')

x = arange(0.1, 4, 0.5)<br /> y = exp(-x)<br /> errorbar(x, y, fmt='b', xerr=0.2, yerr=0.4)<br /> title('Fixed error values example')

变化y误差：

1. x = arange(0.1, 4, 0.5)
   
2. y = exp(-x)
   
3. # example error bar values that vary with x-position
   
4. error = 0.1 + 0.2 * x
   
5. errorbar(x, y, yerr=error, fmt='b-o')
   
6. title('Variable error bar values example')

x = arange(0.1, 4, 0.5)<br /> y = exp(-x)<br /> # example error bar values that vary with x-position<br /> error = 0.1 + 0.2 * x<br /> errorbar(x, y, yerr=error, fmt='b-o')<br /> title('Variable error bar values example')

为直方图（bar）加误差棒：

1. menMeans = [20, 35, 30, 35, 27]std_men = (2, 3, 4, 1, 2)
   
2. n = len(menMeans)
   
3. ind = arange(n)
   
4. width = 0.35
   
5. gap = 0.06
   
6. bar(ind, menMeans, width, yerr=std_men, color='r', ecolor='b', label='Men')
   
7. for j in range(n):
   
8.     text(ind[j] + width / 4, menMeans[j] + 2, str(menMeans[j]))
   
9. 
   
10. womenMeans = [25, 32, 34, 20, 25]
    
11. std_women = (3, 5, 2, 3, 3)
    
12. bar(ind + width + gap, womenMeans, width, yerr=std_women, color='y', ecolor='g', label='Women')
    
13. for j in range(n):
    
14.     text(ind[j] + + width + gap + width / 4, womenMeans[j] + 2, str(womenMeans[j]))
    
15. 
    
16. xlim(-0.2, 5)
    
17. ylim(0, 40)
    
18. ylabel('Scores')
    
19. xticks(ind + width + gap / 2, ['G1','G2','G3','G4','G5'])
    
20. legend()
    
21. title('Scores by group and gender')

menMeans = [20, 35, 30, 35, 27]std_men = (2, 3, 4, 1, 2)<br /> n = len(menMeans)<br /> ind = arange(n)<br /> width = 0.35<br /> gap = 0.06<br /> bar(ind, menMeans, width, yerr=std_men, color='r', ecolor='b', label='Men')<br /> for j in range(n):<br />     text(ind[j] + width / 4, menMeans[j] + 2, str(menMeans[j]))<br /> <br /> womenMeans = [25, 32, 34, 20, 25]<br /> std_women = (3, 5, 2, 3, 3)<br /> bar(ind + width + gap, womenMeans, width, yerr=std_women, color='y', ecolor='g', label='Women')<br /> for j in range(n):<br />     text(ind[j] + + width + gap + width / 4, womenMeans[j] + 2, str(womenMeans[j]))<br /> <br /> xlim(-0.2, 5)<br /> ylim(0, 40)<br /> ylabel('Scores')<br /> xticks(ind + width + gap / 2, ['G1','G2','G3','G4','G5'])<br /> legend()<br /> title('Scores by group and gender')

双Y轴Bar图+误差棒：

1. menMeans = [20, 35, 30, 35, 27]
   
2. std_men = (2, 3, 4, 1, 2)
   
3. n = len(menMeans)
   
4. ind = arange(n)
   
5. width = 0.35
   
6. gap = 0.06
   
7. ax1 = axes(position=[0.11,0.1,0.75,0.8])
   
8. yaxis(ax1, color='r')
   
9. bar1 = bar(ind, menMeans, width, yerr=std_men, color='r', ecolor='b', label='Men')
   
10. for j in range(n):
    
11.     text(ind[j] + width / 4, menMeans[j] + 2, str(menMeans[j]))
    
12. xlim(-0.2, 5)
    
13. ylim(0, 40)
    
14. ylabel('Men Scores', color='r')
    
15. xticks(ind + width + gap / 2, ['G1','G2','G3','G4','G5'])
    
16. title('Scores by group and gender')
    
17. 
    
18. womenMeans = array([25, 32, 34, 20, 25])
    
19. womenMeans = womenMeans * 2
    
20. std_women = (3, 5, 2, 3, 3)
    
21. ax2 = twinx(ax1)
    
22. yaxis(ax2, color='b')
    
23. bar2 = bar(ind + width + gap, womenMeans, width, yerr=std_women, color='b', ecolor='g', label='Women')
    
24. for j in range(n):
    
25.     text(ind[j] + + width + gap + width / 4, womenMeans[j] + 2, str(womenMeans[j]))
    
26. xlim(-0.2, 5)
    
27. ylim(0, 80)
    
28. ylabel('Women Scores', color='b')
    
29. 
    
30. legend(bar1+bar2)

menMeans = [20, 35, 30, 35, 27]<br /> std_men = (2, 3, 4, 1, 2)<br /> n = len(menMeans)<br /> ind = arange(n)<br /> width = 0.35<br /> gap = 0.06<br /> ax1 = axes(position=[0.11,0.1,0.75,0.8])<br /> yaxis(ax1, color='r')<br /> bar1 = bar(ind, menMeans, width, yerr=std_men, color='r', ecolor='b', label='Men')<br /> for j in range(n):<br />     text(ind[j] + width / 4, menMeans[j] + 2, str(menMeans[j]))<br /> xlim(-0.2, 5)<br /> ylim(0, 40)<br /> ylabel('Men Scores', color='r')<br /> xticks(ind + width + gap / 2, ['G1','G2','G3','G4','G5'])<br /> title('Scores by group and gender')<br /> <br /> womenMeans = array([25, 32, 34, 20, 25])<br /> womenMeans = womenMeans * 2<br /> std_women = (3, 5, 2, 3, 3)<br /> ax2 = twinx(ax1)<br /> yaxis(ax2, color='b')<br /> bar2 = bar(ind + width + gap, womenMeans, width, yerr=std_women, color='b', ecolor='g', label='Women')<br /> for j in range(n):<br />     text(ind[j] + + width + gap + width / 4, womenMeans[j] + 2, str(womenMeans[j]))<br /> xlim(-0.2, 5)<br /> ylim(0, 80)<br /> ylabel('Women Scores', color='b')<br /> <br /> legend(bar1+bar2)

rceclx

王老师，第三张图可以改成双Y轴柱状图加误差棒么？谢谢！

MeteoInfo

rceclx 发表于 2015-12-25 13:38
王老师，第三张图可以改成双Y轴柱状图加误差棒么？谢谢！

在1楼加入了这种图的示例。需要重新下载MeteoInfo_Java_1.2.9R2_Files.zip。

rceclx

MeteoInfo 发表于 2015-12-25 15:26
在1楼加入了这种图的示例。需要重新下载MeteoInfo_Java_1.2.9R2_Files.zip。

王老师，棒棒哒！

buzaiyouyu

{:5_213:}{:5_213:}{:5_213:}

0727

{:5_213:}学习了

Linda.W

王老师，您好，我想把error-bar图与plot点线图画在一起，但是legend显示不出来，总是报错，请王老师指点迷津。 r = gca.legend(*args, **kwargs)
AttributeError: 'NoneType' object has no attribute 'legend'
脚本如下：menMeans = [50, 62, 70, 63, 50]
std_men = (41, 26, 39, 25, 19)
n = len(menMeans)
ind = arange(n)
width = 0.5
gap = 0.06
ax1 = axes(position=[0.11,0.1,0.75,0.8])
yaxis(ax1, color='k')
bar1 = bar(ind, menMeans,width, yerr=std_men, color=(160,234,255), linewidth=0,edgecolor=(160,234,255), label='Men')

xlim(-0.5, 5)
ylim(0, 130)

ylabel(u'人数',fontname=u'宋体',fontsize=18,bold='True', color='k')
xticks(ind+width/2, ['(0,0.5]','(0.5,1]','(1,2.5]','(2.5,5]','(5,15]'])

ax2 = twinx(ax1)
yaxis(ax2, color='k')
t=ind+width/2
s2 = [39,21,20,4,4]
line1=plot(t, s2,'.-b')
ylabel(u'天数',fontname=u'宋体',fontsize=18,bold='True', color='k')
xlim(-0.5, 5)
ylim(0, 45)
show()

MeteoInfo

Linda.W 发表于 2019-3-5 21:43
王老师，您好，我想把error-bar图与plot点线图画在一起，但是legend显示不出来，总是报错，请王老师指点迷 ...

legend([bar1[0], line1)

Linda.W

好的，谢谢王老师。

Linda.W

王老师，我在melab官网上看到堆积柱状图的示例脚本，两个堆叠在一起，使用bottom语句，如何实现多个（5~6）柱状图堆积?谢谢老师。
menMeans = [20, 35, 30, 35, 27]
n = len(menMeans)
ind = arange(n)
width = 0.35
bar(ind, menMeans, width, color='r', ecolor='b', label='Men')

womenMeans = [25, 32, 34, 20, 25]
bar(ind, womenMeans, width, bottom=menMeans, color='y', ecolor='g', label='Women')

xlim(-0.2, 5)
ylim(0, 80)
ylabel('Scores')
xticks(ind + width / 2, ['G1','G2','G3','G4','G5'])
legend()
title('Scores by group and gender')