library(igraph)
setwd('~/Documents/Violence Cascades/')
load('Raw Data/lcc.RData')

dim(lcc_verts)[1]
length(vic_ids)
dim(lcc_verts)[1] - length(vic_ids)

birthyears = as.numeric(format(as.Date(lcc_verts$dob),'%Y'))
mean(birthyears)
mean(birthyears[vic_ids])
mean(birthyears[-vic_ids])

mean(lcc_verts$sex=="M")*100
mean(lcc_verts$sex[vic_ids]=="M")*100
mean(lcc_verts$sex[-vic_ids]=="M")*100

mean(lcc_verts$race=="BLK")*100
mean(lcc_verts$race[vic_ids]=="BLK")*100
mean(lcc_verts$race[-vic_ids]=="BLK")*100

mean(lcc_verts$race %in% c("WHI", "WWH"))*100
mean(lcc_verts$race[vic_ids] %in% c("WHI", "WWH"))*100
mean(lcc_verts$race[-vic_ids] %in% c("WHI", "WWH"))*100

mean(lcc_verts$gang.member)*100
mean(lcc_verts$gang.member[vic_ids])*100
mean(lcc_verts$gang.member[-vic_ids])*100

d = degree(lcc)
mean(d)
mean(d[vic_ids])
mean(d[-vic_ids])

is_vic = lcc_verts$vic
vic_nbrs = function(neighbors){ mean(is_vic[neighbors[-1]]) }
nbrs = neighborhood(lcc, order=2)
vic_frac = lapply(nbrs, vic_nbrs)
vic_frac = unlist(vic_frac)

mean(vic_frac)*100
mean(vic_frac[vic_ids])*100
mean(vic_frac[-vic_ids])*100