Ofek - vigilantism

Data preparation

  • Item Q6_4 was reversed.

  • Item Q12_3 was reversed.

  • Following the recommendations of … we removed subjects who hurried through the questionnaire, completing it in less than a third of the average completion time of the entire sample (that is, less than 53 seconds).

  • We remain with \(N = 468\) observations.

Survey data

CFA

  • First, we assessed the reliability of the measurement model for the latent variables using confirmatory factor analysis.

  • The table below details the reliability measures for the latent variables.

  • pj, trust, and efficacy show good \(\alpha\) values.

  • vig is acceptable, just over the common threshold of 0.75/

measure pj trust eff vig
alpha 0.9445 0.9277 0.8562 0.7734
omega 0.9446 0.9276 0.8559 0.7898
omega2 0.9446 0.9276 0.8559 0.7898
omega3 0.9446 0.9275 0.8555 0.7908
avevar 0.7089 0.8103 0.6645 0.5656

Measurement invariance

  • Measurement invariance of past victims and non-victims was assessed by fitting a multigroup SEM model.

  • There were 303 non-victims and 121 victims.

  • The model fitted the data well, CFI = .967, TLI = .956, RMSEA = .059.

  • All item loadings were significant in both groups.

  • We conclude that measurement invariance holds; the grouping variable will be treated as a manifest variable (directly observed, without error).

  • Surprisingly in both groups, higher police efficacy was associated with MORE self-vigilantism.

  • The output below details all the estimates for this model.

0 1
Estimate Std. Err. z p Estimate Std. Err. z p
Factor Loadings
pj
q5.1 1.00+ 1.00+
q5.2 0.94 0.04 25.48 .000 1.01 0.09 11.34 .000
q5.3 0.98 0.04 24.23 .000 1.03 0.11 9.78 .000
q5.4 0.98 0.05 21.37 .000 1.06 0.10 10.98 .000
q5.5 0.91 0.04 21.65 .000 1.04 0.10 10.02 .000
q5.6 1.00 0.03 29.85 .000 1.03 0.09 11.17 .000
q5.7 0.95 0.04 25.47 .000 1.12 0.10 11.49 .000
dj
q8.1 1.00+ 1.00+
trust
q7.1 1.00+ 1.00+
q7.2 1.00 0.03 30.45 .000 1.01 0.06 16.26 .000
q7.3 0.98 0.04 22.63 .000 1.02 0.05 18.94 .000
eff
q6.1 1.00+ 1.00+
q6.2 1.06 0.06 18.15 .000 1.01 0.09 10.84 .000
q6.3 1.07 0.06 17.77 .000 0.89 0.10 9.15 .000
vig
q9.1 1.00+ 1.00+
q9.2 1.32 0.11 11.72 .000 1.52 0.49 3.12 .002
q9.3 0.85 0.10 8.53 .000 0.79 0.14 5.54 .000
Regression Slopes
vig
pj -0.32 0.14 -2.33 .020 -0.32 0.14 -2.33 .020
dj 0.01 0.08 0.17 .868 0.01 0.08 0.17 .868
trust -0.22 0.12 -1.88 .060 -0.22 0.12 -1.88 .060
eff 0.52 0.18 2.89 .004 0.52 0.18 2.89 .004
Intercepts
q5.1 2.24 0.06 36.73 .000 2.02 0.10 19.72 .000
q5.2 2.37 0.06 41.84 .000 1.94 0.09 21.77 .000
q5.3 2.48 0.06 43.27 .000 2.13 0.09 23.53 .000
q5.4 2.84 0.06 47.38 .000 2.58 0.10 26.65 .000
q5.5 2.43 0.06 42.27 .000 2.07 0.10 20.87 .000
q5.6 2.34 0.06 40.19 .000 2.06 0.10 21.41 .000
q5.7 2.40 0.06 43.32 .000 2.13 0.09 22.65 .000
q8.1 2.54 0.06 46.00 .000 2.19 0.10 22.28 .000
q7.1 2.72 0.06 44.69 .000 2.28 0.09 23.98 .000
q7.2 2.75 0.06 44.60 .000 2.40 0.10 23.49 .000
q7.3 2.61 0.06 42.22 .000 2.19 0.10 23.05 .000
q6.1 2.45 0.06 44.38 .000 2.06 0.10 21.41 .000
q6.2 2.40 0.06 41.18 .000 2.23 0.10 22.47 .000
q6.3 2.55 0.06 44.04 .000 2.25 0.09 23.84 .000
q9.1 2.24 0.07 33.09 .000 2.25 0.11 20.83 .000
q9.2 1.98 0.06 31.16 .000 2.02 0.10 20.10 .000
q9.3 1.62 0.06 29.48 .000 1.67 0.09 18.53 .000
Residual Variances
q5.1 0.38 0.04 10.39 .000 0.59 0.12 5.08 .000
q5.2 0.31 0.04 7.34 .000 0.27 0.06 4.80 .000
q5.3 0.28 0.03 8.41 .000 0.27 0.04 6.23 .000
q5.4 0.37 0.04 9.61 .000 0.38 0.06 6.77 .000
q5.5 0.39 0.04 9.65 .000 0.48 0.09 5.11 .000
q5.6 0.27 0.03 8.83 .000 0.41 0.07 6.26 .000
q5.7 0.24 0.03 8.89 .000 0.21 0.04 5.22 .000
q8.1 0.00+ 0.00+
q7.1 0.20 0.04 4.64 .000 0.20 0.04 4.80 .000
q7.2 0.23 0.04 6.30 .000 0.37 0.12 3.03 .002
q7.3 0.29 0.05 5.49 .000 0.17 0.04 3.85 .000
q6.1 0.35 0.04 9.91 .000 0.35 0.10 3.57 .000
q6.2 0.39 0.05 8.44 .000 0.39 0.12 3.16 .002
q6.3 0.37 0.05 7.39 .000 0.46 0.09 5.35 .000
q9.1 0.89 0.11 8.21 .000 0.93 0.21 4.42 .000
q9.2 0.21 0.09 2.43 .015 0.11 0.29 0.37 .711
q9.3 0.54 0.07 7.52 .000 0.72 0.14 5.08 .000
Latent Intercepts
pj 0.00+ 0.00+
dj 0.00+ 0.00+
trust 0.00+ 0.00+
eff 0.00+ 0.00+
vig 0.00+ 0.00+
Latent Variances
pj 0.87 0.07 11.74 .000 0.71 0.14 5.03 .000
dj 1.02 0.07 15.02 .000 1.22 0.13 9.53 .000
trust 1.04 0.08 12.65 .000 0.94 0.11 8.61 .000
eff 0.66 0.07 9.41 .000 0.84 0.13 6.57 .000
vig 0.61 0.10 6.14 .000 0.45 0.24 1.87 .061
Latent Covariances
pj w/dj 0.79 0.06 12.14 .000 0.67 0.11 6.31 .000
pj w/trust 0.80 0.06 12.34 .000 0.68 0.10 7.10 .000
pj w/eff 0.68 0.06 10.92 .000 0.68 0.10 6.52 .000
dj w/trust 0.87 0.07 12.89 .000 0.87 0.11 8.16 .000
dj w/eff 0.69 0.06 10.97 .000 0.79 0.12 6.81 .000
trust w/eff 0.67 0.06 10.44 .000 0.77 0.10 7.33 .000
Fit Indices
χ2 454.45
CFI 0.96
TLI 0.96
RMSEA 0.07
Scaled χ2 376.35(224) .000
+Fixed parameter


SEM fitting

  • The previous encounter with the police was treated as a manifest variable.
    • Including it as a latent variable results in a non-positive-definite covariance matrix.
lavaan 0.6.15 ended normally after 69 iterations

  Estimator                                         ML
  Optimization method                           NLMINB
  Number of model parameters                        49

                                                  Used       Total
  Number of observations                           461         538

Model Test User Model:
                                              Standard      Scaled
  Test Statistic                               277.628     217.037
  Degrees of freedom                               122         122
  P-value (Chi-square)                           0.000       0.000
  Scaling correction factor                                  1.279
    Yuan-Bentler correction (Mplus variant)                       

Model Test Baseline Model:

  Test statistic                              6677.835    4874.889
  Degrees of freedom                               153         153
  P-value                                        0.000       0.000
  Scaling correction factor                                  1.370

User Model versus Baseline Model:

  Comparative Fit Index (CFI)                    0.976       0.980
  Tucker-Lewis Index (TLI)                       0.970       0.975
                                                                  
  Robust Comparative Fit Index (CFI)                         0.981
  Robust Tucker-Lewis Index (TLI)                            0.976

Loglikelihood and Information Criteria:

  Loglikelihood user model (H0)              -9306.387   -9306.387
  Scaling correction factor                                  1.398
      for the MLR correction                                      
  Loglikelihood unrestricted model (H1)      -9167.573   -9167.573
  Scaling correction factor                                  1.313
      for the MLR correction                                      
                                                                  
  Akaike (AIC)                               18710.774   18710.774
  Bayesian (BIC)                             18913.310   18913.310
  Sample-size adjusted Bayesian (SABIC)      18757.798   18757.798

Root Mean Square Error of Approximation:

  RMSEA                                          0.053       0.041
  90 Percent confidence interval - lower         0.044       0.033
  90 Percent confidence interval - upper         0.061       0.049
  P-value H_0: RMSEA <= 0.050                    0.291       0.970
  P-value H_0: RMSEA >= 0.080                    0.000       0.000
                                                                  
  Robust RMSEA                                               0.046
  90 Percent confidence interval - lower                     0.036
  90 Percent confidence interval - upper                     0.056
  P-value H_0: Robust RMSEA <= 0.050                         0.708
  P-value H_0: Robust RMSEA >= 0.080                         0.000

Standardized Root Mean Square Residual:

  SRMR                                           0.027       0.027

Parameter Estimates:

  Standard errors                             Sandwich
  Information bread                           Observed
  Observed information based on                Hessian

Latent Variables:
                   Estimate  Std.Err  z-value  P(>|z|)
  pj =~                                               
    q5_1              1.000                           
    q5_2              0.973    0.036   26.900    0.000
    q5_3              1.000    0.041   24.472    0.000
    q5_4              1.002    0.043   23.297    0.000
    q5_5              0.951    0.040   23.549    0.000
    q5_6              1.014    0.034   29.420    0.000
    q5_7              0.996    0.037   27.089    0.000
  dj =~                                               
    q8_1              1.000                           
  previous =~                                         
    int               1.000                           
  trust =~                                            
    q7_1              1.000                           
    q7_2              0.994    0.029   34.253    0.000
    q7_3              0.984    0.034   28.778    0.000
  eff =~                                              
    q6_1              1.000                           
    q6_2              1.023    0.048   21.230    0.000
    q6_3              1.009    0.050   20.287    0.000
  vig =~                                              
    q9_1              1.000                           
    q9_2              1.342    0.104   12.916    0.000
    q9_3              0.837    0.082   10.157    0.000

Regressions:
                   Estimate  Std.Err  z-value  P(>|z|)
  vig ~                                               
    pj        (h2)   -0.300    0.148   -2.031    0.042
    dj        (h4)   -0.002    0.076   -0.029    0.977
    trust     (h1)   -0.211    0.107   -1.966    0.049
    eff       (h3)    0.499    0.178    2.803    0.005
    previous         -0.023    0.028   -0.820    0.412

Covariances:
                   Estimate  Std.Err  z-value  P(>|z|)
  pj ~~                                               
    dj                0.772    0.056   13.743    0.000
    previous          0.093    0.054    1.736    0.082
    trust             0.789    0.054   14.606    0.000
    eff               0.700    0.053   13.120    0.000
  dj ~~                                               
    previous          0.080    0.061    1.322    0.186
    trust             0.899    0.057   15.776    0.000
    eff               0.742    0.056   13.320    0.000
  previous ~~                                         
    trust             0.049    0.059    0.829    0.407
    eff               0.063    0.053    1.191    0.234
  trust ~~                                            
    eff               0.731    0.055   13.290    0.000

Variances:
                   Estimate  Std.Err  z-value  P(>|z|)
   .q5_1              0.441    0.043   10.194    0.000
   .q5_2              0.303    0.035    8.671    0.000
   .q5_3              0.276    0.027   10.240    0.000
   .q5_4              0.372    0.032   11.767    0.000
   .q5_5              0.413    0.039   10.667    0.000
   .q5_6              0.307    0.029   10.755    0.000
   .q5_7              0.238    0.023   10.456    0.000
   .q8_1              0.000                           
   .int               0.000                           
   .q7_1              0.200    0.033    5.988    0.000
   .q7_2              0.268    0.043    6.225    0.000
   .q7_3              0.260    0.040    6.453    0.000
   .q6_1              0.356    0.038    9.319    0.000
   .q6_2              0.396    0.049    8.137    0.000
   .q6_3              0.391    0.044    8.957    0.000
   .q9_1              0.890    0.093    9.578    0.000
   .q9_2              0.207    0.078    2.640    0.008
   .q9_3              0.583    0.062    9.381    0.000
    pj                0.833    0.069   12.151    0.000
    dj                1.102    0.060   18.406    0.000
    previous          1.282    0.154    8.347    0.000
    trust             1.051    0.067   15.654    0.000
    eff               0.738    0.063   11.727    0.000
   .vig               0.581    0.085    6.873    0.000

Defined Parameters:
                   Estimate  Std.Err  z-value  P(>|z|)
    TRUST            -0.211    0.107   -1.966    0.049
    PJ               -0.300    0.148   -2.031    0.042
    EFFICACY          0.499    0.178    2.803    0.005
    DJ               -0.002    0.076   -0.029    0.977

Hypothesis testing:

  • PJ reduces vigilantism, p = .026.

  • Efficacy INCREASES vigilantism, p = .008.

  • DJ, victim, and previous exp are not associated with vigilantism.

  • Trust is significant (p = .035), and is negatively correlated with vigilantism.

Model
Estimate Std. Err. z p
Constructed
TRUST -0.21 0.11 -1.97 .049
PJ -0.30 0.15 -2.03 .042
EFFICACY 0.50 0.18 2.80 .005
DJ -0.00 0.08 -0.03 .977
Fit Indices
χ2 277.63
CFI 0.98
TLI 0.97
RMSEA 0.05
Scaled χ2 217.04(122) .000
+Fixed parameter


Plot of Structural model only:

Same plot, but with manifest variables:

Experiment data

  • We fitted a parsimonious model to the experimental data (see illustration below).

  • The data are assumed to come from a 2 * 2 between subject design.

    • The control subjects are omitted in this analysis.
    • If you have ideas how to model them, please let me know…
  • In this analysis we use \(N = 358\) observations.

CFA

  • As seen by the following CFA output, the data support our measurement model for vigilantism before and after the experimental manipulation.
measure vig_post vig_pre
alpha 0.8491 0.7687
omega 0.8535 0.7875
omega2 0.8535 0.7875
omega3 0.8555 0.799
avevar 0.5956 0.5578
Model
Estimate Std. Err. z p
Factor Loadings
vig_post
q15.1 1.00+
q15.2 1.05 0.07 15.56 .000
q15.3 0.94 0.09 9.95 .000
q15.4 1.24 0.10 12.76 .000
vig_pre
q9.1 1.00+
q9.2 1.26 0.10 12.42 .000
q9.3 0.90 0.10 8.90 .000
Residual Variances
q15.1 0.68 0.09 7.73 .000
q15.2 0.86 0.10 8.68 .000
q15.3 0.68 0.10 7.12 .000
q15.4 0.38 0.09 4.17 .000
q9.1 0.85 0.10 8.32 .000
q9.2 0.28 0.09 3.23 .001
q9.3 0.58 0.07 8.28 .000
Latent Variances
vig.post 0.85 0.11 7.55 .000
vig.pre 0.64 0.10 6.40 .000
Latent Covariances
vig.post w/vig.pre 0.35 0.07 5.23 .000
Fit Indices
χ2 61.64
CFI 0.95
TLI 0.92
RMSEA 0.10
Scaled χ2 56.24(13) .000
+Fixed parameter


SEM

  • The model fit the data well (see fit measures below).

Hypothesis testing:

  • The high PJ treatment is associated with a significant (p = .024) reduction self-vigilantismץ

  • The high avaןlability condition is not statistically significant, (z = -.68, p = .494).

  • Vigilantism reported before the experiment is significantly positively associated with the vigilantism after p < .001.