diff options
Diffstat (limited to 'R Scripts')
| -rwxr-xr-x | R Scripts/analyze-cascades.R | 1 | ||||
| -rw-r--r-- | R Scripts/fit-background.R | 6 |
2 files changed, 4 insertions, 3 deletions
diff --git a/R Scripts/analyze-cascades.R b/R Scripts/analyze-cascades.R index 04cb0e3..9a47ece 100755 --- a/R Scripts/analyze-cascades.R +++ b/R Scripts/analyze-cascades.R @@ -13,7 +13,6 @@ counts[data$V1] = data$V2 pl = power.law.fit(counts) plot(sizes,counts,log='xy',pch=20,col='#377EB8', xlab='Size of Cascade', ylab='Number of Cascades', main='') -lines(c(1,1)) sizes_pl = displ$new(counts[counts>1]) est = estimate_xmin(sizes_pl) diff --git a/R Scripts/fit-background.R b/R Scripts/fit-background.R index 73ab6b2..506f6e3 100644 --- a/R Scripts/fit-background.R +++ b/R Scripts/fit-background.R @@ -1,6 +1,6 @@ library(igraph) setwd("~/Documents/Violence Cascades/") -load('Results/hyper-lcc.RData') +load('Raw Data/hyper-lcc.RData') # load data hyp_lcc = upgrade_graph(hyp_lcc) @@ -28,5 +28,7 @@ res = nls(formula=fit_form, data=infs, start=list(lambda=3, A=2, phi=4)) co = coef(res); co plot(t) -plot(days,counts,pch=20,cex=0.4,xlab='Day',ylab='Number of Infections') +plot(days,counts,pch=20,cex=0.9,col='#377EB8', + xlab='Day',ylab='Number of Infections') curve(fit(x, lambda=co["lambda"], A=co["A"], phi=co["phi"]), add=TRUE ,lwd=5, col="steelblue") + |