#This program was written for the book chapter
#Use of Monte Carlo Studies in Structural Equation Modeling Research
#code to run the simulation study

setwd("F:/walter/research/simulation/Book_chapter") #set working directory
set.seed(741355) #set random seed

#loop through conditions (in this case, two sample sizes)
for (sample in c(150,500) ) {

#specify all population parameters
lambda =  matrix (c(1,0,1,0,1,0,1,0,0,1,0,1,0,1,0,1),8,2,byrow = T)
phi = matrix(c(0.5,0.25,0.25,0.5),2,2,byrow = T)
theta = diag( 0.5,8)
epsilon = rep(1.39,8)

#read functions
source("function_montecarlo.r")
source("function_run_sem.r")

#simulate 100 datasets
montecarlo(lambda, phi, theta, epsilon,sample,iter=500)
  
#run MPLUS fot fit a CFA to the simulated datasets
run_sem("MPLUS_analysis",number.results=45,conditions=sample)

} #close loop through conditions