blob: 30d458c64892446f7a8bb54bbc0106b1a968800b (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
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')
|
