NYPD Shooting Incident Data

data$OCCUR_DATE[1]

data$DATE <- as.Date(data$OCCUR_DATE, format="%m/%d/%Y")
str(data$DATE)
data$DATE[1]
assert_that(noNA(data$DATE))

str(data$OCCUR_TIME)
assert_that(noNA(data$OCCUR_TIME))

data$Borough <- as.factor(data$BORO)
assert_that(is.factor(data$Borough))
assert_that(noNA(data$Borough))

big_mark <- function(vector, columns){
    kable(table(vector), col.names=columns, format.args=list(big.mark=","))
}
big_mark(data$Borough, c("Borough", "Shootings"))

n_distinct(data$PRECINCT)
str(data$PRECINCT)
data$PRECINCT_FACTOR <- as.factor(data$PRECINCT)
see_if(noNA(data$PRECINCT))

data$JURISDICTION <- as.factor(data$JURISDICTION_CODE)
see_if(noNA(data$JURISDICTION))

big_mark(data$JURISDICTION, c("Jurisdiction", "Shootings"))

data$LOCATION_DESC[1]
see_if(noNA(data$LOCATION_DESC))

n_distinct(data$LOCATION_DESC)

first(na.omit(data$LOCATION_DESC))

big_mark(data$LOCATION_DESC, c("Location", "Shootings"))

see_if(noNA(data$STATISTICAL_MURDER_FLAG))

big_mark(data$STATISTICAL_MURDER_FLAG, c("Victim Died", "Count"))

see_if(noNA(data$PERP_AGE_GROUP))
data$PERP_AGE <- as.factor(data$PERP_AGE_GROUP)

big_mark(data$PERP_AGE, c("Perpetrator's Age Group", "Shootings"))

see_if(noNA(data$PERP_SEX))
data$PERP_SEX <- as.factor(data$PERP_SEX)

big_mark(data$PERP_SEX, c("Perpetrator's Sex", "Shootings"))

see_if(noNA(data$PERP_RACE))
data$PERP_RACE <- as.factor(data$PERP_RACE)

big_mark(data$PERP_RACE, c("Perpetrator's Race", "Shootings"))

see_if(noNA(data$VIC_AGE_GROUP))
data$VIC_AGE_GROUP <- as.factor(data$VIC_AGE_GROUP)

big_mark(data$VIC_AGE_GROUP, c("Victim's Age Group", "Shootings"))

see_if(noNA(data$VIC_SEX))
data$VIC_SEX <- as.factor(data$VIC_SEX)

big_mark(data$VIC_SEX, c("Victims' Sex", "Shootings"))

see_if(noNA(data$VIC_RACE))
data$VIC_RACE <- as.factor(data$VIC_RACE)

big_mark(data$VIC_RACE, c("Victim's Race", "Shootings"))

see_if(noNA(data$X_COORD_CD))
n_distinct(data$X_COORD_CD)

see_if(noNA(data$Y_COORD_CD))
n_distinct(data$Y_COORD_CD)

see_if(noNA(data$Latitude))
n_distinct(data$Latitude)

see_if(noNA(data$Longitude))
n_distinct(data$Longitude)
str(data$Longitude)

see_if(noNA(data$Lon_Lat))
n_distinct(data$Lon_Lat)
str(data$Lon_Lat)

dates <- data
dates$Year <- floor_date(dates$DATE, "year")

annual_shootings <- dates %>% group_by(Year) %>% summarise(Shootings=n())
kable(annual_shootings,
      format.args=list(big.mark=","))

theme_set(theme_minimal())
PATH = "../files/posts/nypd-shooting-incident-data/"
join_path = partial(paste, PATH, sep="")

plot = ggplot(data=annual_shootings, aes(x=Year, y=Shootings)) +
  geom_line() + geom_point() +
  scale_x_date(date_labels="%Y", name="Year", date_breaks="1 year") +
  labs(title="NYPD Annual Shootings Incidents Reported")
ggsave(join_path("annual_shootings.png"),
       plot=plot)

dates$year_month <- floor_date(dates$DATE, "month")
monthly_shootings <- dates %>% group_by(year_month) %>% summarise(Shootings=n())

monthly_shootings$month_number <- as.numeric(
  format(as.Date(monthly_shootings$year_month), "%m"))
monthly_shootings$Month <- months(as.Date(monthly_shootings$year_month),
                                  abbreviate=TRUE)
monthly_shootings$Year <- year(ymd(monthly_shootings$year_month))
monthly_shootings

plot = ggplot(data=monthly_shootings,
              aes(x=month_number,
                  y=Shootings,
                  group=Year, colour=Year)) +
  geom_line() + geom_point() +
  scale_x_continuous(breaks=monthly_shootings$month_number,
                     labels=monthly_shootings$Month, name="Month") +
  labs(title="NYPD Shootings by Month (Year-Over-Year)")

ggsave("../files/posts/nypd-shooting-incident-data/monthly_shootings.png",
       plot=plot)

three_years = monthly_shootings[which(monthly_shootings$Year %in% c(2006, 2018, 2020)),]

plot = ggplot(data=three_years,
              aes(x=month_number,
                  y=Shootings,
                  group=Year, colour=Year)) +
  geom_line() + geom_point() +
  scale_x_continuous(breaks=three_years$month_number,
                     labels=three_years$Month, name="Month") +
  labs(title="NYPD Shootings by Month (2006, 2018, 2020)")

ggsave("../files/posts/nypd-shooting-incident-data/monthly_shootings_three-years.png",
       plot=plot)

dates$Year <- year(ymd(dates$year_month))
last_year <- dates[dates$Year==2020,]
last_year

minimum_latitude = min(last_year$Latitude)
maximum_latitude = max(last_year$Latitude)
minimum_longitude = min(last_year$Longitude)
maximum_longitude = max(last_year$Longitude)


plot = ggplot(data=last_year,
              aes(x=Longitude,
                  y=Latitude, color=BORO)) +
  geom_point(size=1) +
  scale_x_continuous(limits=c(minimum_longitude, maximum_longitude)) +
  scale_y_continuous(limits=c(minimum_latitude, maximum_latitude)) +
  labs(title="NYPD Shootings by Latitude and Longitude (2020)")

ggsave(join_path("monthly_shootings_latitude_longitude.png"),
       plot=plot)