OpenMx Help

Hi,

I'm wanting to add 3 covariates to a saturated model. One of them (age) seems to work fine, but when I try to use either of the other two individually (birthweight or gestation time), I get an error message:

> univTwinSatFit <- mxRun(univTwinSatModel)
Running univTwinSat
Error in mxRun(univTwinSatModel) :
Error NYI: Missing Definition Vars Not Yet Implemented.

I am running a trivariate model using a script from the workshop this year (thanks for that!) but keep getting the same error message. It comes when I try to run the Fit Multivariate ACE Model using Cholesky Decomposition portion of the script. The error message is as follows:

multACEFit <- mxRun(multACEModel)
Running multACE
Error in gzfile(file, mode) : cannot open the connection
In addition: Warning message:
In gzfile(file, mode) :

Hi all,

I appreciate the speedy help that comes with posting on the forums here.

I have two separate problems so thought I'd post them both here:

PROBLEM 1: I'm trying to fit a model using var/cov matrices rather than raw data, as suggested by Steve in an earlier post (in order to get rough ideas of good start values for later FIML fitting). However, I get the following error message (pasted before it is the last bit of the script so you can see how I set up the Objective functions

+ ),
+ mxModel("DZOS",

Hi everybody,
I have a simple question that might be a bit stupid. But I can't find an equivalent to the 'as.matrix' command in R. I understand that 'rvectorize' is the equivalent to 'as.vector' so I'm surprised that I can't find the 'as.matrix' thing anywhere.
Or is it maybe not yet implemented in openmx?

Thanks for your help

Greetings

Jan

Hello everyone,
i am trying to generate a model with several mxAlgebra objects in it. The thing is that they all have the same structure (except for one value), but must be assigned different names. Right now I'm specifying them seperately, so:

mxAlgebra(EVEC%*%(exp(1%x%EVA)-tempa)%*%solve(EVEC), name = "EXPd1"),
mxAlgebra(EVEC%*%(exp(2%x%EVA)-tempa)%*%solve(EVEC), name = "EXPd2"),
mxAlgebra(EVEC%*%(exp(3%x%EVA)-tempa)%*%solve(EVEC), name = "EXPd3"),
mxAlgebra(EVEC%*%(exp(4%x%EVA)-tempa)%*%solve(EVEC), name = "EXPd4"),

Dear all,

I am a new learner of OpenMX and R. I am interested to know if there is alternative estimation (robust ML such as Satorra-Bentler test statisitc, ADF/WLS by Browne) besides normal-theory ML?

Thank you!

Heining Cham

Hi all,

This probably won't be easy to answer given how vague the question is, but...

question1：
when I simulated the model, I found the -2LL value of some nested models (which had less estimate parameters than full model) were smaller than that of full model, without any "warning" in script. I think it is abnormal, and I have no idea why it happened, and could you give me some information or suggestion to solve this problem.
Thank you!