Hi I am generating a plot of control and case sample's diversity indices with the following code:
p <- plot_richness(physeq_rarefied, x="type", measures=c("Shannon", "Observed")) + stat_compare_means(method = "wilcox.test") + geom_boxplot(aes(fill = type), notch = TRUE, outlier.size = -1) + scale_fill_manual(values = c("hotpink", "skyblue"))+ labs(x= "Sample types", y= "Alpha Diversity Measure", fill ="Groups", title = "Alpha diversity of control and case samples") p p$layers <- NULL p
I just get two notch plotes, one for case and another for control samples. Now, I want to get the median value for the case and control groups in a table format and export them. How can I do that? Can anyone please help me?