diff options
Diffstat (limited to 'R Scripts/fit-background.R')
| -rw-r--r-- | R Scripts/fit-background.R | 8 |
1 files changed, 4 insertions, 4 deletions
diff --git a/R Scripts/fit-background.R b/R Scripts/fit-background.R index d1799cd..73ab6b2 100644 --- a/R Scripts/fit-background.R +++ b/R Scripts/fit-background.R @@ -17,14 +17,14 @@ counts = as.vector(t) infs = data.frame(days,counts) # define background function -fit = function(x, lambda, A, phi) {lambda + A*(sin((pi/365)*x+phi))^2} -fit_form = counts ~ lambda + A*(sin((pi/365)*days+phi))^2 +fit = function(x, lambda, A, phi) {lambda + A*(sin((2*pi/365.24)*x+phi))} +fit_form = counts ~ lambda + A*(sin((2*pi/365.24)*days+phi)) # explore data plot(days,counts) -curve(fit(x, lambda=1, A=12, phi=2.8), add=TRUE ,lwd=2, col="steelblue") +curve(fit(x, lambda=3, A=2, phi=4), add=TRUE ,lwd=4, col="steelblue") -res = nls(formula=fit_form, data=infs, start=list(lambda=1, A=12, phi=2.8)) +res = nls(formula=fit_form, data=infs, start=list(lambda=3, A=2, phi=4)) co = coef(res); co plot(t) |