#!/usr/bin/Rscript   
#install.packages("stringr")
library(stringr)
library(tidyverse)
library(ggpubr)
library(rstatix)
library(Hmisc)
library(qwraps2)
source("legend_adjustment.R")

# https://www.datanovia.com/en/blog/how-to-add-p-values-onto-a-grouped-ggplot-using-the-ggpubr-r-package/#perform-all-pairwise-comparisons