OpenMx Help

Hi guys!

I copied the factor analytic example from the OpenMx main website. When I run the model both the original and the resulting model with estimates will have the class MxRAMModel.
However, if I take the same model, add an arbitrary mean structure and mxData of type "raw", the resulting model will only have class "MxModel". To my understanding, it should still be of class MxRAMModel.
Is this a bug or a feature?
I attached my sample code.

cheers,
Andy

Hi All.

Does anyone know how to view the iteration history from the optimization routine? I generally like to examine this information as I consider model fit. I searched the archives and it looks like this issue has received some attention but I could not figure out how to obtain this output.

Thanks in advance,

Rick DeShon

I'm running a multiple-group model with 17 groups (using raw data). However, the number of observations that is reported is wrong. Instead of reporting the total N across all groups, it reports N X the number of groups. In other words, the full sample size across all groups is 20814, but summary reports numObs as 353838 (which is 20814*17). When I ignore the groups and run the model on the full sample, I get the correct number of observations (20814), but if I run it as 17 groups (with each group being a subset of the 20814), I get 353838 observations.

My question probably reflects my inexperience scripting in R, but I thought I'd give it a shot here.

fyi,

got the following error message after mxRun():

Error in runHelper(model, frontendStart, intervals, silent, suppressWarnings, :
BLAS/LAPACK routine 'DPOTRF' gave error code -4

to save y'all some time, turns out that i was running an old script that had an MxAlgebra object with the expression ".5 * VA" where VA was a matrix. to fix it, i just had to change this to ".5 %x% VA."

suggest you trap the DPOTRF error and return a more intelligible error message.

greg

I am running a series of univariate models before I combine my variables into a trivariate cholesky. I am running a full ACE model as well as nested AE, CE and E models and my AE is best fitting (not sure if this is important but I am only requesting CIs for my AE model since it is best-fitting). The following error messages consistently pop up under the output for AE. The output for the full ACE and other nested models does not show this error message.

Hi All:

I just got this strange error message when I ran the script below (trying to use a covariance matrix rather than raw data).

Error in runHelper(model, frontendStart, intervals, silent, suppressWarnings, :
BLAS/LAPACK routine 'DSYMM ' gave error code -9

Here is the script that generated the error:

require(OpenMx)
source("http://www.vipbg.vcu.edu/~vipbg/Tc24/GenEpiHelperFunctions.R")

mzCov <- matrix(c(1, .6, .6, 1), 2,
dimnames=list(c("y1", "y2"), c("y1", "y2")))
dzCov <- matrix(c(1, .4, .4, 1), 2,
dimnames=list(c("y1", "y2"), c("y1", "y2")))

Hi,

I met with the error when I was running a Random Coefficient model:

Error: Illegal label 'Rand[data.ID,1]' detected in matrix 'A'. Square brackets must contain numeric literals when used inside of labels.

So, can we use the definition variable inside of the square brackets? Or maybe there is something wrong with my script. The related script is as below, and the variable ID is numeric (1:100).

-----------------------------------------------------------------------------------------
RandVals <- matrix(c(-0.2, -0.2), numSubjects, 2, byrow=TRUE)
...
...
...

I get the following error messeage from openMX:

Error in mxModel("ACE", mxMatrix(type = "Iden", nrow = 4, ncol = 4, name = "I"), :
argument is missing, with no default

this "ACE" mxModel object contains only mxMatrix and MxAlgebra ellements.

It is one of six submodels in a bigger model this bigger model contains an mx Algebraobjective the other 5 submodels contain data and a mxFIMLObjective each for the different groups.

Which argument am i missing in this ACE mxModel? which are the required model arguments anyway?