run:
MZf <- mxModel(name="MZf", mxData(ObsCovMZf, type="cov", numObs=nMZf),
mxMLObjective(covariance='myModel.PreCovMZf'))
DZf <- mxModel(name="DZf", mxData(ObsCovDZf, type="cov", numObs=nDZf),
mxMLObjective(covariance='myModel.PreCovDZf'))
myModel1 <- mxModel(name="myModel", MZf, DZf)
myModel1
then run and compare:
MZf <- mxModel(name="MZf", mxData(ObsCovMZf, type="cov", numObs=nMZf),
mxMLObjective(covariance='myModel.PreCovMZf'))
DZf <- mxModel(name="DZf", mxData(ObsCovDZf, type="cov", numObs=nDZf),
mxMLObjective(covariance='myModel.PreCovDZf'))
myModel2 <- mxModel(name="myModel")
myModel2 <- mxModel(myModel2, MZf, DZf)
myModel2
myModel1 has only DZf as a submodel while myModel2 has both MZf and DZf as submodels. is this intentional? or should documentation or the MxMOdel constructor be changed?