updated header file with packes for summary stats

Header_TT.R

@@ -6,6 +6,29 @@ library(ggpubr)
 library(rstatix)
 library(Hmisc)
 library(qwraps2)
-source("legend_adjustment.R")
+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
 
-# https://www.datanovia.com/en/blog/how-to-add-p-values-onto-a-grouped-ggplot-using-the-ggpubr-r-package/#perform-all-pairwise-comparisons
+
+########################################################################
+#                       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")