OpenMx General Help

Both models have been run, and I can access a models number of estimated parameters from summary(model).

summary() however, can be costly in terms of execution time: Is there a way to access the EP (or df) property from the model directly without calling summary?

Besten Wünsche,
tim

FYI, the situation I am trying to automate is avoiding this:

|Model        | EP|Delta-2LL |Delta df|p  |      AIC|Compare with Model |

I have been running a model within a loop over several variables, which does not converge well for all variables. Therefore, I applied mxTryHard. There are two things I noticed:

1) for models that run well without calling mxTryHard the best fit does not seem to change using mxTryHard, but the status code changes from 0 to -1

> fit1 <- mxRun(model,intervals =F)
> fit1$output$Minus2LogLikelihood
[1] 4595.342
> fit1$output$status
$code
[1] 0
$status
[1] 0

fit2<-mxTryHard(model,intervals=F)
> fit2$output$Minus2LogLikelihood
[1] 4595.342
fit2$output$status
$code

I'm using a Macbook Air, and I can't work out how to modify the path arrows so they look neat. Can anyone help?

There is nothing in the documentation. The models in the examples all look nice, whereas everytime I add a path or covariance path, the program sticks them in anywhere and I can't move the end-points around to different parts of the observed or latent variable.

Please help!

Jay.

> x$matrices$A$labels
v_1 v_2 v_3 mu
v_1 NA NA NA "b_1"
v_2 NA NA NA "b_2"
v_3 NA NA NA "b_3"
mu NA NA NA NA

I would like to change "b_3" to "b_2":

> x$matrices$A$labels[3,4] <- "b_2"
Error: Right-hand side of assignment operator has illegal value FullMatrix, SymmMatrix, FullMatrix, and FullMatrix

Why does this not work?

This is actually quite a quick question and the title says most of it already.
Anyway, is there anyway to constrain mxAlgebra objects across multiple submodels?

I would like to constrain estimated heritability of two separate twin models to be the same while allowing each submodel to have separate parameter estimations.

Unfortunately I was unable to find a way to do so. I have tried mxConstraint(modelA.algebra1==modelA.algebra2) but that seemed not to work.

I would be thankful for any reply!