Dear all,
I am doing a simulation. Nevertheless, some of its replications have a warning massage of NPSOL none-zero status code, such as: "In model 'Cross Lagged Model3' NPSOL returned a non-zero status code 6. ... "
Since it is important for analyzing the simulation results, how can this status code be written in the simulation output file?
Here is some part of my input list for producing the output file.
5.3 Processing
#
modelFit3 <- mxRun(model3)
sfit3 <- summary(modelFit3)
5.4 Saving Results
#
A1.3 <- sfit3[[1]]$Estimate[[1]]
A2.3 <- sfit3[[1]]$Estimate[[3]]
…
M52 <- sfit3[[1]]$Estimate[[62]]
output3 <- c(k, A1.3, A2.3, … ,M52)
write.table(t(output3), file="outfile3.dat", append=TRUE, sep=" ",
eol="\n", na="NA", dec=".", row.names=FALSE, col.names=FALSE)
Thanks..
Status codes are stored in the fitted model in @output$status.
To get code 6 or code 1 results, you can simply save the first element of the status code. For your example, the code would look something like:
More details of the status field are available in this thread: http://openmx.psyc.virginia.edu/thread/326
Thank you very much...