1 files changed, 43 insertions, 0 deletions
diff --git a/R Scripts/arrest-vs-infection-times.R b/R Scripts/arrest-vs-infection-times.R
new file mode 100644
index 0000000..222be3e
--- /dev/null
+++ b/R Scripts/arrest-vs-infection-times.R
@@ -0,0 +1,43 @@
+library(igraph)
+setwd('~/Documents/Cascade Project/Raw Data/')
+
+##### Load arrest data
+arrests <- read.csv("2006to2014arrests2.csv", header=T, colClass=c("character"))
+arrests$ir2 <- paste("ir", arrests$ir_no)
+sub.arrests <- subset(arrests, select=c(as.character("rd_no"), as.character('ir2'), 
+                                        as.character('arrest_date')))
+colnames(sub.arrests) <- c("events", "individuals", "dates")
+sub.arrests$dates = as.Date(sub.arrests$dates,format='%m/%d/%Y')
+sub.arrests = sub.arrests[order(sub.arrests$dates),]
+# individuals <- unique(sub.arrests$individuals)
+# events <- unique(sub.arrests$events)
+
+##### Look at nonfatal victims
+shootings <- read.csv("shooting-data-withdate2.csv", header = T)
+victims = shootings[shootings$INV_PARTY_TYPE_CD=="VIC",]
+victims = victims[!is.na(victims$IR_NO),]
+victims$ir2 <- paste("ir", victims$IR_NO)
+victims$INCIDENT_DATE = as.Date(victims$INCIDENT_DATE,format='%m/%d/%y')
+victims = victims[order(victims$INCIDENT_DATE),] # sort so match gets first infection
+victims = victims[victims$ir2 %in% arrests$ir2,] # only look at vics who were arrested
+
+vic_dates = victims$INCIDENT_DATE[match(unique(victims$ir2),victims$ir2)]
+arrest_dates = sub.arrests$dates[match(unique(victims$ir2), sub.arrests$individuals)]
+sum(vic_dates==arrest_dates)
+# the number of vics in arrest table doesn't match number of nonfatal vics in person?
+
+##### Look at fatal victims
+murders = read.csv("murder-victims-13nov.csv", header=T)
+murders = murders[!is.na(murders$VICTIM_IR_NO),]
+murders = murders[murders$INJURY_DESCR=="SHOT",]
+murders = murders[match(unique(murders$VICTIM_IR_NO),murders$VICTIM_IR_NO),]
+murders$ir2 = paste("ir", murders$VICTIM_IR_NO)
+start_date = as.Date("2005-12-31")
+murders$INJURY_DATE = as.Date(murders$INJURY_DATE,format='%m/%d/%y')
+murders = murders[murders$INJURY_DATE>=start_date,]
+murders = murders[order(murders$INJURY_DATE),]
+murders = murders[murders$ir2 %in% arrests$ir2,]
+
+mur_dates = murders$INJURY_DATE
+arrest_dates = sub.arrests$dates[match(murders$ir2, sub.arrests$individuals)]
+sum(mur_dates==arrest_dates)