I have tried looking at examples and solutions for this. However, none of them seemed to work. My X-axis labels are overlapping on each other. Thus, I am unable to see which boxplot corresponds to which label. I am comparing gene expression in certain tissue types in different diseases (MND, OND, and NCC).
Here is the R script
my_comparisons <- list(c("MND", "NNC"),c("MND", "OND"),c("NNC", "OND"))
ggboxplot(subset(rips, !is.na(Tissue)), x = "Tissue", y = "ENST00000504942.6",
color = "Tissue", palette = "jco",
ylab = "Expression of ENST00000504942.6", xlab = "Tissue",
add = "jitter")
stat_compare_means(comparisons = my_comparisons)
And this is the output.
I appreciate your help!