OpenSEM Forums

Posted on
No user picture. Yuan Joined: 02/13/2015

Running CholAE

When I was running the CholAE model, it warns that:
In model 'CholAE' Optimizer returned a non-zero status code 1. The final iterate satisfies the optimality conditions to the accuracy requested, but the sequence of iterates has not yet converged. Optimizer was terminated because no further improvement could be made in the merit function (Mx status GREEN).

What does this mean?

Thank you

Posted on
No user picture. slukow Joined: 02/25/2015

Calculation of BIC discrepancy?

Hello - I'm having some trouble with discrepancies in the calculation of BIC in comparing output from Mx with OpenMx on both univariate and bivariate ACE models. I am currently using OpenMx 2.0.1 and R version 3.1.2. Below is the univariate script I used. It seems from the output that the problem may be coming from the fact that mxData objects are pulling number of observations from the total number of rows in mzData and dzData, which includes some rows in which both variables are missing.

Posted on
No user picture. Yuan Joined: 02/13/2015

bivariate genetic model

Hi,
I am new to Twin data analysis and would like to get help on the bivariate genetic analysis. I want to calculate the genetic correlations between two traits in twins. However, I only found the scripts about the univariate twin analysis in documentations. How can I get the bivariate ACE model in OpenMX?

Thank you

Posted on
No user picture. arjenc Joined: 02/17/2015

Factor analysis with weights

I have a question. I looked at the package sem and lavaan.survey for the factor analysis but sem doesn't support weights and lavaan.survey doesn't have FIML. Without FIML, lavaan.survey was still able to do the job of getting the expected graphs (smooth curves) but with OpenMx, it does not work. I get straight lines.

The code is below and the data set is in the attachment. I don't know if the weights are calculated with the data in the analysis.

dsetA <- read.table("dsetA.txt",sep="")
# Possible values of the weighting variable
valm <- seq(21,40,by=1/10)

Posted on
No user picture. GRBS Joined: 02/13/2015

Constrain correlation matrices

I have fit two saturated models and have created correlation matrices from the expected covariance matrices using the cov2cor function. I would like to constrain the correlation matrices to equality between these models without constraining the covariance matrices.

Posted on
No user picture. fife Joined: 07/01/2010

CIs when RMSEA = 0

Forums

Hi all,

I noticed that when RMSEA = 0, the CIs come out as NA. Why is that? Is it a bug? Or are CIs theoretically undefined when RMSEA=0?

Posted on
Picture of user. amullings09 Joined: 12/21/2014

Derivative Estimation Methods for 3 time occasions?

I am planning on using on using multivariate latent differential equation modeling for my MA thesis. I will have the opportunity to plan for the parameters, methods, and measurements required to do so. But currently, I am attempting to utilize one of these derivative estimation methods on some data I currently have. Specifically, I have 3 measurements of PTSD, Cortisol, and Coping Self-Efficacy.

Posted on
Picture of user. ana.martinovici Joined: 12/08/2014

Missing data - OpenMx vs Mplus

Hi,

The data set i use has 214 individuals for which I have different number of observations - varies between 21 and 30.

Posted on
No user picture. K Ram Joined: 08/21/2012

Correlated factors ok to use even when phentypic varaiance not significant?

We have 7 variables, of which one of the variable is significantly heritable - modeled using univariate analyses. Is it theoretically congruent to test for shared genetics between all 7 variables using the correlated factors model despite the fact that 6 of the variables are not significantly heritable?