I have generated these two boxplots. But I am unable to remove the points and the outliers. I already have used outliers.shape = NA which didn't work. I also want to fill the groups like this image. Can anyone please tell me how can I do that?

here's my code:

p <- plot_richness (physeq_rarefied, x="type", color = "type", measures=c("Shannon", "Observed")) +
stat_compare_means(method = "wilcox.test") +
geom_boxplot() +
labs(x= "Sample types", y= "Alpha Diversity Measure", 
title = "Alpha diversity of control and test samples")

Thanks, dpc