DATA
PREPARATION
Load all sheets of an Excel
Matrix we use for IV-UMASEM
## Correlation matrix
## ─────────────────────────────────────────────
## A C ES EX O X Y
## ─────────────────────────────────────────────
## A 1.000 0.290 0.160 0.200 0.190 0.130 0.130
## C 0.290 1.000 0.230 0.190 0.160 0.200 0.210
## ES 0.160 0.230 1.000 0.260 0.070 0.240 0.080
## EX 0.200 0.190 0.260 1.000 0.290 0.190 0.160
## O 0.190 0.160 0.070 0.290 1.000 0.010 0.060
## X 0.130 0.200 0.240 0.190 0.010 1.000 0.180
## Y 0.130 0.210 0.080 0.160 0.060 0.180 1.000
## ─────────────────────────────────────────────
## Sample size matrix
## ────────────────────────────────────────────────────────────────────────
## ES A C EX O X Y
## ────────────────────────────────────────────────────────────────────────
## 1 1.000 -999.000 -999.000 -999.000 -999.000 -999.000 -999.000
## 2 76406.000 1.000 -999.000 -999.000 -999.000 -999.000 -999.000
## 3 106149.000 80305.000 1.000 -999.000 -999.000 -999.000 -999.000
## 4 100801.000 70551.000 96442.000 1.000 -999.000 -999.000 -999.000
## 5 68068.000 67389.000 69753.000 68152.000 1.000 -999.000 -999.000
## 6 24527.000 11856.000 21719.000 20184.000 15196.000 1.000 -999.000
## 7 17274.000 14321.000 41939.000 19868.000 16068.000 54471.000 1.000
## ────────────────────────────────────────────────────────────────────────
MAIN
EFFECTS OF OCEAN ON X AND Y
## X ~ O + C + EX + A + ES
## Y ~ O + C + EX + A + ES
## Y ~~ X
## O ~~ C + EX + A + ES
## C ~~ EX + A + ES
## EX ~~ A + ES
## A ~~ ES
##
## Fit Measures (lavaan):
## χ²(0, N = 475) = 0.000, p = 1.000
## χ²/df = Inf (saturated model)
## AIC = 9219.285 (Akaike Information Criterion)
## BIC = 9335.857 (Bayesian Information Criterion)
## CFI = 1.000 (Comparative Fit Index)
## TLI = 1.000 (Tucker-Lewis Index; Non-Normed Fit Index, NNFI)
## NFI = 1.000 (Normed Fit Index)
## IFI = 1.000 (Incremental Fit Index)
## GFI = 1.000 (Goodness-of-Fit Index)
## AGFI = 1.000 (Adjusted Goodness-of-Fit Index)
## RMSEA = 0.000, 90% CI [0.000, 0.000] (Root Mean Square Error of Approximation)
## SRMR = 0.000 (Standardized Root Mean Square Residual)
##
## Model Estimates (lavaan):
## ────────────────────────────────────────────────────────────────────────
## Estimate S.E. z p LLCI ULCI Beta
## ────────────────────────────────────────────────────────────────────────
## Regression Paths:
## X <- O -0.071 (0.046) -1.533 .125 -0.161 0.020 -0.071
## X <- C 0.132 (0.047) 2.820 .005 ** 0.040 0.224 0.132
## X <- EX 0.130 (0.048) 2.742 .006 ** 0.037 0.223 0.130
## X <- A 0.052 (0.047) 1.108 .268 -0.040 0.143 0.052
## X <- ES 0.172 (0.046) 3.742 <.001 *** 0.082 0.263 0.172
## Y <- O -0.013 (0.047) -0.283 .777 -0.106 0.079 -0.013
## Y <- C 0.172 (0.048) 3.609 <.001 *** 0.079 0.266 0.172
## Y <- EX 0.119 (0.048) 2.459 .014 * 0.024 0.214 0.119
## Y <- A 0.059 (0.047) 1.235 .217 -0.034 0.152 0.059
## Y <- ES 0.001 (0.047) 0.020 .984 -0.091 0.093 0.001
## ────────────────────────────────────────────────────────────────────────
## Note. Raw (Standard) Confidence Interval (CI) and SE.
RESULTS
Longitudinal design
Change the data to CaseID=19
## Statistics of matrix
## ────────────────────────────
## harmonic_mean sum
## ────────────────────────────
## 1 20722.545 1009955.000
## ────────────────────────────
## Descriptive Statistics:
## ──────────────────────────────────────────────────────────────────────
## N Mean SD | Median Min Max Skewness Kurtosis
## ──────────────────────────────────────────────────────────────────────
## 21 48093.10 33825.24 | 41939.00 2987.00 106149.00 0.25 -1.56
## ──────────────────────────────────────────────────────────────────────
## Statistics of matrix
## ─────────────────────────
## harmonic_mean sum
## ─────────────────────────
## 1 69.494 1945.000
## ─────────────────────────
## Descriptive Statistics:
## ─────────────────────────────────────────────────────────
## N Mean SD | Median Min Max Skewness Kurtosis
## ─────────────────────────────────────────────────────────
## 21 92.62 42.25 | 92.00 21.00 163.00 0.00 -1.38
## ─────────────────────────────────────────────────────────
IVs
are OCEAN
## X ~ O + C + EX + A + ES
## Y ~ X
## Y ~~ X
## O ~~ C + EX + A + ES
## C ~~ EX + A + ES
## EX ~~ A + ES
## A ~~ ES
##
## Fit Measures (lavaan):
## χ²(4, N = 519) = 7.868, p = 0.097 .
## χ²/df = 1.967
## AIC = 10073.340 (Akaike Information Criterion)
## BIC = 10175.386 (Bayesian Information Criterion)
## CFI = 0.985 (Comparative Fit Index)
## TLI = 0.923 (Tucker-Lewis Index; Non-Normed Fit Index, NNFI)
## NFI = 0.972 (Normed Fit Index)
## IFI = 0.986 (Incremental Fit Index)
## GFI = 0.996 (Goodness-of-Fit Index)
## AGFI = 0.970 (Adjusted Goodness-of-Fit Index)
## RMSEA = 0.043, 90% CI [0.000, 0.088] (Root Mean Square Error of Approximation)
## SRMR = 0.020 (Standardized Root Mean Square Residual)
##
## Model Estimates (lavaan):
## ────────────────────────────────────────────────────────────────────────
## Estimate S.E. z p LLCI ULCI Beta
## ────────────────────────────────────────────────────────────────────────
## Regression Paths:
## X <- O -0.053 (0.037) -1.422 .155 -0.126 0.020 -0.053
## X <- C 0.165 (0.039) 4.186 <.001 *** 0.088 0.242 0.165
## X <- EX 0.140 (0.039) 3.550 <.001 *** 0.063 0.217 0.140
## X <- A 0.060 (0.037) 1.615 .106 -0.013 0.134 0.060
## X <- ES 0.114 (0.038) 3.024 .002 ** 0.040 0.189 0.114
## Y <- X 0.751 (0.166) 4.518 <.001 *** 0.425 1.077 0.751
## ────────────────────────────────────────────────────────────────────────
## Note. Raw (Standard) Confidence Interval (CI) and SE.
## ─────────────────────────────
## Wald.Statistic Wald.pValue
## ─────────────────────────────
## 1 54.090 0.000
## ─────────────────────────────
## $R2x.z
## [,1]
## [1,] 0.09853
Test for excluding the unqualified IVs
## X ~ C + EX + ES
## Y ~ X
## Y ~~ X
## O ~~ C + EX + A + ES
## C ~~ EX + A + ES
## EX ~~ A + ES
## A ~~ ES
##
## Fit Measures (lavaan):
## χ²(6, N = 519) = 12.049, p = 0.061 .
## χ²/df = 2.008
## AIC = 10073.521 (Akaike Information Criterion)
## BIC = 10167.063 (Bayesian Information Criterion)
## CFI = 0.977 (Comparative Fit Index)
## TLI = 0.920 (Tucker-Lewis Index; Non-Normed Fit Index, NNFI)
## NFI = 0.958 (Normed Fit Index)
## IFI = 0.978 (Incremental Fit Index)
## GFI = 0.993 (Goodness-of-Fit Index)
## AGFI = 0.969 (Adjusted Goodness-of-Fit Index)
## RMSEA = 0.044, 90% CI [0.000, 0.080] (Root Mean Square Error of Approximation)
## SRMR = 0.025 (Standardized Root Mean Square Residual)
##
## Model Estimates (lavaan):
## ─────────────────────────────────────────────────────────────────────
## Estimate S.E. z p LLCI ULCI Beta
## ─────────────────────────────────────────────────────────────────────
## Regression Paths:
## X <- C 0.174 (0.039) 4.498 <.001 *** 0.098 0.249 0.174
## X <- EX 0.133 (0.038) 3.515 <.001 *** 0.059 0.207 0.133
## X <- ES 0.119 (0.038) 3.137 .002 ** 0.045 0.193 0.119
## Y <- X 0.760 (0.172) 4.410 <.001 *** 0.422 1.097 0.760
## ─────────────────────────────────────────────────────────────────────
## Note. Raw (Standard) Confidence Interval (CI) and SE.
## ─────────────────────────────
## Wald.Statistic Wald.pValue
## ─────────────────────────────
## 1 50.449 0.000
## ─────────────────────────────
## $R2x.z
## [,1]
## [1,] 0.09242
## SEM Model fit results
## ──────────────────────────────────────────────────────────────────────────────────────
## Model ChiSqM_Value ChiSqM_DF ChiSqM_PValue CFI SRMR ChiSq_Diff DF_Diff ChiSq_p
## ──────────────────────────────────────────────────────────────────────────────────────
## 1 Model 1 7.868 4.000 0.097 0.985 0.020
## 2 Model 2 12.049 6.000 0.061 0.977 0.025 4.181 2.000 0.124
## ──────────────────────────────────────────────────────────────────────────────────────
IVs
are CEN
## X ~ C + EX + ES
## Y ~ X
## Y ~~ X
## C ~~ EX + ES
## EX ~~ ES
##
## Fit Measures (lavaan):
## χ²(2, N = 519) = 7.106, p = 0.029 *
## χ²/df = 3.553
## AIC = 7232.126 (Akaike Information Criterion)
## BIC = 7287.401 (Bayesian Information Criterion)
## CFI = 0.966 (Comparative Fit Index)
## TLI = 0.830 (Tucker-Lewis Index; Non-Normed Fit Index, NNFI)
## NFI = 0.956 (Normed Fit Index)
## IFI = 0.968 (Incremental Fit Index)
## GFI = 0.995 (Goodness-of-Fit Index)
## AGFI = 0.959 (Adjusted Goodness-of-Fit Index)
## RMSEA = 0.070, 90% CI [0.019, 0.129] (Root Mean Square Error of Approximation)
## SRMR = 0.024 (Standardized Root Mean Square Residual)
##
## Model Estimates (lavaan):
## ─────────────────────────────────────────────────────────────────────
## Estimate S.E. z p LLCI ULCI Beta
## ─────────────────────────────────────────────────────────────────────
## Regression Paths:
## X <- C 0.174 (0.039) 4.498 <.001 *** 0.098 0.249 0.174
## X <- EX 0.133 (0.038) 3.515 <.001 *** 0.059 0.207 0.133
## X <- ES 0.119 (0.038) 3.137 .002 ** 0.045 0.193 0.119
## Y <- X 0.760 (0.172) 4.410 <.001 *** 0.422 1.097 0.760
## ─────────────────────────────────────────────────────────────────────
## Note. Raw (Standard) Confidence Interval (CI) and SE.
## ─────────────────────────────
## Wald.Statistic Wald.pValue
## ─────────────────────────────
## 1 50.449 0.000
## ─────────────────────────────
## $R2x.z
## [,1]
## [1,] 0.09242
IVs
are EN with Y on C
## X ~ C + EX + ES
## Y ~ X + C
## Y ~~ X
## C ~~ EX + ES
## EX ~~ ES
##
## Fit Measures (lavaan):
## χ²(1, N = 519) = 3.805, p = 0.051 .
## χ²/df = 3.805
## AIC = 7230.825 (Akaike Information Criterion)
## BIC = 7290.351 (Bayesian Information Criterion)
## CFI = 0.981 (Comparative Fit Index)
## TLI = 0.813 (Tucker-Lewis Index; Non-Normed Fit Index, NNFI)
## NFI = 0.976 (Normed Fit Index)
## IFI = 0.982 (Incremental Fit Index)
## GFI = 0.997 (Goodness-of-Fit Index)
## AGFI = 0.956 (Adjusted Goodness-of-Fit Index)
## RMSEA = 0.074, 90% CI [0.000, 0.158] (Root Mean Square Error of Approximation)
## SRMR = 0.018 (Standardized Root Mean Square Residual)
##
## Model Estimates (lavaan):
## ─────────────────────────────────────────────────────────────────────
## Estimate S.E. z p LLCI ULCI Beta
## ─────────────────────────────────────────────────────────────────────
## Regression Paths:
## X <- C 0.137 (0.043) 3.162 .002 ** 0.052 0.222 0.137
## X <- EX 0.139 (0.042) 3.295 <.001 *** 0.056 0.222 0.139
## X <- ES 0.158 (0.043) 3.676 <.001 *** 0.074 0.242 0.158
## Y <- X 0.420 (0.197) 2.131 .033 * 0.034 0.806 0.420
## Y <- C 0.126 (0.060) 2.111 .035 * 0.009 0.243 0.126
## ─────────────────────────────────────────────────────────────────────
## Note. Raw (Standard) Confidence Interval (CI) and SE.
## ─────────────────────────────
## Wald.Statistic Wald.pValue
## ─────────────────────────────
## 1 29.583 0.000
## ─────────────────────────────
## $R2x.z
## [,1]
## [1,] 0.07506
## SEM Model fit results
## ──────────────────────────────────────────────────────────────────────────────────────
## Model ChiSqM_Value ChiSqM_DF ChiSqM_PValue CFI SRMR ChiSq_Diff DF_Diff ChiSq_p
## ──────────────────────────────────────────────────────────────────────────────────────
## 1 Model 1 3.805 1.000 0.051 0.981 0.018
## 2 Model 2 7.106 2.000 0.029 0.966 0.024 3.301 1.000 0.069
## ──────────────────────────────────────────────────────────────────────────────────────
IVs
are EN
## X ~ EX + ES
## Y ~ X
## Y ~~ X
## EX ~~ ES
##
## Fit Measures (lavaan):
## χ²(1, N = 519) = 3.542, p = 0.060 .
## χ²/df = 3.542
## AIC = 5811.839 (Akaike Information Criterion)
## BIC = 5850.106 (Bayesian Information Criterion)
## CFI = 0.972 (Comparative Fit Index)
## TLI = 0.833 (Tucker-Lewis Index; Non-Normed Fit Index, NNFI)
## NFI = 0.964 (Normed Fit Index)
## IFI = 0.974 (Incremental Fit Index)
## GFI = 0.997 (Goodness-of-Fit Index)
## AGFI = 0.966 (Adjusted Goodness-of-Fit Index)
## RMSEA = 0.070, 90% CI [0.000, 0.155] (Root Mean Square Error of Approximation)
## SRMR = 0.021 (Standardized Root Mean Square Residual)
##
## Model Estimates (lavaan):
## ─────────────────────────────────────────────────────────────────────
## Estimate S.E. z p LLCI ULCI Beta
## ─────────────────────────────────────────────────────────────────────
## Regression Paths:
## X <- EX 0.163 (0.041) 3.942 <.001 *** 0.082 0.244 0.163
## X <- ES 0.179 (0.042) 4.306 <.001 *** 0.098 0.261 0.179
## Y <- X 0.546 (0.173) 3.165 .002 ** 0.208 0.885 0.546
## ─────────────────────────────────────────────────────────────────────
## Note. Raw (Standard) Confidence Interval (CI) and SE.
## ─────────────────────────────
## Wald.Statistic Wald.pValue
## ─────────────────────────────
## 1 41.524 0.000
## ─────────────────────────────
## $R2x.z
## [,1]
## [1,] 0.07506
## SEM Model fit results
## ──────────────────────────────────────────────────────────────────────────────────────
## Model ChiSqM_Value ChiSqM_DF ChiSqM_PValue CFI SRMR ChiSq_Diff DF_Diff ChiSq_p
## ──────────────────────────────────────────────────────────────────────────────────────
## 1 Model 1 3.542 1.000 0.060 0.972 0.021
## 2 Model 2 3.542 1.000 0.060 0.972 0.021 0.000 0.000 1.000
## ──────────────────────────────────────────────────────────────────────────────────────
RESULTS SUMMARY
## ──────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
## Model ChiSqM_Value ChiSqM_DF ChiSqM_PValue CFI SRMR ChiSq_Diff DF_Diff ChiSq_p
## ──────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
## 1 M1. OCEAN as IVs 7.868 4.000 0.097 0.985 0.020
## 2 M2. CEN as IVs with OA retained 12.049 6.000 0.061 0.977 0.025 4.181 2.000 0.124
## 3 M3. CEN as IVs with OA removed 7.106 2.000 0.029 0.966 0.024 3.301 1.000 0.069
## 4 M4. EN as IVs with the effect of C on Y 3.805 1.000 0.051 0.981 0.018
## 5 M5. EN as IVs 3.542 1.000 0.060 0.972 0.021
## ──────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
## Combined SEM Path Results
## ─────────────────────────────────────────────────────────────────────────────────────────────────────
## label M1.est_sig M1.se M2.est_sig M2.se M3.est_sig M3.se M4.est_sig M4.se M5.est_sig M5.se
## ─────────────────────────────────────────────────────────────────────────────────────────────────────
## 1 A.WITH.ES 0.160*** 0.044 NA NA NA NA
## 2 C.WITH.A 0.289*** 0.046 0.289*** 0.046 NA NA NA
## 3 C.WITH.ES 0.230*** 0.045 0.230*** 0.045 0.230*** 0.045 0.230*** 0.045 NA
## 4 C.WITH.EX 0.190*** 0.045 0.190*** 0.045 0.190*** 0.045 0.190*** 0.045 NA
## 5 C.WITH.O NA 0.160*** 0.044 NA NA NA
## 6 ES.WITH.A NA 0.160*** 0.044 NA NA NA
## 7 ES.WITH.O NA 0.070 0.044 NA NA NA
## 8 EX.WITH.A 0.200*** 0.045 0.200*** 0.045 NA NA NA
## 9 EX.WITH.ES 0.259*** 0.045 0.259*** 0.045 0.259*** 0.045 0.259*** 0.045 0.259*** 0.045
## 10 EX.WITH.O NA 0.289*** 0.046 NA NA NA
## 11 O.WITH.A 0.190*** 0.045 0.190*** 0.045 NA NA NA
## 12 O.WITH.C 0.160*** 0.044 NA NA NA NA
## 13 O.WITH.ES 0.070 0.044 NA NA NA NA
## 14 O.WITH.EX 0.289*** 0.046 NA NA NA NA
## 15 Variances.A 0.998*** 0.062 0.998*** 0.062 NA NA NA
## 16 Variances.C 0.998*** 0.062 0.998*** 0.062 0.998*** 0.062 0.998*** 0.062 NA
## 17 Variances.ES 0.998*** 0.062 0.998*** 0.062 0.998*** 0.062 0.998*** 0.062 0.998*** 0.062
## 18 Variances.EX 0.998*** 0.062 0.998*** 0.062 0.998*** 0.062 0.998*** 0.062 0.998*** 0.062
## 19 Variances.O 0.998*** 0.062 0.998*** 0.062 NA NA NA
## 20 Variances.X 0.904*** 0.056 0.910*** 0.056 0.910*** 0.056 0.906*** 0.056 0.924*** 0.057
## 21 Variances.Y 1.351*** 0.219 1.362*** 0.229 1.362*** 0.229 1.041*** 0.136 1.143*** 0.157
## 22 X.ON.A 0.060 0.037 NA NA NA NA
## 23 X.ON.C 0.165*** 0.039 0.174*** 0.039 0.174*** 0.039 0.137** 0.043 NA
## 24 X.ON.ES 0.114** 0.038 0.119** 0.038 0.119** 0.038 0.158*** 0.043 0.179*** 0.042
## 25 X.ON.EX 0.140*** 0.039 0.133*** 0.038 0.133*** 0.038 0.139*** 0.042 0.163*** 0.041
## 26 X.ON.O -0.053 0.037 NA NA NA NA
## 27 X.WITH.Y -0.610*** 0.160 -0.618*** 0.166 -0.618*** 0.166 -0.305 0.184 -0.406* 0.167
## 28 Y.ON.C NA NA NA 0.126* 0.060 NA
## 29 Y.ON.X 0.751*** 0.166 0.760*** 0.172 0.760*** 0.172 0.420* 0.197 0.546** 0.173
## ─────────────────────────────────────────────────────────────────────────────────────────────────────
## ────────────────────────────────────────────────────────────────────────────────────────────────
## label M1.est_sig M1.se M2.est_sig M2.se M3.est_sig M3.se M4.est_sig M4.se M5.est_sig M5.se
## ────────────────────────────────────────────────────────────────────────────────────────────────
## 1 X.ON.A 0.060 0.037 NA NA NA NA
## 2 X.ON.C 0.165*** 0.039 0.174*** 0.039 0.174*** 0.039 0.137** 0.043 NA
## 3 X.ON.ES 0.114** 0.038 0.119** 0.038 0.119** 0.038 0.158*** 0.043 0.179*** 0.042
## 4 X.ON.EX 0.140*** 0.039 0.133*** 0.038 0.133*** 0.038 0.139*** 0.042 0.163*** 0.041
## 5 X.ON.O -0.053 0.037 NA NA NA NA
## 6 X.WITH.Y -0.610*** 0.160 -0.618*** 0.166 -0.618*** 0.166 -0.305 0.184 -0.406* 0.167
## 7 Y.ON.C NA NA NA 0.126* 0.060 NA
## 8 Y.ON.X 0.751*** 0.166 0.760*** 0.172 0.760*** 0.172 0.420* 0.197 0.546** 0.173
## ────────────────────────────────────────────────────────────────────────────────────────────────