Я сделал ggplot с моими данными.чем я хотел обозначить некоторые взаимодействия между моими данными boxplot (через anova / Tukey), чтобы показать, есть ли существенные различия между специальными группами или нет.Я работаю с одним листом csv под названием «test», который я импортировал, и с одним табелем с именем «final»
[
[
но когда я строю свои данные с меткой из моего TukeyHSD, я получаю каждую метку тройной с пустыми пробелами.
как я могу удалить пустые места.должно быть какое-то решение, чтобы просто показать первые линии дерева с меткой на первом графике, средние линии дерева на среднем графике с их конкретной меткой и последние линии дерева на третьем месте с меткой.Как я могу удалить эти пустые строки с их алфавитными метками?
Вот мой код:
test$int<-interaction(test$geno,test$n)
generate_label_df <- function(TUKEY, variable){
Tukey.levels <- TUKEY[[variable]][,4]
Tukey.labels <- data.frame(multcompLetters(Tukey.levels)['Letters'])
Tukey.labels$int=rownames(Tukey.labels)
Tukey.labels=Tukey.labels[order(Tukey.labels$int) , ]
return(Tukey.labels)
}
model=lm(test$amino~test$int ) ###
ANOVA=aov(model)
TUKEY <- TukeyHSD(x=ANOVA, 'test$int', conf.level=0.95)
labels<-generate_label_df(TUKEY , "test$int")
names(labels)<-c('Letters','int')
yvalue<-aggregate(amino~int, data=test, max)
final<-merge(labels,yvalue)
ggplot(test, aes(x=int, y=amino)) +
stat_boxplot(geom="errorbar") +
geom_boxplot(aes(fill=CO2)) +
geom_text(data = final, aes(label = Letters, vjust = -0.6)) +
labs(title = "Aminosäuren Erdkultur", x = "Genotyp", y = "Aminosäurekonz. [µmol/gFW]") +
scale_x_discrete(labels = c("col.Ammonium" = "col", "sps.Ammonium" = "sps",
"swe.Ammonium" = "swe", "col.Nitrat" = "col",
"sps.Nitrat" = "sps", "swe.Nitrat" = "swe",
"col.NON" = "col", "sps.NON" = "sps","swe.NON" = "swe"),
limits = c("col.Ammonium", "sps.Ammonium", "swe.Ammonium",
"col.Nitrat", "sps.Nitrat","swe.Nitrat", "col.NON",
"sps.NON", "swe.NON") )+
guides(fill = guide_legend(title = "CO2-Behandlung")) +
theme(plot.title = element_text(face = "bold", size="17", hjust = "0.5")) +
scale_fill_manual(values = c("violetred3", "steelblue1")) +
facet_wrap( ~ n) +
theme(legend.title = element_text(face = "bold"))
Мое первое мнение было, чтобы удалить "facet_wrap (~ n)".Проблема решена чем, но, очевидно, мне не хватает подразделения на этом пути