From 15bba91ceee6b5299f8c46c96d021cc1eb238a93 Mon Sep 17 00:00:00 2001
From: Tanushree Tunstall <tanu@tunstall.in>
Date: Fri, 30 Oct 2020 11:53:26 +0000
Subject: [PATCH] updated header file with packes for summary stats

---
 Header_TT.R | 27 +++++++++++++++++++++++++--
 1 file changed, 25 insertions(+), 2 deletions(-)

diff --git a/Header_TT.R b/Header_TT.R
index 787b6c6..3282de3 100755
--- a/Header_TT.R
+++ b/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
\ No newline at end of file
+
+########################################################################
+#                       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")
\ No newline at end of file