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

length(vic_ids)

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

mean(lcc_verts$race=="BLK")

mean(lcc_verts$race %in% c("WHI", "WWH"))

mean(lcc_verts$gang.member)

d = degree(lcc)
mean(d)

vic_frac = rep(NA,vcount(lcc))
for(i in 1:vcount(lcc)){
  if (i %% 10000==0) print(i)
  neighbors = unlist(neighborhood(lcc, nodes=i, order=2, mode='all'))[-1]
  vic_frac[i] = sum(lcc_verts$vic[neighbors])/length(neighbors)
}