OpenSEM Forums
Julia
Joined: 03/29/2012
qingwen
Joined: 08/15/2019
CI is confusing
I fit a nuclear twin family model. AE model is the best. However, the CI of e makes me confused. Here are the output in free parameters and CI details section. The CI details looks strange, ranging from 0 to 2.053589e+06.
- Read more about CI is confusing
- 9 comments
- Log in or register to post comments
liagym1
Joined: 11/06/2019
ACE with moderator-help interpreting results
I am using regular Mx (unfortunatly i cannot find an active forum), and would like to clarify the meanings of the output received after running in an ACE with an ordinal moderator script.
1. Are the estimates of T,U,V the moderated components? (often denoted by betas)
2. If the estimate for a certain parameter is not zero, yet the CI contains zero- does this imply that this component is not significant?
- Read more about ACE with moderator-help interpreting results
- 4 comments
- Log in or register to post comments
liagym1
Joined: 11/06/2019
Error in output
I am trying to run an ACE with moderator script, and this is the error I received in the output. What changes should i make in my data?
Error found in rectangular data; record number 1 :
1004 -999 -999 -999 -999 -999
We have a problem whose error code is 176
and which I ran across at line number 38
of your input script
RE FILE=AGE11DATA.DAT
~~~
An error has occurred while reading a rectangular file
Make sure that you *don't* have a FORMAT at the beginning
- Read more about Error in output
- 7 comments
- Log in or register to post comments
diana
Joined: 04/17/2019
how to avoid NA in confidence interval
Hello, I often got NA in confidence interval when fitting univariate or bivariate SEM and changing optimizers didn't help. Here is the code of a bivariate ACE model, and the output when I use 'SLSQP' optimizer . I wonder how can I aviod NA when estimating CI of parameters and algebra. Thank you so much!
vars <-c('tea_5g','smk_num_4g')
covars <- c("age","region_type")
nv <- 2 # number of variables per twin
ncv <- 2 # number of covariates
- Read more about how to avoid NA in confidence interval
- 36 comments
- Log in or register to post comments
ReddyM
Joined: 10/02/2019
Difference between optimizers
In my twin model i calculate Rmz and Rdz from the saturated model, however when i run it (i run it on several variables) for a few variables CSOLNP fails to calculate confidence intervals, while SLSQP does caculate all. Is there a particular reason for this? Should i be concerned by this, or just use slsqp for the cases csolnp fails?
Here is the sat part, if it helps:
#sat model, rmz, rdz:
PathMZ <- mxMatrix("Lower",nrow = 2,ncol = 2,free = TRUE,values=starting_value_pathcoeff,name = "PathMZ")#cholesky for pos definitness
- Read more about Difference between optimizers
- 3 comments
- Log in or register to post comments
diana
Joined: 04/17/2019
the thresholds in SEM of ordinal variables
Hello, I have got a question about the thresholds in SEM of ordinal variables and really need your help.
- Read more about the thresholds in SEM of ordinal variables
- 20 comments
- Log in or register to post comments
lior abramson
Joined: 07/21/2017
direct comparison between traits' heritability with no raw data
Dear forum
I am trying to conduct a meta-analysis on the heritability of two traits that are different, yet related (I found six studies in which the twins were measured in both traits).
From what I know, in order to directly compare between the two traits heritability, I need to enter them in the same model, and for that, I need either the raw data or the full correlation matrix (and I have neither).
samson91787
Joined: 02/10/2014
Multilevel Logistic Regression
Hi,
- Read more about Multilevel Logistic Regression
- 7 comments
- Log in or register to post comments
Lisa M
Joined: 08/19/2013
Cholesky mx status RED
Hi,
I am trying to run a Cholesky Decomposition with 3 binary variables and I keep getting the error message
“In model 'multiCholACEModel' Optimizer 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)”
- Read more about Cholesky mx status RED
- 8 comments
- Log in or register to post comments