if(!require(ggpubr))install.packages('ggpubr')
library(ggpubr)
ggscatter(iris,x="Sepal.Length",y="Petal.Length",color="Species")

image.png

p <- ggboxplot(iris, x = "Species", y = "Sepal.Length",
          color = "Species", shape = "Species",
          add = "jitter")
p

image.png

my_comparisons <- list( c("setosa", "versicolor"), c("setosa", "virginica"), c("versicolor", "virginica") )
p + stat_compare_means(comparisons = my_comparisons)+ # Add pairwise comparisons p-value
  stat_compare_means(label.y = 9) 

image.png