# NOT RUN {
data(demoOneFactor)
manifests <- names(demoOneFactor)
latents <- c("G")
factorModel <- mxModel(
"One Factor",
type="RAM",
manifestVars = manifests,
latentVars = latents,
mxPath(from=latents, to=manifests,values=0.8),
mxPath(from=manifests, arrows=2,values=1),
mxPath(from=latents, arrows=2,
free=FALSE, values=1.0),
mxData(cov(demoOneFactor), type="cov",
numObs=500)
)
iniPars <- coef(factorModel)
print(iniPars)
pars2 <- imxJiggle(params=iniPars,lbounds=NA,ubounds=NA,dsn="runif",loc=1,scale=0.05)
print(pars2)
mod2 <- mxJiggle(model=factorModel,scale=0.05)
coef(mod2)
mod3 <- mxJiggle(model=factorModel,classic=TRUE)
coef(mod3)
# }
Run the code above in your browser using DataCamp Workspace