******************************************************************************************************************* *Title: JRule-SPSS_for_MPLUS_&_EQS.sps *Description: Implements Power test for MI and EPC and prioritises parameters for modification in * covariance strucutre analysis. Also applies Hancocks type I error adjustment and his sequential * type adjustment. * *Reference: Saris, W, E., Satorra, A., & van der Veld, W.M.(2009). Testing Structural Equation * Models or Detection of Misspecifications? Structural Equation Modeling: A Multidisciplinary * Journal, 16, 561-582. * * Hancock, G. R. (1999). A sequential Scheffe-type respecification procedure for controlling * type I error in exploratory structural equation model modification. Structural Equation Modeling: * A Multidisciplinary Journal, 6, 158-168. * *Software: Currnently set to use Mplus 6.12 MODINDICIES output and EQS output for 6.0. * *Instructions: Cut and paste the Mplus or EQS output as per example below. * *Warning One: Implied under the MI change statistic, althought rarely acknowleged, is that the chi-square test of * perfect fit for the baseline model is not statistically non-significant, that is, that the baseline model * is not misspecified. Therefore the model, theoretically, should only be misspecified by 1 parameter * for r the MI statistic to be reliable (see Yuan & Bentler, 2004). * Use: Free for non-commercial use with acknowledgement.  This spss program is licensed under an Attribution-NonCommercial-ShareAlike 4.0 International. * * Suggested citation: * * Little, J. (2024). JRule-SPSS for MPLUS & EQS (Version 2.1). * https://www.statpsychservices.com.au/resources. Melbourne: Statistical Psychological Services. * * *Writer: Jonathon Little *Date: 5th May 2012 *Revised: 06 October 2024 *Version 2.1 *NOTES: Modify the values where prompted for USER input. ******************************************************************************************************************. *Sample Output from EQS*. DATA LIST LIST /J1(F4) J2(F4) J3(F4) P1(A8) P2(A8) MI(F8.3) Sig(F8.3) Han_sig(F8.3) EPC(F8.3) StdEPC(F8.3) CFI(F8.3) RMSEA(F8.3). BEGIN DATA 1 2 12 V10,F2 367.719 0.000 0.000 0.374 0.374 0.141 0.957 2 2 6 E12,E11 269.071 0.000 0.000 0.169 1.942 0.143 0.956 3 2 6 E5,E3 261.008 0.000 0.000 0.154 0.390 0.143 0.956 4 2 6 E12,E10 226.228 0.000 0.000 -0.104 -0.666 0.144 0.955 5 2 6 E9,E6 226.144 0.000 0.000 0.147 1.192 0.144 0.955 6 2 6 E13,E2 207.803 0.000 0.000 0.118 0.376 0.144 0.955 7 2 12 V10,F4 202.980 0.000 0.000 0.355 0.355 0.144 0.955 8 2 12 V13,F1 202.405 0.000 0.000 0.323 0.323 0.144 0.955 9 2 6 E16,E5 196.470 0.000 0.000 0.167 0.326 0.144 0.955 10 2 6 E9,E2 177.070 0.000 0.000 0.091 0.355 0.145 0.955 11 2 6 E8,E2 176.110 0.000 0.000 0.127 0.567 0.145 0.955 12 2 12 V11,F4 170.457 0.000 0.000 -0.272 -0.272 0.145 0.955 13 2 6 E10,E2 165.253 0.000 0.000 -0.091 -0.337 0.145 0.955 14 2 12 V6,F1 144.830 0.000 0.000 -0.290 -0.290 0.145 0.954 15 2 6 E9,E5 136.905 0.000 0.000 -0.100 -0.318 0.145 0.954 16 2 12 V7,F2 132.735 0.000 0.001 -0.270 -0.270 0.145 0.954 17 2 6 E10,E5 128.872 0.000 0.001 0.090 0.271 0.146 0.954 18 2 6 E13,E9 128.676 0.000 0.001 0.089 0.277 0.146 0.954 19 2 12 V11,F2 125.687 0.000 0.002 -0.169 -0.169 0.146 0.954 20 2 6 E5,E2 123.200 0.000 0.003 -0.089 -0.288 0.146 0.954 21 2 6 E16,E13 118.642 0.000 0.008 -0.132 -0.252 0.146 0.954 22 2 6 E9,E3 113.070 0.000 0.019 -0.086 -0.264 0.146 0.954 23 2 12 V8,F4 108.278 0.000 0.038 -0.212 -0.212 0.146 0.954 24 2 6 E9,E7 108.012 0.000 0.040 -0.067 -0.267 0.146 0.954 25 2 12 V1,F2 105.684 0.000 0.055 0.273 0.273 0.146 0.954 26 2 12 V3,F3 104.751 0.000 0.062 0.227 0.227 0.146 0.954 27 2 6 E13,E5 103.571 0.000 0.073 -0.093 -0.242 0.146 0.954 28 2 6 E16,E8 100.685 0.000 0.104 0.099 0.268 0.146 0.954 29 2 12 V10,F1 100.360 0.000 0.108 0.186 0.186 0.146 0.954 30 2 6 E16,E10 98.864 0.000 0.128 0.104 0.231 0.146 0.954 31 2 12 V11,F1 96.836 0.000 0.160 -0.137 -0.137 0.146 0.954 32 2 6 E11,E1 87.631 0.000 0.372 -0.049 -0.285 0.146 0.954 33 2 6 E13,E3 77.218 0.000 0.687 -0.086 -0.214 0.147 0.954 34 2 6 E6,E1 74.676 0.000 0.757 0.050 0.359 0.147 0.954 35 2 12 V5,F3 74.375 0.000 0.765 0.194 0.194 0.147 0.954 36 2 6 E8,E6 73.832 0.000 0.778 -0.043 -0.402 0.147 0.954 37 2 12 V3,F4 72.618 0.000 0.808 0.201 0.201 0.147 0.954 38 2 12 V16,F1 67.195 0.000 0.910 0.282 0.282 0.147 0.953 39 2 12 V9,F4 64.189 0.000 0.947 -0.202 -0.202 0.147 0.953 40 2 6 E16,E6 59.034 0.000 0.982 -0.067 -0.330 0.147 0.953 41 2 6 E10,E1 55.534 0.000 0.993 0.058 0.187 0.147 0.953 42 2 12 V8,F2 53.317 0.000 0.996 -0.183 -0.183 0.147 0.953 43 2 12 V14,F3 51.724 0.000 0.998 -0.168 -0.168 0.147 0.953 44 2 6 E14,E11 49.753 0.000 0.999 -0.030 -0.227 0.147 0.953 45 2 12 V1,F4 49.713 0.000 0.999 0.156 0.156 0.147 0.953 46 2 12 V5,F4 49.696 0.000 0.999 0.195 0.195 0.147 0.953 47 2 6 E14,E8 48.831 0.000 0.999 -0.041 -0.212 0.147 0.953 48 2 12 V13,F3 48.607 0.000 0.999 0.200 0.200 0.147 0.953 49 2 6 E8,E7 45.127 0.000 1.000 -0.043 -0.195 0.147 0.953 50 2 12 V16,F4 45.045 0.000 1.000 -0.234 -0.234 0.147 0.953 51 2 12 V9,F3 43.737 0.000 1.000 -0.132 -0.132 0.147 0.953 52 2 12 V5,F1 42.697 0.000 1.000 0.177 0.177 0.147 0.953 53 2 6 E3,E2 42.515 0.000 1.000 -0.062 -0.194 0.147 0.953 54 2 6 E15,E10 40.125 0.000 1.000 0.036 0.177 0.147 0.953 55 2 6 E10,E8 39.358 0.000 1.000 0.042 0.174 0.147 0.953 56 2 6 E13,E8 38.800 0.000 1.000 0.048 0.172 0.147 0.953 57 2 12 V7,F1 35.929 0.000 1.000 -0.116 -0.116 0.147 0.953 58 2 6 E14,E5 34.657 0.000 1.000 0.041 0.153 0.147 0.953 59 2 6 E8,E1 33.274 0.000 1.000 -0.054 -0.211 0.147 0.953 60 2 6 E11,E5 33.087 0.000 1.000 -0.031 -0.166 0.147 0.953 61 2 12 V2,F3 32.734 0.000 1.000 -0.112 -0.112 0.147 0.953 62 2 6 E16,E7 31.782 0.000 1.000 -0.056 -0.137 0.147 0.953 63 2 6 E14,E1 28.827 0.000 1.000 0.036 0.147 0.148 0.953 64 2 6 E12,E8 26.257 0.000 1.000 -0.023 -0.179 0.148 0.953 65 2 6 E14,E13 25.786 0.000 1.000 -0.039 -0.144 0.148 0.953 66 2 6 E11,E7 25.457 0.000 1.000 0.023 0.155 0.148 0.953 67 2 12 V2,F2 25.350 0.000 1.000 -0.128 -0.128 0.148 0.953 68 2 6 E16,E15 24.040 0.000 1.000 0.042 0.131 0.148 0.953 69 2 12 V8,F3 23.732 0.000 1.000 -0.092 -0.092 0.148 0.953 70 2 6 E10,E7 23.028 0.000 1.000 -0.032 -0.121 0.148 0.953 71 2 6 E15,E2 22.823 0.000 1.000 -0.028 -0.145 0.148 0.953 72 2 6 E16,E2 22.682 0.000 1.000 -0.050 -0.120 0.148 0.953 73 2 6 E13,E12 22.283 0.000 1.000 0.026 0.143 0.148 0.953 74 2 12 V6,F4 21.816 0.000 1.000 0.119 0.119 0.148 0.953 75 2 12 V9,F1 21.194 0.000 1.000 0.112 0.112 0.148 0.953 76 2 6 E6,E2 21.125 0.000 1.000 -0.024 -0.201 0.148 0.953 77 2 6 E9,E1 20.754 0.000 1.000 -0.034 -0.116 0.148 0.953 78 2 12 V14,F2 20.522 0.000 1.000 0.100 0.100 0.148 0.953 79 2 6 E7,E3 20.473 0.000 1.000 0.037 0.115 0.148 0.953 80 2 6 E12,E6 19.641 0.000 1.000 -0.015 -0.222 0.148 0.953 81 2 6 E13,E10 18.915 0.000 1.000 -0.035 -0.105 0.148 0.953 82 2 6 E14,E7 18.901 0.000 1.000 0.030 0.142 0.148 0.953 83 2 6 E14,E10 17.995 0.000 1.000 0.026 0.112 0.148 0.953 84 2 12 V3,F2 17.697 0.000 1.000 0.119 0.119 0.148 0.953 85 2 6 E12,E2 17.392 0.000 1.000 0.020 0.137 0.148 0.953 86 2 6 E9,E8 16.855 0.000 1.000 0.026 0.116 0.148 0.953 87 2 6 E12,E3 15.397 0.000 1.000 0.022 0.120 0.148 0.953 88 2 6 E15,E6 15.126 0.000 1.000 0.017 0.182 0.148 0.953 89 2 6 E15,E3 14.732 0.000 1.000 0.026 0.108 0.148 0.953 90 2 6 E15,E1 14.013 0.000 1.000 -0.024 -0.108 0.148 0.953 91 2 12 V15,F2 13.464 0.000 1.000 0.079 0.079 0.148 0.953 92 2 6 E6,E5 13.240 0.000 1.000 -0.031 -0.211 0.148 0.953 93 2 6 E11,E6 12.494 0.000 1.000 0.012 0.171 0.148 0.953 94 2 6 E8,E5 12.247 0.000 1.000 0.026 0.096 0.148 0.953 95 2 6 E15,E13 11.619 0.001 1.000 -0.025 -0.106 0.148 0.953 96 2 6 E12,E1 11.604 0.001 1.000 0.018 0.107 0.148 0.953 97 2 6 E13,E6 11.083 0.001 1.000 -0.020 -0.131 0.148 0.953 98 2 12 V7,F3 10.964 0.001 1.000 0.082 0.082 0.148 0.953 99 2 6 E10,E3 10.877 0.001 1.000 -0.028 -0.081 0.148 0.953 100 2 12 V1,F3 10.600 0.001 1.000 0.067 0.067 0.148 0.953 101 2 12 V12,F1 10.173 0.001 1.000 0.044 0.044 0.148 0.953 102 2 6 E16,E11 9.648 0.002 1.000 -0.022 -0.087 0.148 0.953 103 2 12 V13,F2 9.046 0.003 1.000 0.082 0.082 0.148 0.953 104 2 6 E6,E3 8.778 0.003 1.000 -0.018 -0.119 0.148 0.953 105 2 6 E15,E11 8.652 0.003 1.000 -0.012 -0.101 0.148 0.953 106 2 6 E11,E8 8.304 0.004 1.000 0.013 0.097 0.148 0.953 107 2 6 E11,E9 7.958 0.005 1.000 -0.013 -0.084 0.148 0.953 108 2 6 E15,E9 7.946 0.005 1.000 -0.016 -0.080 0.148 0.953 109 2 6 E2,E1 7.242 0.007 1.000 -0.025 -0.087 0.148 0.953 110 2 6 E7,E1 7.025 0.008 1.000 0.020 0.070 0.148 0.953 111 2 12 V14,F1 7.003 0.008 1.000 -0.048 -0.048 0.148 0.953 112 2 6 E11,E10 6.839 0.009 1.000 0.018 0.110 0.148 0.953 113 2 12 V12,F4 6.390 0.011 1.000 0.053 0.053 0.148 0.953 114 2 6 E12,E5 5.226 0.022 1.000 0.012 0.068 0.148 0.953 115 2 6 E7,E5 4.917 0.027 1.000 0.017 0.055 0.148 0.953 116 2 6 E11,E3 4.799 0.028 1.000 0.013 0.065 0.148 0.953 117 2 6 E14,E6 4.732 0.030 1.000 0.010 0.095 0.148 0.953 118 2 6 E15,E14 4.722 0.030 1.000 0.015 0.093 0.148 0.953 119 2 6 E13,E7 4.023 0.045 1.000 0.016 0.053 0.148 0.953 120 2 12 V2,F4 3.878 0.049 1.000 -0.041 -0.041 0.148 0.953 121 2 12 V12,F2 3.418 0.064 1.000 -0.028 -0.028 0.148 0.953 122 2 6 E16,E12 2.999 0.083 1.000 -0.012 -0.050 0.148 0.953 123 2 12 V15,F1 2.631 0.105 1.000 -0.028 -0.028 0.148 0.953 124 2 6 E16,E3 2.376 0.123 1.000 0.019 0.036 0.148 0.953 125 2 6 E15,E8 2.373 0.123 1.000 0.008 0.050 0.148 0.953 126 2 6 E13,E1 2.148 0.143 1.000 -0.013 -0.037 0.148 0.953 127 2 6 E16,E1 1.987 0.159 1.000 0.016 0.034 0.148 0.953 128 2 6 E16,E9 1.968 0.161 1.000 -0.014 -0.034 0.148 0.953 129 2 6 E12,E7 1.781 0.182 1.000 0.006 0.042 0.148 0.953 130 2 6 E15,E12 1.728 0.189 1.000 -0.005 -0.047 0.148 0.953 131 2 12 V16,F3 1.382 0.240 1.000 -0.034 -0.034 0.148 0.953 132 2 6 E14,E2 1.371 0.242 1.000 -0.007 -0.033 0.148 0.953 133 2 6 E5,E1 1.229 0.268 1.000 -0.010 -0.028 0.148 0.953 134 2 6 E16,E14 1.099 0.294 1.000 -0.010 -0.027 0.148 0.953 135 2 6 E15,E7 0.963 0.326 1.000 -0.007 -0.036 0.148 0.953 136 2 6 E10,E6 0.761 0.383 1.000 0.004 0.034 0.148 0.953 137 2 6 E10,E9 0.719 0.396 1.000 0.006 0.021 0.148 0.953 138 2 6 E7,E2 0.630 0.427 1.000 0.005 0.022 0.148 0.953 139 2 12 V15,F3 0.514 0.473 1.000 -0.016 -0.016 0.148 0.953 140 2 6 E15,E5 0.262 0.609 1.000 0.003 0.014 0.148 0.953 141 2 6 E3,E1 0.182 0.670 1.000 0.004 0.012 0.148 0.953 142 2 6 E13,E11 0.118 0.731 1.000 0.002 0.010 0.148 0.953 143 2 6 E14,E9 0.108 0.742 1.000 0.002 0.009 0.148 0.953 144 2 6 E14,E12 0.038 0.846 1.000 0.001 0.006 0.148 0.953 145 2 6 E7,E6 0.035 0.852 1.000 0.001 0.008 0.148 0.953 146 2 6 E8,E3 0.017 0.896 1.000 0.001 0.004 0.148 0.953 147 2 6 E12,E9 0.012 0.914 1.000 0.000 -0.003 0.148 0.953 148 2 12 V6,F3 0.005 0.945 1.000 0.001 0.001 0.148 0.953 149 2 6 E14,E3 0.001 0.970 1.000 0.000 0.001 0.148 0.953 150 2 6 E11,E2 0.001 0.978 1.000 0.000 -0.001 0.148 0.953 151 2 0 F4,F4 0.000 1.000 1.000 0.000 0.000 0.148 0.953 152 2 0 F3,F3 0.000 1.000 1.000 0.000 0.000 0.148 0.953 153 2 0 F2,F2 0.000 1.000 1.000 0.000 0.000 0.148 0.953 154 2 0 F1,F1 0.000 1.000 1.000 0.000 0.000 0.148 0.953 END DATA. DATASET CLOSE ALL. *Sample Output from Mplus*. DATA LIST LIST /P1(A8) P2(A4) P3(A8) MI(F8.3) EPC(F8.3) StdEPC(F8.3) StdYXEPC(F8.3). BEGIN DATA WITHDRAW BY ANSWER5 319.973 0.152 0.152 0.152 WITHDRAW BY ANSWER6 216.292 -0.132 -0.132 -0.132 WITHDRAW BY ANSWER10 98.543 0.054 0.054 0.054 WITHDRAW BY ANSWER11 435.687 -0.129 -0.129 -0.129 WITHDRAW BY ANSWER12 273.390 0.087 0.087 0.087 WITHDRAW BY ANSWER13 319.894 0.139 0.139 0.139 WITHDRAW BY ANSWER14 202.060 -0.095 -0.095 -0.095 WITHDRAW BY ANSWER15 50.762 -0.047 -0.047 -0.047 SELFCARE BY ANSWER1 22.504 -0.038 -0.038 -0.038 SELFCARE BY ANSWER2 413.911 -0.157 -0.157 -0.157 SELFCARE BY ANSWER3 26.247 0.040 0.040 0.040 SELFCARE BY ANSWER7 357.584 -0.152 -0.152 -0.152 SELFCARE BY ANSWER10 437.317 0.124 0.124 0.124 SELFCARE BY ANSWER11 663.755 -0.181 -0.181 -0.181 SELFCARE BY ANSWER12 115.780 0.062 0.062 0.062 SELFCARE BY ANSWER14 145.065 -0.095 -0.095 -0.095 SELFCARE BY ANSWER15 207.783 0.112 0.112 0.112 SELFCARE BY ANSWER16 1037.888 0.303 0.303 0.303 COMPLIAN BY ANSWER1 92.790 -0.063 -0.063 -0.063 COMPLIAN BY ANSWER2 928.467 -0.194 -0.194 -0.194 COMPLIAN BY ANSWER3 885.139 0.194 0.194 0.194 COMPLIAN BY ANSWER5 133.449 0.086 0.086 0.086 COMPLIAN BY ANSWER6 46.815 -0.053 -0.053 -0.053 COMPLIAN BY ANSWER8 97.293 0.063 0.063 0.063 COMPLIAN BY ANSWER9 15.511 -0.029 -0.029 -0.029 COMPLIAN BY ANSWER13 31.156 0.062 0.062 0.062 COMPLIAN BY ANSWER14 141.069 -0.117 -0.117 -0.117 COMPLIAN BY ANSWER15 37.964 0.060 0.060 0.060 ANTISOCI BY ANSWER1 95.434 -0.067 -0.067 -0.067 ANTISOCI BY ANSWER2 1674.506 -0.274 -0.274 -0.274 ANTISOCI BY ANSWER3 821.424 0.196 0.196 0.196 ANTISOCI BY ANSWER5 210.446 0.127 0.127 0.127 ANTISOCI BY ANSWER6 32.275 -0.053 -0.053 -0.053 ANTISOCI BY ANSWER8 302.310 0.116 0.116 0.116 ANTISOCI BY ANSWER9 72.537 -0.075 -0.075 -0.075 ANTISOCI BY ANSWER10 36.573 -0.062 -0.062 -0.062 ANTISOCI BY ANSWER12 36.561 0.063 0.063 0.063 ANTISOCI BY ANSWER16 109.602 0.086 0.086 0.086 ANSWER2 WITH ANSWER1 197.920 0.062 0.062 0.155 ANSWER3 WITH ANSWER1 194.317 0.061 0.061 0.146 ANSWER3 WITH ANSWER2 128.830 -0.050 -0.050 -0.126 ANSWER5 WITH ANSWER1 129.731 0.044 0.044 0.102 ANSWER5 WITH ANSWER2 219.751 -0.056 -0.056 -0.136 ANSWER5 WITH ANSWER3 273.461 0.064 0.064 0.149 ANSWER6 WITH ANSWER1 12.554 -0.012 -0.012 -0.040 ANSWER6 WITH ANSWER3 160.979 -0.044 -0.044 -0.142 ANSWER6 WITH ANSWER5 48.060 -0.048 -0.048 -0.150 ANSWER7 WITH ANSWER2 121.183 -0.040 -0.040 -0.101 ANSWER7 WITH ANSWER3 350.521 0.070 0.070 0.168 ANSWER7 WITH ANSWER6 114.301 -0.037 -0.037 -0.120 ANSWER8 WITH ANSWER1 468.719 -0.095 -0.095 -0.242 ANSWER8 WITH ANSWER2 390.230 0.088 0.088 0.237 ANSWER8 WITH ANSWER3 23.458 -0.021 -0.021 -0.055 ANSWER8 WITH ANSWER5 26.646 -0.019 -0.019 -0.048 ANSWER8 WITH ANSWER6 12.434 -0.012 -0.012 -0.041 ANSWER8 WITH ANSWER7 39.695 0.023 0.023 0.058 ANSWER9 WITH ANSWER1 46.585 -0.027 -0.027 -0.060 ANSWER9 WITH ANSWER2 142.943 0.046 0.046 0.108 ANSWER9 WITH ANSWER3 152.647 -0.049 -0.049 -0.109 ANSWER9 WITH ANSWER5 236.096 -0.082 -0.082 -0.176 ANSWER9 WITH ANSWER6 512.560 0.148 0.148 0.446 ANSWER9 WITH ANSWER7 98.046 -0.039 -0.039 -0.087 ANSWER9 WITH ANSWER8 10.102 0.012 0.012 0.029 ANSWER10 WITH ANSWER1 119.337 0.034 0.034 0.095 ANSWER10 WITH ANSWER3 133.304 -0.035 -0.035 -0.100 ANSWER10 WITH ANSWER5 25.631 0.016 0.016 0.044 ANSWER10 WITH ANSWER6 111.983 0.029 0.029 0.113 ANSWER10 WITH ANSWER7 221.596 -0.046 -0.046 -0.133 ANSWER10 WITH ANSWER8 13.275 -0.011 -0.011 -0.033 ANSWER10 WITH ANSWER9 38.939 0.020 0.020 0.053 ANSWER11 WITH ANSWER1 36.835 -0.017 -0.017 -0.093 ANSWER11 WITH ANSWER3 11.706 0.010 0.010 0.052 ANSWER11 WITH ANSWER6 92.034 -0.026 -0.026 -0.188 ANSWER11 WITH ANSWER7 141.016 0.034 0.034 0.187 ANSWER11 WITH ANSWER10 114.720 0.063 0.063 0.404 ANSWER12 WITH ANSWER1 28.654 -0.016 -0.016 -0.047 ANSWER12 WITH ANSWER3 418.526 0.061 0.061 0.179 ANSWER12 WITH ANSWER10 559.106 -0.240 -0.240 -0.843 ANSWER13 WITH ANSWER1 20.034 -0.020 -0.020 -0.038 ANSWER13 WITH ANSWER2 83.041 0.038 0.038 0.079 ANSWER13 WITH ANSWER5 74.491 -0.039 -0.039 -0.073 ANSWER13 WITH ANSWER7 168.329 0.059 0.059 0.116 ANSWER13 WITH ANSWER8 274.701 0.069 0.069 0.144 ANSWER13 WITH ANSWER10 14.413 -0.014 -0.014 -0.032 ANSWER13 WITH ANSWER12 20.131 0.016 0.016 0.038 ANSWER14 WITH ANSWER2 410.547 -0.070 -0.070 -0.191 ANSWER14 WITH ANSWER3 10.945 0.012 0.012 0.031 ANSWER14 WITH ANSWER5 199.169 0.052 0.052 0.130 ANSWER14 WITH ANSWER6 113.095 -0.035 -0.035 -0.123 ANSWER14 WITH ANSWER7 618.041 0.106 0.106 0.277 ANSWER14 WITH ANSWER8 20.373 0.015 0.015 0.043 ANSWER14 WITH ANSWER9 46.180 -0.026 -0.026 -0.062 ANSWER14 WITH ANSWER10 150.208 -0.037 -0.037 -0.114 ANSWER14 WITH ANSWER11 24.797 0.014 0.014 0.082 ANSWER14 WITH ANSWER12 18.645 -0.013 -0.013 -0.041 ANSWER14 WITH ANSWER13 55.932 -0.033 -0.033 -0.070 ANSWER15 WITH ANSWER1 19.054 -0.015 -0.015 -0.042 ANSWER15 WITH ANSWER2 44.206 -0.022 -0.022 -0.065 ANSWER15 WITH ANSWER3 10.246 -0.011 -0.011 -0.031 ANSWER15 WITH ANSWER6 226.938 0.048 0.048 0.181 ANSWER15 WITH ANSWER7 289.731 -0.073 -0.073 -0.206 ANSWER15 WITH ANSWER9 14.077 0.014 0.014 0.035 ANSWER15 WITH ANSWER10 177.571 0.039 0.039 0.130 ANSWER15 WITH ANSWER11 36.773 -0.017 -0.017 -0.105 ANSWER15 WITH ANSWER13 83.700 -0.040 -0.040 -0.091 ANSWER15 WITH ANSWER14 20.343 0.020 0.020 0.060 ANSWER16 WITH ANSWER3 233.546 -0.074 -0.074 -0.133 ANSWER16 WITH ANSWER5 217.347 0.071 0.071 0.123 ANSWER16 WITH ANSWER6 162.935 0.054 0.054 0.132 ANSWER16 WITH ANSWER7 79.885 -0.041 -0.041 -0.074 ANSWER16 WITH ANSWER9 22.547 0.023 0.023 0.039 ANSWER16 WITH ANSWER10 194.194 0.053 0.053 0.113 ANSWER16 WITH ANSWER11 119.636 -0.038 -0.038 -0.155 ANSWER16 WITH ANSWER12 16.460 -0.015 -0.015 -0.033 ANSWER16 WITH ANSWER13 21.304 -0.025 -0.025 -0.036 ANSWER16 WITH ANSWER14 32.577 0.025 0.025 0.049 END DATA. *************************USER SECTION*START****************************. */minimum parameter size for factor loading - default .40. COMPUTE delta_a = .40. */minimum parameter size for causal effects - default .10. COMPUTE delta_b = .10. */minimum parameter size for correlated errors - default .20. COMPUTE delta_c = .20. */nominal alpha (two-tailed). COMPUTE alpha = .05. */beta. COMPUTE qbeta = .20. */degrees of freedom for your model. COMPUTE df = 98. */chosen power level. COMPUTE power = .80. */Enter 1 = No adjustment to type I error rate - default setting. */Enter 2 = Hancock's determine significance using df for entire model rather than 1 df. */Enter 3 = Sequential Holm-Bonferroni type adjustment.. COMPUTE Type_1_error_adjustment = 1. EXECUTE. *************************USER SECTION*END*******************************. COMPUTE ABS_EPC_ABS = ABS(EPC). COMPUTE ABS_StdYXEPC = ABS(StdYXEPC). FORMATS ABS_EPC_ABS ABS_StdYXEPC (f8.4). COMPUTE alpha_nominal = IDF.CHISQ(1-alpha, 1). EXECUTE. COMPUTE Hancock_sig1tail=sig.chisq(MI,df). COMPUTE sig1tail=sig.chisq(MI,1). EXECUTE. DO IF (( MI <> 0) AND (ABS_StdYXEPC <> 0)). *non-centrality parameter. COMPUTE ncp_a =(MI/(ABS_StdYXEPC**2))*(delta_a**2). COMPUTE ncp_b =(MI/(ABS_StdYXEPC**2))*(delta_b**2). COMPUTE ncp_c =(MI/(ABS_StdYXEPC**2))*(delta_c**2). FORMATS Hancock_sig1tail sig1tail ncp_a ncp_b ncp_c (f8.4). END IF. EXECUTE. DO IF (( MI <> 0) AND (ABS_StdYXEPC <> 0)). COMPUTE power_a = 1 - ncdf.chi(alpha_nominal ,1,ncp_a). COMPUTE power_b = 1 - ncdf.chi(alpha_nominal ,1,ncp_b). COMPUTE power_c = 1 - ncdf.chi(alpha_nominal ,1,ncp_c). END IF. EXECUTE. DO IF (( MI = 0) OR (ABS_StdYXEPC = 0)). COMPUTE power_a = .999. COMPUTE power_b = .999. COMPUTE power_c = .999. END IF. EXECUTE. DO IF (( MI <> 0) AND (ABS_StdYXEPC <> 0)). COMPUTE ABS_StdYXEPC_stdr =ABS_StdYXEPC/(SQRT(MI)). COMPUTE ABS_StdYXEPC_95_lower = ABS_StdYXEPC - (1.96 * ABS_StdYXEPC). COMPUTE ABS_StdYXEPC_95_upper = ABS_StdYXEPC + (1.96 * ABS_StdYXEPC). END IF. EXECUTE. NUMERIC decision_A decision_B (F1.0). VALUE LABELS decision_A 1 'EPC:Misspecified' 2 'EPC: Not Misspecified' 3 'No misspecification' 4 'Misspecified' 5 'Inconclusive'. VALUE LABELS decision_B 1 'EPC CI Overlaps: Inconclusive' 2 'EPC CI Exceeds: Severely Misspecified' 3 'EPC CI Within: Trivally Misspecified'. EXECUTE. DO IF Type_1_error_adjustment = 2. COMPUTE sig1tail = Hancock_sig1tail. END IF. EXECUTE. COMPUTE casetemp = $casenum. EXECUTE. *You could sort these sig1tail values within P2 types if you prefer. SORT CASES BY sig1tail(A). EXECUTE. COMPUTE case = $casenum. EXECUTE. AGGREGATE /OUTFILE=* MODE=ADDVARIABLES /BREAK= /maxi=MAX(case). EXECUTE. DO IF Type_1_error_adjustment = 3. COMPUTE prob = alpha/(maxi-case+1). IF sig1tail < prob sig1tail = 0. IF sig1tail >= prob sig1tail = 1. END IF. EXECUTE. SORT CASES BY casetemp(A). EXECUTE. DELETE VARIABLES casetemp. EXECUTE. *To decide whether a parameter should be freed, inspect its modification index (MI) and expected parameter change (EPC). *These values can be used to evaluate model fit by 2 methods - DECISION A and DECISION B methods. *************************DECISION_A MATRIX************START******************. *Method 1: Saris, Satorra, and van der Veld (2009, pp. 570–573) used power and EPC to decide whether to free a parameter. *This is called the 'JRule' first implemented with LISREL (van der Veld et al., 2008): *If the MI is not significant and the power is low, the test is inconclusive. IF (sig1tail >= alpha) AND (power_a < power) decision_A = 5. *If the MI is not significant and the power is high, there is no misspecification. IF (sig1tail >= alpha) AND (power_a >= power) decision_A = 3. *If the MI is significant and the power is low, the fixed parameter is misspecified. IF (sig1tail < alpha) AND (power_a < power) decision_A = 4. EXECUTE. DO IF P2 = 'BY'. *If the MI is significant and the power is high, the EPC is investigated. If the EPC is large *(greater than the the target EPC), the parameter is misspecified. IF ((sig1tail < alpha) AND (power_a >= power)) AND (ABS_StdYXEPC >= delta_a) decision_A = 1. *If the MI is significant and the power is high, the EPC is investigated. *If the EPC is low (lower than the target EPC), the parameter is not misspecificied. IF ((sig1tail < alpha) AND (power_a >= power)) AND (ABS_StdYXEPC < delta_a) decision_A = 2. END IF. EXECUTE. DO IF P2 = 'ON'. *If the MI is significant and the power is high, the EPC is investigated. If the EPC is large *(greater than the the target EPC), the parameter is misspecified. IF ((sig1tail < alpha) AND (power_a >= power)) AND (ABS_StdYXEPC >= delta_b) decision_A = 1. *If the MI is significant and the power is high, the EPC is investigated. *If the EPC is low (lower than the target EPC), the parameter is not misspecificied. IF ((sig1tail < alpha) AND (power_a >= power)) AND (ABS_StdYXEPC < delta_b) decision_A = 2. END IF. EXECUTE. DO IF P2 = 'WITH'. *If the MI is significant and the power is high, the EPC is investigated. If the EPC is large *(greater than the the target EPC), the parameter is misspecified. IF ((sig1tail < alpha) AND (power_a >= power)) AND (ABS_StdYXEPC >= delta_c) decision_A = 1. *If the MI is significant and the power is high, the EPC is investigated. *If the EPC is low (lower than the target EPC), the parameter is not misspecificied. IF ((sig1tail < alpha) AND (power_a >= power)) AND (ABS_StdYXEPC < delta_c) decision_A = 2. END IF. EXECUTE. *************************DECISION_A MATRIX************END******************. *************************DECISION_B MATRIX************START******************. DO IF P2 = 'BY'. *If a CI overlaps with the range of trivial misspecification, the test is inconclusive. IF (ABS_StdYXEPC_95_lower <= delta_a) AND (ABS_StdYXEPC_95_upper >= delta_a) decision_B = 1. *If a CI completely exceeds the range of trivial misspecification, the fixed parameters are severely misspecified. IF (ABS_StdYXEPC_95_lower > delta_a) decision_B = 2. *IF upper CI is less than delta then the value is trivially misspecified. IF (ABS_StdYXEPC_95_upper < delta_a) decision_B = 3. *If a CI is completely within the range of trivial misspecification, the fixed parameters are trivially misspecified.. IF RANGE(delta_a,ABS_StdYXEPC_95_lower,ABS_StdYXEPC_95_upper) = 1 decision_B = 3. END IF. EXECUTE. DO IF P2 = 'ON'. *If a CI overlaps with the range of trivial misspecification, the test is inconclusive. IF (ABS_StdYXEPC_95_lower <= delta_b) AND (ABS_StdYXEPC_95_upper >= delta_b) decision_B = 1. *If a CI completely exceeds the range of trivial misspecification, the fixed parameters are severely misspecified. IF (ABS_StdYXEPC_95_lower > delta_b) decision_B = 2. *IF upper CI is less than delta then the value is trivially misspecified. IF (ABS_StdYXEPC_95_upper < delta_b) decision_B = 3. *If a CI is completely within the range of trivial misspecification, the fixed parameters are trivially misspecified.. IF RANGE(delta_b,ABS_StdYXEPC_95_lower,ABS_StdYXEPC_95_upper) = 1 decision_B = 3. END IF. EXECUTE. DO IF P2 = 'WITH'. *If a CI overlaps with the range of trivial misspecification, the test is inconclusive. IF (ABS_StdYXEPC_95_lower <= delta_c) AND (ABS_StdYXEPC_95_upper >= delta_c) decision_B = 1. *If a CI completely exceeds the range of trivial misspecification, the fixed parameters are severely misspecified. IF (ABS_StdYXEPC_95_lower > delta_c) decision_B = 2. *IF upper CI is less than delta then the value is trivially misspecified. IF (ABS_StdYXEPC_95_upper < delta_c) decision_B = 3. *If a CI is completely within the range of trivial misspecification, the fixed parameters are trivially misspecified.. IF RANGE(delta_c,ABS_StdYXEPC_95_lower,ABS_StdYXEPC_95_upper) = 1 decision_B = 3. END IF. EXECUTE. *************************DECISION_B MATRIX************END******************. DELETE VARIABLES alpha_nominal Hancock_sig1tail sig1tail ABS_StdYXEPC_stdr case maxi prob ABS_EPC_ABS. EXECUTE. CROSSTABS /TABLES=P2 BY decision_A decision_B /FORMAT=AVALUE TABLES /CELLS=COUNT /COUNT ROUND CELL. FREQUENCIES VARIABLES=decision_A decision_B /ORDER=ANALYSIS.