diff options
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))} |