OpenMx Error Messages

Hi all,

After attending the Workshop in March, we now tried to work with OpenMx for the first time without help... Not completely succesful yet. We only mildly adjusted a script provided at the workshop but it gives an error message: 'The observed data associated with the FIML objective in model 'MZ' does not contain dimnames'. Please find the script below, it seems to us that we did define dimnames (dimnames=selVars). What do we do wrong?

Many thanks in advance for your efforts.

Cathy & Marcel

require(OpenMx)
require(psych)

I'm not sure if this is the right place for this. I'm getting an error message, but I'm fairly confidant that the error is due to the current capabilities of OpenMx, and not to user error.

Given some mxMatrix P and the definition variables foo and bar, the following works:
mxAlgebra(P[data.foo, ], name='Work')

But the following breaks:
mxAlgebra(P[data.foo, data.bar], name='Break')

Here are some simplified test cases:

#-----------------------------
# Works

xdat <- data.frame(foo=rep(1, 10), bar=rep(2, 10))
wmod <- mxModel(
name = 'WorkingModel',

Is there somewhere where error messages are listed, with an interpretation of what they mean and advice about what you might do about the error?

I get
"In model 'multiCholACE' 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). "

Hi Everyone,

I've been running basic, univariate twin models over the last week and have had no troubles, but today (March 14) I began receiving the following error in response to this code. The code is nearly identical to the code used at the Twin Workshop, with the exception of different variables and had worked fine up until now even with the new variables. Any ideas of what might be causing this (it seems like a technical issue rather than a modeling issue, but I may be wrong)?

Thanks,
Jarrod

+ mxAlgebra( MZ.objective + DZ.objective, name="-2sumll" ),

Greetings all,

I finally was able to get OpenMx going on my Mac. I am reading in some data and estimating a very small path model, just to get the hang of things. The syntax is

sciach <- read.table("~/Desktop/sciach.txt",header=TRUE)

require(OpenMx)

#data(sciach)

manifests <- names(sciach)

sciachsimple <- sciach[,c("irtsci", "grades", "probsolv")]

pathmodel <- mxModel("path model",type="RAM",
mxData(observed=sciachsimple,type="raw"),
manifestVars = manifests,

mxPath(from="probsolv", to="grades",arrows=1,free=TRUE),

Hello All,
Based on the generic script (below) from the OpenMx webpage I'm testing a very simple sub-model in which I drop one element of the 'S' matrix to zero and then compare back to the full model. I get an error message which says "Error in parameters[[i]] : subscript out of bounds". Does anyone have an idea as to why this might be happening?
Cheers,
Nathan

---

require(OpenMx)
data(demoOneFactor)
manifests <- names(demoOneFactor)
latents <- c("G")
factorModel <- mxModel("One Factor", type="RAM",
manifestVars = manifests,
latentVars = latents,

I'm trying to fit a very simple model using both the covariance model and the means model. I always get an error message concerning dimnames whether I specify the model via paths or using matrices. If I use matrices, it complains there are no dimnames for the F matrix - yet if I add dimnames for the F matrix it still complains.

All I want to be able to do is add a means model to the following OpenMx script that only has a covariance model:

> one <- mxModel("me",
+ type="RAM",
+ manifestVars = manifests,
+ latentVars = latents,

These errors occur when the mxFIMLObjective has no dimnames for the covariance or means, and cannot find them elsewhere (i.e., in the target algebra). For example:

    model <- mxModel("univSat4", 
        mxMatrix(type="Symm", nrow=1, ncol=1, free=T, values=1, name="expCov"),
        mxMatrix(type="Full", nrow=1, ncol=1, free=T, values=0, name="expMean"),
        mxData(observed=testData, type="raw"),
        mxFIMLObjective(covariance="expCov", means="expMean"))
    fit <- mxRun(model)

Suggested solutions:

Elements of model@output$status

Each of the output codes should report an appropriate error message at the end of evaluation of the mxRun() statement.

@output$status[[1]]

This status element is returned from the optimizer (NPSOL).

A value of 0 means a successful optimization--no error returned.