我在r中制作了一个有两个系列的情节,基本上是时间序列。每个系列都是多个其他系列的集合,因此时间序列中的每个点都有一个均值和一个方差。我想在一个图上汇总系列,我想在每个点显示误差条。误差条代表方差。
以下是我想要添加错误栏的图表示例:
## random generation of series 1 components:
series1a<-rnorm(50,50,25)
series1b<-rnorm(50,49,26)
series1c<-rnorm(50,58,24)
## aggregate components into series 1:
series1.df<-t(data.frame(series1a,series1b,series1c))
series1<-colMeans(series1.df)
## random generation of series 2 components:
series2a<-rnorm(50,45,20)
series2b<-rnorm(50,44,21)
series2c<-rnorm(50,46,19)
## aggregate components into series 1:
series2.df<-t(data.frame(series2a,series2b,series2c))
series2<-colMeans(series2.df)
## plot:
plot(series1,type="o",col="blue",pch=0)
lines(series2,type="o",lty="dashed",col="red",pch=1)
这是我想要添加错误栏的情节图片: plot of two aggregate time series
有没有简单的方法来添加错误栏?
如果可能,我想避免使用ggplot2或其他软件包。
感谢。