Saving warning massage code

2 replies [Last post]
suparman.y's picture
Offline
Joined: 08/19/2011

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
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
Thank you

Thank you very much...