个人多用ggplot2，觉得「图层叠加」的概念用作此图比较合适，很快就能模仿出来
这个页面上的例子里有并在一起的条状图的画法
http://www.cookbook-r.com/Graphs/Bar_and_line_graphs_(ggplot2)/
这个页面上有误差条的画法
http://docs.ggplot2.org/0.9.3.1/geom_errorbar.html
样子：
<image ref="http://imgur.com/czeRFzJ">
</image>
甚至可以做出四张小图片的效果：
code
<br />
> data1 <- data.frame(x=c(0,0,5,5,10,10,20,20), type1=c(0,1,0,1,0,1,0,1), val=c(0,3,15,25,27,27,20,34))<br />
> head(data1)<br />
x type1 val<br />
1 0 0 0<br />
2 0 1 3<br />
3 5 0 15<br />
4 5 1 25<br />
5 10 0 27<br />
6 10 1 27<br />
> data1$x = as.factor(data1$x)<br />
> data1$type1 = as.factor(data1$type1)<br />
> data2<-data.frame(x=c(0,5,10,20), val=c(0,60,70,90))<br />
> data2$x = as.factor(data2$x)<br />
> ggplot() + geom_bar(data=data2, aes(x=x, y=val), stat="identity", width=0.6, fill="white", colour="black") + geom_bar(data=data1, aes(x=x, y=val, fill=type1), position="dodge", stat="identity", colour="black") + geom_errorbar(data=data1, aes(x=x, ymin=val-2, ymax=val+2, fill=type1, width=0.7), position="dodge") + scale_fill_discrete(guide=guide_legend(title=" "), breaks=c(0, 1), labels=c("NT", "eEF2K siRNA")) + ggtitle("Sample graph") + xlab("MK2206 umol/L") + ylab("Percent of cells with\nactivated caspases(%)")</p>
<p># 以下为facet用法（只是将数据复制了三遍）<br />
> data11=data1<br />
> data11$val=data11$val+3<br />
> data12=data1<br />
> data12$val=data12$val/2<br />
> data13=data1<br />
> data13$val=data13$val+5<br />
> data1=rbind(data1,data11,data12,data13)<br />
> data1$typ=c(1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,4,4,4,4,4,4,4,4)</p>
<p>> data21=data2<br />
> data2=rbind(data2,data21,data21,data21)<br />
> data2$typ=c(1,1,1,1,2,2,2,2,3,3,3,3,4,4,4,4)</p>
<p>> data1$typ=as.factor(paste("Type",data1$typ))<br />
> data2$typ=as.factor(paste("Type",data2$typ))</p>
<p>> ggplot() + geom_bar(data=data2, aes(x=x, y=val), stat="identity", width=0.6, fill="white", colour="black") + geom_bar(data=data1, aes(x=x, y=val, fill=type1), position="dodge", stat="identity", colour="black") + geom_errorbar(data=data1, aes(x=x, ymin=val-2, ymax=val+2, fill=type1, width=0.7), position="dodge") + scale_fill_discrete(guide=guide_legend(title=" "), breaks=c(0, 1), labels=c("NT", "eEF2K siRNA")) + ggtitle("Sample graph") + xlab("MK2206 umol/L") + ylab("Percent of cells with\nactivated caspases(%)") + facet_wrap(~typ, ncol=2)</p>
<p>
</p>
