OpenMx Error Messages

Posted on
Picture of user. suzannejak Joined: 01/06/2010

Error: The observed covariance matrix is not a symmetric matrix

Hello all,

I am calculating the observed covariance matrix from the correlation matrix and the standard deviations. However, the resulting matrix apparently is not symmetric anymore (although judged by me it is...). Leading to the error message: The observed covariance matrix is not a symmetric matrix.

I hope somebody can help me out, thanks in advance!

Here is the code:

X = matrix(0,8,8)
val = c(1,.39,.35,.21,.32,.40,.39,.39,
1,.67,.11,.27,.29,.32,.29,1,.16,.29,
.28,.30,.37,1,.38,.30,.31,.42, 1,.47,
.42,.58,1,.41,.51,1,.42,1)
X[lower.tri(X, diag=TRUE)] = val

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

Warning: "Display list redraw incomplete"

what does the warning "Display list redraw incomplete" mean?

In this case it was followed by
"2: In model 'Full G by E' NPSOL returned a non-zero status code 6. The model does not satisfy the first-order optimality conditions to the required accuracy, and no improved point for the merit function could be found during the final linesearch (Mx status RED)", which I understand :-)

Posted on
No user picture. JoopDirksen Joined: 06/28/2011

Expected covariance matrix is not positive-definite in data row ... at iteration ...

Hello,

I'm developing a model, but I run into a couple of problems. Since the OpenMx documentation did not give the answer I was looking for, I want to try it here. I'm having a hard time getting started with SEM and with OpenMx because there is no real good overview of SEM methodology in my opinion.

I'm trying to run the following code:

library(OpenMx)

model <- mxModel(
name="Model",
type="RAM",
manifestVars = c("A", "B", "C", "D", "E","F"),
latentVars = c("var1", "var2", "var3", "var4"),
mxPath(from="var1",to=c("A", "B", "C")),

Posted on
No user picture. mdiaspar Joined: 05/26/2011

BASIC Question about running ACE model

Hello!,
I have prepared a script to compute a factor model (path diagram for Phillips and Fulker model ) of phenotypes of parents and twin children given by Neale.

I get an error message indicating the following sentence:
+ mxAlgebra( expression= (A%*%(S%*%t(C)))+(C%*%(t(S)%*%t(A))), name="J" ),

# MZ Twin
+ mxAlgebra( expression= (A%&%G%)+(C%&%R%)+J+(N%*%t(N)), name="U1" ),

Error: unexpected '&' in:
"mxAlgebra( expression= (A%*%(S%*%t(C)))+(C%*%(t(S)%*%t(A))), name="J" ),
mxAlgebra( expression= (A%&%G%)+(C%&"

Posted on
No user picture. a.ortega Joined: 02/05/2011

Error: "length of 'dimnames' [1] not equal to array extent"

Hi!

I have prepared a script to compute a linear growth model to estimate genetic and environmental influences (ACE) in an intake and a slope with ordinal twin data: 4 variables, 3 categories (2 thresholds) in each of them. The sample is is divided in 6 groups (and so is the script), according to zygosity and sex: MZ men, DZ men, MZ women, DZ women, DZ man-woman and DZ woman-man.

I get an error message indicating the following sentence:
> Error: The algebra 'ACE.expThre' in model 'lgcOrdACE' generated the error message: length of 'dimnames' [1] not equal to array extent

Posted on
No user picture. abhworthington Joined: 06/09/2010

Bivariate ACE model fit causes R front-end shutdown

Hello all,
I'm experiencing a somewhat frustrating error in an attempt to run some bivariate Cholesky ACE models using ordinal data. Rather than use raw data inputs (there are different numbers of thresholds across variables, so coding is tricky), I have taken the "easy" way and simply created polychoric correlation matrices for the variables of interest. However, when I attempt to fit the model, R crashes.

Posted on
No user picture. Mariano Joined: 05/21/2011

Error in readRDS ()

Hi,
I'm a beginner user of OpenMx.
I need help for the following error. It appear after any script, such as univariate regression, regression by path analysis by means of covariance or raw data calls, etc.
I've played with several forms, and check package loadings, including fields, psych, OpenMx, and general packages for R.
Please, if you can give me any tip, I'll be happy.

Thanks,

Mariano Ordano
lillo.org.ar
conicet.gov.ar

Posted on
No user picture. a.ortega Joined: 02/05/2011

Error: "non-conformable arguments"

Hello everyone!
I am trying to fit a latent (linear) growth model (with genetic and environmental influences, ACE) using ordinal twin data (MZ and same-sex + opposite-sex DZ; 6 groups), which includes 3 variables, 7 categories in each of them. The definition of the thresholds (invariant thresholds) seems to be difficult and the execution of the script returns:

"Error: The algebra 'ACE.ThInc' in model 'lgcOrdACE' generated the error message: non-conformable arguments"

Posted on
No user picture. a.ortega Joined: 02/05/2011

Error: NPSOL returned a non-zero status code 3

Hello everyone!
I have prepared a model to estimate genetic and environmental influences (ACE) using a cholesky decomposition with ordinal twin data: 3 phenotypes, 7 categories in each of them. The sample is is divided in 6 groups (and so is the script), according to zygosity and sex: MZ men, DZ men, MZ women, DZ women, DZ man-woman and DZ woman-man. The script including this multivariate threshold model takes about an one and a half hour to be computed, and returns the following error:

Posted on
No user picture. garverapgar Joined: 05/05/2011

Error in substituteOperators(as.list(retval)) : internal error -3 in R_decompress1

Hello all,

I've been running models successfully, but today I started getting the above error message. My code didn't change, but I may have accidentally installed a new version of OpenMx before running my model. Could someone tell me what this error message is about? If there isn't an easy fix, is there a way I can go back to an older version?

Thank you!
Christine