I made a common pathway model in umx. I add CIs through adding umxCI(model, run="yes"). However, the output are NAs and there is a warning that "3: The nonlinear constraints and bounds could not be satisfied. The problem may have no feasible solution." I wonder what it means and what I can do to deal with it.
Here is my information about version and optimizer.
umx version: 3.0.1
OpenMx version: 2.14.11 [GIT v2.14.11]
R version: R version 3.5.3 (2019-03-11)
Platform: x86_64-apple-darwin15.6.0
MacOS: 10.14.6
Default optimizer: SLSQP
Here is my script.
setwd("....")
getwd()
Load Data; return a data frame
data <- read.spss(".....", to.data.frame = TRUE)
Pick the variables
selDVs = c("x","y","z")
mzData <- subset(data, zygosity=="1")
dzData <- subset(data, zygosity=="2")
umx_set_optimizer(opt = "SLSQP")
Common Pathway Models
cp <- umxCP(selDVs = selDVs, dzData = dzData, mzData = mzData, sep = "",tryHard = "yes")
cp <- umxCI(cp, run = "yes")
Here is the output
lbound estimate ubound lbound Code ubound Code
a_cp_r1c1 NA 0.806 NA NA NA
c_cp_r1c1 NA 0.000 NA NA NA
e_cp_r1c1 NA 0.592 NA NA NA
top.as_std[1,1] NA 0.000 NA NA NA
top.as_std[2,1] NA 0.000 NA NA NA
top.as_std[3,1] NA 0.000 NA NA NA
top.as_std[1,2] NA 0.000 NA NA NA
top.as_std[2,2] NA 0.000 NA NA NA
top.as_std[3,2] NA 0.000 NA NA NA
top.as_std[1,3] NA 0.000 NA NA NA
top.as_std[2,3] NA 0.000 NA NA NA
top.as_std[3,3] NA 0.067 NA NA NA
top.cs_std[1,1] NA 0.452 NA NA NA
top.cs_std[2,1] NA 0.000 NA NA NA
top.cs_std[3,1] NA 0.000 NA NA NA
top.cs_std[1,2] 0 0.000 NA 3 NA
top.cs_std[2,2] NA 0.172 NA NA NA
top.cs_std[3,2] NA 0.000 NA NA NA
top.cs_std[1,3] NA 0.000 NA NA NA
top.cs_std[2,3] NA 0.000 NA NA NA
top.cs_std[3,3] NA 0.315 NA NA NA
top.es_std[1,1] NA 0.708 NA NA NA
top.es_std[2,1] NA 0.000 NA NA NA
top.es_std[3,1] NA 0.000 NA NA NA
top.es_std[1,2] NA 0.000 NA NA NA
top.es_std[2,2] NA 0.472 NA NA NA
top.es_std[3,2] NA 0.000 NA NA NA
top.es_std[1,3] NA 0.000 NA NA NA
top.es_std[2,3] NA 0.000 NA NA NA
top.es_std[3,3] NA 0.510 NA NA NA
top.cp_loadings_std[1,1] NA 0.543 NA NA NA
top.cp_loadings_std[2,1] NA 0.864 NA NA NA
top.cp_loadings_std[3,1] NA 0.797 NA NA NA
expMean_Rumination1 NA 24.225 NA NA NA
expMean_Anxiety1 NA 38.829 NA NA NA
expMean_Depression1 NA 41.378 NA NA NA
as_r1c1 NA 0.000 NA NA NA
as_r2c2 NA 0.000 NA NA NA
as_r3c3 NA 0.590 NA NA NA
cs_r1c1 NA 2.677 NA NA NA
cs_r2c2 NA 1.516 NA NA NA
cs_r3c3 NA 2.766 NA NA NA
es_r1c1 NA 4.194 NA NA NA
es_r2c2 NA 4.152 NA NA NA
es_r3c3 NA 4.486 NA NA NA
cp_loadings_r1c1 NA 3.213 NA NA NA
cp_loadings_r2c1 NA 7.598 NA NA NA
cp_loadings_r3c1 NA 7.006 NA NA NA
[1] "NA: "
[1] "3: The nonlinear constraints and bounds could not be satisfied. The problem may have no feasible solution."