Third Party Software
tbates
Joined: 07/31/2009
How do I get a CI on a parameter derived from a model?
umxSummary reports the rA, rC, and rE for models. But this is done in the summary, and doesn't include CIs by default.
It’s computed as
solve(sqrt(I*A)) %*% A %*% solve(sqrt(I*A))
So you can add an mxAlgebra computing that into the model, then do mxSE on that
m2 = mxModel(m1,
mxAlgebra(name="ra", solve(sqrt(top.I*top.A)) %*% top.A %*% solve(sqrt(top.I*top.A)))
)
m2 = mxRun(m2)
mxSE(ra, m2)
That returns SEs, and CI= ± 1.96 * SE
best, t
noamm
Joined: 10/10/2018
Correlated factors solution for umxACE
Hi,
I'm using umxACE for a bivariate model. The current output is in the format of a cholesky model in addition to rG, rC, rE. Is there a way to get an output of the correlated factors solution? I know that I can easily calculate the estimates, but because they are not part of the model's output I didn't succeed in getting CI's for the calculated estimates (or the r's). Is there a way to get these (as these are more intuitive for presenting)?
Thanks for all the help!
- Read more about Correlated factors solution for umxACE
- 1 comment
- Log in or register to post comments
GregoryVerleysen
Joined: 02/26/2020
metaSEM questions
I am learning how to use metaSEM, but I am running into some issues.
Right now I only have 3 studies available to run the analysis (code with data is in attachment).
I am running into the following issues and was wondering if these would be caused by having only 3 studies
or am I doing something wrong here?
- Read more about metaSEM questions
- 1 comment
- Log in or register to post comments
noamm
Joined: 10/10/2018
umxSexLim - heritability doesn't match correlations
Hello,
ucjtcio
Joined: 01/31/2020
autofixtau2
Hi,
I've been using the rerun() function with the autofixtau2 argument, which is leading our models to converge. While I appreciate that it may be choosing suitable starting values for estimated parameters, I’m not entirely sure how the models have converged ( low number of studies and small N). In an effort to ‘fix’ tau2 could it end up running a fixed effects model? Are there conditions where autofixtau2 should/should not be used and how it should be interpreted?
- Read more about autofixtau2
- 5 comments
- Log in or register to post comments
pascofearon
Joined: 05/07/2013
Robert Balazsi
Joined: 10/24/2019
Many missing correlations
Dear Mike,
I am currently working on a metaSEM study. I've got some error, even I tried to use different methods (see below!). I also asked for a pattern of NA. I would like to ask if, in your opinion, the error messages are related to missing correlations? I have 379 lines, but only 1 or 2 correlations for some variables. Is there a way to continue this study? Maybe dropping some variables?
All the best,
Robert
- Read more about Many missing correlations
- 36 comments
- Log in or register to post comments
Leo
Joined: 01/09/2020
Obtain CIs on genetic correlations
I'm interested in obtaining the CIs of the genetic correlations using umx. Running umxCi prior to umxsummary does only provide the CIs for the different paths. Can anyone help me?
Thanks in advance.
- Read more about Obtain CIs on genetic correlations
- 1 comment
- Log in or register to post comments
paulT
Joined: 11/03/2017
default constraints in multi group model?
I am a new OpenMX|umx user and am trying to fit a multigroup model. I have fitted two separate umxRAM models (one for each group) and then tried to repeat using the 'umxSuperModel' command. I seem to get very different model results between the multi group model and the separate group models. In the single groups models, there are fewer significant paths in one group, but the multi group appears to have identical significant paths (only varying between group in path weight magnitude)?
- Read more about default constraints in multi group model?
- 3 comments
- Log in or register to post comments
pascofearon
Joined: 05/07/2013
mediation in OSMASEM
I was wondering whether you could point me to information about how to pass algebraic statements to OpenMx via the osmasem function so that I can test indirect effects – in tssem2 this is done easily, but it’s not so obvious with osmasem. The documentation for osmasem doesn't seem to go into detail about this. I suspect it could be done with the mxModel.args option, but as I don’t know my way around OpenMx all that well I had trouble figuring out how…
Thanks!
Pasco
- Read more about mediation in OSMASEM
- 5 comments
- Log in or register to post comments