# \donttest{
library(gesca)
data(gesca.rick2) # Organizational identification example of Bagozzi
# Model specification
myModel <- "
# Measurement model
OP =~ cei1 + cei2 + cei3
OI =~ ma1 + ma2 + ma3
AC_J =~ orgcmt1 + orgcmt2 + orgcmt3
AC_L =~ orgcmt5 + orgcmt6 + orgcmt8
# Structural model
OI ~ OP
AC_J ~ OI
AC_L ~ OI
"
# Run a multiple-group GSCA with the grouping variable gender:
GSCA.group <- gesca.run(myModel, gesca.rick2, group.name = "gender", nbt=50)
summary(GSCA.group) # default output shows model fit indices as well, but
fitmeasures(GSCA.group) # this gives bootstrapped std errors and CI
# }
Run the code above in your browser using DataLab