##################################################################################### ###### Bivariate Purcell moderation model ######### ##################################################################################### ################# Full model output > summary(ADEmodFit) Summary of ADEmod The final iterate satisfies the optimality conditions to the accuracy requested, but the sequence of iterates has not yet converged. Optimizer was terminated because no further improvement could be made in the merit function (Mx status GREEN). free parameters: name matrix row col Estimate lbound ubound 1 aM MZ.a 1 1 4.409966e-01 0 2 aC MZ.a 2 1 3.128132e-01 -10 3 aU MZ.a 2 2 -6.908945e-04 0! 4 dM MZ.d 1 1 4.314682e-01 0 5 dC MZ.d 2 1 -3.805734e-01 -10 6 dU MZ.d 2 2 1.370199e-01 0 7 eM MZ.e 1 1 7.686737e-01 0 8 eC MZ.e 2 1 -2.368976e-01 -10 9 eU MZ.e 2 2 8.261520e-01 0 10 aMod21 MZ.aMod 2 1 1.045334e-01 11 aMod22 MZ.aMod 2 2 -6.038269e-05 12 dMod21 MZ.dMod 2 1 -2.118522e-01 13 dMod22 MZ.dMod 2 2 8.552178e-02 14 eMod21 MZ.eMod 2 1 4.001176e-02 15 eMod22 MZ.eMod 2 2 3.094684e-01 16 beta_Age_lowsupport_sc_s MZ.bAge 1 1 4.420172e-02 17 beta_Age_ibs_gen_s MZ.bAge 1 2 3.942447e-02 18 beta_Sex_lowsupport_sc_s MZ.bSex 1 1 -1.431220e-01 19 beta_Sex_ibs_gen_s MZ.bSex 1 2 3.755939e-01 20 meanM MZ.intercept 1 1 8.068877e-02 21 ibs_gen_s_thresh MZ.Threshold 1 ibs_gen_s_1 1.381383e+00 Model Statistics: | Parameters | Degrees of Freedom | Fit (-2lnL units) Model: 21 7505 12907.98 Saturated: NA NA NA Independence: NA NA NA Number of observations/statistics: 1881/7526 Constraint 'MZ.Var1' contributes 1 observed statistic. Constraint 'DZ.Var1' contributes 1 observed statistic. Information Criteria: | df Penalty | Parameters Penalty | Sample-Size Adjusted AIC: -2102.015 12949.98 12950.48 BIC: -43676.404 13066.32 12999.60 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-04-17 12:34:33 Wall clock time: 41.89694 secs optimizer: NPSOL OpenMx version number: 2.9.9 Need help? See help(mxSummary) > mxEval(cov2cor(A0)[2,1], ADEmodFit$MZ) [,1] [1,] 0.9999976 > mxEval(cov2cor(D0)[2,1], ADEmodFit$MZ) [,1] [1,] -0.9408768 > mxEval(cov2cor(E0)[2,1], ADEmodFit$MZ) [,1] [1,] -0.2756399 > mxEval(cov2cor(V0)[2,1], ADEmodFit$MZ) [,1] [1,] -0.2113748 ################# Main effects (no moderation) model output > summary(MainEffectsFit) Summary of MainEffects The final iterate satisfies the optimality conditions to the accuracy requested, but the sequence of iterates has not yet converged. Optimizer was terminated because no further improvement could be made in the merit function (Mx status GREEN). free parameters: name matrix row col Estimate lbound ubound 1 aM MZ.a 1 1 4.400746e-01 0 2 aC MZ.a 2 1 5.077977e-01 -10 3 aU MZ.a 2 2 5.000287e-04 0! 4 dM MZ.d 1 1 4.319640e-01 0 5 dC MZ.d 2 1 -2.250430e-01 -10 6 dU MZ.d 2 2 -5.423199e-05 0! 7 eM MZ.e 1 1 7.688411e-01 0 8 eC MZ.e 2 1 5.806634e-02 -10 9 eU MZ.e 2 2 8.295395e-01 0 10 beta_Age_lowsupport_sc_s MZ.bAge 1 1 4.416224e-02 11 beta_Age_ibs_gen_s MZ.bAge 1 2 5.397301e-02 12 beta_Sex_lowsupport_sc_s MZ.bSex 1 1 -1.440010e-01 13 beta_Sex_ibs_gen_s MZ.bSex 1 2 3.323960e-01 14 meanM MZ.intercept 1 1 8.108537e-02 15 ibs_gen_s_thresh MZ.Threshold 1 ibs_gen_s_1 1.430251e+00 Model Statistics: | Parameters | Degrees of Freedom | Fit (-2lnL units) Model: 15 7511 12922.03 Saturated: NA NA NA Independence: NA NA NA Number of observations/statistics: 1881/7526 Constraint 'MZ.Var1' contributes 1 observed statistic. Constraint 'DZ.Var1' contributes 1 observed statistic. Information Criteria: | df Penalty | Parameters Penalty | Sample-Size Adjusted AIC: -2099.97 12952.03 12952.29 BIC: -43707.60 13035.12 12987.47 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-04-17 12:39:31 Wall clock time: 33.07527 secs optimizer: NPSOL OpenMx version number: 2.9.9 Need help? See help(mxSummary) > mxEval(cov2cor(A0)[2,1], MainEffectsFit$MZ) [,1] [1,] 0.9999995 > mxEval(cov2cor(D0)[2,1], MainEffectsFit$MZ) [,1] [1,] -1 > mxEval(cov2cor(E0)[2,1], MainEffectsFit$MZ) [,1] [1,] 0.06982742 > mxEval(cov2cor(V0)[2,1], MainEffectsFit$MZ) [,1] [1,] 0.1734011 ##################################################################################### ###### Bivariate correlated factors moderation model ######### ##################################################################################### ################# Full model output > summary(ADEmodFit) Summary of ADEmod The final iterate satisfies the optimality conditions to the accuracy requested, but the sequence of iterates has not yet converged. Optimizer was terminated because no further improvement could be made in the merit function (Mx status GREEN). free parameters: name matrix row col Estimate lbound ubound 1 aM a 1 1 0.43031292 0 2 aU a 2 2 0.40214526 0 3 dM d 1 1 0.44183711 0 4 dU d 2 2 0.33771102 0 5 eM e 1 1 0.76873979 0 6 eU e 2 2 0.85101736 0 7 aMod22 aMod 2 2 0.04301704 8 dMod22 dMod 2 2 0.15230565 9 eMod22 eMod 2 2 0.15058794 10 ra21 Ra 1 2 1.01417794 -0.9999 0.9999! 11 re21 Re 1 2 -0.07940281 -0.9999 0.9999 12 beta_Age_lowsupport_sc_s MZ.bAge 1 1 0.04393652 13 beta_Age_ibs_gen_s MZ.bAge 1 2 0.05415934 14 beta_Sex_lowsupport_sc_s MZ.bSex 1 1 -0.14383916 15 beta_Sex_ibs_gen_s MZ.bSex 1 2 0.38156152 16 meanM MZ.intercept 1 1 0.08101773 17 ibs_gen_s_thresh MZ.Threshold 1 ibs_gen_s_1 1.41749883 Model Statistics: | Parameters | Degrees of Freedom | Fit (-2lnL units) Model: 17 7509 12910.81 Saturated: NA NA NA Independence: NA NA NA Number of observations/statistics: 1881/7526 Constraint 'MZ.Var1' contributes 1 observed statistic. Constraint 'DZ.Var1' contributes 1 observed statistic. Information Criteria: | df Penalty | Parameters Penalty | Sample-Size Adjusted AIC: -2107.186 12944.81 12945.14 BIC: -43703.734 13038.99 12984.98 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-04-17 12:45:00 Wall clock time: 78.37859 secs optimizer: NPSOL OpenMx version number: 2.9.9 Need help? See help(mxSummary) > mxEval(cov2cor(V0)[2,1], ADEmodFit$MZ) [,1] [1,] 0.01938602 ################# Main effects (no moderation) model output > summary(MainEffectsFit) Summary of MainEffects The final iterate satisfies the optimality conditions to the accuracy requested, but the sequence of iterates has not yet converged. Optimizer was terminated because no further improvement could be made in the merit function (Mx status GREEN). free parameters: name matrix row col Estimate lbound ubound 1 aM a 1 1 0.43189596 0 2 aU a 2 2 0.48728721 0 3 dM d 1 1 0.44051020 0 4 dU d 2 2 0.26929678 0 5 eM e 1 1 0.76861275 0 6 eU e 2 2 0.83068090 0 7 ra21 Ra 1 2 0.99990000 -0.9999 0.9999! 8 re21 Re 1 2 0.06825480 -0.9999 0.9999 9 beta_Age_lowsupport_sc_s MZ.bAge 1 1 0.04417075 10 beta_Age_ibs_gen_s MZ.bAge 1 2 0.05395630 11 beta_Sex_lowsupport_sc_s MZ.bSex 1 1 -0.14401036 12 beta_Sex_ibs_gen_s MZ.bSex 1 2 0.33237975 13 meanM MZ.intercept 1 1 0.08110900 14 ibs_gen_s_thresh MZ.Threshold 1 ibs_gen_s_1 1.43021769 Model Statistics: | Parameters | Degrees of Freedom | Fit (-2lnL units) Model: 14 7512 12922.08 Saturated: NA NA NA Independence: NA NA NA Number of observations/statistics: 1881/7526 Constraint 'MZ.Var1' contributes 1 observed statistic. Constraint 'DZ.Var1' contributes 1 observed statistic. Information Criteria: | df Penalty | Parameters Penalty | Sample-Size Adjusted AIC: -2101.915 12950.08 12950.31 BIC: -43715.081 13027.64 12983.16 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-04-17 12:48:10 Wall clock time: 29.04761 secs optimizer: NPSOL OpenMx version number: 2.9.9 Need help? See help(mxSummary) > mxEval(cov2cor(V0)[2,1], MainEffectsFit$MZ) [,1] [1,] 0.1734787