## You need to import this data into your R workspace using
## the following code
simulated=read.csv("simulated.csv")
#### Next we need to perform Step1 of mediation analysis
s1=glm(simulated[,3]~simulated[,1],family=binomial)
summary(s1)
#### Next we perform Step2 as follows
s2=glm(simulated[,2]~simulated[,1],family=binomial)
summary(s2)
##### And step 3
s3=glm(simulated[,3]~simulated[,2],family=binomial)
summary(s3)
##### Finally, if we find that all the coefficients returned above are statistically
##### significant we need to proceed to step 4
s4=glm(simulated[,3]~simulated[,2] + simulated[,1],family=binomial)
summary(s4)
#### We now compare "s1" to "s4" to verify mediation
summary(s1)
summary(s4)
###### As will be clear, the parameter was statistically significant in one
###### model (step1) but not in the other (step4).