用字母表示 Tukeys 重要性的条形图

如何解决用字母表示 Tukeys 重要性的条形图

我正在尝试使用 ggplot 创建条形图，其中包含 SE 条形和字母以表示 Tukeys 的重要性..

我发现其他代码可以与像这样的箱线图做类似的事情Adding Tukey's significance letters to boxplot

但我无法让该代码为我工作。

我有这个数据：

structure(list(Zone = structure(c(1L,1L,2L,2L),.Label = c("B","D"),class = "factor"),Discharge.Year = structure(c(3L,3L,1L),.Label = c("high","low","medium"),abundance.cpue = c(28.8010906921919,3.72369131138694,11.0539420104989,5.51098996588217,7.37909842836217,10.292812004944,20.3498830724951,15.4479222881813,18.1275906810118,37.2908796546141,30.8621154883972,18.6077701132229,23.7722611980572,20.2087471724977,12.6153743441953,7.0939888237035,13.9233675121425,11.3464651915812,18.4710221930444,20.6208992529895,16.8645439827307,8.97661271793666,10.4020507285483,9.46001714231294,20.8892958275981,10.5984289740361,22.7820776435066,19.9822056271058,33.5941716135996,6.61090124123044,7.20453318942256,9.98381003777658,23.4484619535888,23.7736813701821,34.1442152019594,8.17758295978942,5.42731688171515,9.9128014315335,26.358657632268,26.7710424031483,36.7605417974189,37.1301043353121,5.33924286091559,2.12212178449361,3.40697858236984,1.44220927852318,10.4524668389809,7.31992972867461,7.45517160803601,10.1875612630373,10.6648163852697,14.8268898025492,14.8820996288424,17.8175245633422,16.0924323305093,21.8827741434457,15.4789165852307,10.5586710152749,6.19027967366115,10.9006279742923,24.9283905517041,26.9832703723691,24.7214374146852,13.5266783944604,14.6027110250468,30.4566487926422,33.9020576473356,32.2058388315373,24.8087362305959,21.0498523173697,48.2325957906098,41.1045152005756,59.4034683570813,59.2828874669134,10.5639503507825,4.19892727733295,11.8250214278912,21.1413923367512,12.3954398273071,19.8291532933792,19.851977488243,35.0064874671321,16.2534522713565,3.90547334336922,23.5755893201885,11.3194819212088,25.5932231107622,34.9883072495053,32.4074583962025,19.8731786292499,25.8094981111711,21.0910052094232,16.743119266055,23.8197145356086,37.689535713664,15.4783155322821,7.0649745544859,16.1899622234215,25.2465811622326,21.9852327920326,36.6244838312606,15.734400515836,15.8320208817472,1.85884751454098,17.5226725283982,19.1281405528572,29.2961221185722,25.8063404944254,22.0635932994271,9.67123155920557,4.18376969915016,23.0617017449181,41.8451497571506,14.9890584321237,34.9634290569634,11.2632415093653,4.42525634106854,3.13887022698988,8.32709526545155,19.4220242790129,14.1731713966207,15.6137569706782,15.6392832362303,22.8008634646519,28.7399937039036,66.6527051735387,39.8846465191581,28.0910046297611,20.134799033295,11.9103575138253,26.1348273348874,13.0564211479205,26.5185584937339,15.1106314085267,31.0783784641757,22.1421542173264,29.5453182663112,12.4572764885771,21.9416593291107,30.5268564815116,25.7402208383611,31.6732197466142,25.0515135161654,24.0856350880636,34.8833723091683,35.991321744155,35.730594284291,28.7685271495309)),row.names = c(NA,-148L),class = "data.frame")

并使用过此代码：

ggplot(data.cpue,aes(x = factor(Zone,levels = c("B","D")),y = abundance.cpue,fill = factor(Discharge.Year,levels = c("high","medium","low")))) +  
  stat_summary(fun = mean,geom = "bar",position = position_dodge(0.9),color = "black") +  
  stat_summary(fun.data = mean_se,geom = "errorbar",width = 0.2) + 
  labs(x = "Estuary",y = "CPUE (+/-1SE)",legend = NULL) + 
  scale_y_continuous(expand = c(0,0)) + 
  scale_x_discrete(labels = c("high"="High","medium"="Medium","low"="Low")) + 
  scale_fill_manual(values = c("blue4","royalblue1","lightskyblue1")) + 
  scale_color_hue(l=40,c = 100)+
  coord_cartesian(ylim = c(0,35)) +
  theme(panel.border = element_blank(),panel.grid.major = element_blank(),panel.background = element_blank(),panel.grid.minor = element_blank(),axis.line = element_line(colour = "black"),axis.text.x = element_text(angle = 0))+
  theme(legend.title = element_blank(),legend.position = "bottom",legend.direction = "horizontal")   

这使得：

这是我希望完成的一个例子..

在此先感谢您的帮助！！