#!/usr/bin/Rscript   
#install.packages("stringr")
library(stringr)
library(tidyverse)
library(ggpubr)
library(rstatix)
library(Hmisc)
library(qwraps2)
library(plyr)
library(psych)
library(arsenal)
library(scales)
library(e1071)
# https://www.datanovia.com/en/blog/how-to-add-p-values-onto-a-grouped-ggplot-using-the-ggpubr-r-package/#perform-all-pairwise-comparisons
library(survival)
library(CARS) # for calculating VIF (changed from cars)
#library(BaylorEdPsych)


########################################################################
#                       My functions
########################################################################
#my_kurtosis <- function(x, na.rm = T) {  
#  m4 <- mean((x - mean(x))^4) 
#  kurtosis <- m4/(sd(x)^4) - 3  
#  kurtosis
#}

#my_skewness <-  function(x, na.rm = T) {
#  m3 <- mean((x - mean(x))^3)
#  skewness <- m3/(sd(x)^3)
#  skewness
#}
########################################################################
#                       My local imports
########################################################################
source("legend_adjustment.R")