k-sys / covid-19

A collection of work related to COVID-19

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

Stability in NY estimates

batson opened this issue · comments

I'm trying to get a handle on how to compare estimates from day to day. Here is the delta between yesterday's estimate for NYC and today's from the website rt.live.

Yesterday:
Screen Shot 2020-04-18 at 6 26 29 PM

Today:
Screen Shot 2020-04-19 at 1 14 12 PM

My intuition is that the model is local (modulo the 9 day gaussian smoothing), and forward looking, so adding one day of data shouldn't change the error bars (eg around March 22), nor so dramatically shift the endpoint. A large change in the inferred sigma could cause the expansion of the error bars.

Was there a change in data source?

I think the issue is that the data is for the whole state of NY, not just NYC. I think looking into NYC, we'd see the Rt well above 1. I think it would be more useful if the Rt charts could drill down to cities. JHU data gives the cities in the column "Combined_Key"

The city level data would also be interesting.

My concern here is less about the precise value of Rt in NY, and more about the changes in the model's predictions and error bars between consecutive days/runs. (This may indicate a change of data sources, model, or sensitivity of the whole predicted path to a final day of data.)

commented

Hi guys, actually this was a change in model - we (prematurely) switched to a new method and decided that wasn't right. Instead we decided to keep with the notebook math so that there's a 1:1 relationship.

Separately, the smoothing as it stands today overweights recent observations which is a known issue and I'll be working on a solution over the next day. Sorry for the confusion - you should see lots more stability going forward.

Thanks for the clarification, @k-sys!