以下方案基于MATLAB 2008

方案1：传统法
clear all;clc
d1 = [152 104 111 5 133 85 118 2]';
d2 = [132 90 123 6 123 97 135 4]';
data = [d1,d2 ];
bar(data,'group'); grid on
text(0.5,156,'[0.05,0.1]','FontSize',10);
text(1.5,108,'[0.01,0.05]','FontSize',10);
text(2.5,131,'[0.005,0.01]','FontSize',10);
text(3.5,14,'[0,0.005]','FontSize',10);
text(4.5,141,'[0,-0.005]','FontSize',10);
text(5.5,105,'[-0.005,-0.01]','FontSize',10);
text(6.5,142,'[-0.01,-0.05]','FontSize',10);
text(7.5,12,'[-0.05,-0.1]','FontSize',10);
legend('上海','深圳');
axis([0 9 0 180]);

可以参考Epicalc包中pyramid的代码编写相应的函数
或者参考help.search("histbackback")

原来上面的是matlab，我还想怎么这个格式那么奇怪类，原来是我自己搞错了，以为是R语言