ESEM/exampleScript_ESEMS.R at master · pdparker/ESEM · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
#Install Mplus Automation if Needed and load the package
if("MplusAutomation" %in% rownames(installed.packages())==FALSE) install.packages("MplusAutomation")
require("MplusAutomation")
#Choose a directory for the Mplus Files etc.
cat("Choose a directory that you want to use.
Mplus files will all be sent to there and run from there.
It is best to choose an empty directory or at least one with no existing mplus scripts.")
tempDir<- choose.dir()
#Load data
dataURL <- "http://raw.github.com/pdparker/ESEM/master/ESEM_SIM.Rdata"
load(url(dataURL))
closeAllConnections()
#Run ESEM Invariance Script
scriptURL <- "http://raw.github.com/pdparker/ESEM/master/ESEM_Invariance_Geomin.R"
source(url(scriptURL))
closeAllConnections()
#Prepare Mplus data
prepareMplusData(data, filename=file.path(tempDir, "ESEMdata.dat"))

#---------------------------------------------------#
#Time 2 invariance across treatment group and gender#
#---------------------------------------------------#
#Write Esem Scripts
MplusData <- file.path(tempDir, "ESEMdata.dat")
esemInvaGeomin(2, data, GroupVar = "TbyG", c("treatM", "contM", "treatF", "contF"),
               1:6, FileOut=tempDir, FileIn=MplusData)

#Run Models
runModels(tempDir)

#Readout fit and place in table
FitSummaries <- extractModelSummaries(tempDir)
showSummaryTable(FitSummaries, keepCols=c("Title", "ChiSqM_Value", "ChiSqM_DF",
                                          "CFI", "TLI", "RMSEA_Estimate",
                                          "RMSEA_90CI_LB", "RMSEA_90CI_UB"),
                 sortBy="Title")

#Clean up directory by deleting all mplus files
junk <- dir(path=tempDir,  pattern=".inp$|.out$", full.names = TRUE)
file.remove(junk)

#---------------------------------------------------#
#Time 2 invariance across treatment group#
#---------------------------------------------------#
#Write Esem Scripts
MplusData <- file.path(tempDir, "ESEMdata.dat")
esemInvaGeomin(2, data, GroupVar = "treat", c("treat", "control"),
               1:6, FileOut=tempDir, FileIn=MplusData)

#Run Models
runModels(tempDir)
#Readout fit and place in table
FitSummaries <- extractModelSummaries(tempDir)
showSummaryTable(FitSummaries, keepCols=c("Title", "ChiSqM_Value", "ChiSqM_DF",
                                          "CFI", "TLI", "RMSEA_Estimate",
                                          "RMSEA_90CI_LB", "RMSEA_90CI_UB"),
                 sortBy="Title")
#Clean up directory by deleting all mplus files
junk <- dir(path=tempDir,  pattern=".inp$|.out$", full.names = TRUE)
file.remove(junk)