weight <- (150 - 3) model1 <- mxModel(model1, mxMatrix(name="Va", type="Full", 1, 1, free=TRUE, value=.3), mxMatrix(name="Vc", type="Full", 1, 1, free=TRUE, value=0)) model1 <- mxModel("model1", mxModel("MZ", mxData(matrix(.7, 1, 1), type="cov", numObs=150), mxAlgebra(name="rPreMZ", expression = model1.Va + model1.Vc), mxAlgebra(name="zpreMZ", expression = .5 * log((1 + rPreMZ) / (1 - rpreMZ))))) model1 <- mxModel(model1, mxModel("DZ", mxData(matrix(.3, 1, 1), type="cov", numObs=150), mxAlgebra(name="rPreDZ", expression = .5*model1.Va + model1.Vc), mxAlgebra(name="zpreDZ", expression = .5 * log((1 + rPreDZ) / (1 - rpreDZ))))) model1 <- mxModel(model1, mxAlgebra(name="fvalMZ", expression = weight * (MZ@data@observed - MZ$zpreMZ) * (MZ@data@observed - MZ$zpreMZ))) model1 <- mxModel(model1, mxAlgebra(name="fvalDZ", expression = weight * (DZ@data@observed - DZ$zpreDZ)^2))