OpenMx Error Messages

I realize I am asking for instability using development versions, but thought I would mention that running an OpenMx model kills R.

Here is my system info:

R Under development (unstable) (2012-02-03 r58258)
Platform: x86_64-pc-mingw32/x64 (64-bit)

locale:
[1] LC_COLLATE=English_United States.1252
[2] LC_CTYPE=English_United States.1252
[3] LC_MONETARY=English_United States.1252
[4] LC_NUMERIC=C
[5] LC_TIME=English_United States.1252

attached base packages:
[1] stats graphics grDevices utils datasets methods base

other attached packages:

I have just got this error message. What I have done is

1. define a model say, A, which includes sub-models.
2. define A0<-mxModel(A, name="A0")
3. add a suffix "0" to all parameters in A0, so that A and A0 are exactly the same but with independent parameters
4. define a bigger model M<-mxModel("M", A, A0, some constraints and objective function over A and A0)

Hi all,

I ran a relatively simple path analytic model with four manifest variables, three causal paths, and one unanalyzed relationship. The model has 2 df. I used the code below but got the following error message:

"In model 'Apgar Model 2' NPSOL returned a non-zero status code 1. The final iterate satisfies the optimality conditions to the accuracy requested, but the sequence of iterates has not yet converged. NPSOL was terminated because no further improvement could be made in the merit function (Mx status GREEN). "

Comment from student on getting a helpful error in OpenMX of the form "I'm sorry I can't..."

"Gee: That's the first time I've ever seen R talk in the first person: That's nice!"

:-)

Getting this error
Error: The observed correlation matrix is not 1's along the diagonal

from the statement
mxData(cor, type="cor", numObs=1656)

but looking at the data matrix indicates that its diagonal consists of 1's

> diag(cor)
SOF IWF PWF FSC FSF EF GPA
1 1 1 1 1 1 1

cor
SOF IWF PWF FSC FSF EF GPA
SOF 1.00000000 -0.310734680 0.02337843 -0.01050081 0.324939442 -0.49945847 0.83664278
IWF -0.31073468 1.000000000 0.09947420 -0.03052640 0.001136279 0.49532944 -0.10899843

Hi all, I'm kind of a noob so please forgive me if this is terribly rudimentary.

I get the error message:

Error in runHelper(model, frontendStart, intervals, silent, suppressWarnings, :
BLAS/LAPACK routine 'DPOTRF' gave error code -4

Traceback yields:

3: .Call("callNPSOL", objective, startVals, constraints, matrices,
parameters, algebras, data, intervalList, communication,
options, state, PACKAGE = "OpenMx")
2: runHelper(model, frontendStart, intervals, silent, suppressWarnings,
unsafe, checkpoint, useSocket, onlyFrontend, useOptimizer)

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.

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

I am adapting the attached bivariate script to a 4-variable version. I have expanded the X, Y, and Z matrices to to 4-variables, changed Y to accomodate and AE model. After the follow part of the script:

#### DZ Model Stuff ####
mxModel("DZ",
#Pull out data and covariates for DZs
mxData(data.frame(dzData,dzDefs), type="raw" ),
mxMatrix( type="Full", nrow=1, ncol=ndef, free=FALSE, labels=c("data.age_T1","data.SITE_T1"), name="D1"), #twin1 (N)

There seems to be two issues in using omxSetParameters to reset bounds. Please see below.

rm(list=ls());
require(OpenMx)

A<-mxMatrix(type="Full",nrow=1,ncol=1,free=T,values=.8,label="a",lbound=0,name="A");
model<-mxModel("model",A,mxAlgebra(A^2,name="A2"),mxAlgebraObjective("A2"))

###Now I would like to remove the lower bound

model0<-omxSetParameters(model,"a",lbound=NA)

### This gives an error: " 'lbound' argument must either be NA or a numeric vector. "
### But I was using NA exactly.

model0<-omxSetParameters(model,"a",lbound=NaN)
model0@matrices$A