OpenMx Structural Equation Modeling
P value of 1 in mxCompare
I am working on an analysis of a univariate sex limitation model for depression and I have run into a potential problem. When comparing the homogeneity ACE model to the AE model, I get a p-value of 1 in mxCompare. The fit of the models are very similar:
> mxCompare(HomACEModelFit,HomAEModelFit)
base comparison ep minus2LL df AIC diffLL diffdf p
1 HomACE
- Read more about P value of 1 in mxCompare
- 2 comments
- Log in or register to post comments
Constraining loadings on factor so average of loadings equals 1 and average of intercepts equals 0
oneFactorModel <- mxModel("CFA",
type="RAM",
Simplex model & sex differences
I am working on a simplex script in openMx. There are quantitative and qualitative sex differences for my phenotype, so I estimate separate paths for males and females and I would also like to freely estimate the dos correlations between the latent As. However, I am struggling with the transmission paths.
Probably, someone else has already dealt with this. I would like to know whether there are any simplex scripts available that take into account sex differences?
Regards
Charlotte
- Read more about Simplex model & sex differences
- 1 comment
- Log in or register to post comments
Multilevel SEM with complex sampling
Would you confirm if OpenMx supports multilevel SEM on complex sampling survey data, with categorical outcome? (each datapoint has different weight)
Thank you!
- Read more about Multilevel SEM with complex sampling
- 1 comment
- Log in or register to post comments
Returned matrices (expCov, expMean, A,S, etc) only reflect last row of data
Standardized estimates under equality constraints
I ran a linear latent growth curve model in OpenMx and constrained all residual variances to be equal. Not surprisingly, they all get to be the same unstandardized estimate. However, OpenMx reports only a single standardized estimate for the residual, which happens to be the one associated with the first observation. Even if the unstandardized estimates are constrained to be equal, the unstandardized ones need not be and in reality seldom are the same.
[code]
Adding new algebras in submodels and constraining them
I am running a bivariate moderation model and would like to set some constraints in order to test for some nonlinear effects. I am struggling though with adding new algebras into submodels (with subsequent equating them). Let's say that I specify
pathAm <- mxMatrix(name = "am", type = "Lower", nrow = nv, ncol = nv, free=T, labels = c("amM","amC","amU"), values=pathVal)
pathCm <- mxMatrix(name = "cm", type = "Lower", nrow = nv, ncol = nv, free=T, labels = c("cmM","cmC","cmU"), values=pathVal)
Unable to reproduce MASEM results from a published study
I'd like to reproduce the meta-analytic structural equation modeling (MASEM) results from this study:
soonang[dot]com/wp-content/uploads/2011/04/2007-MISQ-Ang1.pdf
I used the correlation matrix (Table 3, p. 559) as input and specified the paths according to Figure 2 (p. 560).
Additionally, I set the number of observations to 701 (p. 558).
The full openMx code is attached.
The openMx output for the parameter estimates fits the values in Figure 2 quite well.
However, the openMX fit statistics are quite different from the ones in the paper.
Pagination
- Previous page
- Page 5
- Next page