if(interactive()){
data(mstData)
###############
##### SRT #####
###############
outputSRT <- srtFREQ(Posttest~ Intervention + Prettest,
intervention = "Intervention", data = mstData)
outputSRTBoot <- srtFREQ(Posttest~ Intervention + Prettest,
intervention = "Intervention",nBoot=1000, data = mstData)
###############
##### MST #####
###############
outputMST <- mstFREQ(Posttest~ Intervention + Prettest,
random = "School", intervention = "Intervention", data = mstData)
outputMSTBoot <- mstFREQ(Posttest~ Intervention + Prettest,
random = "School", intervention = "Intervention",
nBoot = 1000, data = mstData)
##################
##### Bayesian #####
##################
outputMSTbayes <- mlmBayes(Posttest~ Intervention + Prettest,
random = "School", intervention = "Intervention",
nSim = 10000, data = mstData)
## comparing different results
ComparePlot(list(outputSRT,outputSRTBoot,outputMST,outputMSTBoot,outputMSTbayes),
modelNames =c("ols", "olsBoot","MLM","MLMBoot","MLMBayes"),group=1)
}
Run the code above in your browser using DataLab