# The Holzinger and Swineford (1939) example
HS.model <- ' visual =~ x1 + x2 + x3
textual =~ x4 + x5 + x6
speed =~ x7 + x8 + x9 '
# a data generating function
generateData <- function() simulateData(HS.model, sample.nobs = 100)
fit <- semList(HS.model, dataFunction = generateData, ndat = 10,
store.slots = "partable")
# show parameter estimates, per dataset
coef(fit)
Run the code above in your browser using DataLab