OpenMx General Help

Posted on
Picture of user. neale Joined: 07/31/2009

Problem with making a vector/list of fitted mxModel objects

So I have this function which mxRuns an mxModel, and returns the fitted mxModel object very nicely. I can access the slots of this object just fine. However, I'd like to call this function a number of times in a loop (or using apply variant). I do not seem to be able to build a list of mxModel objects from which I can extract elements.

firstrun <- runACEModelContinuous(TwinPairDTIDemog,varname='FLFAminus58',t1suffix='Twin1',t2suffix='Twin2',idzyg='ZYGOSITYTwin1',extractionMZ=.(ZYGOSITYTwin1=='M'),extractionDZ=.(ZYGOSITYTwin1=='D'))

# that works fine, and I can see, e.g.,

Posted on
No user picture. pgseye Joined: 10/13/2009

Bivariate outliers

Hi,

I've been advised I need to look at bivariate outliers in my twin dataset and that it used to be done with %p in old Mx. Is there an equivalent function in OpenMx?

Thanks,

Paul

Posted on
Picture of user. smedland Joined: 08/04/2009

\stnd cov2cor?

Hi
I'm looking for the OpenMx equivilent to the mx \stnd function. Searching the page the cov2cor() function within mxAlgebra pops up a couple of times but looks to have been discontinued(?) in the more recent releases - any similar function currently availble or implemented in the near future?
thanks
Sarah

Posted on
No user picture. svrieze Joined: 01/03/2010

\m2v operator?

Hello,

What would be the R or OpenMx equivalent of the \m2v operator in Classic Mx?

Related question: if I define in R my own function that performs this operation, can I use that function within the mxAlgebra command?

-Scott

Posted on
Picture of user. mkeller Joined: 08/04/2009

error message in trying to fit a var/covar matrix

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",

Posted on
No user picture. Janosch Joined: 05/10/2010

"as.matrix" in openmx

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

Posted on
No user picture. Janosch Joined: 05/10/2010

naming mxAlgebra objects automatically

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"),

Posted on
No user picture. jam86hk Joined: 05/21/2010

Robust ML estimation in OpenMX

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

Posted on
Picture of user. mkeller Joined: 08/04/2009
Posted on
No user picture. mailyn Joined: 05/19/2010

why is -2LL value of nested model smaller than that of full model, and how to find the best start value

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!