dajmcdon / rtestim

https://dajmcdon.github.io/rtestim/
Other
4 stars 0 forks source link

The degree = 0L case is wrong #48

Closed dajmcdon closed 1 year ago

dajmcdon commented 1 year ago
y <- c(1, rpois(100, dnorm(1:100, 50, 15) * 500 + 1))
out <- estimate_rt(y, nsol = 10)
plot(out)
out0 <- estimate_rt(y, degree = 0L, nsol = 10)
plot(out0)

Seems to be operating on cases rather than Rt.