带有误差条的散点图

Question

我目前拥有大量的数据集，并创建了一个散点图，使用代码查看区域边界与路径复杂度（曲折度）之间的关系：

plot(Tortuosity~Perimeter, data=TR, xlab="Territory perimeter (m)", ylab="Path tortuosity",pch=19,cex.lab=0.85,cex.axis=0.85,cex=0.85)

TR数据框由以下信息组成，其中Prm。代表周长，T代表Tortuosity

Prm.    T.
198.58  0.45
198.58  0.63
198.58  0.24
198.58  0.32
198.58  0.6
198.58  0.49
134.62  0.3
134.62  0.39
134.62  0.12
134.62  0.09
134.62  0.25
134.62  0.45
134.62  0.44
134.62  0.84
128.07  1.99
144.32  0.72
144.32  1.01
144.32  1.09
144.32  1.15
189.96  0.55
189.96  0.76
189.96  0.78
113.74  0.67
113.74  0.73
113.74  0.69
113.74  0.61
134.15  0.39
134.15  0.26
134.15  1.37
156.11  1.1
156.11  0.67
156.11  0.98
88.41   1.09
88.41   0.87
88.41   0.47
88.41   0.8
88.41   0.52
86.71   1.47
94.4    0.92
94.4    1.23
94.4    0.96
94.4    0.56
94.4    0.82
87.53   1.22
10.85   
93.22   1.67
90.13   1.26
90.13   0.67
198.58  1.46
198.58  0.99
198.58  0.68
198.58  0.53
134.62  0.12
134.62  0.32
134.62  0.73
134.62  1.18
144.32  1.03
189.96  0.86
189.96  0.43
189.96  0.42
113.74  0.9
113.74  0.7
134.15  0.97
156.11  0.85
156.11  0.45
156.11  0.84
156.11  0.8
88.41   0.81
88.41   1.74
88.41   0.98
88.41   0.75
86.71   1.1
94.4    0.69
94.4    0.63
94.4    0.67
94.4    0.59
94.4    0.82
87.53   0.64
87.53   0.7
10.85   2.73
10.85   
10.85   2.73
10.85   0.85
10.85   1.01
93.22   1.67
90.13   0.89
90.13   0.74
90.13   0.59
90.13   0.58
90.13   1.14
90.13   0.93
154.33  0.97
154.33  0.53
154.33  0.62

但是我想为每个点添加标准错误'whiskers'（参见下面的图表）。

有没有办法自动执行此操作，类似于交互式效果图的lineplot.CI()功能？对于信息的长度表示道歉，但我不确定将会使用什么......

Answer 1

以下是我将如何使用ggplot2进行操作，但我不得不说，包含错误栏的情节不太可读。无论如何，制作棒材的手动操作并不高。

library(ggplot2)
library(dplyr)

TR %>% 
  group_by(Prm.) %>%
  mutate(sd = sd(T.)) %>%
  ggplot(., aes(x = Prm., y = T.)) +
  geom_point() +
  geom_errorbar(aes(ymin = T. - sd, ymax = T. + sd), width = 2)