带误差区间图

google到很多资源，整理如下：

一.sas

 HILOCTJ  interpolation

1.1 proc gplot gchart等

a.support.sas官网的例子：带误差线的均值图

http://support.sas.com/kb/42/515.html

b.几篇文章

“A Generalized Procedure to Create SAS®/Graph Error Bar Plots”

“How To Plot A Line Chart With A Standard Error”

“How to Get Upper (or Lower) Error Bars with INTERPOL=HILOCTJ?”

1.2proc sgplot等sg语句

a.官网例子：带误差线的均值图

http://support.sas.com/kb/42/542.html

b.两个博客: bar charts

（sgplot）http://blogs.sas.com/content/iml/2011/10/07/creating-bar-charts-with-confidence-intervals/

(gchart and r)http://sas-and-r.blogspot.com/2011/11/example-915-bar-chart-with-error-bars.html

二、R

1.sas-and-r-blogspot博客中提到sciplot包

2.gplots包：plotCI函数

 # plot means and
 data(state)
 tmp   <- split(state.area, state.region)
 means <- sapply(tmp, mean)
 stdev <- sqrt(sapply(tmp, var))
 n     <- sapply(tmp,length)
 ciw   <- qt(0.975, n) * stdev / sqrt(n)

 # plain
 plotCI(x=means, uiw=ciw)

 # prettier
 plotCI(x=means, uiw=ciw, col="black", barcol="blue", lwd=1)

 # give mean values
 plotCI(x=means, uiw=ciw, col="black", barcol="blue",
        labels=round(means,-3), xaxt="n", xlim=c(0,5) )
 axis(side=1, at=1:4, labels=names(tmp), cex=0.7)

3.ggplot2包，

例：http://wiki.stdout.org/rcookbook/Graphs/Plotting%20means%20and%20error%20bars%20(ggplot2)/