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')