From 82bcc4121025f1212b3fcc07224cababa410f32a Mon Sep 17 00:00:00 2001 From: Ben Green Date: Sat, 26 Sep 2015 16:26:19 -0400 Subject: demographic stats --- R Scripts/analyze-cascades.R | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) (limited to 'R Scripts/analyze-cascades.R') 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) -- cgit v1.2.3-70-g09d2