Hi,
I am running Hermine's Cholesky script (Boulder "cdrom" 2016) only change I make is that I Use mxDataWLS to precompute the polychoric correlations and WLS weights:
DZdata <- mxDataWLS(DZ)
This runs without errors.
I then proceed to adjust the mxExpectation to reflect the nature of the data:
expMZ <- mxExpectationNormal( covariance="expCovMZ", means="meanG",thresholds="Tr", dimnames=selVars )
expDZ <- mxExpectationNormal( covariance="expCovDZ", means="meanG",thresholds="Tr", dimnames=selVars )
(Threshodls fixed at 0 and 1, so covariance and means can be freely estimated)
I proceed to specify the correct fit function:
funWLS <- mxFitFunctionWLS()
When I run the model I get the following error:
Running mulACEc with 185 parameters
Error in svd(X) : infinite or missing values in 'x'
Any suggestions?