我在尝试使用矩阵中两列的值将95%CI(下/上)添加到现有绘图时遇到问题。使用此信息添加错误栏的最佳方法是什么?
以下是我的数据示例:
option<-read.table(text=" distance p.move id option mean lower95%CI upper95%CI 1 close 0.05 1 10% 13.682 11.306 15.768 2 close 0.10 2 10% 10.886 9.336 12.270 3 close 0.15 3 10% 8.402 7.262 9.580 4 close 0.20 4 10% 7.240 6.132 8.350 5 close 0.25 5 10% 6.322 5.288 7.370 6 close 0.30 6 10% 5.850 4.920 6.714 7 close 0.35 7 10% 3.838 3.084 4.648 8 close 0.40 8 10% 3.600 2.936 4.200 9 close 0.45 9 10% 3.380 2.702 4.016 10 close 0.50 10 10% 3.152 2.462 3.720 11 close 0.55 11 10% 2.772 2.214 3.286 12 close 0.60 12 10% 3.072 2.458 3.596 13 close 0.65 13 10% 2.670 2.134 3.212 14 close 0.70 14 10% 2.194 1.724 2.634 15 close 0.75 15 10% 1.980 1.612 2.336 16 close 0.80 16 10% 2.028 1.594 2.466 17 close 0.85 17 10% 1.650 1.294 1.974 18 close 0.90 18 10% 1.916 1.564 2.254",header=T) option
这是我的情节:
plot(option$mean~option$p.move,xlim=c(0,1),type="o",ylim=c(0,20), xlab="Probability",ylab="% time",col=1,lwd=1.85)
提前多多感谢,
答案 0 :(得分:3)
您可以使用lines
添加额外的列,例如plot
,但可以使用现有的绘图。 (请参阅?lines
,?points
)。
此外,当您plot
使用数据框时,您可以通过将option$
提供给option
参数(请参阅dat
)来跳过所有?plot
:
# draw original plot
plot(mean ~ p.move, dat=option, xlim=c(0,1), type="o", ylim=c(0,20),
xlab="Probability",ylab="% time",col=1,lwd=1.85)
# draw extra lines (the '%' in the column names gets converted to '.' by R)
# note you can put your usual `plot` arguments into `lines` like lwd, type etc
# if you want
lines(upper95.CI ~ p.move, option)
lines(lower95.CI ~ p.move, option)