demographic stats

1 files changed, 30 insertions, 11 deletions

diff --git a/R Scripts/demographics.R b/R Scripts/demographics.R
index 349eda9..196d485 100644
--- a/R Scripts/demographics.R
+++ b/R Scripts/demographics.R
@@ -2,23 +2,42 @@ 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)
 
-mean(lcc_verts$sex=="M")
-mean(lcc_verts$sex[vic_ids]=="M")
+birthyears = as.numeric(format(as.Date(lcc_verts$dob),'%Y'))
+mean(birthyears)
+mean(birthyears[vic_ids])
+mean(birthyears[-vic_ids])
 
-mean(lcc_verts$race=="BLK")
+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 %in% c("WHI", "WWH"))
+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$gang.member)
+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=1)
+vic_frac = lapply(nbrs, vic_nbrs)
+vic_frac = unlist(vic_frac)
 
-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)
-}
\ No newline at end of file
+mean(vic_frac)
+mean(vic_frac[vic_ids])
+mean(vic_frac[-vic_ids])