library(rio)
library(lavaan)## This is lavaan 0.6-5## lavaan is BETA software! Please report any bugs.FullData <- read.csv("/Users/Lorraine/Desktop/SAQ.csv")SingleF <- '
f =~ q1 + q2 + q3 + q4 + q5 + q6 + q7 + q8 + q9 + q10 + q11 + q12 + q13 + q14 + q15 + q16 + q17 + q18 + q19 + q20 + q21 + q22 + q23
'
fit <- cfa(SingleF, data = FullData)
summary(fit, fit.measures=TRUE)## lavaan 0.6-5 ended normally after 34 iterations
##
## Estimator ML
## Optimization method NLMINB
## Number of free parameters 46
##
## Number of observations 2569
##
## Model Test User Model:
##
## Test statistic 4475.734
## Degrees of freedom 230
## P-value (Chi-square) 0.000
##
## Model Test Baseline Model:
##
## Test statistic 19390.907
## Degrees of freedom 253
## P-value 0.000
##
## User Model versus Baseline Model:
##
## Comparative Fit Index (CFI) 0.778
## Tucker-Lewis Index (TLI) 0.756
##
## Loglikelihood and Information Criteria:
##
## Loglikelihood user model (H0) -75327.499
## Loglikelihood unrestricted model (H1) -73089.632
##
## Akaike (AIC) 150746.998
## Bayesian (BIC) 151016.157
## Sample-size adjusted Bayesian (BIC) 150870.002
##
## Root Mean Square Error of Approximation:
##
## RMSEA 0.085
## 90 Percent confidence interval - lower 0.083
## 90 Percent confidence interval - upper 0.087
## P-value RMSEA <= 0.05 0.000
##
## Standardized Root Mean Square Residual:
##
## SRMR 0.063
##
## Parameter Estimates:
##
## Information Expected
## Information saturated (h1) model Structured
## Standard errors Standard
##
## Latent Variables:
## Estimate Std.Err z-value P(>|z|)
## f =~
## q1 1.000
## q2 -0.492 0.040 -12.422 0.000
## q3 -1.374 0.057 -23.922 0.000
## q4 1.250 0.051 24.421 0.000
## q5 1.093 0.050 21.928 0.000
## q6 1.312 0.058 22.432 0.000
## q7 1.586 0.061 25.893 0.000
## q8 0.985 0.045 21.867 0.000
## q9 -0.675 0.058 -11.548 0.000
## q10 0.771 0.043 17.936 0.000
## q11 1.182 0.048 24.741 0.000
## q12 1.282 0.050 25.425 0.000
## q13 1.340 0.052 25.589 0.000
## q14 1.370 0.055 25.098 0.000
## q15 1.225 0.053 23.060 0.000
## q16 1.292 0.051 25.552 0.000
## q17 1.182 0.048 24.673 0.000
## q18 1.562 0.059 26.396 0.000
## q19 -0.924 0.054 -17.249 0.000
## q20 0.899 0.051 17.759 0.000
## q21 1.345 0.054 25.025 0.000
## q22 -0.598 0.048 -12.344 0.000
## q23 -0.275 0.047 -5.851 0.000
##
## Variances:
## Estimate Std.Err z-value P(>|z|)
## .q1 0.471 0.014 34.286 0.000
## .q2 0.672 0.019 35.576 0.000
## .q3 0.751 0.022 33.999 0.000
## .q4 0.566 0.017 33.818 0.000
## .q5 0.675 0.020 34.543 0.000
## .q6 0.891 0.026 34.427 0.000
## .q7 0.677 0.020 33.116 0.000
## .q8 0.553 0.016 34.557 0.000
## .q9 1.496 0.042 35.617 0.000
## .q10 0.642 0.018 35.163 0.000
## .q11 0.476 0.014 33.689 0.000
## .q12 0.488 0.015 33.373 0.000
## .q13 0.515 0.015 33.287 0.000
## .q14 0.596 0.018 33.532 0.000
## .q15 0.697 0.020 34.264 0.000
## .q16 0.483 0.014 33.307 0.000
## .q17 0.482 0.014 33.717 0.000
## .q18 0.587 0.018 32.795 0.000
## .q19 1.031 0.029 35.234 0.000
## .q20 0.899 0.026 35.182 0.000
## .q21 0.582 0.017 33.565 0.000
## .q22 1.006 0.028 35.580 0.000
## .q23 1.072 0.030 35.788 0.000
## f 0.214 0.014 14.812 0.000Comparative Fit Index (CFI) = 0.778, Tucker-Lewis Index (TLI) = 0.756, and SRMR = 0.063, not a good model.
ThreeF<- '
f1 =~ q1 + q4 + q5 + q8 + q11 + q16 + q17 + q20 + q21
f2 =~ q6 + q7 + q10 + q12 + q13 + q14 + q15 + q18
f3 =~ q2 + q3 + q9 + q19 + q22 + q23
'
fit1 <- cfa(ThreeF, data = FullData)
summary(fit1, fit.measures = TRUE)## lavaan 0.6-5 ended normally after 52 iterations
##
## Estimator ML
## Optimization method NLMINB
## Number of free parameters 49
##
## Number of observations 2569
##
## Model Test User Model:
##
## Test statistic 3299.622
## Degrees of freedom 227
## P-value (Chi-square) 0.000
##
## Model Test Baseline Model:
##
## Test statistic 19390.907
## Degrees of freedom 253
## P-value 0.000
##
## User Model versus Baseline Model:
##
## Comparative Fit Index (CFI) 0.839
## Tucker-Lewis Index (TLI) 0.821
##
## Loglikelihood and Information Criteria:
##
## Loglikelihood user model (H0) -74739.443
## Loglikelihood unrestricted model (H1) -73089.632
##
## Akaike (AIC) 149576.887
## Bayesian (BIC) 149863.599
## Sample-size adjusted Bayesian (BIC) 149707.913
##
## Root Mean Square Error of Approximation:
##
## RMSEA 0.073
## 90 Percent confidence interval - lower 0.070
## 90 Percent confidence interval - upper 0.075
## P-value RMSEA <= 0.05 0.000
##
## Standardized Root Mean Square Residual:
##
## SRMR 0.056
##
## Parameter Estimates:
##
## Information Expected
## Information saturated (h1) model Structured
## Standard errors Standard
##
## Latent Variables:
## Estimate Std.Err z-value P(>|z|)
## f1 =~
## q1 1.000
## q4 1.204 0.047 25.774 0.000
## q5 1.050 0.046 22.905 0.000
## q8 1.060 0.042 24.956 0.000
## q11 1.198 0.044 27.114 0.000
## q16 1.223 0.046 26.744 0.000
## q17 1.206 0.044 27.172 0.000
## q20 0.854 0.047 18.162 0.000
## q21 1.238 0.048 25.598 0.000
## f2 =~
## q6 1.000
## q7 1.102 0.039 28.445 0.000
## q10 0.540 0.028 19.082 0.000
## q12 0.864 0.032 27.220 0.000
## q13 0.952 0.033 28.536 0.000
## q14 0.953 0.035 27.474 0.000
## q15 0.813 0.034 24.028 0.000
## q18 1.119 0.038 29.721 0.000
## f3 =~
## q2 1.000
## q3 2.097 0.115 18.244 0.000
## q9 1.488 0.100 14.849 0.000
## q19 1.478 0.093 15.947 0.000
## q22 1.066 0.078 13.614 0.000
## q23 0.657 0.070 9.443 0.000
##
## Covariances:
## Estimate Std.Err z-value P(>|z|)
## f1 ~~
## f2 0.284 0.014 19.743 0.000
## f3 -0.121 0.008 -14.431 0.000
## f2 ~~
## f3 -0.169 0.012 -14.586 0.000
##
## Variances:
## Estimate Std.Err z-value P(>|z|)
## .q1 0.439 0.013 32.994 0.000
## .q4 0.543 0.017 32.505 0.000
## .q5 0.659 0.020 33.754 0.000
## .q8 0.484 0.015 32.940 0.000
## .q11 0.422 0.013 31.576 0.000
## .q16 0.471 0.015 31.866 0.000
## .q17 0.423 0.013 31.528 0.000
## .q20 0.892 0.026 34.822 0.000
## .q21 0.592 0.018 32.606 0.000
## .q6 0.780 0.024 32.856 0.000
## .q7 0.633 0.020 31.359 0.000
## .q10 0.630 0.018 34.769 0.000
## .q12 0.482 0.015 32.229 0.000
## .q13 0.465 0.015 31.282 0.000
## .q14 0.562 0.018 32.070 0.000
## .q15 0.702 0.021 33.653 0.000
## .q18 0.510 0.017 30.084 0.000
## .q2 0.584 0.018 32.843 0.000
## .q3 0.537 0.025 21.162 0.000
## .q9 1.282 0.039 32.820 0.000
## .q19 0.906 0.029 31.580 0.000
## .q22 0.923 0.027 33.700 0.000
## .q23 1.027 0.029 35.117 0.000
## f1 0.246 0.016 15.811 0.000
## f2 0.479 0.029 16.423 0.000
## f3 0.141 0.014 10.095 0.000parameterEstimates(fit1)## lhs op rhs est se z pvalue ci.lower ci.upper
## 1 f1 =~ q1 1.000 0.000 NA NA 1.000 1.000
## 2 f1 =~ q4 1.204 0.047 25.774 0 1.112 1.295
## 3 f1 =~ q5 1.050 0.046 22.905 0 0.960 1.140
## 4 f1 =~ q8 1.060 0.042 24.956 0 0.977 1.143
## 5 f1 =~ q11 1.198 0.044 27.114 0 1.111 1.285
## 6 f1 =~ q16 1.223 0.046 26.744 0 1.133 1.313
## 7 f1 =~ q17 1.206 0.044 27.172 0 1.119 1.293
## 8 f1 =~ q20 0.854 0.047 18.162 0 0.762 0.947
## 9 f1 =~ q21 1.238 0.048 25.598 0 1.143 1.332
## 10 f2 =~ q6 1.000 0.000 NA NA 1.000 1.000
## 11 f2 =~ q7 1.102 0.039 28.445 0 1.026 1.178
## 12 f2 =~ q10 0.540 0.028 19.082 0 0.484 0.595
## 13 f2 =~ q12 0.864 0.032 27.220 0 0.802 0.926
## 14 f2 =~ q13 0.952 0.033 28.536 0 0.887 1.018
## 15 f2 =~ q14 0.953 0.035 27.474 0 0.885 1.021
## 16 f2 =~ q15 0.813 0.034 24.028 0 0.747 0.879
## 17 f2 =~ q18 1.119 0.038 29.721 0 1.045 1.193
## 18 f3 =~ q2 1.000 0.000 NA NA 1.000 1.000
## 19 f3 =~ q3 2.097 0.115 18.244 0 1.871 2.322
## 20 f3 =~ q9 1.488 0.100 14.849 0 1.291 1.684
## 21 f3 =~ q19 1.478 0.093 15.947 0 1.297 1.660
## 22 f3 =~ q22 1.066 0.078 13.614 0 0.912 1.219
## 23 f3 =~ q23 0.657 0.070 9.443 0 0.520 0.793
## 24 q1 ~~ q1 0.439 0.013 32.994 0 0.413 0.465
## 25 q4 ~~ q4 0.543 0.017 32.505 0 0.510 0.576
## 26 q5 ~~ q5 0.659 0.020 33.754 0 0.621 0.697
## 27 q8 ~~ q8 0.484 0.015 32.940 0 0.456 0.513
## 28 q11 ~~ q11 0.422 0.013 31.576 0 0.396 0.448
## 29 q16 ~~ q16 0.471 0.015 31.866 0 0.442 0.500
## 30 q17 ~~ q17 0.423 0.013 31.528 0 0.396 0.449
## 31 q20 ~~ q20 0.892 0.026 34.822 0 0.842 0.943
## 32 q21 ~~ q21 0.592 0.018 32.606 0 0.557 0.628
## 33 q6 ~~ q6 0.780 0.024 32.856 0 0.734 0.827
## 34 q7 ~~ q7 0.633 0.020 31.359 0 0.594 0.673
## 35 q10 ~~ q10 0.630 0.018 34.769 0 0.595 0.666
## 36 q12 ~~ q12 0.482 0.015 32.229 0 0.453 0.511
## 37 q13 ~~ q13 0.465 0.015 31.282 0 0.436 0.494
## 38 q14 ~~ q14 0.562 0.018 32.070 0 0.528 0.596
## 39 q15 ~~ q15 0.702 0.021 33.653 0 0.661 0.743
## 40 q18 ~~ q18 0.510 0.017 30.084 0 0.476 0.543
## 41 q2 ~~ q2 0.584 0.018 32.843 0 0.549 0.619
## 42 q3 ~~ q3 0.537 0.025 21.162 0 0.487 0.587
## 43 q9 ~~ q9 1.282 0.039 32.820 0 1.205 1.359
## 44 q19 ~~ q19 0.906 0.029 31.580 0 0.850 0.962
## 45 q22 ~~ q22 0.923 0.027 33.700 0 0.869 0.977
## 46 q23 ~~ q23 1.027 0.029 35.117 0 0.970 1.085
## 47 f1 ~~ f1 0.246 0.016 15.811 0 0.216 0.277
## 48 f2 ~~ f2 0.479 0.029 16.423 0 0.422 0.536
## 49 f3 ~~ f3 0.141 0.014 10.095 0 0.113 0.168
## 50 f1 ~~ f2 0.284 0.014 19.743 0 0.256 0.312
## 51 f1 ~~ f3 -0.121 0.008 -14.431 0 -0.137 -0.104
## 52 f2 ~~ f3 -0.169 0.012 -14.586 0 -0.192 -0.147modificationIndices(fit1, sort.= TRUE, minimum.value = 3.84)## lhs op rhs mi epc sepc.lv sepc.all sepc.nox
## 162 q8 ~~ q11 468.641 0.219 0.219 0.484 0.484
## 164 q8 ~~ q17 315.539 0.180 0.180 0.397 0.397
## 232 q20 ~~ q21 259.000 0.247 0.247 0.339 0.339
## 62 f1 =~ q3 235.566 -1.361 -0.675 -0.628 -0.628
## 182 q11 ~~ q17 199.518 0.138 0.138 0.326 0.326
## 53 f1 =~ q6 124.411 -1.054 -0.523 -0.466 -0.466
## 70 f2 =~ q8 97.598 -0.526 -0.364 -0.417 -0.417
## 63 f1 =~ q9 96.200 0.827 0.410 0.325 0.325
## 85 f3 =~ q8 94.856 0.692 0.259 0.297 0.297
## 103 q1 ~~ q16 94.066 0.098 0.098 0.215 0.215
## 77 f2 =~ q3 91.989 -0.615 -0.425 -0.396 -0.396
## 208 q16 ~~ q15 85.610 0.114 0.114 0.198 0.198
## 344 q3 ~~ q22 78.907 -0.177 -0.177 -0.251 -0.251
## 261 q6 ~~ q7 77.998 0.140 0.140 0.199 0.199
## 248 q21 ~~ q7 74.721 0.116 0.116 0.189 0.189
## 89 f3 =~ q20 74.364 -0.797 -0.299 -0.288 -0.288
## 75 f2 =~ q21 73.449 0.509 0.352 0.358 0.358
## 263 q6 ~~ q12 73.231 -0.116 -0.116 -0.189 -0.189
## 338 q2 ~~ q9 70.694 0.159 0.159 0.184 0.184
## 351 q22 ~~ q23 69.243 0.167 0.167 0.171 0.171
## 181 q11 ~~ q16 68.589 -0.085 -0.085 -0.190 -0.190
## 91 f3 =~ q6 67.464 0.745 0.279 0.249 0.249
## 90 f3 =~ q21 58.538 -0.605 -0.227 -0.231 -0.231
## 184 q11 ~~ q21 56.542 -0.085 -0.085 -0.169 -0.169
## 59 f1 =~ q15 55.118 0.650 0.323 0.320 0.320
## 163 q8 ~~ q16 54.881 -0.079 -0.079 -0.165 -0.165
## 177 q8 ~~ q9 54.428 0.123 0.123 0.156 0.156
## 166 q8 ~~ q21 52.916 -0.085 -0.085 -0.159 -0.159
## 88 f3 =~ q17 49.630 0.483 0.181 0.205 0.205
## 324 q15 ~~ q18 47.755 -0.094 -0.094 -0.157 -0.157
## 78 f2 =~ q9 47.520 0.418 0.289 0.229 0.229
## 131 q4 ~~ q12 47.105 0.076 0.076 0.148 0.148
## 267 q6 ~~ q18 44.367 0.097 0.097 0.154 0.154
## 61 f1 =~ q2 43.610 0.375 0.186 0.219 0.219
## 87 f3 =~ q16 41.517 -0.463 -0.174 -0.189 -0.189
## 100 q1 ~~ q5 40.865 0.073 0.073 0.137 0.137
## 202 q16 ~~ q6 38.661 -0.082 -0.082 -0.135 -0.135
## 94 f3 =~ q12 37.813 -0.445 -0.167 -0.182 -0.182
## 143 q5 ~~ q11 37.097 -0.070 -0.070 -0.134 -0.134
## 228 q17 ~~ q9 35.865 0.095 0.095 0.129 0.129
## 216 q17 ~~ q20 35.687 -0.079 -0.079 -0.128 -0.128
## 217 q17 ~~ q21 35.414 -0.067 -0.067 -0.134 -0.134
## 233 q20 ~~ q6 33.202 -0.100 -0.100 -0.120 -0.120
## 347 q9 ~~ q22 33.132 0.135 0.135 0.124 0.124
## 123 q4 ~~ q11 32.763 -0.062 -0.062 -0.129 -0.129
## 76 f2 =~ q2 31.783 0.231 0.160 0.188 0.188
## 269 q6 ~~ q3 29.176 0.086 0.086 0.133 0.133
## 236 q20 ~~ q12 28.946 0.074 0.074 0.113 0.113
## 165 q8 ~~ q20 28.668 -0.074 -0.074 -0.112 -0.112
## 99 q1 ~~ q4 28.323 0.057 0.057 0.117 0.117
## 56 f1 =~ q12 27.579 0.397 0.197 0.215 0.215
## 121 q4 ~~ q5 26.444 0.066 0.066 0.111 0.111
## 250 q21 ~~ q12 25.806 0.059 0.059 0.110 0.110
## 145 q5 ~~ q17 25.714 -0.059 -0.059 -0.111 -0.111
## 170 q8 ~~ q12 23.710 -0.051 -0.051 -0.105 -0.105
## 189 q11 ~~ q13 23.518 0.048 0.048 0.108 0.108
## 241 q20 ~~ q2 22.872 -0.071 -0.071 -0.099 -0.099
## 244 q20 ~~ q19 22.662 -0.090 -0.090 -0.100 -0.100
## 142 q5 ~~ q8 21.998 -0.057 -0.057 -0.100 -0.100
## 72 f2 =~ q16 21.749 0.253 0.175 0.191 0.191
## 102 q1 ~~ q11 21.737 -0.045 -0.045 -0.104 -0.104
## 262 q6 ~~ q10 20.421 0.067 0.067 0.095 0.095
## 81 f2 =~ q23 20.361 0.225 0.156 0.149 0.149
## 276 q7 ~~ q13 20.336 -0.057 -0.057 -0.105 -0.105
## 302 q12 ~~ q3 20.326 -0.057 -0.057 -0.112 -0.112
## 66 f1 =~ q23 19.951 0.310 0.154 0.147 0.147
## 125 q4 ~~ q17 19.729 -0.048 -0.048 -0.100 -0.100
## 221 q17 ~~ q12 18.426 -0.043 -0.043 -0.094 -0.094
## 137 q4 ~~ q3 18.400 -0.057 -0.057 -0.106 -0.106
## 308 q13 ~~ q15 17.832 -0.053 -0.053 -0.093 -0.093
## 348 q9 ~~ q23 17.823 0.101 0.101 0.088 0.088
## 186 q11 ~~ q7 17.794 -0.048 -0.048 -0.094 -0.094
## 97 f3 =~ q15 17.547 -0.354 -0.133 -0.131 -0.131
## 150 q5 ~~ q10 17.544 0.055 0.055 0.086 0.086
## 256 q21 ~~ q3 17.287 -0.058 -0.058 -0.103 -0.103
## 340 q2 ~~ q22 17.010 0.065 0.065 0.089 0.089
## 289 q10 ~~ q15 16.966 0.057 0.057 0.085 0.085
## 199 q16 ~~ q17 16.773 -0.042 -0.042 -0.094 -0.094
## 188 q11 ~~ q12 16.731 -0.040 -0.040 -0.090 -0.090
## 242 q20 ~~ q3 16.626 -0.067 -0.067 -0.097 -0.097
## 206 q16 ~~ q13 16.086 -0.042 -0.042 -0.089 -0.089
## 190 q11 ~~ q14 15.827 -0.043 -0.043 -0.087 -0.087
## 106 q1 ~~ q21 15.643 -0.044 -0.044 -0.086 -0.086
## 113 q1 ~~ q15 15.493 -0.046 -0.046 -0.083 -0.083
## 311 q13 ~~ q3 15.451 0.049 0.049 0.099 0.099
## 264 q6 ~~ q13 14.998 0.053 0.053 0.087 0.087
## 297 q12 ~~ q13 14.439 0.041 0.041 0.087 0.087
## 175 q8 ~~ q2 14.190 0.042 0.042 0.080 0.080
## 343 q3 ~~ q19 13.973 -0.083 -0.083 -0.119 -0.119
## 174 q8 ~~ q18 13.830 -0.041 -0.041 -0.082 -0.082
## 253 q21 ~~ q15 13.792 -0.051 -0.051 -0.079 -0.079
## 207 q16 ~~ q14 13.572 0.041 0.041 0.081 0.081
## 249 q21 ~~ q10 13.555 -0.047 -0.047 -0.077 -0.077
## 235 q20 ~~ q10 13.542 -0.056 -0.056 -0.075 -0.075
## 120 q1 ~~ q23 13.184 0.050 0.050 0.075 0.075
## 204 q16 ~~ q10 13.125 0.042 0.042 0.076 0.076
## 107 q1 ~~ q6 12.288 -0.044 -0.044 -0.075 -0.075
## 224 q17 ~~ q15 12.091 0.041 0.041 0.075 0.075
## 104 q1 ~~ q17 12.018 -0.033 -0.033 -0.077 -0.077
## 327 q15 ~~ q9 11.867 -0.069 -0.069 -0.072 -0.072
## 226 q17 ~~ q2 11.836 0.037 0.037 0.074 0.074
## 185 q11 ~~ q6 11.338 0.042 0.042 0.073 0.073
## 222 q17 ~~ q13 11.250 0.033 0.033 0.075 0.075
## 229 q17 ~~ q19 11.170 0.045 0.045 0.073 0.073
## 144 q5 ~~ q16 11.099 0.041 0.041 0.073 0.073
## 60 f1 =~ q18 10.952 -0.275 -0.137 -0.130 -0.130
## 212 q16 ~~ q9 10.499 -0.054 -0.054 -0.070 -0.070
## 268 q6 ~~ q2 10.167 0.046 0.046 0.068 0.068
## 247 q21 ~~ q6 10.081 -0.046 -0.046 -0.068 -0.068
## 108 q1 ~~ q7 9.960 -0.036 -0.036 -0.069 -0.069
## 211 q16 ~~ q3 9.706 -0.039 -0.039 -0.078 -0.078
## 65 f1 =~ q22 9.637 0.215 0.106 0.102 0.102
## 129 q4 ~~ q7 9.486 0.040 0.040 0.067 0.067
## 187 q11 ~~ q10 9.293 0.033 0.033 0.064 0.064
## 275 q7 ~~ q12 9.273 -0.038 -0.038 -0.070 -0.070
## 255 q21 ~~ q2 9.212 -0.038 -0.038 -0.065 -0.065
## 140 q4 ~~ q22 9.123 0.045 0.045 0.064 0.064
## 68 f2 =~ q4 9.063 0.172 0.119 0.125 0.125
## 73 f2 =~ q17 8.955 -0.155 -0.107 -0.122 -0.122
## 258 q21 ~~ q19 8.952 -0.047 -0.047 -0.064 -0.064
## 349 q19 ~~ q22 8.363 0.058 0.058 0.064 0.064
## 342 q3 ~~ q9 8.195 -0.070 -0.070 -0.085 -0.085
## 299 q12 ~~ q15 8.076 -0.036 -0.036 -0.062 -0.062
## 243 q20 ~~ q9 8.023 -0.063 -0.063 -0.058 -0.058
## 272 q6 ~~ q22 7.980 -0.050 -0.050 -0.059 -0.059
## 132 q4 ~~ q13 7.933 -0.031 -0.031 -0.062 -0.062
## 309 q13 ~~ q18 7.670 0.033 0.033 0.067 0.067
## 122 q4 ~~ q8 7.624 -0.031 -0.031 -0.061 -0.061
## 95 f3 =~ q13 7.585 0.200 0.075 0.079 0.079
## 293 q10 ~~ q9 7.581 -0.051 -0.051 -0.057 -0.057
## 67 f2 =~ q1 7.503 -0.139 -0.096 -0.116 -0.116
## 282 q7 ~~ q9 7.301 0.053 0.053 0.058 0.058
## 101 q1 ~~ q8 7.241 -0.027 -0.027 -0.058 -0.058
## 254 q21 ~~ q18 7.080 0.033 0.033 0.059 0.059
## 191 q11 ~~ q15 7.032 0.031 0.031 0.057 0.057
## 151 q5 ~~ q12 6.623 0.031 0.031 0.055 0.055
## 314 q13 ~~ q22 6.589 -0.036 -0.036 -0.055 -0.055
## 286 q10 ~~ q12 6.523 -0.030 -0.030 -0.054 -0.054
## 320 q14 ~~ q9 6.479 0.046 0.046 0.055 0.055
## 231 q17 ~~ q23 6.289 -0.035 -0.035 -0.053 -0.053
## 116 q1 ~~ q3 6.246 -0.030 -0.030 -0.061 -0.061
## 71 f2 =~ q11 5.784 -0.124 -0.086 -0.098 -0.098
## 306 q12 ~~ q23 5.732 0.035 0.035 0.050 0.050
## 136 q4 ~~ q2 5.670 0.029 0.029 0.051 0.051
## 69 f2 =~ q5 5.434 0.141 0.098 0.101 0.101
## 141 q4 ~~ q23 5.368 0.036 0.036 0.048 0.048
## 105 q1 ~~ q20 5.324 -0.030 -0.030 -0.048 -0.048
## 292 q10 ~~ q3 5.248 0.032 0.032 0.055 0.055
## 345 q3 ~~ q23 5.119 -0.044 -0.044 -0.059 -0.059
## 169 q8 ~~ q10 5.079 -0.026 -0.026 -0.047 -0.047
## 319 q14 ~~ q3 5.019 -0.031 -0.031 -0.056 -0.056
## 79 f2 =~ q19 4.884 -0.119 -0.082 -0.075 -0.075
## 259 q21 ~~ q22 4.787 0.034 0.034 0.046 0.046
## 252 q21 ~~ q14 4.736 0.027 0.027 0.047 0.047
## 139 q4 ~~ q19 4.731 0.033 0.033 0.047 0.047
## 237 q20 ~~ q13 4.665 -0.030 -0.030 -0.046 -0.046
## 288 q10 ~~ q14 4.591 -0.027 -0.027 -0.046 -0.046
## 173 q8 ~~ q15 4.580 0.026 0.026 0.045 0.045
## 111 q1 ~~ q13 4.485 0.021 0.021 0.046 0.046
## 115 q1 ~~ q2 4.445 0.023 0.023 0.045 0.045
## 283 q7 ~~ q19 4.407 -0.035 -0.035 -0.046 -0.046
## 192 q11 ~~ q18 4.281 -0.022 -0.022 -0.047 -0.047
## 346 q9 ~~ q19 4.122 0.049 0.049 0.046 0.046lavTestLRT(fit, fit1)## Chi-Squared Difference Test
##
## Df AIC BIC Chisq Chisq diff Df diff Pr(>Chisq)
## fit1 227 149577 149864 3299.6
## fit 230 150747 151016 4475.7 1176.1 3 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1Test<- '
f1 =~ q1 + q4 + q5 + q8 + q11 + q16 + q17 + q20 + q21
f2 =~ q6 + q7 + q10 + q12 + q13 + q14 + q15 + q18
f3 =~ q2 + q3 + q9 + q19 + q22 + q23
q8 ~~ q11
'
fit2 <- cfa(Test, data = FullData)
summary(fit2, fit.measures = TRUE)## lavaan 0.6-5 ended normally after 52 iterations
##
## Estimator ML
## Optimization method NLMINB
## Number of free parameters 50
##
## Number of observations 2569
##
## Model Test User Model:
##
## Test statistic 2815.403
## Degrees of freedom 226
## P-value (Chi-square) 0.000
##
## Model Test Baseline Model:
##
## Test statistic 19390.907
## Degrees of freedom 253
## P-value 0.000
##
## User Model versus Baseline Model:
##
## Comparative Fit Index (CFI) 0.865
## Tucker-Lewis Index (TLI) 0.849
##
## Loglikelihood and Information Criteria:
##
## Loglikelihood user model (H0) -74497.334
## Loglikelihood unrestricted model (H1) -73089.632
##
## Akaike (AIC) 149094.668
## Bayesian (BIC) 149387.231
## Sample-size adjusted Bayesian (BIC) 149228.367
##
## Root Mean Square Error of Approximation:
##
## RMSEA 0.067
## 90 Percent confidence interval - lower 0.065
## 90 Percent confidence interval - upper 0.069
## P-value RMSEA <= 0.05 0.000
##
## Standardized Root Mean Square Residual:
##
## SRMR 0.052
##
## Parameter Estimates:
##
## Information Expected
## Information saturated (h1) model Structured
## Standard errors Standard
##
## Latent Variables:
## Estimate Std.Err z-value P(>|z|)
## f1 =~
## q1 1.000
## q4 1.207 0.046 26.166 0.000
## q5 1.061 0.045 23.386 0.000
## q8 0.934 0.041 22.822 0.000
## q11 1.095 0.043 25.699 0.000
## q16 1.241 0.045 27.388 0.000
## q17 1.132 0.043 26.295 0.000
## q20 0.870 0.047 18.651 0.000
## q21 1.261 0.048 26.289 0.000
## f2 =~
## q6 1.000
## q7 1.110 0.039 28.345 0.000
## q10 0.542 0.028 19.035 0.000
## q12 0.873 0.032 27.215 0.000
## q13 0.954 0.034 28.328 0.000
## q14 0.960 0.035 27.393 0.000
## q15 0.816 0.034 23.927 0.000
## q18 1.126 0.038 29.585 0.000
## f3 =~
## q2 1.000
## q3 2.122 0.116 18.288 0.000
## q9 1.486 0.100 14.814 0.000
## q19 1.480 0.093 15.928 0.000
## q22 1.057 0.078 13.516 0.000
## q23 0.650 0.070 9.347 0.000
##
## Covariances:
## Estimate Std.Err z-value P(>|z|)
## .q8 ~~
## .q11 0.226 0.012 18.670 0.000
## f1 ~~
## f2 0.292 0.015 19.898 0.000
## f3 -0.127 0.009 -14.618 0.000
## f2 ~~
## f3 -0.168 0.012 -14.559 0.000
##
## Variances:
## Estimate Std.Err z-value P(>|z|)
## .q1 0.433 0.013 32.738 0.000
## .q4 0.533 0.017 32.159 0.000
## .q5 0.647 0.019 33.504 0.000
## .q8 0.541 0.016 33.566 0.000
## .q11 0.473 0.015 32.382 0.000
## .q16 0.452 0.014 31.236 0.000
## .q17 0.458 0.014 32.075 0.000
## .q20 0.882 0.025 34.692 0.000
## .q21 0.569 0.018 32.079 0.000
## .q6 0.786 0.024 32.936 0.000
## .q7 0.632 0.020 31.373 0.000
## .q10 0.630 0.018 34.780 0.000
## .q12 0.478 0.015 32.194 0.000
## .q13 0.468 0.015 31.387 0.000
## .q14 0.561 0.017 32.080 0.000
## .q15 0.703 0.021 33.680 0.000
## .q18 0.509 0.017 30.119 0.000
## .q2 0.585 0.018 32.936 0.000
## .q3 0.527 0.025 20.931 0.000
## .q9 1.285 0.039 32.921 0.000
## .q19 0.908 0.029 31.707 0.000
## .q22 0.927 0.027 33.805 0.000
## .q23 1.029 0.029 35.153 0.000
## f1 0.252 0.016 15.976 0.000
## f2 0.474 0.029 16.316 0.000
## f3 0.139 0.014 10.076 0.000parameterEstimates(fit2)## lhs op rhs est se z pvalue ci.lower ci.upper
## 1 f1 =~ q1 1.000 0.000 NA NA 1.000 1.000
## 2 f1 =~ q4 1.207 0.046 26.166 0 1.117 1.297
## 3 f1 =~ q5 1.061 0.045 23.386 0 0.973 1.150
## 4 f1 =~ q8 0.934 0.041 22.822 0 0.854 1.015
## 5 f1 =~ q11 1.095 0.043 25.699 0 1.011 1.178
## 6 f1 =~ q16 1.241 0.045 27.388 0 1.152 1.330
## 7 f1 =~ q17 1.132 0.043 26.295 0 1.047 1.216
## 8 f1 =~ q20 0.870 0.047 18.651 0 0.779 0.962
## 9 f1 =~ q21 1.261 0.048 26.289 0 1.167 1.355
## 10 f2 =~ q6 1.000 0.000 NA NA 1.000 1.000
## 11 f2 =~ q7 1.110 0.039 28.345 0 1.033 1.187
## 12 f2 =~ q10 0.542 0.028 19.035 0 0.486 0.598
## 13 f2 =~ q12 0.873 0.032 27.215 0 0.811 0.936
## 14 f2 =~ q13 0.954 0.034 28.328 0 0.888 1.020
## 15 f2 =~ q14 0.960 0.035 27.393 0 0.891 1.029
## 16 f2 =~ q15 0.816 0.034 23.927 0 0.750 0.883
## 17 f2 =~ q18 1.126 0.038 29.585 0 1.051 1.201
## 18 f3 =~ q2 1.000 0.000 NA NA 1.000 1.000
## 19 f3 =~ q3 2.122 0.116 18.288 0 1.895 2.350
## 20 f3 =~ q9 1.486 0.100 14.814 0 1.289 1.683
## 21 f3 =~ q19 1.480 0.093 15.928 0 1.298 1.662
## 22 f3 =~ q22 1.057 0.078 13.516 0 0.904 1.210
## 23 f3 =~ q23 0.650 0.070 9.347 0 0.514 0.787
## 24 q8 ~~ q11 0.226 0.012 18.670 0 0.202 0.249
## 25 q1 ~~ q1 0.433 0.013 32.738 0 0.407 0.459
## 26 q4 ~~ q4 0.533 0.017 32.159 0 0.501 0.566
## 27 q5 ~~ q5 0.647 0.019 33.504 0 0.609 0.685
## 28 q8 ~~ q8 0.541 0.016 33.566 0 0.510 0.573
## 29 q11 ~~ q11 0.473 0.015 32.382 0 0.445 0.502
## 30 q16 ~~ q16 0.452 0.014 31.236 0 0.423 0.480
## 31 q17 ~~ q17 0.458 0.014 32.075 0 0.430 0.486
## 32 q20 ~~ q20 0.882 0.025 34.692 0 0.832 0.931
## 33 q21 ~~ q21 0.569 0.018 32.079 0 0.534 0.604
## 34 q6 ~~ q6 0.786 0.024 32.936 0 0.739 0.833
## 35 q7 ~~ q7 0.632 0.020 31.373 0 0.592 0.671
## 36 q10 ~~ q10 0.630 0.018 34.780 0 0.595 0.666
## 37 q12 ~~ q12 0.478 0.015 32.194 0 0.449 0.508
## 38 q13 ~~ q13 0.468 0.015 31.387 0 0.439 0.497
## 39 q14 ~~ q14 0.561 0.017 32.080 0 0.526 0.595
## 40 q15 ~~ q15 0.703 0.021 33.680 0 0.662 0.744
## 41 q18 ~~ q18 0.509 0.017 30.119 0 0.476 0.542
## 42 q2 ~~ q2 0.585 0.018 32.936 0 0.550 0.620
## 43 q3 ~~ q3 0.527 0.025 20.931 0 0.478 0.577
## 44 q9 ~~ q9 1.285 0.039 32.921 0 1.209 1.362
## 45 q19 ~~ q19 0.908 0.029 31.707 0 0.852 0.964
## 46 q22 ~~ q22 0.927 0.027 33.805 0 0.873 0.981
## 47 q23 ~~ q23 1.029 0.029 35.153 0 0.972 1.087
## 48 f1 ~~ f1 0.252 0.016 15.976 0 0.221 0.283
## 49 f2 ~~ f2 0.474 0.029 16.316 0 0.417 0.530
## 50 f3 ~~ f3 0.139 0.014 10.076 0 0.112 0.166
## 51 f1 ~~ f2 0.292 0.015 19.898 0 0.263 0.321
## 52 f1 ~~ f3 -0.127 0.009 -14.618 0 -0.144 -0.110
## 53 f2 ~~ f3 -0.168 0.012 -14.559 0 -0.190 -0.145modificationIndices(fit2, sort.= TRUE, minimum.value = 3.84)## lhs op rhs mi epc sepc.lv sepc.all sepc.nox
## 232 q20 ~~ q21 239.850 0.234 0.234 0.330 0.330
## 63 f1 =~ q3 238.361 -1.493 -0.749 -0.697 -0.697
## 164 q8 ~~ q17 213.441 0.138 0.138 0.277 0.277
## 54 f1 =~ q6 151.651 -1.276 -0.640 -0.571 -0.571
## 208 q16 ~~ q15 88.944 0.115 0.115 0.204 0.204
## 64 f1 =~ q9 88.617 0.837 0.420 0.333 0.333
## 78 f2 =~ q3 84.046 -0.596 -0.410 -0.382 -0.382
## 261 q6 ~~ q7 80.611 0.142 0.142 0.202 0.202
## 182 q11 ~~ q17 77.252 0.079 0.079 0.171 0.171
## 344 q3 ~~ q22 77.146 -0.174 -0.174 -0.248 -0.248
## 104 q1 ~~ q16 75.351 0.087 0.087 0.196 0.196
## 263 q6 ~~ q12 72.850 -0.115 -0.115 -0.188 -0.188
## 338 q2 ~~ q9 72.054 0.160 0.160 0.185 0.185
## 351 q22 ~~ q23 70.889 0.169 0.169 0.173 0.173
## 92 f3 =~ q6 69.813 0.763 0.285 0.254 0.254
## 248 q21 ~~ q7 67.728 0.109 0.109 0.182 0.182
## 90 f3 =~ q20 60.721 -0.762 -0.284 -0.275 -0.275
## 177 q8 ~~ q9 56.702 0.115 0.115 0.138 0.138
## 86 f3 =~ q8 54.359 0.501 0.187 0.214 0.214
## 60 f1 =~ q15 51.644 0.689 0.346 0.343 0.343
## 267 q6 ~~ q18 47.344 0.100 0.100 0.159 0.159
## 324 q15 ~~ q18 47.136 -0.093 -0.093 -0.155 -0.155
## 79 f2 =~ q9 44.967 0.410 0.282 0.223 0.223
## 62 f1 =~ q2 43.677 0.396 0.199 0.234 0.234
## 57 f1 =~ q12 43.180 0.544 0.273 0.298 0.298
## 202 q16 ~~ q6 43.067 -0.085 -0.085 -0.143 -0.143
## 228 q17 ~~ q9 40.729 0.105 0.105 0.137 0.137
## 89 f3 =~ q17 39.527 0.473 0.176 0.200 0.200
## 95 f3 =~ q12 38.759 -0.451 -0.168 -0.184 -0.184
## 132 q4 ~~ q12 38.127 0.068 0.068 0.134 0.134
## 347 q9 ~~ q22 35.003 0.138 0.138 0.127 0.127
## 233 q20 ~~ q6 34.656 -0.102 -0.102 -0.122 -0.122
## 76 f2 =~ q21 34.333 0.389 0.268 0.272 0.272
## 71 f2 =~ q8 33.079 -0.305 -0.210 -0.241 -0.241
## 101 q1 ~~ q5 31.897 0.065 0.065 0.122 0.122
## 107 q1 ~~ q21 31.129 -0.062 -0.062 -0.124 -0.124
## 91 f3 =~ q21 30.653 -0.464 -0.173 -0.176 -0.176
## 77 f2 =~ q2 30.612 0.228 0.157 0.184 0.184
## 216 q17 ~~ q20 27.410 -0.071 -0.071 -0.112 -0.112
## 189 q11 ~~ q13 27.214 0.047 0.047 0.100 0.100
## 269 q6 ~~ q3 26.612 0.082 0.082 0.127 0.127
## 68 f2 =~ q1 25.552 -0.286 -0.197 -0.238 -0.238
## 236 q20 ~~ q12 24.295 0.067 0.067 0.104 0.104
## 67 f1 =~ q23 24.008 0.357 0.179 0.172 0.172
## 217 q17 ~~ q21 23.931 -0.056 -0.056 -0.110 -0.110
## 262 q6 ~~ q10 21.373 0.068 0.068 0.097 0.097
## 241 q20 ~~ q2 21.232 -0.068 -0.068 -0.095 -0.095
## 100 q1 ~~ q4 21.023 0.049 0.049 0.102 0.102
## 244 q20 ~~ q19 20.891 -0.086 -0.086 -0.096 -0.096
## 82 f2 =~ q23 20.175 0.225 0.155 0.149 0.149
## 187 q11 ~~ q10 19.191 0.044 0.044 0.080 0.080
## 276 q7 ~~ q13 18.946 -0.055 -0.055 -0.101 -0.101
## 348 q9 ~~ q23 18.736 0.104 0.104 0.090 0.090
## 250 q21 ~~ q12 18.537 0.049 0.049 0.094 0.094
## 340 q2 ~~ q22 18.372 0.068 0.068 0.092 0.092
## 264 q6 ~~ q13 17.969 0.058 0.058 0.095 0.095
## 98 f3 =~ q15 17.938 -0.359 -0.134 -0.133 -0.133
## 206 q16 ~~ q13 17.871 -0.043 -0.043 -0.094 -0.094
## 122 q4 ~~ q5 17.758 0.054 0.054 0.092 0.092
## 88 f3 =~ q16 17.591 -0.320 -0.119 -0.130 -0.130
## 302 q12 ~~ q3 17.553 -0.052 -0.052 -0.104 -0.104
## 289 q10 ~~ q15 17.217 0.057 0.057 0.086 0.086
## 185 q11 ~~ q6 16.946 0.047 0.047 0.077 0.077
## 221 q17 ~~ q12 16.845 -0.042 -0.042 -0.089 -0.089
## 151 q5 ~~ q10 16.632 0.054 0.054 0.084 0.084
## 181 q11 ~~ q16 16.393 -0.037 -0.037 -0.080 -0.080
## 224 q17 ~~ q15 16.293 0.049 0.049 0.086 0.086
## 308 q13 ~~ q15 16.039 -0.051 -0.051 -0.088 -0.088
## 114 q1 ~~ q15 15.896 -0.047 -0.047 -0.084 -0.084
## 249 q21 ~~ q10 15.881 -0.050 -0.050 -0.084 -0.084
## 343 q3 ~~ q19 15.844 -0.088 -0.088 -0.128 -0.128
## 222 q17 ~~ q13 15.816 0.041 0.041 0.088 0.088
## 175 q8 ~~ q2 15.695 0.041 0.041 0.073 0.073
## 146 q5 ~~ q17 15.450 -0.047 -0.047 -0.086 -0.086
## 253 q21 ~~ q15 15.201 -0.053 -0.053 -0.083 -0.083
## 72 f2 =~ q11 15.175 0.200 0.137 0.156 0.156
## 311 q13 ~~ q3 14.875 0.048 0.048 0.097 0.097
## 235 q20 ~~ q10 14.684 -0.058 -0.058 -0.078 -0.078
## 138 q4 ~~ q3 14.186 -0.050 -0.050 -0.094 -0.094
## 109 q1 ~~ q7 14.049 -0.043 -0.043 -0.082 -0.082
## 297 q12 ~~ q13 13.901 0.040 0.040 0.085 0.085
## 108 q1 ~~ q6 13.637 -0.046 -0.046 -0.079 -0.079
## 242 q20 ~~ q3 13.546 -0.060 -0.060 -0.088 -0.088
## 165 q8 ~~ q20 13.338 -0.046 -0.046 -0.066 -0.066
## 121 q1 ~~ q23 13.234 0.050 0.050 0.075 0.075
## 226 q17 ~~ q2 13.080 0.040 0.040 0.077 0.077
## 256 q21 ~~ q3 12.721 -0.049 -0.049 -0.089 -0.089
## 184 q11 ~~ q21 12.522 -0.036 -0.036 -0.069 -0.069
## 204 q16 ~~ q10 12.222 0.040 0.040 0.074 0.074
## 83 f3 =~ q1 11.753 0.246 0.092 0.111 0.111
## 327 q15 ~~ q9 11.693 -0.068 -0.068 -0.072 -0.072
## 169 q8 ~~ q10 11.504 -0.036 -0.036 -0.061 -0.061
## 275 q7 ~~ q12 11.449 -0.042 -0.042 -0.077 -0.077
## 247 q21 ~~ q6 11.422 -0.049 -0.049 -0.073 -0.073
## 186 q11 ~~ q7 11.232 -0.035 -0.035 -0.064 -0.064
## 124 q4 ~~ q11 11.182 -0.033 -0.033 -0.065 -0.065
## 87 f3 =~ q11 10.904 -0.216 -0.080 -0.091 -0.091
## 144 q5 ~~ q11 10.683 -0.034 -0.034 -0.062 -0.062
## 141 q4 ~~ q22 10.332 0.048 0.048 0.068 0.068
## 229 q17 ~~ q19 10.132 0.045 0.045 0.069 0.069
## 190 q11 ~~ q14 10.100 -0.031 -0.031 -0.060 -0.060
## 61 f1 =~ q18 10.007 -0.289 -0.145 -0.137 -0.137
## 66 f1 =~ q22 9.934 0.229 0.115 0.111 0.111
## 106 q1 ~~ q20 9.682 -0.041 -0.041 -0.066 -0.066
## 349 q19 ~~ q22 9.547 0.062 0.062 0.068 0.068
## 207 q16 ~~ q14 9.415 0.034 0.034 0.068 0.068
## 268 q6 ~~ q2 9.412 0.044 0.044 0.065 0.065
## 342 q3 ~~ q9 8.859 -0.073 -0.073 -0.089 -0.089
## 75 f2 =~ q20 8.855 -0.226 -0.155 -0.150 -0.150
## 299 q12 ~~ q15 8.853 -0.037 -0.037 -0.065 -0.065
## 272 q6 ~~ q22 8.794 -0.053 -0.053 -0.062 -0.062
## 309 q13 ~~ q18 8.779 0.035 0.035 0.071 0.071
## 133 q4 ~~ q13 8.124 -0.031 -0.031 -0.063 -0.063
## 137 q4 ~~ q2 7.944 0.034 0.034 0.060 0.060
## 199 q16 ~~ q17 7.796 -0.029 -0.029 -0.064 -0.064
## 96 f3 =~ q13 7.760 0.204 0.076 0.080 0.080
## 255 q21 ~~ q2 7.610 -0.034 -0.034 -0.059 -0.059
## 282 q7 ~~ q9 7.544 0.053 0.053 0.059 0.059
## 293 q10 ~~ q9 7.531 -0.051 -0.051 -0.057 -0.057
## 126 q4 ~~ q17 7.379 -0.030 -0.030 -0.061 -0.061
## 170 q8 ~~ q12 7.341 -0.026 -0.026 -0.051 -0.051
## 314 q13 ~~ q22 7.319 -0.038 -0.038 -0.058 -0.058
## 258 q21 ~~ q19 7.267 -0.042 -0.042 -0.059 -0.059
## 286 q10 ~~ q12 7.087 -0.031 -0.031 -0.057 -0.057
## 212 q16 ~~ q9 7.070 -0.044 -0.044 -0.058 -0.058
## 320 q14 ~~ q9 6.674 0.047 0.047 0.055 0.055
## 140 q4 ~~ q19 6.623 0.039 0.039 0.056 0.056
## 130 q4 ~~ q7 6.553 0.033 0.033 0.056 0.056
## 116 q1 ~~ q2 6.449 0.027 0.027 0.054 0.054
## 243 q20 ~~ q9 6.415 -0.056 -0.056 -0.052 -0.052
## 197 q11 ~~ q22 6.002 -0.030 -0.030 -0.045 -0.045
## 231 q17 ~~ q23 5.848 -0.035 -0.035 -0.050 -0.050
## 174 q8 ~~ q18 5.678 -0.024 -0.024 -0.046 -0.046
## 211 q16 ~~ q3 5.670 -0.030 -0.030 -0.061 -0.061
## 259 q21 ~~ q22 5.570 0.036 0.036 0.050 0.050
## 55 f1 =~ q7 5.549 0.230 0.115 0.105 0.105
## 147 q5 ~~ q20 5.534 -0.037 -0.037 -0.049 -0.049
## 292 q10 ~~ q3 5.465 0.032 0.032 0.056 0.056
## 306 q12 ~~ q23 5.459 0.034 0.034 0.049 0.049
## 191 q11 ~~ q15 5.385 0.025 0.025 0.043 0.043
## 103 q1 ~~ q11 5.381 -0.020 -0.020 -0.044 -0.044
## 179 q8 ~~ q22 5.338 0.030 0.030 0.042 0.042
## 142 q4 ~~ q23 5.257 0.035 0.035 0.048 0.048
## 203 q16 ~~ q7 5.196 -0.027 -0.027 -0.051 -0.051
## 80 f2 =~ q19 5.195 -0.123 -0.085 -0.077 -0.077
## 195 q11 ~~ q9 5.176 -0.033 -0.033 -0.042 -0.042
## 237 q20 ~~ q13 5.051 -0.031 -0.031 -0.048 -0.048
## 166 q8 ~~ q21 4.892 -0.023 -0.023 -0.042 -0.042
## 345 q3 ~~ q23 4.700 -0.042 -0.042 -0.057 -0.057
## 346 q9 ~~ q19 4.688 0.052 0.052 0.048 0.048
## 288 q10 ~~ q14 4.677 -0.027 -0.027 -0.046 -0.046
## 167 q8 ~~ q6 4.565 -0.026 -0.026 -0.039 -0.039
## 127 q4 ~~ q20 4.341 -0.030 -0.030 -0.044 -0.044
## 112 q1 ~~ q13 4.324 0.020 0.020 0.045 0.045
## 319 q14 ~~ q3 4.323 -0.028 -0.028 -0.052 -0.052
## 283 q7 ~~ q19 4.293 -0.034 -0.034 -0.045 -0.045
## 118 q1 ~~ q9 4.041 0.032 0.032 0.043 0.043
## 254 q21 ~~ q18 3.879 0.024 0.024 0.044 0.044lavTestLRT(fit1, fit2)## Chi-Squared Difference Test
##
## Df AIC BIC Chisq Chisq diff Df diff Pr(>Chisq)
## fit2 226 149095 149387 2815.4
## fit1 227 149577 149864 3299.6 484.22 1 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1Test1<- '
f1 =~ q1 + q4 + q5 + q8 + q11 + q16 + q17 + q20 + q21
f2 =~ q6 + q7 + q10 + q12 + q13 + q14 + q15 + q18
f3 =~ q2 + q3 + q9 + q19 + q22 + q23
q8 ~~ q11
q8 ~~ q17
q20 ~~ q21
f1 =~ q3
'
fit3 <- cfa(Test1, data = FullData)
summary(fit3, fit.measures = TRUE)## lavaan 0.6-5 ended normally after 53 iterations
##
## Estimator ML
## Optimization method NLMINB
## Number of free parameters 53
##
## Number of observations 2569
##
## Model Test User Model:
##
## Test statistic 2141.610
## Degrees of freedom 223
## P-value (Chi-square) 0.000
##
## Model Test Baseline Model:
##
## Test statistic 19390.907
## Degrees of freedom 253
## P-value 0.000
##
## User Model versus Baseline Model:
##
## Comparative Fit Index (CFI) 0.900
## Tucker-Lewis Index (TLI) 0.886
##
## Loglikelihood and Information Criteria:
##
## Loglikelihood user model (H0) -74160.438
## Loglikelihood unrestricted model (H1) -73089.632
##
## Akaike (AIC) 148426.875
## Bayesian (BIC) 148736.992
## Sample-size adjusted Bayesian (BIC) 148568.597
##
## Root Mean Square Error of Approximation:
##
## RMSEA 0.058
## 90 Percent confidence interval - lower 0.056
## 90 Percent confidence interval - upper 0.060
## P-value RMSEA <= 0.05 0.000
##
## Standardized Root Mean Square Residual:
##
## SRMR 0.045
##
## Parameter Estimates:
##
## Information Expected
## Information saturated (h1) model Structured
## Standard errors Standard
##
## Latent Variables:
## Estimate Std.Err z-value P(>|z|)
## f1 =~
## q1 1.000
## q4 1.205 0.045 26.488 0.000
## q5 1.059 0.045 23.643 0.000
## q8 0.894 0.039 22.671 0.000
## q11 1.084 0.042 25.876 0.000
## q16 1.234 0.045 27.661 0.000
## q17 1.111 0.042 26.300 0.000
## q20 0.785 0.046 17.142 0.000
## q21 1.214 0.047 25.896 0.000
## f2 =~
## q6 1.000
## q7 1.108 0.039 28.359 0.000
## q10 0.544 0.028 19.124 0.000
## q12 0.872 0.032 27.235 0.000
## q13 0.953 0.034 28.368 0.000
## q14 0.958 0.035 27.414 0.000
## q15 0.816 0.034 23.962 0.000
## q18 1.125 0.038 29.623 0.000
## f3 =~
## q2 1.000
## q3 0.993 0.078 12.693 0.000
## q9 1.592 0.098 16.322 0.000
## q19 1.341 0.084 16.044 0.000
## q22 1.128 0.075 15.052 0.000
## q23 0.700 0.066 10.691 0.000
## f1 =~
## q3 -0.838 0.052 -16.114 0.000
##
## Covariances:
## Estimate Std.Err z-value P(>|z|)
## .q8 ~~
## .q11 0.190 0.011 17.356 0.000
## .q17 0.141 0.010 13.872 0.000
## .q20 ~~
## .q21 0.231 0.016 14.111 0.000
## f1 ~~
## f2 0.296 0.015 20.017 0.000
## f3 -0.100 0.008 -12.588 0.000
## f2 ~~
## f3 -0.148 0.011 -13.414 0.000
##
## Variances:
## Estimate Std.Err z-value P(>|z|)
## .q1 0.428 0.013 32.549 0.000
## .q4 0.527 0.016 31.954 0.000
## .q5 0.642 0.019 33.387 0.000
## .q8 0.520 0.015 34.121 0.000
## .q11 0.473 0.015 32.332 0.000
## .q16 0.448 0.014 31.036 0.000
## .q17 0.463 0.014 32.071 0.000
## .q20 0.914 0.026 34.817 0.000
## .q21 0.591 0.018 32.307 0.000
## .q6 0.785 0.024 32.932 0.000
## .q7 0.633 0.020 31.403 0.000
## .q10 0.629 0.018 34.771 0.000
## .q12 0.479 0.015 32.211 0.000
## .q13 0.468 0.015 31.395 0.000
## .q14 0.561 0.017 32.097 0.000
## .q15 0.702 0.021 33.681 0.000
## .q18 0.509 0.017 30.137 0.000
## .q2 0.544 0.018 29.939 0.000
## .q3 0.630 0.022 29.270 0.000
## .q9 1.136 0.040 28.647 0.000
## .q19 0.889 0.030 29.332 0.000
## .q22 0.853 0.027 31.073 0.000
## .q23 1.000 0.029 34.292 0.000
## f1 0.257 0.016 16.154 0.000
## f2 0.475 0.029 16.342 0.000
## f3 0.180 0.017 10.816 0.000parameterEstimates(fit3)## lhs op rhs est se z pvalue ci.lower ci.upper
## 1 f1 =~ q1 1.000 0.000 NA NA 1.000 1.000
## 2 f1 =~ q4 1.205 0.045 26.488 0 1.115 1.294
## 3 f1 =~ q5 1.059 0.045 23.643 0 0.971 1.147
## 4 f1 =~ q8 0.894 0.039 22.671 0 0.817 0.971
## 5 f1 =~ q11 1.084 0.042 25.876 0 1.002 1.166
## 6 f1 =~ q16 1.234 0.045 27.661 0 1.146 1.321
## 7 f1 =~ q17 1.111 0.042 26.300 0 1.029 1.194
## 8 f1 =~ q20 0.785 0.046 17.142 0 0.696 0.875
## 9 f1 =~ q21 1.214 0.047 25.896 0 1.122 1.306
## 10 f2 =~ q6 1.000 0.000 NA NA 1.000 1.000
## 11 f2 =~ q7 1.108 0.039 28.359 0 1.031 1.184
## 12 f2 =~ q10 0.544 0.028 19.124 0 0.488 0.600
## 13 f2 =~ q12 0.872 0.032 27.235 0 0.809 0.935
## 14 f2 =~ q13 0.953 0.034 28.368 0 0.887 1.019
## 15 f2 =~ q14 0.958 0.035 27.414 0 0.890 1.027
## 16 f2 =~ q15 0.816 0.034 23.962 0 0.749 0.883
## 17 f2 =~ q18 1.125 0.038 29.623 0 1.050 1.199
## 18 f3 =~ q2 1.000 0.000 NA NA 1.000 1.000
## 19 f3 =~ q3 0.993 0.078 12.693 0 0.840 1.146
## 20 f3 =~ q9 1.592 0.098 16.322 0 1.401 1.783
## 21 f3 =~ q19 1.341 0.084 16.044 0 1.177 1.505
## 22 f3 =~ q22 1.128 0.075 15.052 0 0.981 1.275
## 23 f3 =~ q23 0.700 0.066 10.691 0 0.572 0.829
## 24 q8 ~~ q11 0.190 0.011 17.356 0 0.168 0.211
## 25 q8 ~~ q17 0.141 0.010 13.872 0 0.121 0.160
## 26 q20 ~~ q21 0.231 0.016 14.111 0 0.199 0.264
## 27 f1 =~ q3 -0.838 0.052 -16.114 0 -0.939 -0.736
## 28 q1 ~~ q1 0.428 0.013 32.549 0 0.402 0.454
## 29 q4 ~~ q4 0.527 0.016 31.954 0 0.495 0.559
## 30 q5 ~~ q5 0.642 0.019 33.387 0 0.605 0.680
## 31 q8 ~~ q8 0.520 0.015 34.121 0 0.490 0.550
## 32 q11 ~~ q11 0.473 0.015 32.332 0 0.444 0.502
## 33 q16 ~~ q16 0.448 0.014 31.036 0 0.420 0.477
## 34 q17 ~~ q17 0.463 0.014 32.071 0 0.435 0.492
## 35 q20 ~~ q20 0.914 0.026 34.817 0 0.862 0.965
## 36 q21 ~~ q21 0.591 0.018 32.307 0 0.555 0.627
## 37 q6 ~~ q6 0.785 0.024 32.932 0 0.738 0.831
## 38 q7 ~~ q7 0.633 0.020 31.403 0 0.594 0.673
## 39 q10 ~~ q10 0.629 0.018 34.771 0 0.594 0.664
## 40 q12 ~~ q12 0.479 0.015 32.211 0 0.450 0.508
## 41 q13 ~~ q13 0.468 0.015 31.395 0 0.439 0.497
## 42 q14 ~~ q14 0.561 0.017 32.097 0 0.527 0.595
## 43 q15 ~~ q15 0.702 0.021 33.681 0 0.662 0.743
## 44 q18 ~~ q18 0.509 0.017 30.137 0 0.476 0.542
## 45 q2 ~~ q2 0.544 0.018 29.939 0 0.508 0.580
## 46 q3 ~~ q3 0.630 0.022 29.270 0 0.588 0.672
## 47 q9 ~~ q9 1.136 0.040 28.647 0 1.058 1.214
## 48 q19 ~~ q19 0.889 0.030 29.332 0 0.830 0.948
## 49 q22 ~~ q22 0.853 0.027 31.073 0 0.800 0.907
## 50 q23 ~~ q23 1.000 0.029 34.292 0 0.943 1.057
## 51 f1 ~~ f1 0.257 0.016 16.154 0 0.226 0.288
## 52 f2 ~~ f2 0.475 0.029 16.342 0 0.418 0.532
## 53 f3 ~~ f3 0.180 0.017 10.816 0 0.148 0.213
## 54 f1 ~~ f2 0.296 0.015 20.017 0 0.267 0.325
## 55 f1 ~~ f3 -0.100 0.008 -12.588 0 -0.115 -0.084
## 56 f2 ~~ f3 -0.148 0.011 -13.414 0 -0.170 -0.127modificationIndices(fit3, sort.= TRUE, minimum.value = 3.84)## lhs op rhs mi epc sepc.lv sepc.all sepc.nox
## 183 q11 ~~ q17 325.387 0.189 0.189 0.405 0.405
## 57 f1 =~ q6 144.061 -1.254 -0.636 -0.566 -0.566
## 209 q16 ~~ q15 90.658 0.116 0.116 0.206 0.206
## 248 q21 ~~ q7 82.439 0.114 0.114 0.187 0.187
## 261 q6 ~~ q7 80.447 0.142 0.142 0.201 0.201
## 263 q6 ~~ q12 73.195 -0.116 -0.116 -0.189 -0.189
## 67 f1 =~ q19 67.824 -0.495 -0.251 -0.228 -0.228
## 106 q1 ~~ q16 67.642 0.082 0.082 0.188 0.188
## 82 f2 =~ q19 66.696 -0.366 -0.253 -0.229 -0.229
## 78 f2 =~ q21 57.620 0.481 0.332 0.337 0.337
## 92 f3 =~ q20 52.279 -0.470 -0.199 -0.193 -0.193
## 203 q16 ~~ q6 48.480 -0.090 -0.090 -0.152 -0.152
## 324 q15 ~~ q18 47.373 -0.093 -0.093 -0.156 -0.156
## 267 q6 ~~ q18 46.697 0.100 0.100 0.158 0.158
## 63 f1 =~ q15 46.607 0.660 0.335 0.332 0.332
## 351 q22 ~~ q23 44.157 0.134 0.134 0.146 0.146
## 60 f1 =~ q12 41.750 0.539 0.273 0.298 0.298
## 134 q4 ~~ q12 38.150 0.068 0.068 0.135 0.135
## 70 f2 =~ q1 37.654 -0.354 -0.244 -0.295 -0.295
## 344 q3 ~~ q22 33.619 -0.106 -0.106 -0.144 -0.144
## 190 q11 ~~ q13 30.441 0.052 0.052 0.110 0.110
## 178 q8 ~~ q9 29.284 0.077 0.077 0.101 0.101
## 88 f3 =~ q8 29.000 0.245 0.104 0.122 0.122
## 103 q1 ~~ q5 27.169 0.059 0.059 0.113 0.113
## 182 q11 ~~ q16 26.914 -0.049 -0.049 -0.107 -0.107
## 94 f3 =~ q6 26.795 0.356 0.151 0.135 0.135
## 66 f1 =~ q9 26.153 0.356 0.180 0.143 0.143
## 269 q6 ~~ q3 25.214 0.078 0.078 0.111 0.111
## 339 q2 ~~ q19 23.384 -0.087 -0.087 -0.125 -0.125
## 233 q20 ~~ q6 22.831 -0.079 -0.079 -0.093 -0.093
## 81 f2 =~ q9 20.947 0.237 0.164 0.130 0.130
## 109 q1 ~~ q21 20.505 -0.047 -0.047 -0.094 -0.094
## 262 q6 ~~ q10 20.420 0.067 0.067 0.095 0.095
## 90 f3 =~ q16 20.291 -0.241 -0.102 -0.112 -0.112
## 338 q2 ~~ q9 19.652 0.092 0.092 0.117 0.117
## 126 q4 ~~ q11 19.472 -0.045 -0.045 -0.090 -0.090
## 185 q11 ~~ q21 19.289 -0.044 -0.044 -0.083 -0.083
## 207 q16 ~~ q13 19.121 -0.045 -0.045 -0.098 -0.098
## 146 q5 ~~ q11 19.045 -0.047 -0.047 -0.086 -0.086
## 276 q7 ~~ q13 18.406 -0.054 -0.054 -0.100 -0.100
## 236 q20 ~~ q12 17.890 0.055 0.055 0.083 0.083
## 242 q20 ~~ q3 17.729 -0.065 -0.065 -0.085 -0.085
## 264 q6 ~~ q13 17.529 0.057 0.057 0.094 0.094
## 223 q17 ~~ q13 17.345 0.040 0.040 0.087 0.087
## 311 q13 ~~ q3 17.046 0.051 0.051 0.093 0.093
## 73 f2 =~ q8 17.044 -0.217 -0.150 -0.176 -0.176
## 97 f3 =~ q12 16.759 -0.223 -0.095 -0.103 -0.103
## 148 q5 ~~ q17 16.663 -0.046 -0.046 -0.084 -0.084
## 302 q12 ~~ q3 16.660 -0.050 -0.050 -0.091 -0.091
## 289 q10 ~~ q15 16.545 0.056 0.056 0.084 0.084
## 116 q1 ~~ q15 16.331 -0.047 -0.047 -0.086 -0.086
## 308 q13 ~~ q15 16.274 -0.051 -0.051 -0.089 -0.089
## 110 q1 ~~ q6 16.170 -0.050 -0.050 -0.086 -0.086
## 102 q1 ~~ q4 15.921 0.042 0.042 0.089 0.089
## 244 q20 ~~ q19 15.652 -0.071 -0.071 -0.079 -0.079
## 188 q11 ~~ q10 15.297 0.040 0.040 0.074 0.074
## 153 q5 ~~ q10 15.203 0.051 0.051 0.080 0.080
## 111 q1 ~~ q7 15.173 -0.044 -0.044 -0.085 -0.085
## 85 f3 =~ q1 15.082 0.195 0.083 0.100 0.100
## 337 q2 ~~ q3 14.384 0.057 0.057 0.097 0.097
## 124 q4 ~~ q5 14.222 0.048 0.048 0.083 0.083
## 128 q4 ~~ q17 14.177 -0.040 -0.040 -0.080 -0.080
## 253 q21 ~~ q15 14.110 -0.048 -0.048 -0.075 -0.075
## 297 q12 ~~ q13 14.053 0.040 0.040 0.086 0.086
## 268 q6 ~~ q2 13.943 0.053 0.053 0.081 0.081
## 241 q20 ~~ q2 13.555 -0.052 -0.052 -0.073 -0.073
## 186 q11 ~~ q6 13.442 0.043 0.043 0.071 0.071
## 225 q17 ~~ q15 13.180 0.042 0.042 0.073 0.073
## 100 f3 =~ q15 13.136 -0.232 -0.099 -0.098 -0.098
## 125 q4 ~~ q8 13.120 0.035 0.035 0.068 0.068
## 91 f3 =~ q17 12.431 0.175 0.074 0.084 0.084
## 234 q20 ~~ q7 12.150 -0.053 -0.053 -0.069 -0.069
## 191 q11 ~~ q14 12.006 -0.035 -0.035 -0.068 -0.068
## 327 q15 ~~ q9 11.475 -0.066 -0.066 -0.074 -0.074
## 187 q11 ~~ q7 11.346 -0.037 -0.037 -0.067 -0.067
## 229 q17 ~~ q9 11.016 0.051 0.051 0.070 0.070
## 275 q7 ~~ q12 10.824 -0.041 -0.041 -0.075 -0.075
## 105 q1 ~~ q11 10.659 -0.029 -0.029 -0.065 -0.065
## 205 q16 ~~ q10 10.492 0.037 0.037 0.069 0.069
## 222 q17 ~~ q12 10.436 -0.031 -0.031 -0.066 -0.066
## 232 q17 ~~ q23 10.387 -0.043 -0.043 -0.064 -0.064
## 213 q16 ~~ q9 10.150 -0.051 -0.051 -0.072 -0.072
## 84 f2 =~ q23 9.855 0.129 0.089 0.085 0.085
## 93 f3 =~ q21 9.403 -0.170 -0.072 -0.073 -0.073
## 123 q1 ~~ q23 9.261 0.042 0.042 0.064 0.064
## 64 f1 =~ q18 8.932 -0.275 -0.139 -0.132 -0.132
## 250 q21 ~~ q12 8.866 0.032 0.032 0.061 0.061
## 299 q12 ~~ q15 8.848 -0.037 -0.037 -0.065 -0.065
## 309 q13 ~~ q18 8.784 0.035 0.035 0.071 0.071
## 135 q4 ~~ q13 8.749 -0.032 -0.032 -0.065 -0.065
## 208 q16 ~~ q14 8.618 0.033 0.033 0.065 0.065
## 341 q2 ~~ q23 8.607 -0.048 -0.048 -0.066 -0.066
## 282 q7 ~~ q9 8.377 0.055 0.055 0.065 0.065
## 89 f3 =~ q11 8.328 -0.138 -0.059 -0.066 -0.066
## 249 q21 ~~ q10 8.270 -0.034 -0.034 -0.057 -0.057
## 74 f2 =~ q11 8.265 0.158 0.109 0.124 0.124
## 272 q6 ~~ q22 8.059 -0.050 -0.050 -0.061 -0.061
## 184 q11 ~~ q20 7.964 0.033 0.033 0.051 0.051
## 192 q11 ~~ q15 7.852 0.031 0.031 0.054 0.054
## 320 q14 ~~ q9 7.786 0.050 0.050 0.062 0.062
## 170 q8 ~~ q10 7.782 -0.028 -0.028 -0.049 -0.049
## 286 q10 ~~ q12 7.572 -0.032 -0.032 -0.058 -0.058
## 104 q1 ~~ q8 7.450 0.024 0.024 0.050 0.050
## 69 f1 =~ q23 7.447 0.149 0.076 0.073 0.073
## 179 q8 ~~ q19 7.421 -0.034 -0.034 -0.050 -0.050
## 293 q10 ~~ q9 7.240 -0.049 -0.049 -0.058 -0.058
## 107 q1 ~~ q17 7.147 -0.025 -0.025 -0.056 -0.056
## 283 q7 ~~ q19 7.132 -0.045 -0.045 -0.060 -0.060
## 198 q11 ~~ q22 7.124 -0.033 -0.033 -0.052 -0.052
## 292 q10 ~~ q3 7.093 0.036 0.036 0.057 0.057
## 314 q13 ~~ q22 6.918 -0.036 -0.036 -0.058 -0.058
## 346 q9 ~~ q19 6.830 -0.071 -0.071 -0.070 -0.070
## 243 q20 ~~ q9 6.806 -0.054 -0.054 -0.053 -0.053
## 132 q4 ~~ q7 6.327 0.032 0.032 0.056 0.056
## 189 q11 ~~ q12 6.320 -0.024 -0.024 -0.049 -0.049
## 127 q4 ~~ q16 6.179 -0.028 -0.028 -0.058 -0.058
## 235 q20 ~~ q10 6.176 -0.036 -0.036 -0.047 -0.047
## 218 q17 ~~ q21 5.863 -0.025 -0.025 -0.048 -0.048
## 200 q16 ~~ q17 5.815 -0.024 -0.024 -0.052 -0.052
## 306 q12 ~~ q23 5.693 0.035 0.035 0.050 0.050
## 214 q16 ~~ q19 5.689 -0.034 -0.034 -0.054 -0.054
## 204 q16 ~~ q7 5.672 -0.028 -0.028 -0.053 -0.053
## 80 f2 =~ q3 5.651 0.178 0.123 0.114 0.114
## 254 q21 ~~ q18 5.545 0.027 0.027 0.049 0.049
## 230 q17 ~~ q19 5.532 0.031 0.031 0.049 0.049
## 58 f1 =~ q7 5.183 0.224 0.114 0.103 0.103
## 288 q10 ~~ q14 5.082 -0.029 -0.029 -0.048 -0.048
## 176 q8 ~~ q2 4.892 0.022 0.022 0.041 0.041
## 143 q4 ~~ q22 4.885 0.032 0.032 0.048 0.048
## 86 f3 =~ q4 4.837 0.125 0.053 0.056 0.056
## 258 q21 ~~ q19 4.807 -0.033 -0.033 -0.045 -0.045
## 257 q21 ~~ q9 4.416 0.036 0.036 0.044 0.044
## 77 f2 =~ q20 4.217 -0.151 -0.104 -0.101 -0.101
## 290 q10 ~~ q18 4.116 -0.025 -0.025 -0.045 -0.045
## 301 q12 ~~ q2 4.092 -0.023 -0.023 -0.044 -0.044
## 321 q14 ~~ q19 4.020 -0.031 -0.031 -0.044 -0.044
## 319 q14 ~~ q3 3.926 -0.026 -0.026 -0.044 -0.044
## 114 q1 ~~ q13 3.873 0.019 0.019 0.043 0.043lavTestLRT(fit2, fit3)## Chi-Squared Difference Test
##
## Df AIC BIC Chisq Chisq diff Df diff Pr(>Chisq)
## fit3 223 148427 148737 2141.6
## fit2 226 149095 149387 2815.4 673.79 3 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1