diff options
| author | Ben Green <bgreen@g.harvard.edu> | 2015-09-18 23:49:02 -0400 |
|---|---|---|
| committer | Ben Green <bgreen@g.harvard.edu> | 2015-09-18 23:49:02 -0400 |
| commit | 8ac291db5fa8bd20df1ae6fc74970ca223e1c810 (patch) | |
| tree | 3ebe9717048f1d7a748d514dfd63468c181a37ef /R Scripts/fit-background.R | |
| parent | bb68c49552f69b90ce840dbd7346fc76eb93e3e3 (diff) | |
| download | criminal_cascades-8ac291db5fa8bd20df1ae6fc74970ca223e1c810.tar.gz | |
updated images
Diffstat (limited to 'R Scripts/fit-background.R')
| -rw-r--r-- | R Scripts/fit-background.R | 6 |
1 files changed, 4 insertions, 2 deletions
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") + |