Hello,
I'm developing a model, but I run into a couple of problems. Since the OpenMx documentation did not give the answer I was looking for, I want to try it here. I'm having a hard time getting started with SEM and with OpenMx because there is no real good overview of SEM methodology in my opinion.
I'm trying to run the following code:
library(OpenMx)
model <- mxModel(
name="Model",
type="RAM",
manifestVars = c("A", "B", "C", "D", "E","F"),
latentVars = c("var1", "var2", "var3", "var4"),
mxPath(from="var1",to=c("A", "B", "C")),
mxPath(from="var2",to=c("B", "D", "E")),
mxPath(from="var3",to=c("F", "E")),
mxPath(from="var4",to=c("C", "D")),
mxPath(from="one",to=c("A", "B", "C", "D", "E","F", "var1", "var2", "var3", "var4")),
mxData(observed=data,type="raw")
)
fit <- mxRun(model)
summary(fit)
I'm getting the error message "Expected covariance matrix is not positive-definite in data row 126 at iteration 0"
Can anybody tell me what this error message means? I read on the errors page that it means that you should change your "starting values", but where should I do that? And besides that, what are "starting values" anyway? Somewhere it also mentions "free variables", but I cannot find a clear explanation what a "free variable" is exactly.
Who can give me the exact code I need to add to modify these starting values? I also read somewhere that it has to do with matrices for which no inverse can be calculated and that a Cholesky decomposition could help. Is anybody familiar with this method? And if so, what exactly do I need to change in the code above to make this work?
Thanks,
Steven