diff options
| author | Ben Green <bgreen@g.harvard.edu> | 2015-09-22 08:28:13 -0400 |
|---|---|---|
| committer | Ben Green <bgreen@g.harvard.edu> | 2015-09-22 08:28:16 -0400 |
| commit | 28945245f49fafed3d05579166a46530845d2d8a (patch) | |
| tree | 9703409e142334339d5c145f9415501be8597ec7 /R Scripts/fit-background.R | |
| parent | 8c2d3d070a3db1b469e3e32e3c20ef67bc274b3b (diff) | |
| download | criminal_cascades-28945245f49fafed3d05579166a46530845d2d8a.tar.gz | |
script to plot entire lcc
Diffstat (limited to 'R Scripts/fit-background.R')
| -rw-r--r-- | R Scripts/fit-background.R | 8 |
1 files changed, 5 insertions, 3 deletions
diff --git a/R Scripts/fit-background.R b/R Scripts/fit-background.R index 06e91c9..5bc1ce0 100644 --- a/R Scripts/fit-background.R +++ b/R Scripts/fit-background.R @@ -7,14 +7,16 @@ load('Raw Data/lcc.RData') vic_days = as.numeric(unlist(lcc_verts[,16:21])) vic_days = vic_days[!is.na(vic_days)] -hist(as.numeric(vic_days,2000,col='lightblue', +hist(vic_days,2000,col='lightblue', xlab='Day of Study Period',main='Infections During the Study Period') # get infection counts per day days = 1:3012 -t = table(factor(hyp_lcc_verts$vic.day[vic_ids],levels=days)) +t = table(factor(vic_days,levels=days)) counts = as.vector(t) -infs = data.frame(days,counts) +start_date = as.Date("2005-12-31") +dates = start_date + days +infs = data.frame(days,counts,dates) # define background function fit = function(x, lambda, A, phi) {lambda + A*(sin((2*pi/365.24)*x+phi))} |