I think this can be done, but I'm having problems seeing a clean solution. Is there a way with OpenMx to predict an outcome using the interaction of two latent variables?

The standard citation on latent variable interactions is a Psycometrika article by Klein & Moosbrugger (2000). Mplus implements some form of this, though I recall it does some variation on the EM algorithm. Could be wrong, though. I'd imagine implementation is possible in OpenMx, but no one has done it yet to my knowledge.

There are some general approaches to latent variable interactions which make the specification of the non-linear constraints much simpler (basically two matrix algebra expressions are required). I wrote about them here:

Neale, M.C. (1998) Modeling interaction and nonlinear effects with Mx: A general approach. In: G. Marcoulides & R. Schumacker Interaction and Non-linear Effects in Structural Equation Modeling pp. 43-61. New York: Lawrence Erlbaum Associates

Among many latent interaction methods have been and being developed, the Kenny-Judd (1984) method (also see Dr Neale's (1998) chapter) is available for OpenMX. Marsh, Wen, and Hau (2004) unconstrained approach helps us to have less non-linear constraint commands.

For the LMS approach, if OpenMX provides numerical integration, Monte Carlo integration, or Gauss-Hermite quadrature integration, it is possible to be implemented by OpenMX.

I think the methods with product indicators (methods mentioned by Heining Cham, and others) are outdated now. There are other methods available in other packages, like Mplus and EQS.

However, for me it is still unclear how the ML method (with all kinds of its variations) can be implemented in openMx. Is there someone who has done this? And is the dev. team planning to implement this method?

Furthermore, in "Structural Equation Modeling" 2010, pages 357- 373, Mooijaart and Bentler propose a method using higher order moments in stead of using only means and covariances. This method is compared to the method of Mplus and seems to be promising also in cases with nonnormally distributed predictors.

I'm not sure if analyzing higher order moments can be done easily in openMx.

In classic Mx, I had some success with implementing Gauss-Hermite quadrature using a mixture distribution (weighted with fixed weights). That should be possible with OpenMx also.

It is possible to specify user-defined fit functions using OpenMx. mxAlgebraObjective() would be the way to go. I recall Peter Bentler publishing a paper on higher moments early in the 80's (BJMSP if I recall correctly) and the matrix algebra was a bit heavy. However, I doubt whether it's anything that OpenMx can't handle.

Hey Pascal,

The standard citation on latent variable interactions is a Psycometrika article by Klein & Moosbrugger (2000). Mplus implements some form of this, though I recall it does some variation on the EM algorithm. Could be wrong, though. I'd imagine implementation is possible in OpenMx, but no one has done it yet to my knowledge.

There are some general approaches to latent variable interactions which make the specification of the non-linear constraints much simpler (basically two matrix algebra expressions are required). I wrote about them here:

Neale, M.C. (1998) Modeling interaction and nonlinear effects with Mx: A general approach. In: G. Marcoulides & R. Schumacker Interaction and Non-linear Effects in Structural Equation Modeling pp. 43-61. New York: Lawrence Erlbaum Associates

see:

http://www.vipbg.vcu.edu/vipbg/Articles/Chap3-ModelingInteraction-43.pdf

Thanks for the pointers! I'll look into both articles. Thanks!

Dear all,

Among many latent interaction methods have been and being developed, the Kenny-Judd (1984) method (also see Dr Neale's (1998) chapter) is available for OpenMX. Marsh, Wen, and Hau (2004) unconstrained approach helps us to have less non-linear constraint commands.

For the LMS approach, if OpenMX provides numerical integration, Monte Carlo integration, or Gauss-Hermite quadrature integration, it is possible to be implemented by OpenMX.

Heining Cham

I think the methods with product indicators (methods mentioned by Heining Cham, and others) are outdated now. There are other methods available in other packages, like Mplus and EQS.

However, for me it is still unclear how the ML method (with all kinds of its variations) can be implemented in openMx. Is there someone who has done this? And is the dev. team planning to implement this method?

Furthermore, in "Structural Equation Modeling" 2010, pages 357- 373, Mooijaart and Bentler propose a method using higher order moments in stead of using only means and covariances. This method is compared to the method of Mplus and seems to be promising also in cases with nonnormally distributed predictors.

I'm not sure if analyzing higher order moments can be done easily in openMx.

Ab Mooijaart

A couple of points:

In classic Mx, I had some success with implementing Gauss-Hermite quadrature using a mixture distribution (weighted with fixed weights). That should be possible with OpenMx also.

It is possible to specify user-defined fit functions using OpenMx. mxAlgebraObjective() would be the way to go. I recall Peter Bentler publishing a paper on higher moments early in the 80's (BJMSP if I recall correctly) and the matrix algebra was a bit heavy. However, I doubt whether it's anything that OpenMx can't handle.

HTH

Mike