You are here

Saving warning massage code

3 posts / 0 new
Last post
suparman.y's picture
Offline
Joined: 08/19/2011 - 05:49
Saving warning massage code

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..

tbrick's picture
Offline
Joined: 07/31/2009 - 15:10
Status codes are stored in the fitted model in @output$status.

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:

...
Status <- modelFit3@output$status[[1]]
...
output3 <- c(k, A1, ..., M52, Status)

More details of the status field are available in this thread: http://openmx.psyc.virginia.edu/thread/326

suparman.y's picture
Offline
Joined: 08/19/2011 - 05:49
Thank you

Thank you very much...