!!グラフにエラーバーを入れる

! データ

 ID <-  c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82,83,84,85,86,87,88,89,90,91,92,93,94,95,96,97,98,99,100)
 status <- c("A","A","A","A","A","A","A","A","A","A","A","A","A","A","A","A","A","A","A","A","A","A","A","A","A","B","B","B","B","B","B","B","B","B","B","B","B","B","B","B","B","B","B","B","B","B","B","B","B","B","C","C","C","C","C","C","C","C","C","C","C","C","C","C","C","C","C","C","C","C","C","C","C","C","C","D","D","D","D","D","D","D","D","D","D","D","D","D","D","D","D","D","D","D","D","D","D","D","D","D")
 point <- c(8,5,6,3,7,8,3,2,5,3,2,2,3,5,8,1,7,10,10,2,3,6,8,9,6,4,5,9,2,3,9,5,2,3,7,1,9,6,6,8,5,2,3,3,4,1,4,6,1,1,9,6,3,2,4,9,1,7,7,9,2,10,2,7,8,10,10,6,3,7,3,7,4,9,6,3,4,5,7,4,2,1,1,3,6,10,5,2,2,7,4,2,2,1,7,8,5,7,3,3)

 test <- tibble(ID = as.integer(ID), status, point = as.integer(point))

!! 自分のやってみた方法

 summary.test <-
  test |> group_by(status) |>
  summarise( mean = mean(point), sd = sd(point), lower = mean - sd, upper = mean + sd) 
 ggplot(test) +  
  ggbeeswarm::geom_beeswarm(aes(x = status, y = point, color = status)) +
  geom_errorbar(data = summary.test, aes(x = status, ymin = lower, ymax = upper, width = 0.2)) +
  geom_point(aes(x = status, y = mean), size = 2)

!! 教わった方法
stat_summary()を使うと事前の計算が要らないようだ。う&#12316;ん、便利すぎる。


 test %>% 
  ggplot2::ggplot(aes(x=status,y=point))+
  ggbeeswarm::geom_beeswarm()+
  ggplot2::stat_summary(fun.data = mean_cl_normal,  
                        geom = "errorbar")+
  ggplot2::stat_summary(fun = mean,  geom = "point",color="red")