summaryrefslogtreecommitdiffstats
path: root/R Scripts/demographics.R
blob: 196d485c5259cdbdaf3b4bda7ca687e839a4ef70 (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
32
33
34
35
36
37
38
39
40
41
42
43
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=1)
vic_frac = lapply(nbrs, vic_nbrs)
vic_frac = unlist(vic_frac)

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