Data Analysis Code
ANOVA
library(ggplot2)
library(ggpubr)
library(tidyverse)
library(broom)
library(AICcmodavg)

data <- Analysis
data1_omit <- data %>% select("SubjectID", "Group", "stroop", "antisaccadelatency", "antisaccadepercentcorrect")

ANTILATENCY <- aov(antisaccadelatency ~ Group,data = data1_omit)
summary(ANTILATENCY)

TukeyHSD(ANTILATENCY)


ANTICORRECT <- aov(antisaccadepercentcorrect ~ Group,data = data1_omit)
summary(ANTICORRECT)
TukeyHSD(ANTICORRECT)

STROOP <- aov(stroop ~ Group,data = data1_omit)
summary(STROOP)
TukeyHSD(STROOP)

Linear Mixed Effects Modelling

library(broom)
library(gridExtra)
library(here)
library(lme4)
library(patchwork)
library(sjPlot)
library(tidyverse)
library("lmerTest")

data <- Analysis
View (data)

data_omit <- data %>% select("SubjectID", "Group", "product", "differencescore", "stroop", "antisaccadelatency", "antisaccadepercentcorrect", "ACEscore")

data1.lmer <-lmer(differencescore ~ Group + stroop + antisaccadelatency + antisaccadepercentcorrect + ACEscore + (1|subjectID) + (1|product),data = data_omit)
summary(data1.lmer)