Failure to obtain CI's for a umx model
Hi all,
I have used umx for a bivariate sex-limitation model, using the umxSexLim function, which have been a great help. However, I've failed in trying to compute CI's for my estimates, or any other indication of significance. Attached below are all the options I've tried. I would appreciate any help in how to compute any indication of significance.
Thanks!
umxCI(m2, run = "yes")
Running Scalar with 19 parameters
CSOLNP 64815751 -1.00592 -1.93e-09 Error in runHelper(model, frontendStart, intervals, silent, suppressWarnings, :
BLAS/LAPACK routine 'DGEBAL' gave error code -3
umxCI_boot(m2, type = "par.observed")
Error in if (model$data$type == "raw") { : argument is of length zero
umxConfint(m2, parm = "all",run = TRUE)
Running Scalar with 19 parameters
CSOLNP 50077024 503487 5.035e+05 Error in runHelper(model, frontendStart, intervals, silent, suppressWarnings, :
BLAS/LAPACK routine 'DGEBAL' gave error code -3
x <- umxRun(m2, intervals = TRUE)
Running Scalar with 19 parameters
summary(x)
Summary of Scalar
free parameters:
name matrix row col Estimate lbound ubound
1 am_r1c1 top.am 1 1 0.7895632830 1e-04
2 am_r2c2 top.am 2 2 0.0556255100 1e-04
3 cm_r1c1 top.cm 1 1 0.0001000004 1e-04!
4 cm_r2c2 top.cm 2 2 0.5058710146 1e-04
5 em_r1c1 top.em 1 1 0.6265102953 1e-04
6 em_r2c2 top.em 2 2 0.7983627259 1e-04
7 af_r1c1 top.af 1 1 1.2427580390 1e-04
8 af_r2c2 top.af 2 2 0.0001000354 1e-04!
9 cf_r1c1 top.cf 1 1 0.0001000000 1e-04!
10 cf_r2c2 top.cf 2 2 0.1961368709 1e-04
11 ef_r1c1 top.ef 1 1 0.3057593395 1e-04
12 ef_r2c2 top.ef 2 2 0.9729683133 1e-04
13 Ra_r2c1 top.Ram 1 2 -0.2442409485 -1 1
14 Rc_r2c1 top.Rcm 1 2 0.0629173480 -1 1
15 Re_r2c1 top.Rem 1 2 0.0883929800 -1 1
16 x_mean_m top.expMeanGm 1 parents_1 -0.0529413841
17 y_mean_m top.expMeanGm 1 peers_1 0.1286657693
18 x_mean_f top.expMeanGf 1 parents_1 0.0457339934
19 y_mean_f top.expMeanGf 1 peers_1 -0.0864138755
Model Statistics:
| Parameters | Degrees of Freedom | Fit (-2lnL units)
Model: 19 2259 7060.65
Saturated: NA NA NA
Independence: NA NA NA
Number of observations/statistics: 618/2278
Constraint 'top.Keep_it_Positive' contributes 0 observed statistics.
Information Criteria:
| df Penalty | Parameters Penalty | Sample-Size Adjusted
AIC: 2542.650 7098.650 7099.921
BIC: -7456.787 7182.754 7122.432
CFI: NA
TLI: 1 (also known as NNFI)
RMSEA: 0 [95% CI (NA, NA)]
Prob(RMSEA <= 0.05): NA
To get additional fit indices, see help(mxRefModels)
timestamp: 2019-06-27 15:12:35
Wall clock time: 0.2890499 secs
optimizer: CSOLNP
OpenMx version number: 2.12.2
Need help? See help(mxSummary)
umxSummary(m2, SE = TRUE)
umxSummarySexLim is a beta feature. Some things are broken. If any desired stats are not presented, let me know what's missing
Scalar -2 × log(Likelihood) = 7060.65
Scalar: 2-variable sex limitation analysis
Genetic Factor Correlations (male (Ram) lower triangle, female (Raf) upper)
|
x |
y |
x |
1.00 |
-0.24 |
y |
-0.24 |
1.00 |
Opposite sex A Correlations Rao
|
x |
y |
x |
1.00 |
-0.24 |
y |
-0.24 |
1.00 |
C Factor Correlations (male (Rcm) lower triangle, female (Rcf) upper)
|
x |
y |
x |
1.00 |
0.06 |
y |
0.06 |
1.00 |
Opposite sex C Correlations Rco
|
x |
y |
x |
1.00 |
0.06 |
y |
0.06 |
1.00 |
E Factor Correlations (male (Rem) lower triangle, female (Ref) upper)
|
x |
y |
x |
1.00 |
0.09 |
y |
0.09 |
1.00 |
Standardized solution (top.[ACE][mf]Std + R[ac]o matrices). Use std=F for raw variance.
|
am |
af |
cm |
cf |
em |
ef |
Rao |
Rco |
x |
0.56 |
0.80 |
. |
. |
0.44 |
0.20 |
1 |
1 |
y |
0.07 |
0.06 |
0.36 |
0.16 |
0.57 |
0.78 |
1 |
1 |
note: no plots for umxPlotSexLim
as yet.