Добавление значимых писем Тьюки в коробку

Question

Я пытаюсь создать boxplot с ggplot подсчета (MedMean) на yaxis и различных независимых образцов (Site_Name) на xaxis.

ggplot (медианный список,aes(x= изменить порядок (Site_Name,MedMean,FUN = median),y=MedMean))+geom_boxplot()

Я хочу добавить буквы значимости Тьюки в коробки.

Спасибо

Answer 1

С помощью agricolae::HSD.test ты можешь сделать

library(dplyr)
library(agricolae)
library(ggplot2)
iris.summarized = iris %>% group_by(Species) %>% summarize(Max.Petal.Length=max(Petal.Length))
hsd=HSD.test(aov(Petal.Length~Species,data=iris), "Species", group=T)
ggplot(iris,aes(x=reorder(Species,Petal.Length,FUN = median),y=Petal.Length))+geom_boxplot()+geom_text(data=iris.summarized,aes(x=Species,y=0.2+Max.Petal.Length,label=hsd$groups$groups),vjust=0)