From 28945245f49fafed3d05579166a46530845d2d8a Mon Sep 17 00:00:00 2001 From: Ben Green Date: Tue, 22 Sep 2015 08:28:13 -0400 Subject: script to plot entire lcc --- R Scripts/fit-background.R | 8 +++++--- 1 file changed, 5 insertions(+), 3 deletions(-) (limited to 'R Scripts/fit-background.R') 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))} -- cgit v1.2.3-70-g09d2