NYC Data Science Academy| Blog
Bootcamps
Lifetime Job Support Available Financing Available
Bootcamps
Data Science with Machine Learning Flagship ๐Ÿ† Data Analytics Bootcamp Artificial Intelligence Bootcamp New Release ๐ŸŽ‰
Free Lesson
Intro to Data Science New Release ๐ŸŽ‰
Find Inspiration
Find Alumni with Similar Background
Job Outlook
Occupational Outlook Graduate Outcomes Must See ๐Ÿ”ฅ
Alumni
Success Stories Testimonials Alumni Directory Alumni Exclusive Study Program
Courses
View Bundled Courses
Financing Available
Bootcamp Prep Popular ๐Ÿ”ฅ Data Science Mastery Data Science Launchpad with Python View AI Courses Generative AI for Everyone New ๐ŸŽ‰ Generative AI for Finance New ๐ŸŽ‰ Generative AI for Marketing New ๐ŸŽ‰
Bundle Up
Learn More and Save More
Combination of data science courses.
View Data Science Courses
Beginner
Introductory Python
Intermediate
Data Science Python: Data Analysis and Visualization Popular ๐Ÿ”ฅ Data Science R: Data Analysis and Visualization
Advanced
Data Science Python: Machine Learning Popular ๐Ÿ”ฅ Data Science R: Machine Learning Designing and Implementing Production MLOps New ๐ŸŽ‰ Natural Language Processing for Production (NLP) New ๐ŸŽ‰
Find Inspiration
Get Course Recommendation Must Try ๐Ÿ’Ž An Ultimate Guide to Become a Data Scientist
For Companies
For Companies
Corporate Offerings Hiring Partners Candidate Portfolio Hire Our Graduates
Students Work
Students Work
All Posts Capstone Data Visualization Machine Learning Python Projects R Projects
Tutorials
About
About
About Us Accreditation Contact Us Join Us FAQ Webinars Subscription An Ultimate Guide to
Become a Data Scientist
    Login
NYC Data Science Acedemy
Bootcamps
Courses
Students Work
About
Bootcamps
Bootcamps
Data Science with Machine Learning Flagship
Data Analytics Bootcamp
Artificial Intelligence Bootcamp New Release ๐ŸŽ‰
Free Lessons
Intro to Data Science New Release ๐ŸŽ‰
Find Inspiration
Find Alumni with Similar Background
Job Outlook
Occupational Outlook
Graduate Outcomes Must See ๐Ÿ”ฅ
Alumni
Success Stories
Testimonials
Alumni Directory
Alumni Exclusive Study Program
Courses
Bundles
financing available
View All Bundles
Bootcamp Prep
Data Science Mastery
Data Science Launchpad with Python NEW!
View AI Courses
Generative AI for Everyone
Generative AI for Finance
Generative AI for Marketing
View Data Science Courses
View All Professional Development Courses
Beginner
Introductory Python
Intermediate
Python: Data Analysis and Visualization
R: Data Analysis and Visualization
Advanced
Python: Machine Learning
R: Machine Learning
Designing and Implementing Production MLOps
Natural Language Processing for Production (NLP)
For Companies
Corporate Offerings
Hiring Partners
Candidate Portfolio
Hire Our Graduates
Students Work
All Posts
Capstone
Data Visualization
Machine Learning
Python Projects
R Projects
About
Accreditation
About Us
Contact Us
Join Us
FAQ
Webinars
Subscription
An Ultimate Guide to Become a Data Scientist
Tutorials
Data Analytics
  • Learn Pandas
  • Learn NumPy
  • Learn SciPy
  • Learn Matplotlib
Machine Learning
  • Boosting
  • Random Forest
  • Linear Regression
  • Decision Tree
  • PCA
Interview by Companies
  • JPMC
  • Google
  • Facebook
Artificial Intelligence
  • Learn Generative AI
  • Learn ChatGPT-3.5
  • Learn ChatGPT-4
  • Learn Google Bard
Coding
  • Learn Python
  • Learn SQL
  • Learn MySQL
  • Learn NoSQL
  • Learn PySpark
  • Learn PyTorch
Interview Questions
  • Python Hard
  • R Easy
  • R Hard
  • SQL Easy
  • SQL Hard
  • Python Easy
Data Science Blog > R > Data Study on The Global Crime Landscape

Data Study on The Global Crime Landscape

John Montroy
Posted on Nov 19, 2015
The skills the author demoed here can be learned through taking Data Science with Machine Learning bootcamp with NYC Data Science Academy.
Contributed by John Montroy. John took NYC Data Science Academy 12 week full time Data Science Bootcamp program between Sept 23 to Dec 18, 2015. The post was based on his second class project(due at 4th week of the program).
Try out the app!   ||   View the code   ||   Download the data

Given today's constant barrage of negative news media, harboring concerns about public safety is understandable. For many, data shows the specter of crime hangs much heavier than the reality of it - the world is, on average, much safer today than in the past. This is, however, certainly not always the case across the globe. Things many of us might take for granted, such as a low likelihood of a burglary, or a responsive police force, are near non-existent in other countries, some of them similarly industrialized.

Data Study on The Global Crime Landscape

Artwork by Gary Halvorson, Oregon State Archives

The International Crime Victimization Survey (ICVS) is a quadrennial survey performed across approximately 90 countries world-wide. Every four years, it takes the temperature of countries on issues such as crime, safety, and policing, serving as a standardized lense of comparison for these statistics. Although the dataset does not yet provide full country coverage, it offers a more well-rounded story on global crime and policing.

To aid in understanding the ICVS, this app helps the user to interactively explore the results of all surveys conducted over the past 25 years. The app was developed using R, Shiny, and ggplot2, and offers visualizations of world results, a country comparison tool, and a look at question correlations using GoogleViz.

Disastrous Data

In terms of complexity, the dataset here is massive - a modestly large number of records, but a huge number of variables. My starting point in a nutshell:

  1. Two main datasets (1989 - 2005, 2010) with 700+ variables each and 350,000+ records in total.
  2. One codebook, a 56-page PDF.
  3. Multiple questionnaires, spanning all iterations of the survey.

The first steps, then, are clear: limit scope, clean the data, translate the data to useful categories, and aggregate scores based on question and country. All of this data wrangling was performed in R. The full wrangling code (+ all other code) is available on GitHub; I discuss some highlights here.

Aggregation for (Actual) Days

The basic workflow for aggregation was as follows:

Data Study on The Global Crime Landscape

Each question varies by year and survey.

  1. Identify response labels using the provided codebook ("Yes", "No", "Unknown")
  2. Modify aggregation function to handle new labels
  3. Parse labels to numbers ("Yes" -> 1, "No" -> 0)
  4. Generate aggregate score per country per question, based on response values
  5. Bind aggregated scores together into master dataframe

Curveballs abounded. Some questions had four possible respondent answers, some three, some two. A handful of questions actually changed the number of possible answers over time. Ambiguous language was abundant - "Do Not Know" vs "Unknown"? Variation also existed with survey years per country; the survey was conducted mostly on a quadrennial basis, but some were conducted off-year. Three countries might've had a question asked in 1995, 40 in 1996, and two more in 1997. And lastly, per above, two datasets were provided (1989 - 2005, 2010); these datasets occasionally contained slight variations in structure.

A number of strong assumptions had to be made to make aggregation feasible and productive.

  • Six main survey sweeps were conducted every four years, starting from 1988. Surveys conducted on off-years are associated with one of these six sweeps.
  • Respondent answers without a strong sentiment attached ("Don't know", "Can't remember", "Unknown") are considered to contain no meaningful information, and are discarded as missingness.
  • Variation in the number of possible respondent answers was ignored - only the intersection of answer sets was considered for aggregation.

In short, some information loss occurred at this step, but gains in clarity far outweighed the minimal losses in specificity.

Functions were generated to handle this aggregation process on the R side; the below function dynamically generates country scores by year, given an input column to group on. The function itself was tweaked on-the-fly to handle fundamental differences (like response labels) - the 2.0 version of this function will lean heavily towards the functional programming side (and would preferably be much cleaner...):

library(lazyeval)

Clean.Reporting <- function(df, df2, scorecol) {
 df.rest <- df %>% 
  select_(~Year, ~Country, scorecat = interp(~x, x = as.name(scorecol))) %>% 
  group_by(Year, Country, scorecat) %>% 
  summarise(count = n())
 
 df.2010 <- df2 %>%
  select_(~Year, ~Country, scorecat = interp(~x, x = as.name(scorecol))) %>% 
  group_by(Year, Country, scorecat) %>% 
  summarise(count = n())
 
 df.net <- data.frame(rbind(df.rest, df.2010))
 df.net <- df.net[!is.na(df.net$Country),]
 df.net$scorecat <- as.character(df.net$scorecat)
 
 df.net$scorecat[df.net$scorecat == "yes"] <- 1
 df.net$scorecat[df.net$scorecat == "no"] <- 0
 df.net$scorecat[!(df.net$scorecat %in% 0:1)] <- NA
 
 df.final <- data.frame(GetPercents(df.net[!is.na(df.net$scorecat),], "scorecat"))
 
 df.final <- df.final %>% 
  filter(scorecat == 1) %>%
  rename(Score = Perc) %>%
  mutate(Question = qtrans$Name[qtrans$Id == scorecol]) %>%
  select(Year, Country, Question, Score)
 
 return(df.final)
}

This function is called by the above, and generates percentages per category for the aggregations.

GetPercents <- function(df, perccol) {
 df.net <- df %>%
 group_by(Country, Year) %>%
 summarise(count.net = sum(count))
 
 df.group <- df %>%
 group_by_(~Country, ~Year, interp(~x, x = as.name(perccol))) %>%
 summarise(count.d = sum(count))
 
 df.final <- inner_join(df.group,df.net, by = c("Country","Year"))
 
 df.final <- df.final %>%
 mutate(Perc = count.d/count.net) %>%
 select_(~Country,~Year,interp(~x, x = as.name(perccol)),~Perc)
 
 return(df.final)
}

I have great admiration for whoever originally compiled this dataset - the fact that it's as clean, standardized, and well-documented as it is is a marker of obvious labor and scrupulousness. That said, generating aggregated statistics on top of survey-level data was difficult and time-consuming. My personal toolkit of R functions for such tasks has greatly expanded as a result, but this was an ambitious dataset to work with.

Data on Shiny Adventures

Data Study on The Global Crime Landscape

The app allows for country comparisons per question.

Shiny is R's interactive UI module - it allows for easy creation of flexible data applications, and supports almost any R functionality you might need. Of course, Shiny isn't without its quirks; the core functionality is a bit unintuitive, and the syntax is mercilessly specific. You can read more about Shiny here, and view a gallery of applications here.

This application makes use of the Shiny Dashboard layout, and takes advantage of reactive programming styles to generate a handful of ggplot2 visualizations. Maps, question text, and the list of selectable countries are all generated dynamically based on other inputs.

Google's Motion Chart allows for an easy peek at question correlations.

Google's Motion Chart allows for an easy peek at question correlations.

The application itself is broken into four tabs: a world map result visualization, a tool for comparing results between countries, a GoogleViz Motion Chart, and an About tab. The world map visualizer gives a broad overview perspective of results by question and year; the country compare tab allows the user to drill down on particular countries of interest. Google's Motion Chart is good at examining correlations between particular questions over time - for a great demonstration, see the popular TED talk given by Hans Rosling.

Code That Worked

The full code can be found on GitHub - a few choice bits are reviewed here.

Error handling between reactive calls was configurable through a combo of Shiny's validate + need functions, seen here:

validate(
   need(dim(ICVS.country.df1())[1] != 0, "Loading data set...")
 )

This calls validate + need on a reactive function - every time the reactive is called, this validate condition needs to be fulfilled. This is only a problem due to wait times between loading data sets.

The app needs a way to dynamically render a number of faceted country plots based on the user input of how many countries to compare. Moreover, these plots need to be rendered with only one legend, rather than 1 - 3 separate ones. This was achieved with the below code, which creates a list of plots generated by CountryMapViz based on how many country comparisons were requested:

glist <- list()
for(x in 1:input$ncountries) {
 
  ICVS.country.looper <- paste0("ICVS.country.df",x,"()")
  ICVS.map.looper <- paste0("ICVS.country.map",x,"()")
 
  glist[[length(glist)+1]] <- CountryMapViz(
   eval(parse(text=ICVS.country.looper)), 
   eval(parse(text=ICVS.map.looper)), 
   qinfo.country()
  ) 
} 
 
grid_arrange_shared_legend(glist)

The grid_arrange_shared_legend functions comes courtesy of the great Hadley Wickham; it extracts the legend (via "guide-box") from one plot of the several passed, and grid.arranges the plots with the legend.

grid_arrange_shared_legend <- function(plots) {
 g <- ggplotGrob(plots[[1]] + theme(legend.position="top"))$grobs
 legend <- g[[which(sapply(g, function(x) x$name) == "guide-box")]]
 lheight <- sum(legend$height)
 grid.arrange(
 do.call(arrangeGrob, lapply(plots, function(x)
 x + theme(legend.position="none"))),
 legend,
 ncol = 1,
 heights = unit.c(unit(1, "npc") - lheight, lheight))
}

One of the last big hang-ups was getting the "softer" data into the app - question texts, question names, legend scales, and so on. This came down to the unfortunate task of diving into the survey questionnaires, finding the question text, and piecing that together into an R data frame. Note the use of the I() function to create a column of lists - subsetting this data frame allows for the extraction of legend scales and labels for each question.

C01B400.Info <- data.frame(Id = "C01B400", Title = "Car Theft", 
 FullText = "REPORT TO POLICE - Over the past five years have you or other members of your household had any of their cars/vans/trucks stolen? If so, did you or anyone else report the incident to the police?", 
 Min = 0, Max = 1, marks = I(list(c(0,0.5,1))), 
 labels = I(list(c("0% Reported","50% Reported", "100% Reported"))))

Thoughts

And so what? What new perspectives on crime and safety can this app provide us?

Sexual offenses are woefully under-reported.

Sexual offenses are woefully under-reported.

The world map paints a not-unexpected picture: despite sabre-rattling, Westernized countries in Europe and beyond have, on average, a positive outlook on public safety. Countries in eastern Europe or South America typically don't fare as well - all three sentiment categories (police approval, safety at night, chance of break-in) are similarly negative. This correlation in the sentiment categories is quickly confirmed by the GViz Correlations tab - any combination of those three variables seem to be roughly linearly correlated with one another.

With police reporting rates, however, the story is more complex. Some reporting rates vary across the globe similar to the sentiment questions, but other rates are jarringly static. Percentages of sexual offenses reported, for example, are abysmal world-wide. Equally low are percentages of consumer fraud reported, but these low rates are likely derived from separate causes. Consumer fraud may be more typically handled in the private sphere; sexual offenses are complex and difficult to pin down, as exemplified by recent controversies in the United States.

Car theft is likely the best-reported of any crime.

Car theft is likely the best-reported of any crime.

Conclusion

What is clear is that high-cost crimes are exceptionally well-reported. Car theft is likely the most frequently reported crime across the world; burglary and robbery fair well too. Lesser-impact crimes like bicycle theft and car vandalism make up the other side of spectrum, likely due to lesser cost and increased difficulty of a resolution.

I encourage you to try out the app and find your own patterns. The app allows for compelling visualizations at the click of a button, and there is surely more to find. Perhaps in the future, this app can be expanded to handle other datasets, based on interest.

About Author

John Montroy

John Montroy is a graduate of Middlebury College with a B.A. in Physics. After a summer of particle physics at CERN with the Harvard ATLAS team, he began his career as a data analyst in the auto industry....
View all posts by John Montroy >

Leave a Comment

Cancel reply

You must be logged in to post a comment.

No comments found.

View Posts by Categories

All Posts 2399 posts
AI 7 posts
AI Agent 2 posts
AI-based hotel recommendation 1 posts
AIForGood 1 posts
Alumni 60 posts
Animated Maps 1 posts
APIs 41 posts
Artificial Intelligence 2 posts
Artificial Intelligence 2 posts
AWS 13 posts
Banking 1 posts
Big Data 50 posts
Branch Analysis 1 posts
Capstone 206 posts
Career Education 7 posts
CLIP 1 posts
Community 72 posts
Congestion Zone 1 posts
Content Recommendation 1 posts
Cosine SImilarity 1 posts
Data Analysis 5 posts
Data Engineering 1 posts
Data Engineering 3 posts
Data Science 7 posts
Data Science News and Sharing 73 posts
Data Visualization 324 posts
Events 5 posts
Featured 37 posts
Function calling 1 posts
FutureTech 1 posts
Generative AI 5 posts
Hadoop 13 posts
Image Classification 1 posts
Innovation 2 posts
Kmeans Cluster 1 posts
LLM 6 posts
Machine Learning 364 posts
Marketing 1 posts
Meetup 144 posts
MLOPs 1 posts
Model Deployment 1 posts
Nagamas69 1 posts
NLP 1 posts
OpenAI 5 posts
OpenNYC Data 1 posts
pySpark 1 posts
Python 16 posts
Python 458 posts
Python data analysis 4 posts
Python Shiny 2 posts
R 404 posts
R Data Analysis 1 posts
R Shiny 560 posts
R Visualization 445 posts
RAG 1 posts
RoBERTa 1 posts
semantic rearch 2 posts
Spark 17 posts
SQL 1 posts
Streamlit 2 posts
Student Works 1687 posts
Tableau 12 posts
TensorFlow 3 posts
Traffic 1 posts
User Preference Modeling 1 posts
Vector database 2 posts
Web Scraping 483 posts
wukong138 1 posts

Our Recent Popular Posts

AI 4 AI: ChatGPT Unifies My Blog Posts
by Vinod Chugani
Dec 18, 2022
Meet Your Machine Learning Mentors: Kyle Gallatin
by Vivian Zhang
Nov 4, 2020
NICU Admissions and CCHD: Predicting Based on Data Analysis
by Paul Lee, Aron Berke, Bee Kim, Bettina Meier and Ira Villar
Jan 7, 2020

View Posts by Tags

#python #trainwithnycdsa 2019 2020 Revenue 3-points agriculture air quality airbnb airline alcohol Alex Baransky algorithm alumni Alumni Interview Alumni Reviews Alumni Spotlight alumni story Alumnus ames dataset ames housing dataset apartment rent API Application artist aws bank loans beautiful soup Best Bootcamp Best Data Science 2019 Best Data Science Bootcamp Best Data Science Bootcamp 2020 Best Ranked Big Data Book Launch Book-Signing bootcamp Bootcamp Alumni Bootcamp Prep boston safety Bundles cake recipe California Cancer Research capstone car price Career Career Day ChatGPT citibike classic cars classpass clustering Coding Course Demo Course Report covid 19 credit credit card crime frequency crops D3.js data data analysis Data Analyst data analytics data for tripadvisor reviews data science Data Science Academy Data Science Bootcamp Data science jobs Data Science Reviews Data Scientist Data Scientist Jobs data visualization database Deep Learning Demo Day Discount disney dplyr drug data e-commerce economy employee employee burnout employer networking environment feature engineering Finance Financial Data Science fitness studio Flask flight delay football gbm Get Hired ggplot2 googleVis H20 Hadoop hallmark holiday movie happiness healthcare frauds higgs boson Hiring hiring partner events Hiring Partners hotels housing housing data housing predictions housing price hy-vee Income industry Industry Experts Injuries Instructor Blog Instructor Interview insurance italki Job Job Placement Jobs Jon Krohn JP Morgan Chase Kaggle Kickstarter las vegas airport lasso regression Lead Data Scienctist Lead Data Scientist leaflet league linear regression Logistic Regression machine learning Maps market matplotlib Medical Research Meet the team meetup methal health miami beach movie music Napoli NBA netflix Networking neural network Neural networks New Courses NHL nlp NYC NYC Data Science nyc data science academy NYC Open Data nyc property NYCDSA NYCDSA Alumni Online Online Bootcamp Online Training Open Data painter pandas Part-time performance phoenix pollutants Portfolio Development precision measurement prediction Prework Programming public safety PwC python Python Data Analysis python machine learning python scrapy python web scraping python webscraping Python Workshop R R Data Analysis R language R Programming R Shiny r studio R Visualization R Workshop R-bloggers random forest Ranking recommendation recommendation system regression Remote remote data science bootcamp Scrapy scrapy visualization seaborn seafood type Selenium sentiment analysis sentiment classification Shiny Shiny Dashboard Spark Special Special Summer Sports statistics streaming Student Interview Student Showcase SVM Switchup Tableau teachers team team performance TensorFlow Testimonial tf-idf Top Data Science Bootcamp Top manufacturing companies Transfers tweets twitter videos visualization wallstreet wallstreetbets web scraping Weekend Course What to expect whiskey whiskeyadvocate wildfire word cloud word2vec XGBoost yelp youtube trending ZORI

NYC Data Science Academy

NYC Data Science Academy teaches data science, trains companies and their employees to better profit from data, excels at big data project consulting, and connects trained Data Scientists to our industry.

NYC Data Science Academy is licensed by New York State Education Department.

Get detailed curriculum information about our
amazing bootcamp!

Please enter a valid email address
Sign up completed. Thank you!

Offerings

  • HOME
  • DATA SCIENCE BOOTCAMP
  • ONLINE DATA SCIENCE BOOTCAMP
  • Professional Development Courses
  • CORPORATE OFFERINGS
  • HIRING PARTNERS
  • About

  • About Us
  • Alumni
  • Blog
  • FAQ
  • Contact Us
  • Refund Policy
  • Join Us
  • SOCIAL MEDIA

    ยฉ 2025 NYC Data Science Academy
    All rights reserved. | Site Map
    Privacy Policy | Terms of Service
    Bootcamp Application