Guidali_et_el_2023_EJN_RR/Dataset and analysis/Robust rmANOVA_rMT_Biphasic_explorative.R

setwd("/Users/giacomoguidali/Library/CloudStorage/GoogleDrive-g.guidali@campus.unimib.it/Il mio Drive/Post-Doc BS/MoBi P15/4 - Analisi/MEP") #set directory

rm(list = ls()) # pulisco il workspace, male non fa
library(WRS2)

MEPamplitude_raw <- read.csv("Database rMT_robustANOVABiphasic.csv",header=T,sep=";", dec=",") #importo database

#rmANOVA + post-hoc
rmanova(MEPamplitude_raw$rMT, MEPamplitude_raw$Condition, MEPamplitude_raw$Subject , tr = 0.2)
rmmcp(MEPamplitude_raw$rMT, MEPamplitude_raw$Condition, MEPamplitude_raw$Subject, tr = 0.2, alpha = 0.05)