
library(metafor)

year2attdata <- read.csv("C:/Users/genni/OneDrive/Documents/Dissertation/Paper 1 - Systematic Review/Meta-Analysis/year2att.csv")


#Year 2 Disengagement of Attention Meta-Analysis
year2att <- rma(yi=z, 
                vi=var,
                method="HE", 
                data=year2attdata)

summary(year2att)
forest(year2att,
       order=order(year2attdata$Author),
       xlab="Correlation Coefficient",
       slab=paste(year2attdata$Author, year2attdata$Year, sep=" "),
       transf=transf.ztor
)

#Year 2 Attention with MODERATOR TASK (Frustration vs Fear)
atty2task <- rma(yi=z,vi=var, method="HE", mods=~factor(task),data=year2attdata)

summary(atty2task)

#####

#Year 2 Parent Focused Meta-Analysis

year2pardata <- read.csv("C:/Users/genni/OneDrive/Documents/Dissertation/Paper 1 - Systematic Review/Meta-Analysis/year2par.csv")

year2par <- rma(yi=z, 
                vi=var,
                method="HE", 
                data=year2pardata)

summary(year2par)
forest(year2par,
       order=order(year2pardata$Author),
       xlab="Correlation Coefficient",
       slab=paste(year2pardata$Author, year2pardata$Year, sep=" "),
       transf=transf.ztor
)

#Year 2 Parent Focused with MODERATOR TASK (Frustration vs Fear)
pary2task <- rma(yi=z,vi=var, method="HE", mods=~factor(task),data=year2pardata)

summary(pary2task)


#####

#Year 2 Self Focused Meta-Analysis

year2selfdata <- read.csv("C:/Users/genni/OneDrive/Documents/Dissertation/Paper 1 - Systematic Review/Meta-Analysis/year2self.csv")

year2self <- rma(yi=z, 
                 vi=var,
                 method="HE", 
                 data=year2selfdata)

summary(year2self)
forest(year2self,
       order=order(year2selfdata$Author),
       xlab="Correlation Coefficient",
       slab=paste(year2selfdata$Author, year2selfdata$Year, sep=" "),
       transf=transf.ztor
)


#Year 2 Self Focused with MODERATOR TASK (Frustration vs Fear)

selfy2task <- rma(yi=z,vi=var, method="HE", mods=~factor(task),data=year2selfdata)

summary(selfy2task)


