R 绘制连续数据的平均值 + SD

Question

我有一些连续的数据：

• mean_force - 2700 个值反映了我的相关变量在 2700 毫秒内的变化
• sd_force - 另一个 2700 个值，因变量随时间变化的标准差

我在图上绘制 mean_force ，如下所示：

    x=c(1:2700)
    plot(x, mean_force, ty="l", col="blue", ylim=c(-15, 15),
    ylab="force (mN)", 
    xlab='time (ms)',
    lty=1,lwd=3)

问题是我如何也绘制 SD，如下图所示（平均值周围的阴影区域）？ 在 R 中有什么简单的方法可以做到吗？

Answer 1

如果您正在寻找基本图形解决方案，也许这对您有用：

#create data
x<-1:100
mean_force<-0.5*x+rnorm(100)
#assume constant standard deviation across the 
sd<-5
#determine error band
psd<-mean_force+sd
nsd<-mean_force-sd

plot(x, mean_force, ty="l", col="blue", 
     ylab="force (mN)", 
     xlab='time (ms)',
     lty=1,lwd=3)
#draw boundary and fill
lines(x, psd)
lines(x, nsd)
polygon(x=c(x, rev(x)), y=c(psd, rev(nsd)), col="lightblue", density = 40, angle=90)
#redraw line on top
lines(x, mean_force, col="blue",lwd=3)