R
TukeyHSD
stats
multicomp
ggboxplot(clean_twoway_data, x = "Fertilizer", y = "Yield", fill = "Manure")