1 files changed, 31 insertions, 0 deletions
diff --git a/R Scripts/prior-arrests-by-day.R b/R Scripts/prior-arrests-by-day.R
new file mode 100644
index 0000000..30d458c
--- /dev/null
+++ b/R Scripts/prior-arrests-by-day.R
@@ -0,0 +1,31 @@
+library(igraph)
+setwd('~/Documents/Violence Cascades/Raw Data/')
+
+load('arrests.RData')
+load('lcc.RData')
+start_date = as.Date("2005-12-31")
+
+arrests = arrests[arrests$ir2 %in% lcc_verts$ir_no,c('ir2','arrest_date')]
+arrests$arrest_date = as.Date(arrests$arrest_date,format='%m/%d/%Y')
+arrests$arrest_day = as.numeric(arrests$arrest_date - start_date)
+arrests$id = match(arrests$ir2, lcc_verts$ir_no)
+
+arrest_days = function(arr,i){return(arr$arrest_day[arr$id==i)}
+
+prior_arrests = sapply(1:vcount(lcc),function(x) NULL)
+for(i in 1:vcount(lcc)){
+  if(i%%10000==0) print(i)
+  prior_arrests[[i]] = sort(arrests$arrest_day[arrests$id==i])
+}
+save(prior_arrests,file='prior-arrests.RData')
+
+
+#### turn vic_times_lcc into a list
+vic_times = sapply(1:vcount(lcc),function(x) NULL)
+for(i in 1:vcount(lcc)){
+  if(i%%10000==0) print(i)
+  days = as.numeric(vic_times_lcc[i,2:7])
+  days = days[!is.na(days)]
+  vic_times[[i]] = days
+}
+save(vic_times,file='vic_times_list.RData')