OpenMx Help

Hi.
Is it possible to extract values of upper and lower bounds of CI separately? Let's say I would like to estimate CI's for path estimates for males by introducing the following command:

ciVm = mxCI(c('am','cm','em'))


The only thing I am able to get is:

> QualACEFit@output$confidenceIntervals
lbound ubound
MZM.am[1,1] 6.128870 10.699917
MZM.cm[1,1] -7.901978 7.901978
MZM.em[1,1] 3.775610 5.115599


What I am interested in is to get, e.g., a value for the lower boud of CI for am. Is that feasible somehow?

Dear all,
I have data from a survey, e.g, several Items answered by pupils within several schools. And further I want to confirm that there are 4 latent with assumed loadings on the Items (CFA assuming a normal distribution). Is it possible to analyze this multilevel structure (school/pupil) with OpenMx (some way to get a Mplus-like handling 'Type=complex cluster=school')?
Thanks in advance,
René

First I ran this model using my own data, and met the error message:

Running Growth Mixture Model
Error in substituteOperators(as.list(retval), algebra@name) :
cannot allocate memory block of size 2.9 Gb
In addition: Warning message:
In substituteOperators(as.list(retval), algebra@name) :
restarting interrupted promise evaluation

So I just tried to run the example data and example R code "GrowthMixtureModel_PathRaw.R", but still the same error.
Anyone knows what happened here?

Thank you.

HI,
I am using this script to run a multivariate-Cholesky AE model

multAEModel <- mxModel(choleskyACEFit, name="multAE",
mxModel(choleskyACEFit$ACE,
mxMatrix( type="Lower", nrow=nv, ncol=nv, free=F, values=0, name="c" ),
mxCI(c("ACE.sta","ACE.stc","ACE.ste")))
)

This is the error message I get:

Error: The following error occurred while evaluating the subexpression 'solve(sqrt(AE.I * AE.C))' during the evaluation of 'AE.Rc' in model 'univAE' : Lapack routine dgesv: system is exactly singular

Sorry for the beginner's question...
But please help me understand: When I specify a path like

mxPath(
from=latents, # names of latent variables
arrows=2,
connect='single',
free=T,
values=NA,
labels=('lab1','lab2','lab3')
)

and run a model (RAM) with this and other path for the indicators with mxRun(), are 'lab1', 'lab2', 'lab3' the estimates of variances of the latent variables or are they estimates of variances of the latent disturbance?

Thanks in advance!!

Hello, I have a problem when reading the data, for example, when I type "data(mv.csv)", there always returns a warning message: data set 'mv.csv' not found. Even when I add the document path "data('/Users/Documents/mv.csv'), the data are also not be able to find. Actually, this problem is not all about OpenMx, it relates to the whole R system. But I can only think of here to ask this question. Sorry, I forgot to say my computer system, it is MAC OS X 10.8.
Your help would be appreciated.
Ting

Hi everyone,

I have a problem running multilevel model using a SEM framework to get a likelihood-based CI of the effect size. I have two factors representing (random) intercept and (fixed) slope. I attach my data file here. Use the following code to get the dataset:

datawide <- read.csv("datawide.csv")

The dataset has been transformed into a wide format. Y is the response variable. X is the upper-level predictor. Z is the lower-level predictor without any centering. My target model is

L1: Y = B_0 + B_1 * Z + r
L2: B_0 = X + u

Hello again :)

I need some assistance:
When running my model (factor analysis with four indicators) with mxRun() I get a negative residual variance. So I used lbound = 0 for the estimate of this variance. The output gives me 0 for the estimate. I wonder whether OpenMx really estimates there or does it try to estimate in the direction of the old, negative estimate and stop at the lower bound of 0.

If this is not detailed enough I could provide my script and full output.

Thanks in advance!
Sophia

I would like to be able to handle Poisson counts. The only lead to doing this in OpenMx seems to involve using mxRowObjective. I was wondering if any one could share some examples.