diff options
| author | Ben Green <bgreen@g.harvard.edu> | 2015-09-09 20:27:12 -0400 |
|---|---|---|
| committer | Ben Green <bgreen@g.harvard.edu> | 2015-09-09 20:27:15 -0400 |
| commit | 717e7f2edda52100406d17826bd617b7915c5daf (patch) | |
| tree | a3501a0132e2d9024d8a0eb5afe0771311404b34 /R Scripts/sim-analysis.R | |
| parent | 26c65bfeeb22a46e6b688ea3765b9a6e0479f748 (diff) | |
| download | criminal_cascades-717e7f2edda52100406d17826bd617b7915c5daf.tar.gz | |
updated R scripts
Diffstat (limited to 'R Scripts/sim-analysis.R')
| -rwxr-xr-x | R Scripts/sim-analysis.R | 11 |
1 files changed, 11 insertions, 0 deletions
diff --git a/R Scripts/sim-analysis.R b/R Scripts/sim-analysis.R index af280cc..92c3d26 100755 --- a/R Scripts/sim-analysis.R +++ b/R Scripts/sim-analysis.R @@ -85,3 +85,14 @@ dates = c(V(lcc)$fatal_date[!is.na(V(lcc)$fatal_date)], dates = as.Date(dates) hist(dates, breaks='months',freq=T,col='lightblue',format='%m/%Y', xlab='Month',ylab='Number of Shootings',main='') + +### data vs sim hist comparison +hist(vic.time.data,breaks=seq(0,3000,by=10),col=rgb(0,0,1,1/2),freq=F, + xlim=c(0,max(vic.time,vic.time.data)), + xlab='Days between infections',main='') +hist(vic.time,breaks=seq(0,3000,by=10),col=rgb(1,0,0,1/2),add=T,freq=F) +legend("topright", inset=0.05, + c("Data", "Simulation"), + fill=c(rgb(0,0,1,1/2),rgb(1,0,0,1/2))) +box(which='plot') + |