Hi,
I'm wanting to add 3 covariates to a saturated model. One of them (age) seems to work fine, but when I try to use either of the other two individually (birthweight or gestation time), I get an error message:
> univTwinSatFit <- mxRun(univTwinSatModel)
Running univTwinSat
Error in mxRun(univTwinSatModel) :
Error NYI: Missing Definition Vars Not Yet Implemented.
I have full data for age, but a lot of missing data with the other two variables. I just want to confirm what the error message says - does OpenMx choke if your definition variable contains NA's? Will this be resolved anytime soon?
The definition variable script format is as follows:
univTwinSatModel <- mxModel("univTwinSat",
mxModel("MZ",
mxMatrix( type="Lower", nrow=ntv, ncol=ntv, free=TRUE, values=1, name="CholMZ" ),
mxAlgebra( expression=CholMZ %*% t(CholMZ), name="expCovMZ" ),
mxData( observed=mzData, type="raw" ),
# Adjust for age
mxMatrix( type="Full", nrow=1, ncol=nv, free=TRUE, values=0, name="Mean" ),
mxMatrix( type="Full", nrow=1, ncol=nv, free=TRUE, values=0, labels=c("beta1","beta2"), name="b" ),
mxMatrix( type="Full", nrow=1, ncol=1, free=FALSE, labels=c("data.Age_1"), name="A1"), #twin1
mxMatrix( type="Full", nrow=1, ncol=1, free=FALSE, labels=c("data.Age_2"), name="A2"), #twin2
mxAlgebra( expression= A1 %x% MZ.b, name="A1R"),
mxAlgebra( expression= A2 %x% MZ.b, name="A2R"),
mxAlgebra( expression= cbind((MZ.Mean + A1R),(MZ.Mean + A2R)), name="expMeanMZ"),
mxFIMLObjective( covariance="expCovMZ", means="expMeanMZ", dimnames=selVars),
...
Thanks,
Paul