37 lines
No EOL
1.1 KiB
R
Executable file
37 lines
No EOL
1.1 KiB
R
Executable file
#!/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") |