Los Angeles, California 😎. The City of Angels. Tinseltown. The Entertainment Capital of the World!
Known for its warm weather, palm trees, sprawling coastline, and Hollywood, along with producing some of the most iconic films and songs. However, as with any highly populated city, it isn't always glamorous and there can be a large volume of crime. That's where you can help!
You have been asked to support the Los Angeles Police Department (LAPD) by analyzing crime data to identify patterns in criminal behavior. They plan to use your insights to allocate resources effectively to tackle various crimes in different areas.
The Data
They have provided you with a single dataset to use. A summary and preview are provided below.
It is a modified version of the original data, which is publicly available from Los Angeles Open Data.
crimes.csv
| Column | Description |
|---|---|
'DR_NO' | Division of Records Number: Official file number made up of a 2-digit year, area ID, and 5 digits. |
'Date Rptd' | Date reported - MM/DD/YYYY. |
'DATE OCC' | Date of occurrence - MM/DD/YYYY. |
'TIME OCC' | In 24-hour military time. |
'AREA NAME' | The 21 Geographic Areas or Patrol Divisions are also given a name designation that references a landmark or the surrounding community that it is responsible for. For example, the 77th Street Division is located at the intersection of South Broadway and 77th Street, serving neighborhoods in South Los Angeles. |
'Crm Cd Desc' | Indicates the crime committed. |
'Vict Age' | Victim's age in years. |
'Vict Sex' | Victim's sex: F: Female, M: Male, X: Unknown. |
'Vict Descent' | Victim's descent:
|
'Weapon Desc' | Description of the weapon used (if applicable). |
'Status Desc' | Crime status. |
'LOCATION' | Street address of the crime. |
# Import required libraries
import pandas as pd
import numpy as np
import matplotlib.pyplot as plt
import seaborn as sns
crimes = pd.read_csv("crimes.csv", parse_dates=["Date Rptd", "DATE OCC"], dtype={"TIME OCC": str})
crimes.head()crimes['TIME OCC'] = pd.to_datetime(crimes['TIME OCC'], format='%H%M', errors='coerce')
crimes['hour'] = crimes['TIME OCC'].dt.hour
peak_crime_hour1 = crimes['hour'].value_counts().idxmax()
peak_crime_hour = int(peak_crime_hour1)
# Plot a histogram showing crime occurrences over time
plt.figure(figsize=(10, 6))
plt.hist(crimes['hour'], bins=24, color='skyblue', edgecolor='black')
plt.title('Crime Occurrences by Hour')
plt.xlabel('Hour of the Day')
plt.ylabel('Frequency')
plt.xticks(np.arange(0, 24))
plt.grid(axis='y')
plt.show()
# Print the hour with the highest frequency of crimes
print(f"The hour with the highest frequency of crimes is: {peak_crime_hour1}")crimes['hour'] = crimes['TIME OCC'].dt.hour
crimes['night_crime'] = ((crimes['hour'] >= 22) | (crimes['hour'] <= 3))
# Filter for night crimes and find the area with the largest frequency
night_crimes = crimes[crimes['night_crime']]
peak_night_crime_location = night_crimes['AREA NAME'].value_counts().idxmax()
# Print the area with the largest frequency of night crimes
print(f"The area with the largest frequency of night crimes is: {peak_night_crime_location}")
# Create a pivot table with 'AREA NAME' and 'hour' columns
heatmap_data = night_crimes.pivot_table(index='AREA NAME', columns='hour', aggfunc='size', fill_value=0)
# Define the order of hours
hour_order = list(range(22, 24)) + list(range(0, 4))
# Create a heatmap with specified order
plt.figure(figsize=(15, 8))
sns.heatmap(heatmap_data, cmap='YlOrRd', cbar_kws={'label': 'Number of Night Crimes'}, xticklabels=hour_order)
plt.title('Night Crimes Heatmap by Area and Hour')
plt.show()age_bins = [0, 17, 25, 34, 44, 54, 64, np.inf] # Define the age group bins
age_labels = ['0-17', '18-25', '26-34', '35-44', '45-54', '55-64', '65+'] # Labels for age groups
# Categorize victim ages into age groups
crimes['Victim Age Group'] = pd.cut(crimes['Vict Age'], bins=age_bins, labels=age_labels)
# Count the number of crimes per age group
victim_ages = crimes['Victim Age Group'].value_counts()
# Plotting the number of crimes per age group
plt.figure(figsize=(8, 6))
victim_ages.sort_index().plot(kind='bar', color='skyblue')
plt.title('Number of Crimes Committed against Victims by Age Group')
plt.xlabel('Age Group')
plt.ylabel('Number of Crimes')
plt.xticks(rotation=45)
plt.tight_layout()
plt.show()
# Display the series of number of crimes per age group
print("Number of crimes committed against victims by age group:")
print(victim_ages)