Hi guys,
I'm trying to fit a saturated model where the variable, 'manifests', includes all of the variables in the model. The non-saturated structural model runs fine, but I get an error when I fit the saturated model:
Error: The job for model 'Saturated Model2' exited abnormally with the error message: Expected covariance matrix is not positive-definite in data row 79 at iteration 0.
I have tried fiddling with the starting values, but to no avail. By the way, I know that all=TRUE is being deprecated, and I have written the same code in a longer form without all=TRUE, but this was shorter, so I thought it'd be easier to post. Below is the model:
saturated_model <- mxModel(paste("Saturated Model"),
type="RAM",
manifestVars=manifests,
#Variances and Covariances
mxPath(
from=manifests,
arrows=2,
free=TRUE,
all=TRUE,
values=.5,
#Means
mxPath(
from="one",
to=manifests,
arrows=1,
free=TRUE,
mxData(observed=modeldata, type="raw")
)
saturatedModel_run <- mxRun(saturated_model)
Any help would be greatly appreciated. Thanks!