faulty warning message?

When I load the workspace, there is no cccFit object. Evaluating the cccModel object yields NA.

mxEval(objective, cccModel, compute=TRUE)


Similarly for each of the submodels objectives.

> mxEval(MZf.objective, cccModel, compute=TRUE)
[,1]
[1,] NA
> mxEval(DZf.objective, cccModel, compute=TRUE)
[,1]
[1,] NA
> mxEval(MZm.objective, cccModel, compute=TRUE)
[,1]
[1,] NA
> mxEval(DZm.objective, cccModel, compute=TRUE)
[,1]
[1,] NA
> mxEval(DZo.objective, cccModel, compute=TRUE)
[,1]
[1,] NA


cccFit could be left over from a different model run where the objective was well not NA.

Ah yes, I think one is supposed to do this:

cccFit <- mxRun(cccModel, unsafe=T)

in order to generate the error code and to make the mxEval's work.

I haven't been able to figure out much about this issue. I do know that it is not a faulty warning message, rather it is a faulty error.

Making some changes only leads to more mysteries. The m2LL algebra seems to compute okay, and the model runs fine without an objective function.

dddModel <- mxModel("Container",
cccModel[["MZf"]],
cccModel[["MZm"]],
cccModel[["DZm"]],
cccModel[["DZf"]],
cccModel[["DZo"]],
mxAlgebra(MZf.objective + MZm.objective + DZm.objective + DZf.objective + DZo.objective, name="m2LL")
)
dddRun <- mxRun(dddModel)
dddRun$m2LL@result


Moreover, replacing the algebra with a simpler one and using an algebra objective gives the same error as before.

eeeModel <- mxModel("Container",
cccModel[["MZf"]],
cccModel[["MZm"]],
cccModel[["DZm"]],
cccModel[["DZf"]],
cccModel[["DZo"]],
mxAlgebra(DZf.objective, name="m2LL"),
mxAlgebraObjective("m2LL")
)
eeeRun <- mxRun(eeeModel)


Could this be a corner case of badness for the algebra objective when there are submodels, an algebra objective, and no free parameters? Do the same errors persist when there is a free parameter?

Not with this example:

testMZf <- cccModel[["MZf"]]
testMZf$af@free[1,1] <- TRUE
testDZf <- cccModel[["DZf"]]
testDZf$af@free[1,1] <- TRUE
testDZo <- cccModel[["DZo"]]
testDZo$af@free[1,1] <- TRUE
fffModel <- mxModel("Container",
testMZf,
cccModel[["MZm"]],
cccModel[["DZm"]],
testDZf,
testDZo,
mxAlgebra(MZf.objective + MZm.objective + DZm.objective + DZf.objective + DZo.objective, name="m2LL"),
mxAlgebraObjective("m2LL")
)
fffRun <- mxRun(fffModel) #give warning
fffRun$m2LL@result


Thoughts?

I bet this is fixed in SVN trunk. Can somebody test it again?

If it's still broken, can we have an R script to reproduce the problem? Thanks.

Current thought is this was related to

http://openmx.psyc.virginia.edu/thread/3871

where setting independent=TRUE fixed the problem. The code in the linked thread runs fine unmodified in the r3750 Beta, so I'm closing this.