library(tidyverse)
library(data.table)
library(readxl)
library(emmeans)
library(lm.beta)
EF <- read.csv("MEANS_EFs.csv", header = TRUE, na.strings=c(""," ","NA"))
SI_SigTrans<- read.csv("SI_SigTrans.csv", header = TRUE, na.strings=c(""," ","NA"))
Data <- Reduce(inner_join, list (EF,SI_SigTrans))
ALLdata <- Data %>% dplyr::select( 'ParticipantID','MeasureTime','Condition', 'FT_EFF','FM_EFF','CB_EFF_backward','SigTran','SigTran_InfoCom', 'SigTran_FluDel','SigTran_TLQual', 'SI', 'SI_InfoCom','SI_FluDel','SI_TLQual')
summary(ALLdata)
ALLdata_pre <- ALLdata[ALLdata$MeasureTime == "pretest",]
ALLdata_post <- ALLdata[ALLdata$MeasureTime == "posttest",]
ALLdata_pre <-  rename (ALLdata_pre, CB_EFF_pre = CB_EFF_backward)%>%
rename (FT_EFF_pre = FT_EFF) %>%
rename (FM_EFF_pre = FM_EFF) %>%
rename (SigTran_pre = SigTran) %>%
rename (SigTran_InfoCom_pre = SigTran_InfoCom)%>%
rename (SigTran_FluDel_pre = SigTran_FluDel)%>%
rename (SigTran_TLQual_pre = SigTran_TLQual)%>%
rename (SI_pre = SI)%>%
rename (SI_InfoCom_pre = SI_InfoCom) %>%
rename (SI_FluDel_pre = SI_FluDel)%>%
rename (SI_TLQual_pre = SI_TLQual)%>%
rename (MeasureTime_pre = MeasureTime) %>%
rename (Condition_pre = Condition)
ALLdata_post <-  rename (ALLdata_post, CB_EFF_post = CB_EFF_backward)%>%
rename (FT_EFF_post = FT_EFF) %>%
rename (FM_EFF_post = FM_EFF) %>%
rename (SigTran_post = SigTran) %>%
rename (SigTran_InfoCom_post = SigTran_InfoCom)%>%
rename (SigTran_FluDel_post = SigTran_FluDel)%>%
rename (SigTran_TLQual_post = SigTran_TLQual)%>%
rename (SI_post = SI)%>%
rename (SI_InfoCom_post = SI_InfoCom) %>%
rename (SI_FluDel_post = SI_FluDel)%>%
rename (SI_TLQual_post = SI_TLQual)%>%
rename (MeasureTime_post = MeasureTime) %>%
rename (Condition_post = Condition)
ALLdata_difference <- Reduce(inner_join, list(ALLdata_pre, ALLdata_post))%>%
mutate(FT_EFF_diff = FT_EFF_post - FT_EFF_pre)%>%
mutate(FM_EFF_diff = FM_EFF_post - FM_EFF_pre)%>%
mutate(CB_EFF_diff = CB_EFF_post - CB_EFF_pre)%>%
mutate (SigTran_diff = SigTran_post - SigTran_pre) %>%
mutate (SigTran_InfoCom_diff = SigTran_InfoCom_post - SigTran_InfoCom_pre)%>%
mutate (SigTran_FluDel_diff = SigTran_FluDel_post - SigTran_FluDel_pre) %>%
mutate (SigTran_TLQual_diff = SigTran_TLQual_post - SigTran_TLQual_pre)%>%
mutate (SI_diff = SI_post - SI_pre)%>%
mutate (SI_InfoCom_diff = SI_InfoCom_post - SI_InfoCom_pre) %>%
mutate (SI_FluDel_diff = SI_FluDel_post - SI_FluDel_pre)%>%
mutate (SI_TLQual_diff = SI_TLQual_post -SI_TLQual_pre)%>%
rename (Condition = Condition_post)
control <- ALLdata_difference %>% filter(Condition == "control")
experimental <- ALLdata_difference %>% filter(Condition == "experiment")
#pre-test
shapiro.test(control$FT_EFF_pre)
shapiro.test(control$FM_EFF_pre)
shapiro.test(control$CB_EFF_pre)
shapiro.test(control$SigTran_pre)
shapiro.test(control$SigTran_InfoCom_pre)
shapiro.test(control$SigTran_FluDel_pre)
shapiro.test(control$SigTran_TLQual_pre)
shapiro.test(control$SI_pre)
shapiro.test(control$SI_InfoCom_pre)
shapiro.test(control$SI_FluDel_pre)
shapiro.test(control$SI_TLQual_pre)
shapiro.test(experimental$FT_EFF_pre)
shapiro.test(experimental$FM_EFF_pre)
shapiro.test(experimental$CB_EFF_pre)
shapiro.test(experimental$SigTran_pre)
shapiro.test(experimental$SigTran_InfoCom_pre)
shapiro.test(experimental$SigTran_FluDel_pre)
shapiro.test(experimental$SigTran_TLQual_pre)
shapiro.test(experimental$SI_pre)
shapiro.test(experimental$SI_InfoCom_pre)
shapiro.test(experimental$SI_FluDel_pre)
shapiro.test(experimental$SI_TLQual_pre)
#post-test
shapiro.test(control$FT_EFF_post)
shapiro.test(control$FM_EFF_post)
shapiro.test(control$CB_EFF_post)
shapiro.test(control$SigTran_post)
shapiro.test(control$SigTran_InfoCom_post)
shapiro.test(control$SigTran_FluDel_post)
shapiro.test(control$SigTran_TLQual_post)
shapiro.test(control$SI_post)
shapiro.test(control$SI_InfoCom_post)
shapiro.test(control$SI_FluDel_post)
shapiro.test(control$SI_TLQual_post)
shapiro.test(experimental$FT_EFF_post)
shapiro.test(experimental$FM_EFF_post)
shapiro.test(experimental$CB_EFF_post)
shapiro.test(experimental$SigTran_post)
shapiro.test(experimental$SigTran_InfoCom_post)
shapiro.test(experimental$SigTran_FluDel_post)
shapiro.test(experimental$SigTran_TLQual_post)
shapiro.test(experimental$SI_post)
shapiro.test(experimental$SI_InfoCom_post)
shapiro.test(experimental$SI_FluDel_post)
shapiro.test(experimental$SI_TLQual_post)
#differences
shapiro.test(ALLdata_difference$FT_EFF_diff)
shapiro.test(ALLdata_difference$FM_EFF_diff)
shapiro.test(ALLdata_difference$CB_EFF_diff)
shapiro.test(ALLdata_difference$SigTran_diff)
shapiro.test(ALLdata_difference$SigTran_InfoCom_diff)
shapiro.test(ALLdata_difference$SigTran_FluDel_diff)
shapiro.test(ALLdata_difference$SigTran_TLQual_diff)
shapiro.test(ALLdata_difference$SI_diff)
shapiro.test(ALLdata_difference$SI_InfoCom_diff)
shapiro.test(ALLdata_difference$SI_FluDel_diff)
shapiro.test(ALLdata_difference$SI_TLQual_diff)
#EF
aov_FM <- aov(
rank(FM_EFF) ~ MeasureTime*Condition+ Error(ParticipantID), data = ALLdata)
print(summary(aov_FM),digits=10)
# posthoc
print(summary(pairs(emmeans(aov_FM, ~ MeasureTime),adjust = "tukey")),digits=10)
print(summary(pairs(emmeans(aov_FM, ~ Condition),adjust = "tukey")),digits=10)
print(summary(pairs(emmeans(aov_FM, pairwise ~ MeasureTime | Condition),adjust = "tukey")),digits=10)
aov_CB <- aov(rank(CB_EFF_backward) ~ MeasureTime*Condition+ Error(ParticipantID), data = ALLdata)
print(summary(aov_CB),digits=10)
# posthoc
print(summary(pairs(emmeans(aov_CB, ~ MeasureTime),adjust = "tukey")),digits=10)
aov_FT <- aov(
rank(FT_EFF) ~ MeasureTime*Condition+ Error(ParticipantID), data = ALLdata)
print(summary(aov_FT),digits=10)
# posthoc
print(summary(pairs(emmeans(aov_FT, ~ MeasureTime),adjust = "tukey")),digits=10)
#SigTran
aov_SigTran <- aov(
rank(SigTran) ~ MeasureTime*Condition+ Error(ParticipantID), data = ALLdata)
print(summary(aov_SigTran),digits=10)
aov_SigTran_InfoCom <- aov(
rank(SigTran_InfoCom) ~ MeasureTime*Condition+ Error(ParticipantID), data = ALLdata)
print(summary(aov_SigTran_InfoCom),digits=10)
aov_SigTran_TLQual <- aov(
rank(SigTran_TLQual) ~ MeasureTime*Condition+ Error(ParticipantID), data = ALLdata)
print(summary(aov_SigTran_TLQual),digits=10)
aov_SigTran_FluDel <- aov(
rank(SigTran_FluDel) ~ MeasureTime*Condition+ Error(ParticipantID), data = ALLdata)
print(summary(aov_SigTran_FluDel),digits=10)
# posthoc
print(summary(pairs(emmeans(aov_SigTran_FluDel, ~ MeasureTime),adjust = "tukey")),digits=10)
print(summary(pairs(emmeans(aov_SigTran_FluDel, pairwise ~ MeasureTime | Condition),adjust = "tukey")),digits=10)
#SI
aov_SI <- aov(
rank(SI) ~ MeasureTime*Condition+ Error(ParticipantID), data = ALLdata)
print(summary(aov_SI),digits=10)
# posthoc
print(summary(pairs(emmeans(aov_SI, ~ MeasureTime),adjust = "tukey")),digits=10)
print(summary(pairs(emmeans(aov_SI, pairwise ~ MeasureTime | Condition),adjust = "tukey")),digits=10)
aov_SI_InfoCom <- aov(
rank(SI_InfoCom) ~ MeasureTime*Condition+ Error(ParticipantID), data = ALLdata)
print(summary(aov_SI_InfoCom),digits=10)
# posthoc
print(summary(pairs(emmeans(aov_SI_InfoCom, ~ MeasureTime),adjust = "tukey")),digits=10)
print(summary(pairs(emmeans(aov_SI_InfoCom, pairwise ~ MeasureTime | Condition),adjust = "tukey")),digits=10)
aov_SI_TLQual <- aov(
rank(SI_TLQual) ~ MeasureTime*Condition+ Error(ParticipantID), data = ALLdata)
print(summary(aov_SI_TLQual),digits=10)
# posthoc
print(summary(pairs(emmeans(aov_SI_TLQual, ~ MeasureTime),adjust = "tukey")),digits=10)
print(summary(pairs(emmeans(aov_SI_TLQual, pairwise ~ MeasureTime | Condition),adjust = "tukey")),digits=10)
aov_SI_FluDel <- aov(
rank(SI_FluDel) ~ MeasureTime*Condition+ Error(ParticipantID), data = ALLdata)
print(summary(aov_SI_FluDel),digits=10)
# posthoc
print(summary(pairs(emmeans(aov_SI_FluDel, ~ MeasureTime),adjust = "tukey")),digits=10)
print(summary(pairs(emmeans(aov_SI_FluDel, pairwise ~ MeasureTime | Condition),adjust = "tukey")),digits=10)
pre_FT_Regression = lm(experimental$FT_EFF_diff ~ experimental$FT_EFF_pre)
beta <- lm.beta(pre_FT_Regression)
summary(beta)
confint(beta)
pre_FM_Regression = lm(experimental$FM_EFF_diff ~ experimental$FM_EFF_pre)
beta <- lm.beta(pre_FM_Regression)
summary(beta)
confint(beta)
pre_CB_Regression = lm(experimental$CB_EFF_diff ~ experimental$CB_EFF_pre)
beta <- lm.beta(pre_CB_Regression)
summary(beta)
confint(beta)
SigTran_Regression_diff = lm(experimental$SigTran_diff ~ experimental$FT_EFF_pre +experimental$FM_EFF_pre +experimental$CB_EFF_pre )
beta <- lm.beta(SigTran_Regression_diff)
summary(beta)
confint(beta)
SigTran_InfoCom_Regression_diff = lm(experimental$SigTran_InfoCom_diff ~ experimental$FT_EFF_pre +experimental$FM_EFF_pre +experimental$CB_EFF_pre )
beta <- lm.beta(SigTran_InfoCom_Regression_diff)
summary(beta)
confint(beta)
SigTran_FluDel_Regression_diff = lm(experimental$SigTran_FluDel_diff ~ experimental$FT_EFF_pre +experimental$FM_EFF_pre +experimental$CB_EFF_pre )
beta <- lm.beta(SigTran_FluDel_Regression_diff)
summary(beta)
confint(beta)
SigTran_TLQual_Regression_diff = lm(experimental$SigTran_TLQual_diff ~ experimental$FT_EFF_pre +experimental$FM_EFF_pre +experimental$CB_EFF_pre )
beta <- lm.beta(SigTran_TLQual_Regression_diff)
summary(beta)
confint(beta)
SI_Regression_diff = lm(experimental$SI_diff ~ experimental$FT_EFF_pre +experimental$FM_EFF_pre +experimental$CB_EFF_pre )
beta <- lm.beta(SI_Regression_diff)
summary(beta)
confint(beta)
SI_InfoCom_Regression_diff = lm(experimental$SI_InfoCom_diff ~ experimental$FT_EFF_pre +experimental$FM_EFF_pre +experimental$CB_EFF_pre )
beta <- lm.beta(SI_InfoCom_Regression_diff)
summary(beta)
confint(beta)
SI_FluDel_Regression_diff = lm(experimental$SI_FluDel_diff ~ experimental$FT_EFF_pre +experimental$FM_EFF_pre +experimental$CB_EFF_pre )
beta <- lm.beta(SI_FluDel_Regression_diff)
summary(beta)
confint(beta)
SI_TLQual_Regression_diff = lm(experimental$SI_TLQual_diff ~ experimental$FT_EFF_pre +experimental$FM_EFF_pre +experimental$CB_EFF_pre )
beta <- lm.beta(SI_TLQual_Regression_diff)
summary(beta)
confint(beta)
pre_FT_Regression = lm(experimental$FT_EFF_diff ~ experimental$FT_EFF_pre)
beta_FT <- lm.beta(pre_FT_Regression)
summary(beta_FT)
confint(beta_FT)
pre_FM_Regression = lm(experimental$FM_EFF_diff ~ experimental$FM_EFF_pre)
beta_FM <- lm.beta(pre_FM_Regression)
summary(beta_FM)
confint(beta_FM)
pre_CB_Regression = lm(experimental$CB_EFF_diff ~ experimental$CB_EFF_pre)
beta_CB <- lm.beta(pre_CB_Regression)
summary(beta_CB)
confint(beta_CB)
SigTran_Regression_diff = lm(experimental$SigTran_diff ~ experimental$FT_EFF_pre +experimental$FM_EFF_pre +experimental$CB_EFF_pre )
beta_SigTran <- lm.beta(SigTran_Regression_diff)
summary(beta_SigTran)
confint(beta_SigTran)
SigTran_InfoCom_Regression_diff = lm(experimental$SigTran_InfoCom_diff ~ experimental$FT_EFF_pre +experimental$FM_EFF_pre +experimental$CB_EFF_pre )
beta_SigTran_InfoCom <- lm.beta(SigTran_InfoCom_Regression_diff)
summary(beta_SigTran_InfoCom)
confint(beta_SigTran_InfoCom)
#FluDel
SigTran_FluDel_Regression_diff = lm(experimental$SigTran_FluDel_diff ~ experimental$FT_EFF_pre +experimental$FM_EFF_pre +experimental$CB_EFF_pre )
beta_SigTran_FluDel <- lm.beta(SigTran_FluDel_Regression_diff)
summary(beta_SigTran_FluDel)
confint(beta_SigTran_FluDel)
SigTran_TLQual_Regression_diff = lm(experimental$SigTran_TLQual_diff ~ experimental$FT_EFF_pre +experimental$FM_EFF_pre +experimental$CB_EFF_pre )
beta_SigTran_TLQual <- lm.beta(SigTran_TLQual_Regression_diff)
summary(beta_SigTran_TLQual)
confint(beta_SigTran_TLQual)
SI_Regression_diff = lm(experimental$SI_diff ~ experimental$FT_EFF_pre +experimental$FM_EFF_pre +experimental$CB_EFF_pre )
beta_SI <- lm.beta(SI_Regression_diff)
summary(beta_SI)
confint(beta_SI)
#InfoCom
SI_InfoCom_Regression_diff = lm(experimental$SI_InfoCom_diff ~ experimental$FT_EFF_pre +experimental$FM_EFF_pre +experimental$CB_EFF_pre )
beta_SI_InfoCom <- lm.beta(SI_InfoCom_Regression_diff)
summary(beta_SI_InfoCom)
confint(beta_SI_InfoCom)
#FluDel
SI_FluDel_Regression_diff = lm(experimental$SI_FluDel_diff ~ experimental$FT_EFF_pre +experimental$FM_EFF_pre +experimental$CB_EFF_pre )
beta_SI_FluDel <- lm.beta(SI_FluDel_Regression_diff)
summary(beta_SI_FluDel)
confint(beta_SI_FluDel)
#TLQual
SI_TLQual_Regression_diff = lm(experimental$SI_TLQual_diff ~ experimental$FT_EFF_pre +experimental$FM_EFF_pre +experimental$CB_EFF_pre )
beta_SI_TLQual <- lm.beta(SI_TLQual_Regression_diff)
summary(beta_SI_TLQual)
confint(beta_SI_TLQual)
library(tidyverse)
library(data.table)
library(readxl)
library(emmeans)
library(lm.beta)
EF <- read.csv("MEANS_EFs.csv", header = TRUE, na.strings=c(""," ","NA"))
SI_SigTrans<- read.csv("SI_STR.csv", header = TRUE, na.strings=c(""," ","NA"))
Data <- Reduce(inner_join, list (EF,SI_STR))
SI_STR<- read.csv("SI_STR.csv", header = TRUE, na.strings=c(""," ","NA"))
EF <- read.csv("MEANS_EFs.csv", header = TRUE, na.strings=c(""," ","NA"))
SI_STR<- read.csv("SI_STR.csv", header = TRUE, na.strings=c(""," ","NA"))
Data <- Reduce(inner_join, list (EF,SI_STR))
ALLdata <- Data %>% dplyr::select( 'ParticipantID','MeasureTime','Condition', 'FT_EFF','FM_EFF','CB_EFF_backward','SigTran','SigTran_InfoCom', 'SigTran_FluDel','SigTran_TLQual', 'SI', 'SI_InfoCom','SI_FluDel','SI_TLQual')
summary(ALLdata)
library(tidyverse)
library(data.table)
library(readxl)
library(emmeans)
library(lm.beta)
EF <- read.csv("MEANS_EFs.csv", header = TRUE, na.strings=c(""," ","NA"))
SI_SigTrans<- read.csv("SI_SigTrans.csv", header = TRUE, na.strings=c(""," ","NA"))
Data <- Reduce(inner_join, list (EF,SI_SigTrans))
ALLdata <- Data %>% dplyr::select( 'ParticipantID','MeasureTime','Condition', 'FT_EFF','FM_EFF','CB_EFF_backward','SigTran','SigTran_InfoCom', 'SigTran_FluDel','SigTran_TLQual', 'SI', 'SI_InfoCom','SI_FluDel','SI_TLQual')
summary(ALLdata)
