绘制R中的错误栏

时间:2017-02-04 12:25:43

标签: r statistics errorbar

我有以下数据集:

            conf    variable        value
1           C1      T1      0.6578166
2           C2      T1      0.6729989
3           C3      T1      0.6629129
4           C4      T1      0.6602419
5           C5      T1      0.6951054
6           C1      T2      0.6764889
7           C2      T2      0.7008580
8           C3      T2      0.6921985
9           C4      T2      0.6868215
10          C5      T2      0.7104877
11          C1      T3      0.6834098
12          C2      T3      0.7120974
13          C3      T3      0.7045760
14          C4      T3      0.6986261
15          C5      T3      0.7180474
16          C1      T4      0.6875082
17          C2      T4      0.7185950
18          C3      T4      0.7122448
19          C4      T4      0.7051376
20          C5      T4      0.7229975
21          C1      T5      0.6902167
22          C2      T5      0.7228775
23          C3      T5      0.7167748
24          C4      T5      0.7093825
25          C5      T5      0.7261777
....

我想将这些数据绘制为带有误差线的线图。为此,我使用了website中定义的summarySE函数。我希望绘图是一个线图,其中x轴代表variable,y轴代表value,组代表conf。考虑到从T1到当前变量点的值,每个变量都应显示错误条。例如,T6中错误条的考虑值应为T1T6的值。

基于此,我将summarySE函数称为:

m2 <- summarySE(m, measurevar = "value", groupvars = c("variable","conf") , na.rm=TRUE)

我得到的结果是:

     variable    conf N  value     sd se ci
1    T1          C1   1  0.6578166 NA NA NA
2    T1          C2   1  0.6729989 NA NA NA
3    T1          C3   1  0.6951054 NA NA NA
4    T1          C4   1  0.6629129 NA NA NA
5    T1          C5   1  0.6602419 NA NA NA
6    T2          C1   1  0.6764889 NA NA NA
......

我也尝试了以下内容:

m2 <- summarySE(m, measurevar = "value", groupvars = c("conf") , na.rm=TRUE)

我得到了以下结果:

             conf  N     value          sd          se          ci
1              C1 60 0.7038542 0.009609437 0.001240573 0.002482381
2              C2 60 0.7395447 0.012869318 0.001661422 0.003324497
3              C3 60 0.7418830 0.010397609 0.001342326 0.002685987
4              C4 60 0.7365407 0.014438260 0.001863971 0.003729798
5              C5 60 0.7301954 0.014389284 0.001857649 0.003717146

但这并不能用每个变量绘制错误栏。

绘制情节:

coverage_plot <- ggplot(m2, aes(x=m$variable, y=m$value, group=conf, color=conf)) +
      geom_errorbar(aes(ymin=value-ci, ymax=value+ci), width=0.1, position=pd) +
      geom_line(position=pd) + geom_point(position=pd) +
      scale_x_discrete(labels = c(1:60))

当我运行代码时,我也会收到以下错误:

Error: Aesthetics must be either length 1 or the same as the data (5): ymin, ymax, x, y, group, colour

有人可以帮我解决这个问题吗?

1 个答案:

答案 0 :(得分:-1)

尝试删除此{/ 1}}中的$

aes

此外,您的数据中没有名为变量的列。确保你写的是正确的名字