diff options
| author | Ben Green <bgreen@g.harvard.edu> | 2015-09-26 16:26:19 -0400 |
|---|---|---|
| committer | Ben Green <bgreen@g.harvard.edu> | 2015-09-26 16:26:21 -0400 |
| commit | 82bcc4121025f1212b3fcc07224cababa410f32a (patch) | |
| tree | fd16ecc292a6f0ae9dd11ab8ec6ad25693648935 /R Scripts/analyze-cascades.R | |
| parent | 786714e2a9acdbd650092bcaf351128dbebd160d (diff) | |
| download | criminal_cascades-82bcc4121025f1212b3fcc07224cababa410f32a.tar.gz | |
demographic stats
Diffstat (limited to 'R Scripts/analyze-cascades.R')
| -rwxr-xr-x | R Scripts/analyze-cascades.R | 4 |
1 files changed, 3 insertions, 1 deletions
diff --git a/R Scripts/analyze-cascades.R b/R Scripts/analyze-cascades.R index f03af36..02f9349 100755 --- a/R Scripts/analyze-cascades.R +++ b/R Scripts/analyze-cascades.R @@ -11,8 +11,10 @@ sizes = 1:max(data$V1) counts = rep(0,max(sizes)) counts[data$V1] = data$V2 pl = power.law.fit(counts) -plot(sizes,counts,log='xy',pch=20,col='#377EB8', +plot(sizes,counts,log='xy',pch=20,col='#1f78b4',xaxt='n',yaxt='n', xlab='Size of Cascade', ylab='Number of Cascades', main='') +axis(1,at=c(1,5,10,50,100,500)) +axis(2,at=c(1,10,100,1000)) sizes_pl = displ$new(counts[counts>1]) est = estimate_xmin(sizes_pl) |