Who Will Save You Now?

Evan Alexander Peters, Samra McCullin, Sabir Hashim

College of Science, George Mason University, Fairfax, VA, USA

Address: College of Science, George Mason University, 4400 University Drive, Fairfax, VA

This study investigates the spatial and demographic factors influencing gun violence in the U.S. by first creating weighted Kernel Density Estimation (KDE) maps to normalize centroid and coordinate-specific incidents and rates where available. These weighted coefficients were then used in regression analyses and hot/cold spot evaluations to identify significant predictors. A comprehensive "All Violence Rate" map was generated as the primary visualization, complemented by four subset maps representing individual violence types, each highlighting their significant relationships and spatial patterns. Finally, an additional map was created to illustrate the impact of economic factors, particularly median income, on the overall violence rate

Brief Presentation on Firearm Incident Risk Analysis

Introduction: General Results Overview

  1. Purpose of Analysis:

    • This study examines the spatial dynamics of firearm violence across various regions in the U.S., analyzing how factors like population density, income, and proximity to churches and schools influence violence rates.

    • Key types of violence examined include mass shootings, homicides, suicides, and injuries.

  2. Main Predictors:

    • Population Density: Areas with higher density tend to have elevated violence rates, likely due to urban pressures and social proximity.

    • Income Density: Surprisingly, areas with higher income density also showed a positive relationship with violence rates, particularly suicides, highlighting socio-economic stressors.

    • Area Per Person: Larger area per person, often found in rural settings, was negatively associated with violence rates, suggesting that lower density may act as a protective factor.

  3. Churches and Schools:

    • Church Proximity: Found to have a weak, statistically significant negative effect on violence, suggesting localized dampening of violence near religious sites.

    • School Proximity: Showed a small but significant positive relationship, reflecting higher violence rates in urban areas where schools are more densely located.

  4. Model Fit and Diagnostics:

    • The model explains a substantial portion of the variance (R² ~51%) but highlights geographic disparities and unmodeled factors influencing violence.

    • Diagnostics show heteroskedasticity and non-normality, indicating regional and cultural dynamics not fully captured.

Key Findings by Violence Type

  1. Education Density:

    • One of the most consistent and significant mitigators of violence. Higher density of schools correlates with reduced homicide and suicide rates but is positively associated with urban-related violence due to population clustering.

  2. Religious Density:

    • A mixed relationship, with slight protective effects against homicides but positive associations with injury and suicide rates, suggesting complex socio-cultural factors at play.

  3. Spatial Patterns:

    • Hot spots in urban areas with high population and income density.

    • Cold spots in rural regions with lower density and larger area per person.

Conclusion: Church vs. School Proximity

  1. General Safety Dynamics:

    • The analysis reveals collinearity between school and church density with other variables like population and income. These factors heavily influence violence rates, making it difficult to isolate their unique effects.

    • Near schools, the increased violence rate is tied to urban settings, while near churches, the slight decrease suggests localized community cohesion or deterrence.

  2. Critical Assumptions:

    • Proximity alone does not define safety; broader socio-economic and cultural dynamics drive these patterns.

    • For schools, urban clustering might mask their protective societal role.

    • For churches, varying regional cultural attitudes toward religion could shape the observed dampening effect.

  3. Takeaway:

    • While churches may show a slight advantage for localized violence reduction, the broader context indicates that population density and socio-economic inequalities are more significant determinants of safety.

    • This nuanced relationship emphasizes the need for deeper, targeted strategies in urban planning and community resource allocation.

Key Findings by Violence Type

  1. Education Density:

    • One of the most consistent and significant mitigators of violence. Higher density of schools correlates with reduced homicide and suicide rates but is positively associated with urban-related violence due to population clustering.

  2. Religious Density:

    • A mixed relationship, with slight protective effects against homicides but positive associations with injury and suicide rates, suggesting complex socio-cultural factors at play.

  3. Spatial Patterns:

    • Hot spots in urban areas with high population and income density.

    • Cold spots in rural regions with lower density and larger area per person.

Conclusion: Church vs. School Proximity

  1. General Safety Dynamics:

    • The analysis reveals collinearity between school and church density with other variables like population and income. These factors heavily influence violence rates, making it difficult to isolate their unique effects.

    • Near schools, the increased violence rate is tied to urban settings, while near churches, the slight decrease suggests localized community cohesion or deterrence.

  2. Critical Assumptions:

    • Proximity alone does not define safety; broader socio-economic and cultural dynamics drive these patterns.

    • For schools, urban clustering might mask their protective societal role.

    • For churches, varying regional cultural attitudes toward religion could shape the observed dampening effect.

  3. Takeaway:

    • While churches may show a slight advantage for localized violence reduction, the broader context indicates that population density and socio-economic inequalities are more significant determinants of safety.

    • This nuanced relationship emphasizes the need for deeper, targeted strategies in urban planning and community resource allocation.

This presentation reflects the complexity of firearm violence and underscores the interplay between spatial, social, and economic factors in shaping safety near schools and churches.

Summary of Getis Ord Gi* Overall Firearm Incident Risk Analysis


Key Takeaways:

  1. Intercept:

    • Coefficient: 0.002499

    • Indicates the baseline weighted violence rate when all explanatory variables are at their average values. It is positive and statistically significant.

  2. Key Predictors:

    • Area Per Person (AREAPP_KDEZONALMEAN):

      • Coefficient: -5.804769

      • Strong negative relationship with overall violence rate, suggesting less densely populated areas have lower violence rates.

    • Population Density (POPULATION_KDEZONALMEAN):

      • Coefficient: 0.000015

      • Positive and highly significant, indicating that higher population density increases the weighted violence rate.

    • Income Density (INCOME_KDEZONALMEAN):

      • Coefficient: 0.000153

      • Strong positive relationship, suggesting that areas with higher income density have a slightly increased violence rate.

    • ReligiousAssembly_point:

      • Coefficient: -0.000000

      • Weak but statistically significant negative effect, suggesting localized dampening of violence near religious points.

    • EducationCenters_point:

      • Coefficient: 0.000002

      • Small but statistically significant positive effect, reflecting increased violence in denser urban educational zones.

  3. Model Diagnostics:

    • R-Squared: 0.511538

      • Indicates that the model explains about 51% of the variance in the weighted violence rate. This is a solid fit given the complexity and variability of the data.

    • Koenker Statistic: Highly significant, suggesting heteroskedasticity. Robust probabilities should be relied on.

    • Jarque-Bera Statistic: Highly significant, indicating that residuals are not normally distributed. Additional geographic or unmodeled factors may exist.


      Summary of Results for the Income-Violence Model

      Map 3.1

      {

      Response Variable: KDE_MeanWeighted_Violence

      Explanatory Variable: Income_KDEzonalMEAN

      }

      Key Takeaways

      1. Intercept (0.000976):

        • Positive and statistically significant (p < 0.01).

        • Represents the baseline violence rate when income density is at its mean value.

      2. Income Density (INCOME_KDEZONALMEAN):

        • Coefficient (0.000204): Positive and highly significant (p < 0.01).

        • Indicates a strong positive relationship between income density and violence rates. This suggests that areas with higher income density tend to experience slightly elevated violence rates, possibly reflecting socio-economic inequalities or stressors.

      Model Diagnostics

      • R-Squared (0.459528):

        • Indicates that the model explains ~46% of the variance in violence rates using income density.

        • A reasonable fit for a single-variable regression.

      • Adjusted R-Squared (0.459356):

        • Adjusts for the number of predictors; consistent with the R-Squared.

      • Akaike's Information Criterion (AICc):

        • Value of -22,190.17 indicates model fit; lower values signify better fit.

      • Joint F-Statistic (2662.087469):

        • Highly significant (p < 0.01), indicating that income density as a predictor contributes significantly to the model.

      • Koenker (BP) Statistic (234.994106):

        • Significant (p < 0.01), indicating heteroskedasticity (non-constant variance). Use robust probabilities (Robust_Pr) for interpreting variable significance.

      • Jarque-Bera Statistic (45,676.423636):

        • Significant (p < 0.01), suggesting non-normality in residuals. Geographic or unmodeled factors may influence violence rates.

      Key Insights

      • Income Density's Influence:

        • Higher income areas exhibit slightly elevated violence rates, potentially reflecting complex socio-economic dynamics. While income may indicate affluence, it could also highlight disparities within those areas that contribute to violence.

      • Model Fit:

        • The model provides a meaningful explanation of the relationship between income and violence but leaves significant variance unexplained, likely due to other contributing factors like education, population density, or cultural dynamics.


          Getis Ord Gi* = Weighted Gun Violence Analysis

          Map 1.1

          • WHERE TARGET VARIABLE IS

            KDE_MeanWeighted_Violence = {

            (0.35 * !MassShooting_KDEzonalMEAN!) +

            (0.25 * !Homicide_KDEzonalMEAN!) +

            (0.25 * !Suicide_KDEzonalMEAN!) +

            (0.15 * !Injury_KDEzonalMEAN!)

            }

        • AND AND EXPLANATORY VARIABLES ARE

          AreaPP_KDEzonalMEAN

          Population_KDEzonalMEAN

          Income_KDEzonalMEAN

          ReligiousAssembly_point

          EducationCenters_point

So, where is safer? This study's Overall Analysis (Map 1.1 and OLS Regression Results) is based on county-level centroid data derived from KDE means for the variables in our study. The results reveal contrasting impacts of education and religious density on violence rates. Educational density generally correlates negatively with all violence types, significantly reducing suicide, homicide, and injury incidents, but remains a contributor to mass casualty events. Conversely, religious density shows mixed results, positively associated with injuries and negatively linked to suicide and homicide rates, while remaining insignificant in mass casualty incidents. These findings suggest that statistically, educational institutions broadly mitigate violence risks, while religious buildings exhibit more variable impacts depending on the violence type, highlighting their complex community role.

It should be noted that sub-chart analysis is driven by the result of individual Kernel Density Estimations prioritizing the precise locations of educational buildings, religious congregations, and mass firearm violence with increased scrutiny for generalized public safety, showing that proximity to religious congregations, on average, offers a safer environment compared to educational institutions. These factors are significant, but of weak correlation.

Summary of Kernel Density Results for Firearm Violence Type


Key Takeaways:

  1. Significance of Variables:

    • Across all four target variables (Mass Shooting, Homicide, Suicide, and Injury KDE densities), certain explanatory variables consistently show statistically significant relationships, while others are less impactful or redundant:

      • Area Per Person (AREAPP):

        • Shows a significant positive or negative relationship depending on the target variable, indicating that population density plays a role in influencing violence rates.

        • A larger area per person is negatively associated with homicide and suicide rates, suggesting that areas with lower population density may experience lower violence rates.

      • Education (EDUCATION_KDE):

        • One of the strongest predictors across all target variables, with a consistently negative relationship for Suicide, Homicide, and Injury rates.

        • Higher density of educational buildings correlates with reduced violence, likely due to socio-economic or educational factors that mitigate crime and self-harm.

      • Income (INCOME_KDE):

        • Positively correlated with Suicide rates but less significant for Mass Shooting and Homicide rates. This suggests that higher income areas might experience specific socio-psychological stressors related to suicide.

        • Its lack of strong correlation with Homicide and Mass Shootings might indicate income inequalities or other socio-economic factors not captured by the model.

      • Religious Building Density (RELIGIOUS_KDE):

        • A significant positive predictor for some outcomes, such as Suicide and Injury rates, potentially reflecting complex societal or cultural dynamics.

        • Conversely, it has weaker significance for Homicide rates, suggesting religion may not always act as a mitigating factor for violence.

      • Population Density (POPULATION_KDE):

        • Often a secondary predictor, it shows mixed significance. In areas with high residuals, it may reflect either under- or over-prediction by the model.

      • Homicide, Suicide, and Injury KDE Means:

        • These variables frequently interact with each other as predictors, showing that one form of violence is often strongly correlated with another (e.g., suicide and homicide tend to co-occur in regions).

  2. Model Fit:

    • All models (Mass Shooting, Suicide, Homicide, Injury) show moderate-to-strong R² values (e.g., ~0.45–0.72), indicating the explanatory variables capture a substantial portion of the variance in the target variables.

    • Education density and income density consistently contribute to the model’s explanatory power, with education density emerging as the most robust and impactful predictor overall.

  3. Variance Inflation Factor (VIF):

    • VIF values highlight possible collinearity between variables:

      • Education and Religious Density (VIF > 7.5) show redundancy in some cases, suggesting overlapping socio-environmental contexts.

      • Income and Population Density (VIF ~3-5) also indicate minor collinearity but are still interpretable within the model.

  4. Model Diagnostics:

    • Koenker (BP) statistic significance: Suggests heteroskedasticity (non-constant variance), so you should rely on robust probabilities for interpretation.

    • Jarque-Bera statistic significance: Indicates residuals are not normally distributed, suggesting potential for geographic or unmodeled influences.

Key Relationships to Mention

  1. Education as a Mitigator of Violence:

    • Higher density of educational institutions is a consistent and statistically significant negative predictor for multiple violence types (Suicide, Homicide, Injury).

    • This indicates education access or proximity may have a protective influence against violence and self-harm, likely through socio-economic upliftment, community resources, or awareness.

  2. Area Per Person and Violence Dynamics:

    • Low population density areas (large area per person) are often associated with lower homicide and suicide rates, while denser urban regions may amplify risk factors, especially for homicides.

    • This suggests urban planning and population density play important roles in violence mitigation strategies.

  3. Income and Suicide Dynamics:

    • Income density is positively associated with suicide rates, highlighting the potential socio-economic and psychological stressors unique to higher-income communities.

    • Addressing mental health resources in affluent regions could be an actionable recommendation.

  4. Religious Density as a Complex Factor:

    • Religious building density shows mixed relationships:

      • Positively associated with Suicide and Injury rates, potentially reflecting social stressors or cultural influences.

      • Not significantly protective for Mass Shootings or Homicides, suggesting it may not act as a universal mitigator of violence.

  5. Mass Shootings:

    • Relationships are less clear for Mass Shootings, with mixed significance for most predictors.

    • Localized factors, such as access to firearms, political/cultural tensions, or media influence, may require additional data for better modeling.

  6. Spatial and Regional Disparities:

    • Residual maps indicate geographic regions where the model over- or under-performs:

      • Hot spots: Areas with higher violence rates than predicted suggest missing variables or localized factors (e.g., systemic inequality, gun access laws, or cultural norms).

      • Cold spots: Areas with lower violence rates than predicted may highlight protective community features (e.g., tight-knit communities, better law enforcement).

        Kernel Density Maps and Sub-Charts

        Maps (2.1, 2.2, 2.3, 2.4)

  • WHERE TARGET VARIABLE(s) ( non-redundant )*

    MassShooting_KDEzonalMEAN

    Homicide_KDEzonalMEAN

    Suicide_KDEzonalMEAN

    Injury_KDEzonalMEAN

    • AND EXPLANATORY VARIABLES ARE

      {

      Population_KDEzonalMEAN

      AreaPP_KDEzonalMEAN

      Income_KDEzonalMEAN

      Income_KDEzonalMEAN

      Religious_KDEzonalMEAN

      Education_KDEzonalMEAN

      ReligiousAssembly_point

      EducationCenters_point
      INCLUSION* for collinear analysis
      MassShooting_KDEzonalMEAN

      Homicide_KDEzonalMEAN

      Suicide_KDEzonalMEAN

      Injury_KDEzonalMEAN

      }

Source Data

Gun Violence Analysis in ArcGIS Pro - Data Sources

  1. Religious Data

    • Source: 2020 U.S. Religion Census, Association of Statisticians of American Religious Bodies.

    • Data Type: Tabular data with summary statistics by county.

    • Attributes: Includes the number of congregations and adherents by county.

    • Link: 2020 U.S. Religion Census Summaries

  2. Gun Violence Data

    • Source: Centers for Disease Control and Prevention (CDC), National Center for Health Statistics (via WONDER Online Database); Gun Violence Archive (street address level).

    • Data Type: Tabular data by county, with classifications for homicide, suicide, injury/other, and mass shooting events.

    • Temporal Coverage: 2018–2022.

    • Links:

  3. County Health Rankings

    • Source: Robert Wood Johnson Foundation and University of Wisconsin Population Health Institute.

    • Data Type: Polygon data for health rankings by county.

    • Temporal Coverage: 2015–2019.

    • Link: County Health Rankings Data and Documentation

  4. Structure Data for Educational and Religious Sites

    • Source: Oak Ridge National Laboratory (ORNL); FEMA Geospatial Response Office.

    • Data Type: Polygon data converted to points for ease of analysis.

    • Data Currency: November 2023.

    • Access: This data is available through FEMA’s USA Structures State GDB on the Esri Portal. (Direct link not available; accessible via portal navigation.)

  5. County Boundaries

    • Source: U.S. Census Bureau.

    • Data Type: Polygon data for U.S. county boundaries.

    • Attributes: Geographic boundaries for all U.S. counties.

    • Link: 2023 U.S. County Boundaries

  6. American Community Survey (ACS) Data

    • Source: U.S. Census Bureau, ACS.

    • Data Type: Tabular data by county.

    • Attributes: Includes median household income for the years 2018–2022.

    • Link: County 2018-2022 ACS Data

References

[[1]        J. Gramlich, “What the data says about gun deaths in the U.S.,” Pew Research Center. Accessed: Nov. 11, 2024. [Online]. Available: https://www.pewresearch.org/short-reads/2023/04/26/what-the-data-says-about-gun-deaths-in-the-u-s/

[2]        E. Grinshteyn and D. Hemenway, “Violent Death Rates: The US Compared with Other High-income OECD Countries, 2010,” Am. J. Med., vol. 129, no. 3, pp. 266–273, Mar. 2016, doi: 10.1016/j.amjmed.2015.10.025.

[3]        W. Stroebe, N. P. Leander, and A. W. Kruglanski, “Is It a Dangerous World Out There? The Motivational Bases of American Gun Ownership,” Pers. Soc. Psychol. Bull., vol. 43, no. 8, pp. 1071–1085, Aug. 2017, doi: 10.1177/0146167217703952.

[4]        J. Kreienkamp, M. Agostini, N. P. Leander, and W. Stroebe, “How news exposure and trust in law enforcement relate to defensive gun ownership.,” Psychol. Violence, vol. 11, no. 4, pp. 417–427, Jul. 2021, doi: 10.1037/vio0000375.

[5]        J. Beardslee, M. Docherty, E. Mulvey, and D. Pardini, “The Direct and Indirect Associations between Childhood Socioeconomic Disadvantage and Adolescent Gun Violence,” J. Clin. Child Adolesc. Psychol., vol. 50, no. 3, pp. 326–336, May 2021, doi: 10.1080/15374416.2019.1644646.

[6]        D. Bergen-Cico et al., “Community Gun Violence as a Social Determinant of Elementary School Achievement,” Soc. Work Public Health, vol. 33, no. 7–8, pp. 439–448, Nov. 2018, doi: 10.1080/19371918.2018.1543627.

[7]        P. M. Reeping, A. Mak, C. C. Branas, A. N. Gobaud, and M. L. Nance, “Firearm Death Rates in Rural vs Urban US Counties,” JAMA Surg., vol. 158, no. 7, pp. 771–772, Jul. 2023, doi: 10.1001/jamasurg.2023.0265.

[8]        A. M. Banks, “Attack part of pattern against black churches,” National Catholic Reporter, vol. 51, no. 19, National Catholic Reporter Publishing Company, Kansas City, United States, p. 13, Jul. 03, 2015.

[9]        Andrew L. Whitehead, “Gun Control in the Crosshairs: Christian Nationalism and Opposition to Stricter Gun Laws.” Accessed: Nov. 11, 2024. [Online]. Available: https://journals.sagepub.com/doi/epub/10.1177/2378023118790189

[10]      H. Schramm, “Operations Research Responses to Mass Shootings,” Phalanx, vol. 46, no. 1, pp. 30–35, 2013.

[11]      B. K. Roberts, C. P. Nofi, E. Cornell, S. Kapoor, L. Harrison, and C. Sathya, “Trends and Disparities in Firearm Deaths Among Children,” Pediatrics, vol. 152, no. 3, p. e2023061296, Sep. 2023, doi: 10.1542/peds.2023-061296.

[12]      J. Tinnes, “Bibliography: Internet-Driven Right-Wing Terrorism,” Perspect. Terror., vol. 14, no. 3, pp. 168–189, 2020.

[13]      S. Hofer, “Lockout: Spacing Trauma and Recovery in the Aftermath of the Virginia Tech Shootings,” Am. Imago, vol. 72, no. 3, pp. 231–283, 2015.

[14]      J. Schildkraut and L. B. Geller, “Mass Shootings in the United States: Prevalence, Policy, and a Way Forward,” Ann. Am. Acad. Pol. Soc. Sci., vol. 704, no. 1, pp. 181–203, Nov. 2022, doi: 10.1177/00027162231164484.

[15]      H. Schramm, “Operations Research Responses to Mass Shootings,” Phalanx, vol. 46, no. 1, pp. 30–35, 2013.

[16]      D. D. McDaniel, “Overcoming our displaced fears: THE NEW RELIGIOUS INTOLERANCE: OVERCOMING THE POLITICS OF FEAR IN AN ANXIOUS AGE,” National Catholic Reporter, vol. 49, no. 8, National Catholic Reporter Publishing Company, Kansas City, United States, p. 16, Feb. 01, 2013.

[17]      A. Mirpuri, “Racial Violence, Mass Shootings, and the U.S. Neoliberal State,” Crit. Ethn. Stud., vol. 2, no. 1, pp. 73–106, 2016, doi: 10.5749/jcritethnstud.2.1.0073.

[18]      M. Gius, “Relationship Between Red Flag Laws and Mass Shootings,” Atl. Econ. J., vol. 52, no. 1, pp. 31–38, Mar. 2024, doi: 10.1007/s11293-024-09791-2.

[19]      R. Wamser-Nanney, “Understanding gun violence: Factors associated with beliefs regarding guns, gun policies, and gun violence.,” Psychol. Violence, vol. 11, no. 4, pp. 349–353, Jul. 2021, doi: 10.1037/vio0000392.