如何将 tukey 测试的结果放在箱线图上

问题描述

我正在尝试为箱线图上的显着差异添加标签。我想将 Tukey 检验的结果放在箱线图上（带字母）。我尝试了以下代码，但没有出现字母。有没有人有这方面的经验并且可以帮助我？ ????

My treatment = zone
My variable = pH 
My data = thesis_data_D1

TukeyHSD(res_anov_pH)
plot(TukeyHSD(res_anov_pH))
TUKEY <- TukeyHSD(res_anov_pH)
plot(TUKEY,las=1,col="brown")

generate_label_df <- function(TUKEY,variable){
Tukey.levels <- TUKEY[[variable]][,4]
Tukey.labels <- data.frame(multcompLetters(Tukey.levels)['Letters'])
Tukey.labels$treatment=rownames(Tukey.labels)
Tukey.labels=Tukey.labels[order(Tukey.labels$treatment),]
return(Tukey.labels)}


LABELS <- generate_label_df(TUKEY,"zone")
my_colors <- c(rgb(143,199,74,maxColorValue = 255),rgb(242,104,34,maxColorValue = 
255),rgb(111,145,202,maxColorValue = 255))

boxplot_pH_labels<-boxplot(thesis_data_D1$pH ~ thesis_data_D1$zone,ylim=c(min(thesis_data_D1$pH),1.1*max(thesis_data_D1$pH)),col=my_colors[as.numeric(LABELS[,1])],ylab="pH",main="")

over=0.1*max( boxplot_pH_labels$stats[nrow(boxplot_pH_labels$stats),] )

text( c(1:nlevels(thesis_data_D1$zone)),boxplot_pH_labels$stats[nrow(boxplot_pH_labels$stats),]+over,LABELS[,1],1])]

代码运行但没有出现字母（见图）。

如果有人能在我完全发疯之前提供帮助，那就太好了！

谢谢

摩根

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）