diff options
| author | Ben Green <bgreen@g.harvard.edu> | 2015-08-01 12:01:19 -0400 |
|---|---|---|
| committer | Ben Green <bgreen@g.harvard.edu> | 2015-08-01 12:01:19 -0400 |
| commit | 9350ee2c6359562a23cf8efdefdd7de80b2a682e (patch) | |
| tree | b1217fcc2098bebaa6fed641bb759ac71eae139f /R Scripts/demographics.R | |
| parent | 37c23e0d1fcbdee116856d59c0dab98686e265ae (diff) | |
| download | criminal_cascades-9350ee2c6359562a23cf8efdefdd7de80b2a682e.tar.gz | |
minor updates
Diffstat (limited to 'R Scripts/demographics.R')
| -rw-r--r-- | R Scripts/demographics.R | 24 |
1 files changed, 24 insertions, 0 deletions
diff --git a/R Scripts/demographics.R b/R Scripts/demographics.R new file mode 100644 index 0000000..349eda9 --- /dev/null +++ b/R Scripts/demographics.R @@ -0,0 +1,24 @@ +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) +}
\ No newline at end of file |