From f23e9e5da58d1d19913703e11f7ff96a7adff96b Mon Sep 17 00:00:00 2001
From: Ben Green <bgreen@g.harvard.edu>
Date: Sun, 13 Sep 2015 23:52:40 -0400
Subject: updated prediction code for new model and data format

---
 R Scripts/prior-arrests-by-day.R | 31 +++++++++++++++++++++++++++++++
 1 file changed, 31 insertions(+)
 create mode 100644 R Scripts/prior-arrests-by-day.R

(limited to 'R Scripts/prior-arrests-by-day.R')

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')
-- 
cgit v1.2.3-70-g09d2