library(readr)
library(lubridate)
library(dplyr)
nypd <- read_csv("~/ShinyApps/madlan/nypd.csv")
nypd$year<-year(mdy(nypd$CMPLNT_FR_DT))
nypd<-nypd%>%filter(year>2005)
crime<-subset(nypd, !is.na(nypd$Latitude))
crime<-subset(crime, !is.na(crime$Longitude))
save(crime, file="crime.rda")
library(lubridate)
library(ggplot2)
library(dplyr)
load("crime.rda")
crime_table<-crime%>%group_by(OFNS_DESC)%>%summarize(n=sum(!is.na(year)))%>%arrange(desc(n))
crime_table
boro_table<-crime%>%group_by(BORO_NM)%>%summarize(n=sum(!is.na(year)))%>%arrange(desc(n))
#add population size
boro_table
crime%>%filter(LAW_CAT_CD=="FELONY")%>%group_by(BORO_NM)%>%summarize(n=sum(!is.na(year)))%>%arrange(desc(n))

man_fal%>%group_by(OFNS_DESC)%>%summarize(n=sum(!is.na(year)))%>%arrange(desc(n))
LS0tCnRpdGxlOiAiUiBOb3RlYm9vayIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQoKYGBge3IsIGV2YWw9RkFMU0UsIGluY2x1ZGU9VFJVRX0KbGlicmFyeShyZWFkcikKbGlicmFyeShsdWJyaWRhdGUpCmxpYnJhcnkoZHBseXIpCm55cGQgPC0gcmVhZF9jc3YoIn4vU2hpbnlBcHBzL21hZGxhbi9ueXBkLmNzdiIpCm55cGQkeWVhcjwteWVhcihtZHkobnlwZCRDTVBMTlRfRlJfRFQpKQpueXBkPC1ueXBkJT4lZmlsdGVyKHllYXI+MjAwNSkKY3JpbWU8LXN1YnNldChueXBkLCAhaXMubmEobnlwZCRMYXRpdHVkZSkpCmNyaW1lPC1zdWJzZXQoY3JpbWUsICFpcy5uYShjcmltZSRMb25naXR1ZGUpKQpzYXZlKGNyaW1lLCBmaWxlPSJjcmltZS5yZGEiKQpgYGAKCmBgYHtyLCBlY2hvPVRSVUUsIGluY2x1ZGU9VFJVRX0KbGlicmFyeShsdWJyaWRhdGUpCmxpYnJhcnkoZ2dwbG90MikKbGlicmFyeShkcGx5cikKCmxvYWQoImNyaW1lLnJkYSIpCgpjcmltZV90YWJsZTwtY3JpbWUlPiVncm91cF9ieShPRk5TX0RFU0MpJT4lc3VtbWFyaXplKG49c3VtKCFpcy5uYSh5ZWFyKSkpJT4lYXJyYW5nZShkZXNjKG4pKQpjcmltZV90YWJsZQoKYm9yb190YWJsZTwtY3JpbWUlPiVncm91cF9ieShCT1JPX05NKSU+JXN1bW1hcml6ZShuPXN1bSghaXMubmEoeWVhcikpKSU+JWFycmFuZ2UoZGVzYyhuKSkKI2FkZCBwb3B1bGF0aW9uIHNpemUKYm9yb190YWJsZQoKY3JpbWUlPiVmaWx0ZXIoTEFXX0NBVF9DRD09IkZFTE9OWSIpJT4lZ3JvdXBfYnkoQk9ST19OTSklPiVzdW1tYXJpemUobj1zdW0oIWlzLm5hKHllYXIpKSklPiVhcnJhbmdlKGRlc2MobikpCmBgYAoKCmBgYHtyLCBlY2hvPVRSVUV9Cm55Y19iYXNlID0gcW1hcCgiTmV3IFlvcmsgQ2l0eSIsIHpvb20gPSAxMSwgCiAgICBzb3VyY2U9InN0YW1lbiIsIG1hcHR5cGU9InRvbmVyIixkYXJrZW4gPSBjKC4zLCIjQkJCQkJCIikpCgptYW5fZmFsPC1jcmltZSU+JWZpbHRlcihCT1JPX05NPT0iTUFOSEFUVEFOIiwgTEFXX0NBVF9DRD09IkZFTE9OWSIpCgpueXBkX21hcDwtbnljX2Jhc2UgKyBnZW9tX3BvaW50KGRhdGE9bWFuX2ZhbCwgYWVzKHg9TG9uZ2l0dWRlLCB5PUxhdGl0dWRlKSwgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGNvbG9yPSJkYXJrIGdyZWVuIiwgYWxwaGE9LjE1LCBzaXplPTIpCm55cGRfbWFwCiMgbnlwZF9kZW5zaXR5X21hcApgYGAKCgpgYGB7ciwgZWNobz1UUlVFfQptYW5fZmFsJT4lZ3JvdXBfYnkoT0ZOU19ERVNDKSU+JXN1bW1hcml6ZShuPXN1bSghaXMubmEoeWVhcikpKSU+JWFycmFuZ2UoZGVzYyhuKSkKCmBgYAoK