nature . human behaviour LETTERS PUBLISHED: 30 JANUARY 2017 | VOLUME: 1| ARTICLE NUMBER: 0040 Economic insecurity and the rise in gun violence at US schools A. R. Pah’, J. Hagan?4, A. L. Jennings?, A. Jain°®, K. Albrecht3, A. J. Hockenberry® and L. A. N. Amaral'*7”* Frequent school shootings are a unique US phenomenon that has defied understanding"?. Uncovering the aetiology of this problem is hampered by the lack of an established dataset?“. Here we assemble a carefully curated dataset for the period 1990-2013 that is built upon an exhaustive review of exist- ing data and original sources. Using this dataset, we find that the rate of gun violence is time-dependent and that this rate is heightened from 2007 to 2013. We further find that periods of increased shooting rates are significantly correlated with increases in the unemployment rate across different geographic aggregation levels (national, regional and city). Consistent with the hypothesis that increasing uncertainty in the school-to-work transition contributes to school shoot- ings, we find that multiple indicators of economic distress significantly correlate with increases in the rate of gun vio- lence when events at both K12 and post-secondary schools are considered. Although there is extensive work on school shootings'**"", the result is a patchwork of contradictory claims. Some studies report an insignificant increase in the rate of school shooting incidents or homicides over time?*°, while others point to copycat effects and an increasing frequency in certain population segments and settings*'®. Going beyond these basic insights muddies the waters further, as the identification of sociological and risk factors related to being a shooter!" are of doubtful predictive value because of their lack of specificity to individuals that commit gun violence at schools*!*'8, Resolving the differing perspectives on gun violence at schools is challenging for several reasons. Not least is the lack of a definitive data source with clear event inclusion criteria. Event counts of gun violence at schools vary dramatically across sources (Fig. 1a), which creates concern about the reliability of previous quantitative analy- ses. Even simple questions—such as whether the rate of shootings is increasing—are impossible to answer without valid data. Multiple datasets with different inclusion criteria were used in previous research on mass killings, shootings and gun violence at schools (Fig. 1). To create a consistent dataset to investigate the phenomenon of gun violence at schools, we advance the following criteria for event inclusion: (1) the shooting must involve a firearm being discharged, even if by accident; (2) it must occur on a school campus; and (3) it must involve students or school employees, either as perpetrators, bystanders or victims. As an example, gang violence on a playground at night during the summer months would not be included since it violates the last criterion, while a student being shot at the school’s baseball field after a game is included. To build our dataset, we merged events from six original datasets pertaining to school violence (Fig. 1), resulting in 535 events. We cleaned the merged dataset and corrected dates using primary news sources, resulting in 529 events for potential inclusion. Three coders then independently evaluated each individual event against the defined criteria. If at least two coders agreed that an event should be included, then the event was added to the final consensus dataset. This process yielded 379 events meeting our strict criteria and two additional events found during the discovery process that were not present in any of the original six datasets. Our consensus dataset is positively correlated with all of the original datasets—demonstrating that our criteria do not exclude events from any one dataset at a noticeably higher rate than any of the others (Supplementary Fig. 2). We categorized the events within the consensus dataset to gain a more concrete understanding of what constitutes gun violence at schools (Fig. 2). Consistent with previous reports*”’, we find that most events are targeted, that is, the shooter intends to harm a specific person. In our dataset, gang-related violence constitutes 6.6% of all incidents; this is a much smaller fraction compared with what is observed for urban violence outside of schools”. The average number of fatalities per event is one and the number of incidents with three or more deaths constitutes 6.3% of included events. Moreover, gun violence at schools has not become more deadly over time (Fig. 3a and Supplementary Table 3). It is also important to note that that this dataset is focused on all gun violence at schools and is not limited to mass shootings. Although there are notable mass shooting events on school cam- puses, most mass shootings happen at locations other than schools. Similarly, this dataset includes all instances of gun usage, whether someone dies in the course of the event or not, since the discharge of a gun is not permitted on school campuses (despite allowances ‘to carry on some state college campuses”'). The inclusion of attempted violence distinguishes this dataset from other measures of violence, such as the homicide rate, since that rate is by definition only concerned with acts of violence resulting in a death. We next evaluate the timing of these events to determine whether they follow a Poisson process. Since our dates cover an extended period of time, we allow for the possibility that the rate parame- ter, A, varies over time in a stepwise manner. We fit models with an increasing number of change-points to the monthly time series of events and find that the best fit has four distinct periods sepa- rated by three change-points (see Supplementary Information for methods, parameter values and information criteria scores; Fig. 3a for model fit). ‘Northwestern Institute on Complex Systems (NICO), Northwestern University, Evanston, Illinois 60208, USA. *Management and Organizations, Kellogg School of Management, Northwestern University, Evanston, Illinois 60208, USA. ?Department of Sociology, Northwestern University, Evanston, Illinois 60208, USA. “American Bar Foundation, Chicago, Illinois 60611, USA. °Department of Chemical and Biological Engineering, Northwestern University, Evanston, Illinois 60208, USA. Department of Physics and Astronomy, Northwestern University, Evanston, Illinois 60208, USA. Department of Medicine, Northwestern University, Evanston, Illinois 60208, USA. *e-mail: amaral@northwestern.edu NATURE HUMAN BEHAVIOUR 1, 0040 (2017) | DOI: 10.1038/s41562-016-0040 | www.nature.com/nathumbehav 1 © 2017 Macmillan Publishers Limited, part of Springer Nature. All rights reserved. LETTERS NATURE HUMAN BEHAVIOUR a ~ 404 / Brady Campaign 40 4 SAVD Ss 2 events 236 events £2 301 /- 304 o 3 £ 204 Lo 20 4 8 = 104 Lo 10 4 < 0-4 0+ 1990 1995 2000 2005 2010 2015 1990 1995 2000 2005 2010 2015 ~ 407 Shultz 40 + Slate SF 215 events 166 events. wn 7 € 304 304 3 VA. VA © 204 20 4 3 a y = 104 10 4 < o- o+ 1990 1995 2000 2005 2010 2015 1990 1995 2000 2005 2010 2015 4 Virginia Tech 4 Wikipedia S 40 aan 40 174 events 2 304 / 304 vo 3 € 204 Lo 20 4 mi) a = 104 10 4 and that the attitudes of youths have a significant impact on their future employment prospects and earnings**°. We posit that gun violence at schools is a response, in part, to the breakdown in the expectation that sustained participation in the educational system will improve economic opportunities and outcomes. 2 NATURE HUMAN BEHAVIOUR 1, 0040 (2017) | DOI: 10.1038/s41562-016-0040 | www.nature.com/nathumbehav © 2017 Macmillan Publishers Limited, part of Springer Nature. All rights reserved. NATURE HUMAN BEHAVIOUR LETTERS a a _— 70- 50- nv Change-point model 60 + 40 | Z 10- % 07 & 5 8 £ 404 #30 : & 30 8 594 £ 64 oO o = 2 204 a 10 = 44 104 = 21 0s O° > 0 ‘S < x &< & 4. . Le 7 se sx SF M&M & < 1990 1995 2000 2005 2010 S we ne <2 Year ; ; b 12 Figure 2 | Breakdown of urbanity and event type in the consensus 3 - dataset. a, The urbanity of the events mirrors the population distribution Ss 10 wn oO in the United States (6% rural, 26% suburban, 67% urban), with 5 6-4 3 3 most events occurring in urban areas. b, Events categorized based 3 cy Cc on the reported outcome and shooter intent (when available). 344 6 a Targeted events (an attacker intending to shoot one or two others) = > a . . fo} 4 x is the primary event type. = 4s 0-5 2 National level 1990 1995 2000 2005 2010 Since we hypothesize that increased school shootings are a response Year to increasing unemployment, we fit the data using Poisson regression: fc 4) 95% Cl E[Sypltms Ms] = eF0tP itt Poms (1) < 6. 5 we ; where S,, is the number of shootings per month, u,, isthe monthly % 44 ese P e208 unemployment rate, m, is a dummy variable that accounts for the = @| eee eee @| ee emmme @o ww summer months, E is the expected value, and f,, 6,, and 8, arethe § 27 Cmnen G00 WHOMOS GO C8 © O@ OEOw parameters being estimated. We find a significant (P<10~, pseudo- = aero © neem oman ¢ Ene © eae : : + © SER STS 1S Tre ao~ R’?=0.074) relationship between the unemployment rate and num- 0 | —_ oe ber of incidents per month (95% confidence intervals (CIs) shown in Fig. 3c and Supplementary Table 8). Although the pseudo-R? value is low, this is largely due to the inherent noise in Poisson pro- cesses—the fit captures 53% of the maximum variance that would be expected for this number of observations even if there was a perfect correlation between unemployment and the number of inci- dents per month (Supplementary Fig. 5). The unemployment rate is still a significant predictor if we control for the change in student population over time (Supplementary Table 9). To further confirm the robustness of this finding, we test our hypothesis in two additional ways. First, we model the relationship of the average time between events to unemployment. Using this formulation, we find again that there is a significant relationship between increasing unemployment and decreasing time between event incidences (P=0.011, R’=0.10; Supplementary Table 10). Second, we normalize the unemployment rate into the range [0,1] during the time period studied and categorize the months based on the number of shootings within each month. Since the period (1994-2007) with a lowered rate of shootings has an aver- age of approximately one shooting per month, we use that number of shootings per month as a threshold to separate the two groups. If unemployment is a factor in school shootings then we would expect that months with more shootings would have a significantly larger mean normalized unemployment rate. This is indeed what we observe. We find that the two distributions are significantly dif- ferent and that months with two or more events have a larger mean normalized unemployment rate (0.43 versus 0.35; Kolmogorov— Smirnov (K-S) two-sample test, P=0.006; Fig. 3d). Next, we test our hypothesis at different levels of spatial aggrega- tion to assess whether this relationship is conditional on location or might arise from an ecological fallacy*’. Regional level We partition the continental United States into seven regions according to geography and socioeconomic similarity (Fig. 4a and Supplementary T T T 2000 2005 2010 Year T T 1990 1995 a 2.5 = @<1event @ >1 event Probability density 0.0 O1 02 03 Normalized unemployment rate 04 05 06 07 O08 09 Figure 3 | The rate of school shootings is time-dependent and correlated with increasing unemployment above ‘normal’ levels. a, The monthly number of events categorized based on number of fatalities (green O-1, orange 2-5 and red >5). We fit Poisson process change-point models to the monthly incident time series and find that the best fit model has four distinct periods (see Supplementary Table 7). b, National unemployment rate peaks (black line) qualitatively align with periods of elevated rates of school shootings (blue bars). ¢, Confidence intervals for the fit of national unemployment rate to monthly shootings (blue dots) (see Equation (1)). d, The distributions of normalized unemployment in months with <1 and >1 event differ significantly and months with >1 shooting have a larger mean normalized unemployment. Fig. 7). We examine the distribution of normalized unemployment rates, with each region having its unemployment normalized into the range [0,1] individually. Due to the lower frequency of events at a regional scale, we partition months into those with no shootings and those with shootings. As before, we find that months with one or more shootings have a normalized unemployment rate distribution that sig- nificantly differs with a larger mean normalized unemployment rate (0.41 versus 0.37; K-S two-sample test, P=0.017; Fig. 4b). NATURE HUMAN BEHAVIOUR 1, 0040 (2017) | DOI: 10.1038/s41562-016-0040 | www.nature.com/nathumbehav 3 © 2017 Macmillan Publishers Limited, part of Springer Nature. All rights reserved. a S 12 5 Northeast 9 | S 12 5 Southeast We ¥ 104 8 ® 4 104 9 38 5 84 73 5 a 8 3 3s o 3 co b 4 < Vv 4 < 2 6 6 3 2 6 7 3 B No events = 44 5 3 2 44 6 3 » 2556 Events 2 24 43 © 24 Ee 8 iS} x [e} & 4 2.0 = 04 3° = 04 4 ~- o 1990 1995 2000 2005 2010 1990 1995 2000 2005 2010 3 15 = 10 Se 125 Great lakes ae Se 125 Pacific Ze g S ns + nz 9 05 4» 104 10 2 4 10-4 ic) a 5 3 92 5 a1 pe 9.0 3 6 8 g 3 6 4 3 Sg 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 = a Z 3 = a4 7 a Normalized unemployment rate s 5 6% e 24 48 5§ 27 5 § = 04 3“ = 04 4 ~ 1990 1995 2000 2005 2010 1990 1995 2000 2005 2010 Year Year Figure 4 | Elevated unemployment is correlated with an increased rate of school shootings at the regional level. a, Average unemployment (black lines) and school shooting incidents (blue bars) for the four most active regions (all eight regions shown in Supplementary Figure 7). b, When the distributions of these normalized unemployment rates are compared, we find that there is a significant difference and months with shootings have an over-representation of months with a high normalized unemployment rate. City level We analyse the six cities with the most gun violence at schools: New York City, Detroit, Chicago, Memphis, Los Angeles and Houston (Fig. 5a). As for the national and regional levels, we find that months with one or more shootings have a normalized unemployment rate distribution that significantly differs and has a larger mean normalized unemployment rate (0.51 versus 0.41; K-S two-sample test, P=0.005; Fig. 5b). Educational attainment Our results strongly support the hypothesis that a breakdown in the school-to-work transition contributes to an increase in gun violence at schools. Taking this hypothesis a step further, we would expect that there would be a shift in the temporal location of these shootings during the period when post-secondary education has increasingly supplanted high school in determining successful school-to-work transitions*’*”. When we analyse the post-secondary event series separately, we do find that the rate of gun violence is elevated from November 2005 to December 2013 (Supplementary Fig. 8 and Supplementary Table 15). When these individual time series are fit against correspond- ing unemployment metrics (‘less than high school’ unemploy- ment levels for K12 schools and ‘some college’ unemployment for 4a 6, Chicago r13 4 ~ 64 Los Angeles -14 7 s 12 & s & 2°] bi wS 2°] b 12 WS g4t plas B44 lio a3 S +9 50 3 ie} £37 lg So = 34 lg #2 2 24 17 23 224 <3 z 1 Fé sa = 44 Le 38 L 5 0 09 ele Li bss 27) Uh Lie 1990 1995 2000 2005 2010 1990 1995 2000 2005 2010 2 65 Houston ro Ss 2 65 Detroit [ #0 > b 2.5 Noevents So] Leg Sec] r & 2 B Events @ 45 t7 33 @ 45 p14 38 S45 s 30S -f12 50 > £37 ;&§ Sg 2 34 Lio 32 24 L ge 21 mimirsé 23) 4 HL fs 8" 20, 0 T T T 3 0 T T 2~ 0.0 1990 1995 2000 2005 2010 1990 1995 2000 2005 2010 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 6 New York 10 6 M hi n Normalized unemployment rate ~ 67 ew Yor r 10 ~ 69 emphis rc A s 5 J Lo & Ss 54 L10 & 2 lg wS 2 Lg wS o 44 32 6 44 3 3 L782 8 re 82 23) Po es 234 b? 38 wo bt og = 27 i5 ag £71 83 ae Laas 514 A at 2o/1 Jl {1 13% =2,| ll | IIs ¢ 1990 1995 2000 2005 2010 1990 1995 2000 2005 2010 Year Year Figure 5 | Elevated unemployment is significantly correlated with an increased rate of school shootings in the most active cities. a, The city-wide unemployment rate (black lines) and school shooting incidents (blue bars) for the six cities with the most events. b, The distributions of normalized unemployment rates between months with and without a shooting differ significantly and months with a shooting have a larger mean normalized unemployment rate. 4 NATURE HUMAN BEHAVIOUR 1, 0040 (2017) | DOI: 10.1038/s41562-016-0040 | www.nature.com/nathumbehav © 2017 Macmillan Publishers Limited, part of Springer Nature. All rights reserved. NATURE HUMAN BEHAVIOUR LETTERS 85 K12 25 0.20 5 e774 20 3 # 2 64 15 3 § 0154 s 5 4 1.0 g ‘S £44 05 8 2 010+ > 34 00 2 = ny E € 24 054 60.054 fo} fal S 214 -10 § & o4 15 ~ 0.004 1990 1995 2000 2005 2010 44 Post-secondary p25 _, 0.85 fo} S L 2.0 3 = 0.7 5 ¥% 34 lis 9 @& 067 c . Pat 2 Lio a § 057 £27 os 5 044 = p