> ### stage 2 > require("OpenMx") > library(OpenMx) > A <- create.mxMatrix( + c( 0,"0.1*CEO2RND",0,0,0, "0.1*ENV2RND",0,0, + 0,0,0,0,0,0,0,0, + 0,"0.1*CEO2ABS",0,0,0,0,"0.1*REL2ABS",0, + 0,"0.1*CEO2COM",0,0,"0.1*TEC2COM","0.1*ENV2COM","0.1*REL2COM",0, + "0.1*RND2TEC",0,"0.1*ABS2TEC",0,0,0,0,0, + 0,0,0,0,0,0,0,0, + 0,0,0,0,0,0,0,0, + 0,0,0,"0.1*COM2PER",0,0,0,0), + type = "Full", + nrow = 8, + ncol = 8, + byrow = TRUE, + name = "A", + dimnames = list(varnames,varnames)) > dimnames(A)[[1]] <-dimnames(A)[[2]] <- c("RND","CEO","ABS","COM","TEC","ENV","REL","PER") > A FullMatrix 'A' $labels RND CEO ABS COM TEC ENV REL PER RND NA "CEO2RND" NA NA NA "ENV2RND" NA NA CEO NA NA NA NA NA NA NA NA ABS NA "CEO2ABS" NA NA NA NA "REL2ABS" NA COM NA "CEO2COM" NA NA "TEC2COM" "ENV2COM" "REL2COM" NA TEC "RND2TEC" NA "ABS2TEC" NA NA NA NA NA ENV NA NA NA NA NA NA NA NA REL NA NA NA NA NA NA NA NA PER NA NA NA "COM2PER" NA NA NA NA $values RND CEO ABS COM TEC ENV REL PER RND 0.0 0.1 0.0 0.0 0.0 0.1 0.0 0 CEO 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0 ABS 0.0 0.1 0.0 0.0 0.0 0.0 0.1 0 COM 0.0 0.1 0.0 0.0 0.1 0.1 0.1 0 TEC 0.1 0.0 0.1 0.0 0.0 0.0 0.0 0 ENV 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0 REL 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0 PER 0.0 0.0 0.0 0.1 0.0 0.0 0.0 0 $free RND CEO ABS COM TEC ENV REL PER RND FALSE TRUE FALSE FALSE FALSE TRUE FALSE FALSE CEO FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE ABS FALSE TRUE FALSE FALSE FALSE FALSE TRUE FALSE COM FALSE TRUE FALSE FALSE TRUE TRUE TRUE FALSE TEC TRUE FALSE TRUE FALSE FALSE FALSE FALSE FALSE ENV FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE REL FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE PER FALSE FALSE FALSE TRUE FALSE FALSE FALSE FALSE $lbound: No lower bounds assigned. $ubound: No upper bounds assigned. > S <- create.mxMatrix( + c("1*var_RND", + 0,1, + 0,0,"1*var_ABS", + 0,0,0,"1*var_COM", + 0,0,0,0,"1*var_TEC", + 0,"0.1*cor_CEOENV",0,0,0,1, + 0,"0.1*cor_CEOREL",0,0,0,"0.1*cor_ENVREL",1, + 0,0,0,0,0,0,0,"1*var_PER"), + type = "Symm", + byrow = TRUE, + name="S", + dimnames = list(varnames,varnames)) > S SymmMatrix 'S' $labels RND CEO ABS COM TEC ENV REL PER RND "var_RND" NA NA NA NA NA NA NA CEO NA NA NA NA NA "cor_CEOENV" "cor_CEOREL" NA ABS NA NA "var_ABS" NA NA NA NA NA COM NA NA NA "var_COM" NA NA NA NA TEC NA NA NA NA "var_TEC" NA NA NA ENV NA "cor_CEOENV" NA NA NA NA "cor_ENVREL" NA REL NA "cor_CEOREL" NA NA NA "cor_ENVREL" NA NA PER NA NA NA NA NA NA NA "var_PER" $values RND CEO ABS COM TEC ENV REL PER RND 1 0.0 0 0 0 0.0 0.0 0 CEO 0 1.0 0 0 0 0.1 0.1 0 ABS 0 0.0 1 0 0 0.0 0.0 0 COM 0 0.0 0 1 0 0.0 0.0 0 TEC 0 0.0 0 0 1 0.0 0.0 0 ENV 0 0.1 0 0 0 1.0 0.1 0 REL 0 0.1 0 0 0 0.1 1.0 0 PER 0 0.0 0 0 0 0.0 0.0 1 $free RND CEO ABS COM TEC ENV REL PER RND TRUE FALSE FALSE FALSE FALSE FALSE FALSE FALSE CEO FALSE FALSE FALSE FALSE FALSE TRUE TRUE FALSE ABS FALSE FALSE TRUE FALSE FALSE FALSE FALSE FALSE COM FALSE FALSE FALSE TRUE FALSE FALSE FALSE FALSE TEC FALSE FALSE FALSE FALSE TRUE FALSE FALSE FALSE ENV FALSE TRUE FALSE FALSE FALSE FALSE TRUE FALSE REL FALSE TRUE FALSE FALSE FALSE TRUE FALSE FALSE PER FALSE FALSE FALSE FALSE FALSE FALSE FALSE TRUE $lbound: No lower bounds assigned. $ubound: No upper bounds assigned. > stage2 <- tssem2(stage1random, Amatrix=A, Smatrix=S, + diag.constraints = TRUE, intervals.type = "LB") Warning message: In .solve(x = object$mx.fit@output$calculatedHessian, parameters = my.name) : Error in solving the Hessian matrix. Generalized inverse is used. The standard errors may not be trustworthy. > summary(stage2) Call: wls(Cov = pooledS, aCov = aCov, n = tssem1.obj$total.n, Amatrix = Amatrix, Smatrix = Smatrix, Fmatrix = Fmatrix, diag.constraints = diag.constraints, cor.analysis = cor.analysis, intervals.type = intervals.type, mx.algebras = mx.algebras, model.name = model.name, suppressWarnings = suppressWarnings, silent = silent, run = run) 95% confidence intervals: Likelihood-based statistic Coefficients: Estimate Std.Error lbound ubound z value Pr(>|z|) CEO2ABS 0.474405 NA 0.385012 0.560683 NA NA REL2ABS 0.269626 NA 0.143822 0.385825 NA NA CEO2COM 0.287259 NA 0.089165 0.474484 NA NA ENV2COM 0.327082 NA NA 0.483294 NA NA REL2COM 0.435905 NA NA 0.614947 NA NA TEC2COM 0.100417 NA -0.067508 0.261062 NA NA COM2PER 0.519724 NA 0.455970 0.585385 NA NA CEO2RND 0.629701 NA 0.533918 0.726958 NA NA ENV2RND 0.014355 NA -0.093106 0.117244 NA NA ABS2TEC 0.561304 NA 0.370425 0.753955 NA NA RND2TEC -0.046011 NA -0.236604 0.138792 NA NA var_ABS 0.582054 NA 0.506237 0.649113 NA NA cor_CEOENV 0.075803 NA -0.026694 0.178967 NA NA cor_CEOREL 0.469808 NA 0.374715 0.565709 NA NA var_COM 0.395267 NA 0.228393 0.522889 NA NA cor_ENVREL 0.139974 NA -0.068350 0.351312 NA NA var_PER 0.729887 NA 0.657316 0.792091 NA NA var_RND 0.601900 NA 0.472553 0.711500 NA NA var_TEC 0.702425 NA 0.508399 0.831507 NA NA Goodness-of-fit indices: Value Sample size 18902.0000 Chi-square of target model 144.9502 DF of target model 14.0000 p value of target model 0.0000 Number of constraints imposed on "Smatrix" 5.0000 DF manually adjusted 0.0000 Chi-square of independence model 1338.1849 DF of independence model 28.0000 RMSEA 0.0222 RMSEA lower 95% CI 0.0190 RMSEA upper 95% CI 0.0256 SRMR 0.1490 TLI 0.8001 CFI 0.9001 AIC 116.9502 BIC 7.0919 OpenMx status1: 0 ("0" or "1": The optimization is considered fine. Other values indicate problems.)