现在的位置: 首页 > 综合 > 正文

带误差区间图

2018年10月21日 ⁄ 综合 ⁄ 共 1175字 ⁄ 字号 评论关闭
文章目录

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语句

LIMITS= BOTH | LOWER | UPPER ,
LIMITSTAT= CLM | STDDEV | STDERR

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)/

 

【上篇】
【下篇】

抱歉!评论已关闭.