Skip to Main Content

COVID-19 projections for Connecticut

May 28, 2020

Associate Professor of Biostatistics, Statistics & Data Science, Operation, and Ecology & Evolutionary Biology

5.25.2020 Biostatistics Seminar

ID
5250

Transcript

  • 00:00- From a colleague
  • 00:05asking for help with planning for the intensive care unit
  • 00:10and floor bed capacity at the Yale New Haven Hospital
  • 00:14Health System and Yale New Haven in particular.
  • 00:18Margret and Sohei had previously, or around the same time,
  • 00:21been working with the statistics policy
  • 00:25modeling an epidemiology collective on a queuing model
  • 00:28or discussing the parameters of the queuing model
  • 00:32for the dynamics of Covid-19 patient flow through hospitals.
  • 00:37So we decided to use this model setup to make a concrete
  • 00:41software product in the form of a web application
  • 00:44that Yale New Haven Health System and other hospital systems
  • 00:47could use for capacity planning.
  • 00:50We wanted to respond to their very immediate need
  • 00:54to know how full the hospital would get if Covid patients
  • 01:00kept coming at the rates that they were seeing
  • 01:02and how they might expand capacity to accommodate
  • 01:06these new patients.
  • 01:09So we created a Slack channel,
  • 01:12a way of communicating directly in real time
  • 01:15with the team members, who created a GitHub repository.
  • 01:17Within, I think, only about two hours,
  • 01:20we had a web application written in R,
  • 01:22using Shiny framework,
  • 01:25where you could sort of dial in the
  • 01:30current bed capacity at a hospital system.
  • 01:32You could enter parameters that govern the length
  • 01:35of stay of Covid patients and how they move through
  • 01:37the hospital from the emergency department
  • 01:40to the floor to the ICU
  • 01:41and then toward discharge or possibly death.
  • 01:47So that product went live very, very quickly.
  • 01:53There are many other collaborators and contributors
  • 01:54to the application beyond just our group.
  • 01:59Our goal here was to produce
  • 02:03something very quickly and immediately useful.
  • 02:08The structure of this model is shown
  • 02:10in this very complicated diagram.
  • 02:12It's not as complicated as it looks.
  • 02:14The basic idea is that patients enter through
  • 02:18to the emergency department.
  • 02:21They move to the floor then to the ICU.
  • 02:27There are many things that could happen
  • 02:28if those places are full.
  • 02:31Each of those parts of the hospital is treated as a queue.
  • 02:34That is, it's essentially a pool of patients
  • 02:37who are waiting to exit.
  • 02:39One of the ways they can exit is to step up
  • 02:41from the floor to the ICU.
  • 02:43One of the ways they can exit is to die.
  • 02:45Another is to be discharged
  • 02:47if they are no longer acutely ill.
  • 02:50So sort of taking into account all of this schematic, this
  • 02:56stylized depiction of the way Covid patients
  • 03:00would flow through a hospital,
  • 03:02we wrote a system of ordinary differential equations,
  • 03:05which describe formally, the dynamics of this system.
  • 03:10It's a very simple type of modeling
  • 03:13that is very useful when the number of patients is large
  • 03:17and when you want sort of aggregate dynamics over time.
  • 03:21So we're not modeling, it's not an agent-based model.
  • 03:24We're not modeling individual patients trajectories
  • 03:26through the hospital.
  • 03:27Rather, this idea of patient flow through the hospital.
  • 03:31So this model, depicted schematically here,
  • 03:34is formalized in a system about ordinary differential
  • 03:37equations with many parameters.
  • 03:39Those parameters are calibrated to data that we have
  • 03:41from the Yale New Haven Health System
  • 03:46and to values from the literature.
  • 03:49We wrote this web application, which is now live
  • 03:53at the Shiny apps URL that you can see below.
  • 03:56You can interact with it if you like.
  • 04:01It basically allows the user to specify time horizon,
  • 04:04how quickly or slowly they think new Covid patients
  • 04:07will present to the emergency department
  • 04:10and then on subsequent tabs, you can dial in the current
  • 04:15hospital capacity at your institution.
  • 04:17You can dial in capacity increases that you anticipate
  • 04:20being able to implement into the future
  • 04:23to see how dynamics would change if say,
  • 04:25you could add 100 new ICU beds over the course of two weeks
  • 04:30a month from now, for example.
  • 04:32Then there are many, many input parameters.
  • 04:34Things like the age-specific rates
  • 04:39of death or of stepping up from the floor in the ICU,
  • 04:42to the average length of stay in each of those compartments
  • 04:47for patients who come to the hospital.
  • 04:50You can generate reports, downloadable PDF reports.
  • 04:54We sort of envisioned this tool being responsive
  • 04:57to the needs of hospital decision makers
  • 05:00who wanted to be able to add this planning capability
  • 05:04to their existing bed management software applications
  • 05:09and then to be able to generate reports for say,
  • 05:12supervisors and higher up decision makers
  • 05:14that would describe the scenario that the analysts
  • 05:17was most interested in.
  • 05:19The reports would also describe the consequences
  • 05:22of a capacity expansion strategy that might be implemented
  • 05:25by the system.
  • 05:27So I think this tool was very useful
  • 05:30to the Yale New Haven Health System.
  • 05:33It was publicized kind of broadly and we got some interest
  • 05:36from hospital systems throughout the U.S..
  • 05:39I had spoke to some of them about the ways
  • 05:41that they were making decisions, planning capacity increases
  • 05:45and using this application and others
  • 05:47that are also publicly available online,
  • 05:50to help guide their decision making.
  • 05:54This is an open source project.
  • 05:56You can get all of the source code for the Shiny application
  • 05:59on our GitHub repository here, shown below.
  • 06:07So what are the next steps for this project?
  • 06:10Fortunately, hospitalization in Connecticut is declining.
  • 06:15This figure that I've shown here is kind of compressed.
  • 06:17It's declining slowly.
  • 06:21But it has been declining for I think,
  • 06:23more than three weeks now.
  • 06:25Yale New Haven Health System, along with hospitals
  • 06:28heath systems throughout the state,
  • 06:30are doing much better than they were in mid-April.
  • 06:33They have enough bed capacity to accommodate
  • 06:35all the Covid patients and many more who may arrive
  • 06:39in the coming months.
  • 06:40So this is very good news for the hospitals
  • 06:42and for the state.
  • 06:43It's one of the reasons that the governor initiated
  • 06:48the first phase of the reopening plan on May 20 this week.
  • 06:53However, a lot of the projections and some that I'll show
  • 06:56in a few minutes, indicate a substantial risk of resurgence
  • 06:58of new cases, hospitalizations and deaths
  • 07:01following reopening the state.
  • 07:03This resurgence is anticipated to occur in July,
  • 07:06August, maybe September, depending on how things go
  • 07:09with reopening.
  • 07:11So I think that
  • 07:13the model, the web application,
  • 07:16and this work in general will unfortunately,
  • 07:20become useful again and very relevant again
  • 07:23later on in the summer if hospitalization
  • 07:25of Covid patients increases again.
  • 07:28So we want to maintain our capacity to continue developing
  • 07:31this model and responding to the needs
  • 07:34of decision makers within hospital systems.
  • 07:37We're taking this down time though,
  • 07:38to write a technical report and a lessons learned paper
  • 07:42about the way that we interact with health systems
  • 07:46and how we might improve the way
  • 07:48that we do that in the future.
  • 07:50This work is of course also gotten us very interested
  • 07:52in the ways that hospitals manage Covid patients.
  • 07:55We're very interested in comparative evaluation
  • 07:57and comparative effectiveness in the evaluation
  • 08:00of Covid-19 medical interventions.
  • 08:03That's something that Margret Erlensdottir
  • 08:06an MD PhD student in biostat is working on.
  • 08:09- Forrest? - Yes?
  • 08:10- Can you take a question?
  • 08:12- Yes, please go ahead.
  • 08:16- Have you only applied this to Yale New Haven?
  • 08:19- We have, the model itself is generic.
  • 08:22This is a good question, but we have calibrated many
  • 08:26of the length of stay and probability parameters
  • 08:30based on data that we received from Yale New Haven.
  • 08:34So in that sense, the dynamics that we present by default
  • 08:39are specific to Yale New Haven.
  • 08:41The user has the ability to change all of those parameters,
  • 08:45so we anticipate that this could be useful
  • 08:47for hospital systems of any size
  • 08:50with different patient demographics,
  • 08:51different age distributions for example.
  • 08:54So we want it to be as useful as possible,
  • 08:57but having said all this, the customer in this case,
  • 08:59was very clearly for us, Yale New Haven
  • 09:02and they had a very specific need and--
  • 09:04- Have you had a reaction?
  • 09:06Did you have an ongoing reaction with the people
  • 09:08at Yale New Haven who were using this product,
  • 09:10whether or not it was helping them or was it accurate
  • 09:13or did they have any complaints about it?
  • 09:16I'm sure they did.
  • 09:16Can you tell me about that interaction?
  • 09:19- Sure, sure.
  • 09:22I think that they made a few requests of us.
  • 09:25Some of them were very qualitative.
  • 09:27They wanted very early to be able to generate reports.
  • 09:30A lot of the requests were for additional functionality
  • 09:33rather than additional structure in the OD model
  • 09:37but they really,
  • 09:39I think many of the requests were about flexibility
  • 09:41and granularity in the predictions.
  • 09:44They wanted to be able to dial in the exact patient
  • 09:46demographics and the care parameters that were actually
  • 09:50being implemented at Yale New Haven.
  • 09:52So we tried to give them that ability and that control.
  • 09:58I think mostly, successfully.
  • 09:59We retained some of the generality of the model,
  • 10:02while allowing users to input the parameters
  • 10:05that they felt were right for their system.
  • 10:07In terms of the way it was used at Yale New Haven,
  • 10:11I think that by the time they asked us for help,
  • 10:15many of the actual capacity expansions
  • 10:18had already been implemented.
  • 10:19I'm talking about taking over high school gymnasia,
  • 10:22changing the configuration of parking lots
  • 10:25to provide drive through testing and turning,
  • 10:30I guess parts of the hospital into ICUs.
  • 10:33Many of those--
  • 10:36- You might say that they over expanded a little bit
  • 10:38since they quickly came not needed capacity.
  • 10:41So did you help them, saying hey,
  • 10:43you guys don't need to do that much?
  • 10:47- I think that based on the projections for population
  • 10:50level incidence that they were receiving in early
  • 10:52to mid-April, the capacity expansion was appropriate.
  • 10:59This model here did not provide
  • 11:03population level projections,
  • 11:05which I'll show in a few minutes.
  • 11:06So we were not telling them
  • 11:07that they had over expanded capacity.
  • 11:11I think that at the state level,
  • 11:13the total hospitalization in the state came very close
  • 11:17to the preexisting capacity, as it was in early March.
  • 11:22So I think that there was a big concern that it was unclear
  • 11:26what the doubling rate of new cases would be.
  • 11:29We had not yet seen some of the benefits of state lock down
  • 11:31and closure of schools.
  • 11:33So the hospital systems were expanding very aggressively,
  • 11:37I think for good reason.
  • 11:40- Okay, but they were just doing that by looking at
  • 11:44the daily or maybe the weekly case counts right,
  • 11:47and seeing what the doubling rate was and things like that.
  • 11:50They were doing anything more subtle than that?
  • 11:52- That is what they were doing when they called us on.
  • 11:55We tried to give them projections under their own
  • 12:00in-house assumed doubling rates.
  • 12:03So we were very interested in showing them when the hospital
  • 12:06would fill up and under what circumstance
  • 12:08and how different parts of the hospital would fill up.
  • 12:13- Okay.
  • 12:14I'll let you go on.
  • 12:16- In a few minutes I'll show state level projections
  • 12:18that might answer some of your questions.
  • 12:24All right.
  • 12:28- By the way, I'm John Hardigen, by the way.
  • 12:31Used to be in the statistics department.
  • 12:33- Yes, I know, good to see you.
  • 12:38All right, second project.
  • 12:40On April 14, so just as we
  • 12:47finished the most fundamental software development
  • 12:50on the application that I just showed you,
  • 12:52on April 14 we were asked to start producing projections
  • 12:55for the governor's Reopen Connecticut Advisory Panel,
  • 12:58which was charged with
  • 13:01making recommendations to the state, to the government,
  • 13:04to the Department of Public Health
  • 13:06on how reopening should proceed and what the timeline
  • 13:09should be and what business sectors could safely reopen
  • 13:12at which times.
  • 13:14The panel consisted of public health researchers,
  • 13:17including Albert Ko and several other people from Yale
  • 13:20and many business leaders in Connecticut.
  • 13:23It was a mixed group.
  • 13:27The panel needed projections at that time of Covid-19
  • 13:31incidence, hospitalizations and deaths
  • 13:33under future reopening scenarios, to plan testing expansion,
  • 13:37seroprevalence studies and most importantly,
  • 13:39to assess the risk of a second wave of infections.
  • 13:43So this was in mid-April,
  • 13:45around the time when hospitalization was peaking.
  • 13:48Of course, nobody knew exactly at that time
  • 13:51that the peak was occurring and there was a lot of concern
  • 13:55that things would continue to get much worse,
  • 13:58in terms of hospitalization in Connecticut.
  • 14:03The work of that committee
  • 14:07advised the governor in his reopening strategy,
  • 14:09which we've all probably heard about,
  • 14:11if you're following press releases from the state.
  • 14:14The state began reopening on May 20th and there's now,
  • 14:17I think, although the work
  • 14:19of the advisory panel may be wrapping up,
  • 14:22there's now an ongoing need for projections
  • 14:24to inform decision making and epidemiological study design,
  • 14:28that further informs decision making
  • 14:32for the Connecticut response and reopening.
  • 14:36This part that I'll talk about now is joint work
  • 14:37with Olga Morozova and Richard Li.
  • 14:42So at the beginning of this project, we had to explain
  • 14:46to decision makers and members of the advisory panel,
  • 14:50how data are different from model projections
  • 14:54and what sort of...
  • 14:59What the differences between these two products were.
  • 15:01But I think there is a recognition at that time
  • 15:03on the part of policy makers and committee members
  • 15:05that the policy makers have access
  • 15:07to a real-time data stream, which is very high quality.
  • 15:11They have access to all sorts of state dashboards
  • 15:14describing the current state of the Connecticut pandemic.
  • 15:18They know about hospitalization and bed capacity
  • 15:20information from the Connecticut Hospital Association.
  • 15:23They know about test counts and nearly real-time
  • 15:25case counts, number of tests positive at hospitals
  • 15:28and in the community.
  • 15:30They know how many deaths have occurred to attributable
  • 15:33to Covid-19 or that are suspicious,
  • 15:37that are possibly related.
  • 15:39They might have information about excess deaths
  • 15:42that are not attributed to Covid-19 but are above
  • 15:45and beyond what you might normally expect in a typical year.
  • 15:49They have access to all this information.
  • 15:51They have access to very responsive staff
  • 15:53and many very smart people working for the state
  • 15:57Department of Public Health and other state agencies.
  • 16:02So there might be a sense that policy makers have access
  • 16:04to all the information and the most timely information
  • 16:07they could possibly need to make good decisions
  • 16:08for the state.
  • 16:10We tried to argue that there was more information
  • 16:12that they might be able to use constructively
  • 16:16to guide reopening, and that was information that was not
  • 16:19directly derived from contemporaneous data streams,
  • 16:24but rather these would be projections from transmission
  • 16:27models about possible futures.
  • 16:30So projections here can tell us about what might
  • 16:33happen in the future, possible hypothetical
  • 16:34or counterfactual scenarios to be defined
  • 16:39by the governor and the outcomes that would occur
  • 16:42under those reopening scenarios.
  • 16:44So I'm talking about phases, business sectors,
  • 16:47reopening back-to-school, what might happen in late August,
  • 16:50early September, as children go back to school
  • 16:52or back to summer camp in June and July.
  • 16:56What might happen under expanded testing
  • 17:00and contact tracing or continue to modified
  • 17:03social distancing guidelines.
  • 17:04Things like wearing masks or keeping six feet apart
  • 17:10and all of those things.
  • 17:11So we tried to explain how projections from these types
  • 17:15of models might be very different from simple plots
  • 17:19of the data streams that policy makers have access to.
  • 17:23This is a figure I showed them at the very beginning.
  • 17:26On the left, we have the number of death,
  • 17:28I think by early May, that had accumulated in Connecticut.
  • 17:33These are the red dots on the left hand side.
  • 17:36On the right hand side, we have a projection of what might
  • 17:38occur in the future on this day and I think it was
  • 17:40first week of May.
  • 17:42Right, and I think this may seem silly as a projection
  • 17:47exercise or it seems silly to
  • 17:52make a distinction between data and predictions,
  • 17:56but it may have useful in the setting to emphasize
  • 17:59that the real-time data that policy makers were using
  • 18:02was just the stuff on the left
  • 18:04and that if one believed the assumptions underlying
  • 18:07some of these dynamic transmission models,
  • 18:10that they could be provided with the stuff on the right,
  • 18:13which would be a projection
  • 18:14of what might happen in the future.
  • 18:16Here, I happen to have shown projections starting
  • 18:19on March 1, just to emphasize sort of how the line follows
  • 18:24the data points in the projection.
  • 18:27But the idea is that these projections would come
  • 18:29with some sort of uncertainty windows or sets
  • 18:34that would represent, in some sense,
  • 18:38the most likely possible futures under what we know today
  • 18:41and what we believe may happen about the future.
  • 18:44So the-- - May I stop you for a second?
  • 18:47- Of course. - Forrest.
  • 18:50First of all, I think you'll agree that the points
  • 18:53on the line at the left are extremely highly correlated
  • 18:55with each other, since they're just cumulatives.
  • 18:59And that's not a good way to show what's happening,
  • 19:01is to look at cumulatives.
  • 19:02You have to kind of guess what the derivatives are
  • 19:06and people aren't so good at that.
  • 19:07You would be much better off trying to project
  • 19:10and look at say, the weekly values.
  • 19:13Certainly can't look at daily values because God knows
  • 19:15what the daily values goes from,
  • 19:17but you know, you see in a week they kind of catch up
  • 19:19with the truth.
  • 19:20So if you looked at weekly values, you would tell on
  • 19:22what the present situation was.
  • 19:24Surely, that's what the hospitals need to know.
  • 19:27They don't need to know how many people they had
  • 19:29a long time ago or what the total was.
  • 19:31They want to know that the present charge is.
  • 19:33So I would just suggest that the thing you should be
  • 19:36working on is something closer.
  • 19:38Can't use daily values, it's too small,
  • 19:40but a weekly value and then that's what really matters.
  • 19:43That's the present situation.
  • 19:46- Certainly and have access to all that information.
  • 19:48The State Department of Public Health produces
  • 19:50weekly smoothed and unsmoothed count.
  • 19:54In fact, daily counts as well.
  • 19:56They're very volatile.
  • 19:58They jump up--
  • 19:59- The daily counts have a huge weekly effect.
  • 20:02You just don't want to rely on them at all.
  • 20:05The docs aren't bothered to do things on the weekends
  • 20:07is my interpretation of it.
  • 20:08But maybe it's someone not bothering, but whatever it is,
  • 20:11it's a big weekly effect.
  • 20:12It's something you don't want to have.
  • 20:14But if you take a weekly value, that's always averaged out.
  • 20:17I just think that projecting the future
  • 20:20and I think you would find there's quite a lot
  • 20:22more error in that.
  • 20:23You're getting the benefit of the fact that all this
  • 20:26is highly correlated but if you were trying to project
  • 20:28the future, these things would be whoa, of stuff.
  • 20:32- Yes, totally agree.
  • 20:33This figure was generated in response to a very specific
  • 20:36question, which is how many deaths will the state expect
  • 20:39to have accumulated on a future date.
  • 20:42- Okay.
  • 20:45Thanks.
  • 20:48- Okay, so I wanted to answer this question
  • 20:50because I hope that you're all wondering about it.
  • 20:54Does the world need another Covid-19 projection model?
  • 20:58There are lots of them out there.
  • 21:00Vary in quality, some from very experienced research groups
  • 21:05and experienced epidemiologists, some from
  • 21:09Silicon Valley software developers
  • 21:11who just learned about regression.
  • 21:13I don't think that the world needs another Covid-19 model
  • 21:17at the national or international level.
  • 21:19But I think Connecticut does for several reasons
  • 21:23that I wanted to describe briefly here.
  • 21:27We wanted to develop a scenario analysis tool
  • 21:29that was responsive to specific questions
  • 21:31from the Connecticut leadership,
  • 21:33who were planning to reopen the state.
  • 21:38We thought there were several reasons that we could add
  • 21:39some value here, beyond what is provided by some of the more
  • 21:42generic models that are available
  • 21:47for national, state and also local projections.
  • 21:50The first thing is access to epidemiologists
  • 21:53at the School of Public Health
  • 21:54and in the Public Health Modeling Unit.
  • 21:57We have pretty unique access to data
  • 21:59from the Connecticut Hospital Association
  • 22:01on the bed capacity and bed occupancy throughout the state.
  • 22:06We can use information on individual patient trajectories
  • 22:10through the healthcare system from using data
  • 22:12from Yale New Haven.
  • 22:14We have access to empirical epidemiological studies
  • 22:18from Yale emerging infections program and data streams
  • 22:21from the Department of Public Health through Yale EIP.
  • 22:25We have connection to the people who are running
  • 22:28the testing and seroprevalence studies
  • 22:30to be conducted in the future and the model projections
  • 22:34that we produce will be very closely tied
  • 22:38to the conduct of those studies.
  • 22:39Some of them can give information that we can use
  • 22:42for calibrating the model, and in turn,
  • 22:44we can use model projections
  • 22:48to provide preliminary estimates of say,
  • 22:51cumulative incidence of Covid-19 for study planning,
  • 22:56in order to do sample size calculations.
  • 23:00And of course, we are hoping to be able to help
  • 23:04with the Department of Public Health's
  • 23:07implementation of optimal testing and sampling strategies
  • 23:10as they look for new cases and try to control outbreaks
  • 23:13that may occur in the future in Connecticut.
  • 23:20So the modeling principle here,
  • 23:22this is an infections disease model that I'm gonna show you.
  • 23:25It's not a model for hospital
  • 23:29patient flow through hospitals.
  • 23:31But I think in introducing this to people who have not seen
  • 23:35these models before, the operating principle
  • 23:38is that of mass action.
  • 23:41I think if mathematical infectious disease epidemiology
  • 23:44has a central dogma or a single principle that governs
  • 23:47the structure of quantitative models for infections,
  • 23:51it's something like the Law of Mass Action,
  • 23:54that in a small time interval, the number of new cases
  • 23:57that accrue is proportional to the number of ways
  • 24:01that susceptible individuals and infectious individuals
  • 24:05can come together.
  • 24:07This means that new cases or incidences
  • 24:09is driven by the product of--
  • 24:11Or sorry, I should have said the product of susceptibles
  • 24:14and infectives or the number of ways that people
  • 24:18susceptible individual can come into contact
  • 24:20with an infected person.
  • 24:23This general principle is what underlies all transmission
  • 24:27models and many transmission models are compartmentalized
  • 24:29or they are separated in space and geography
  • 24:33or by age group or by different risk categories,
  • 24:36but this is the essential principle.
  • 24:38That new cases of a certain type and a certain place
  • 24:41arise at a rate that is proportional to the product
  • 24:45of the number of susceptibles and infectives.
  • 24:47The number of ways that disease can be transmitted.
  • 24:53So we have divided the population of Connecticut
  • 24:56into many compartments.
  • 24:59Those who have not had the disease,
  • 25:01those who are susceptible,
  • 25:03those who have been infected,
  • 25:06they are exposed but not yet infectious.
  • 25:08So they don't have symptoms.
  • 25:11Those who are infectious but remain asymptomatic,
  • 25:14those with mild symptoms.
  • 25:16They know they're sick
  • 25:17but they do not require hospitalization.
  • 25:19Those with severe symptoms who do require hospitalization.
  • 25:24Those who have mild symptoms but are successfully isolated
  • 25:28because they realized they have symptoms or they got
  • 25:33a viral test that told them that they are infected.
  • 25:36So they successful isolate themselves.
  • 25:38Those people with severe disease who are hospitalized,
  • 25:42those who have severe disease but remain unhospitalized
  • 25:45because there's no space for them.
  • 25:47This is very important in projecting deaths in the future
  • 25:50scenario, in which we run out of hospital capacity.
  • 25:55Then we have severe institutionalized populations,
  • 25:58who are not in the hospital,
  • 25:59such as people in nursing homes, correctional institutions
  • 26:03and other long-term care facilities.
  • 26:06Those who have been infected but did not die
  • 26:09and are now recovered or successfully isolated
  • 26:12and recovering and those who have died.
  • 26:15So the idea here is to divide up the population
  • 26:18of Connecticut into a number of people
  • 26:22in each of these compartments.
  • 26:28The model that we put together is a variation
  • 26:30on the susceptible exposed infected and removed model.
  • 26:38We divide up the infectious individuals into three
  • 26:42categories that I told you about, severe, mild
  • 26:44and asymptomatic infections.
  • 26:47We have two different types of patients
  • 26:49who need hospitalization.
  • 26:51We have unhospitalized patients.
  • 26:55We can remove patients by isolating them
  • 26:57and they can recover after some amount of time,
  • 27:01if they do not die.
  • 27:04This is the basic structure of the SEIR model.
  • 27:09The usual model structure is just this linear part,
  • 27:12SEI and then R.
  • 27:15We divided up into these additional components,
  • 27:19not because we believed that these components
  • 27:21cover every possible scenario or every possible type
  • 27:25of illness or state of the world or state of patients,
  • 27:28but because this is the most parsimonious model
  • 27:31that we can think of that captures the dynamics
  • 27:34of infection that are most likely to lead to the outcomes
  • 27:38that a state government cares most about.
  • 27:40Those are state-level hospitalizations and deaths
  • 27:44and possibly cumulative incidents.
  • 27:48Right, so this model is not intended to capture
  • 27:51every biological or epidemiological feature of Covid-19
  • 27:55transmission in Connecticut.
  • 27:57Rather, it is the simplest model that captures the features
  • 28:00that policy makers care most about.
  • 28:06It's also structured by geography.
  • 28:10We found that the...
  • 28:14We looked at information about travel and commuting patterns
  • 28:17throughout the state to look at where people
  • 28:18might be mixing, where they live, where they work,
  • 28:21others things like that.
  • 28:22But we found that that information did not give us
  • 28:26much more information than simple adjacency matrix
  • 28:28of counties in the state.
  • 28:33We're well aware that many people in Connecticut
  • 28:38work or commute or travel often to New York City area.
  • 28:41We'll try to accommodate that in the model
  • 28:45or in our interpretation of the model.
  • 28:48Rather, the adjacency matrix of counties in Connecticut
  • 28:54gives us much of the information that we use
  • 28:55for the geographically dependent nature of transmission.
  • 29:04Basic idea--
  • 29:05- Rhode Island and Massachusetts
  • 29:07aren't doing too good either.
  • 29:09- They're not doing well, I agree.
  • 29:12To avoid turning this into a very granular or national model
  • 29:19we are going to treat the exogenous force of infection
  • 29:22experience by Connecticut residents as something else.
  • 29:28So we sort of imagined that it is
  • 29:31subsumed into the force of infection experience
  • 29:33by everyone in Connecticut.
  • 29:35I agree, both a lot of infections
  • 29:38and a lot of heterogeneity outside of Connecticut
  • 29:42in bordering states.
  • 29:45So most of this we don't specifically take into account.
  • 29:50The basic idea here, I'm just showing two compartments
  • 29:54of the ODE system, the basic idea is that in county I,
  • 29:59the number of susceptibles or the rate of new infections
  • 30:02is governed by the number of infectious individuals
  • 30:07in that county and the number of infectious individuals
  • 30:10in neighboring counties.
  • 30:13Here in beta is the transmission rate of infection.
  • 30:17So individuals who are susceptible transition to the exposed
  • 30:22infectious state and then to other states down the road.
  • 30:25But these are sort of the mass action equations
  • 30:28for a heterogeneous population in which the force
  • 30:30of infection is coming from outside and within
  • 30:34individual counties.
  • 30:39I'm not going to go into a great deal of detail
  • 30:42about the system of ODs that is most useful here.
  • 30:48I'll just say that we solve in numerically.
  • 30:50It's a system of 11 differential equations
  • 30:53given the parameters, which I'm just gonna bundle into
  • 30:56a vector theta.
  • 30:58Let Y of T given theta, be the solution to the OD system
  • 31:01at time T with parameters theta.
  • 31:03You can solve this system with pretty good accuracy
  • 31:07using modern OD solvers.
  • 31:10This solution--
  • 31:11- They're just linear equations, are they?
  • 31:13Linear right?
  • 31:15- They're non-linear in the right hand side
  • 31:18is non-linear in the other model compartments.
  • 31:21Right, that's what mass action is.
  • 31:24It's proportional to the product.
  • 31:26So OD is proportional to the product of S and I.
  • 31:31So it's--
  • 31:33- On the other hand, you agree that S doesn't change much
  • 31:35because unless you've got a very fully infected population,
  • 31:39S doesn't change that much.
  • 31:41You've got--
  • 31:42- S is most quickly when infections are increasing
  • 31:45most quickly and Connecticut right now,
  • 31:48S is still pretty large.
  • 31:50I think cumulative incidence is between 5% and 15%.
  • 31:55So S has not changed.
  • 31:57- S is 85% and is gonna change.
  • 31:59I'm just making it linear for myself, that's all.
  • 32:02- Sure, yeah.
  • 32:04So right now S has not decreased that much.
  • 32:08You know, between, it's still at 95% to maybe 85%,
  • 32:13something like that.
  • 32:15As the pandemic progresses and into the fall,
  • 32:17if there's another resurgence of infections,
  • 32:19we will expect S to change quite a lot more.
  • 32:22If it changes a lot, then we'll be in herd immunity
  • 32:25territory where depletion of susceptibles
  • 32:28plays a prominent role in altering the dynamics
  • 32:32of the pandemic, but we're not there yet.
  • 32:34- But I am right in thinking that this is linear,
  • 32:36it's really just a matrix problem isn't it,
  • 32:38that we have to solve.
  • 32:39- If it were linear it would be a matrix problem.
  • 32:41- Yeah, okay.
  • 32:44- Right. - Yeah.
  • 32:47- So this system is a deterministic system.
  • 32:52Engineers, mostly and some epidemiologists,
  • 32:55have been thinking for a very long time about principled
  • 32:58ways of estimating parameters for deterministic system.
  • 33:02Unfortunately, for models of this type,
  • 33:04which is generally the case in infectious disease
  • 33:06epidemiology, there are some serious
  • 33:11identifiability problems.
  • 33:11Not all parameters can be uniquely estimated from the data
  • 33:16or infinitely many combinations of parameters
  • 33:18that appear to fit equally well.
  • 33:23We only observe in this case,
  • 33:25the hospitalization and death compartments.
  • 33:27There's some information from PCR testing
  • 33:30about the prevalence of infection at different times,
  • 33:34but because the testing strategy in Connecticut
  • 33:36and elsewhere has varied so dramatically
  • 33:40over the last few months, we didn't feel like we could use
  • 33:42any information from testing alone to inform the sizes
  • 33:47of the currently infected compartments.
  • 33:51So basically, we're trying to estimate many parameters
  • 33:54for a system with 11 components using only the time series
  • 33:59of hospitalizations and deaths.
  • 34:02So it's quite challenging and in practice,
  • 34:04this necessitates taking parameter values
  • 34:08from the literature, from clinical studies,
  • 34:10from our knowledge of how hospitals treat patients
  • 34:14and also using a statistical estimation scheme
  • 34:18to learn about elements of theta, of the unknown parameters.
  • 34:22I wish that I could give you a more coherent statistical
  • 34:26inference strategy in which all of the parameters
  • 34:31were learned from the data and I could tell you
  • 34:34that they were being consistently estimated
  • 34:36and that as the epidemic went on, we would get more and more
  • 34:39precise estimates of each of those parameters.
  • 34:41Unfortunately, it's just not true.
  • 34:43That the model structure that we need here to be able
  • 34:46to accommodate
  • 34:51the structure of the pandemic is more complicated
  • 34:54than the model structure that we could possibly identify
  • 34:59non-parametrically or semi-parametrically
  • 35:01or even in this parametric model.
  • 35:05So I just wanted to give you some examples
  • 35:06of how people do this in practice.
  • 35:08These are not exactly endorsements
  • 35:10of statistical frameworks.
  • 35:12The basic idea is that given theta,
  • 35:14we can solve the ODE system, it gives us deterministic
  • 35:17solutions at time points where we have an observation
  • 35:21and then calibration or statistical inference
  • 35:23essentially amounts to minimizing a loss criteria
  • 35:26and are comparing the observed values to model predictions.
  • 35:31The two frameworks that are most frequently used here
  • 35:34are imposing a normal errors or gaussian errors,
  • 35:37almost normal gaussian errors
  • 35:41or equivalently minimizing at least squares type
  • 35:44of loss function or doing this plus on maximum likelihood
  • 35:49estimation for elements of theta
  • 35:51that you can identify in this way.
  • 35:53I think in this project we used
  • 35:58the Poisson maximum likelihood.
  • 36:01There are many things about this, one of which is that
  • 36:03a Poisson random variable could take values
  • 36:06that are larger than the size of the population.
  • 36:08In practice here, that's not what occurs
  • 36:12because the number of infections here is small,
  • 36:15but this is basically a framework for doing a type
  • 36:18of statistical inference or learning about
  • 36:22a posterior distribution on parameters
  • 36:25from a model, which gives deterministic predictions
  • 36:27and which doesn't have any inherent stochasticity.
  • 36:31The procedure that we used here, which I'm not gonna
  • 36:33talk about in great detail here, was developed
  • 36:36by postdoc Olga, is a hybrid approach that fixes
  • 36:39some parameters and imposes uncertainty distributions
  • 36:42on them from our prior knowledge and the literature
  • 36:46and conducts Bayesian posterior inference
  • 36:49on known parameters and initial conditions.
  • 36:51So we try to learn jointly about parameter--
  • 36:55Yes go.
  • 36:56- Forrest, there's a question people always ask of this
  • 36:58whenever I give a talk like this,
  • 37:00how do you determine your prior distribution?
  • 37:03- In this case, I would say we're in a very good position
  • 37:06to interpret priors as literally being prior beliefs.
  • 37:11We have for example, point estimates and confidence
  • 37:14intervals from published studies.
  • 37:19We also have parameters which are intrinsic to the model
  • 37:22but for which we have very little information.
  • 37:24So we assign to them, what we believe qualitatively,
  • 37:27to be an appropriate representation of our uncertainty
  • 37:31or ignorance about those parameters
  • 37:33under the parametrization.
  • 37:36But to your question--
  • 37:38- What you believe to be true then, is that right?
  • 37:42- Oh certainly.
  • 37:43It is a mixture of what other people believe to be true
  • 37:46and what we believe to be true as well.
  • 37:48So I would take a subject of interpretation
  • 37:49to the priors here.
  • 37:52They are subjective in the sense that we believe
  • 37:56these uncertainty distributions.
  • 37:59They are quantitative in the sense in that
  • 38:02some of them come from published studies.
  • 38:06- Okay, I'm sorry, I do just a little bit longer.
  • 38:11You know, I know you've got a lot of parameters in here,
  • 38:13many of which I don't know anything about,
  • 38:15but I suspect the very important one is parameter
  • 38:18which says what is the ratio of new cases,
  • 38:24assuming that susceptibility isn't changing
  • 38:26to the infection rate, right.
  • 38:28What's the...
  • 38:30That's the number,
  • 38:32that ratio is an important ratio.
  • 38:34New cases against the number that are infected
  • 38:38and that number out to extract is an important number
  • 38:42because it changes a lot, according to the conditions
  • 38:44that the government sets.
  • 38:46Changes all the time because you're trying to reduce
  • 38:48contacts and effectively reducing that contacts
  • 38:50is to change that ratio.
  • 38:52I assume that that's built into the model somehow,
  • 38:55but I would think you probably don't know very much
  • 38:57about how the government's policies and whatever
  • 39:00are gonna change that ratio.
  • 39:02So if you said you know, I know it's gonna be a month
  • 39:04from now, I'd say no you don't.
  • 39:06- Oh sure. - Yeah.
  • 39:08So how do you handle it?
  • 39:10- We certainly do parametrize that rate,
  • 39:14that is the transmission rate that you were talking about.
  • 39:16It's the parameter that multiples the product
  • 39:19of the number of susceptibles
  • 39:20and the number of infectious individuals.
  • 39:23That's called beta in the model.
  • 39:25Beta does change over time.
  • 39:27It's parametrized as a sum of step functions.
  • 39:33Those step functions change in their value
  • 39:36around when the governor closes schools, which happened,
  • 39:42I think on March 25th and when the governor--
  • 39:45Or sorry, a little bit earlier, maybe March 20th,
  • 39:48I can't remember.
  • 39:49Then when the governor issued the stay at home order,
  • 39:52the stay safe stay at home order,
  • 39:56which I think took effect on the 23rd.
  • 39:59So those step functions are in the model for historical
  • 40:02interventions that were implemented by the state.
  • 40:04For future interventions which are implemented by the state,
  • 40:07we are guessing.
  • 40:09Fortunately, we are guessing using information
  • 40:11from the people who will actually make those decisions.
  • 40:14So I will show how we assume that that transmission rate
  • 40:20or contact rate might change in the future
  • 40:23under guidelines expressed by the governor
  • 40:26and policy makers.
  • 40:28Right, so in the future of course,
  • 40:29I don't know what going to actually occur.
  • 40:31The best I can do is ask the people
  • 40:33who will implement the change.
  • 40:35- All right, well.
  • 40:37I'm sorry, this is my last remark.
  • 40:39I won't keep on doing this, but I would think that
  • 40:42these rates that we're talking about,
  • 40:43which seems to be are really critical to what happens
  • 40:45in the model, that you and find invasion inferency
  • 40:49you have to give a plausible, defensible probability
  • 40:52for them, which I would find hard to do,
  • 40:55and I also find it hard to do because I know that those
  • 40:58rates differ huge amount in Connecticut
  • 41:00between the different counties, that you can just see
  • 41:03if you look at what's happening in different counties.
  • 41:05Those rates are different because different
  • 41:10amount of separation and different amount
  • 41:12of personal contact.
  • 41:13- Sure.
  • 41:14- I think so kind of do that on an average way
  • 41:16of all the counties, seven or eight of them,
  • 41:19you'd think you at least got a vary among the counties
  • 41:23and have some number among the counties.
  • 41:25Then if there's a change of policy from the governor,
  • 41:27there'd be a change in sum or expected you need
  • 41:29to have that built in somehow here.
  • 41:31- Certainly.
  • 41:32In this work, I guess in all policy-relevant work,
  • 41:37there is a constant tension between the need
  • 41:41for parsimony and parametrization
  • 41:45and the need for these rich ways
  • 41:48of accommodating heterogeneity.
  • 41:52What we have found in this setting is that we lacked
  • 41:55the information or data to be able to separately
  • 41:58parametrize transmission rates at the county level
  • 42:04but that we can capture the aggregate number of cases,
  • 42:07hospitalizations and other relevant outcomes
  • 42:09at the state level by averaging over them.
  • 42:13The reason is because the counties themselves
  • 42:15have very different incidence, which actually does explain
  • 42:19quite a lot in the differing trajectories
  • 42:22of case counts and hospitalizations and deaths
  • 42:25within the counties.
  • 42:30- Hi Forrest thank you, this is very interesting.
  • 42:32This is Donna.
  • 42:33I have a question.
  • 42:35Do you have, the para--
  • 42:37Hi.
  • 42:38Are the parameters identifiable without Bayesian priors
  • 42:41or just from the data that we have or do you need
  • 42:46the priors in order to estimate the parameters?
  • 42:49- A subset of parameters is uniquely identifiable
  • 42:52by maximum likelihood or is point identified.
  • 42:55But really speaking, the answer to your question is no.
  • 42:59There are infinitely many combinations of parameters,
  • 43:04which fit any given loss function criteria equally well.
  • 43:07So we do need parameters here.
  • 43:09It is unfortunate and I think--
  • 43:12Yeah, go ahead.
  • 43:15- Priors you mean.
  • 43:17Do you know like what's the simplest possible model
  • 43:21that's just identifiable from the data
  • 43:23and is that model useful at all or is it so simple
  • 43:26that it's not even helpful?
  • 43:29- Two parts to that question, the simplest model
  • 43:31that is identifiable from the data is probably one in which
  • 43:34there is no heterogeneity in types of infection,
  • 43:39no asymptomatic infection.
  • 43:40We just lump all those people together
  • 43:42and there's only one kind of hospitalization
  • 43:45and people just transition, a certain proportion
  • 43:47of people transition to hospitalization.
  • 43:49That model is probably, has all the parameters identified.
  • 43:56And no, it's not useful.
  • 44:00That seems to be what we have found.
  • 44:03But I would say, I think there are two kinds
  • 44:05of usefulness, right.
  • 44:06One is answering the questions that policy makers have
  • 44:09and the other one is what Charles Manski calls credibility,
  • 44:12that there is a need to take into account
  • 44:16known heterogeneity and known mechanisms
  • 44:20when we construct these models.
  • 44:21So if I produce a useful projection that a policy maker
  • 44:24likes but I have not separated out asymptomatic infections,
  • 44:28then the numbers that I'm producing may become less,
  • 44:32regarded as less credible, right.
  • 44:34There's always this rhetorical function of modeling
  • 44:38beyond the numbers that are being produced,
  • 44:40to being able to accommodate or capture known mechanisms
  • 44:44by which data are generated
  • 44:48is one way that we can produce more believable
  • 44:50and actionable projections, right.
  • 44:53So I think there's this balance right,
  • 44:55between parsimony and richness and also this balance between
  • 45:03simplicity and believability of the assumptions.
  • 45:07So here we tried to you know, strike that balance.
  • 45:09If you think we've done it wrong, then please let us know.
  • 45:14- No, I definitely don't think you did it wrong,
  • 45:16but it would be interesting to see how much you lose
  • 45:20and sort of,
  • 45:22sort of cross validated predictability
  • 45:26by adding in priors, as opposed to just using the data
  • 45:30itself in a very simple model.
  • 45:32- Right, so--
  • 45:33- I don't know if you know the answer to that or not
  • 45:34but you should probably go on and I know other people
  • 45:37are wanting you to go on and not spend time answering
  • 45:40a lot of individual questions and we can always
  • 45:42talk another time.
  • 45:44- Okay, sounds good.
  • 45:46The model fits pretty well, fits observe data pretty well.
  • 45:49Here, I'm showing projections that start on March 1st,
  • 45:51rather than at the current day or any intermediate day,
  • 45:54just to emphasize that
  • 45:58model projections and uncertainty intervals here,
  • 46:00which are point-wise 95%, I call it--
  • 46:04They are not proper confidence intervals.
  • 46:07They're point-wise projections from draws,
  • 46:10using draws of parameters and initial conditions
  • 46:12from the posterior distribution over those quantities.
  • 46:15They're not confidence intervals in the strict sense.
  • 46:19But they do appear to
  • 46:22match observed data quite well.
  • 46:26So I think we're capturing dynamics that govern
  • 46:29what has occurred already.
  • 46:30We can learn quite a lot about the transmission rate
  • 46:34and under historical circumstances because we know
  • 46:37when those circumstances changed.
  • 46:39So we can estimate for example,
  • 46:41the percent decrease in transmission in Connecticut
  • 46:45following closure of schools and implementation
  • 46:48of the stay at home order.
  • 46:50That is what causes actually,
  • 46:51this downturn in hospitalizations and flattening
  • 46:54of cumulative deaths in the state.
  • 46:58So here, just to get a little bit more concrete,
  • 47:01on the upper left-hand corner, we see what we call
  • 47:05the contact intervention.
  • 47:06This is a function that multiples that transmission rate
  • 47:10parameter that we were discussing.
  • 47:12So in early March, schools are closed,
  • 47:15people start staying home
  • 47:16and so this intervention drops down.
  • 47:19The level to which it drops is a little more,
  • 47:23it drops more than 85%, I think, or somewhere around 85%.
  • 47:28That is an estimated quantity.
  • 47:30So the drops in historical contact are estimated
  • 47:35based on the changes in hospitalizations and deaths
  • 47:39and the implied changes in new infections.
  • 47:42Then what happens after the dotted line,
  • 47:44that is after May 20th, this is just a scenario
  • 47:49in which the amount of contact between individuals
  • 47:53increases at, I think here, monthly intervals by 10%
  • 47:57of the suppressed latent contact.
  • 48:00Under this historical and hypothetical future scenario,
  • 48:07we see cumulative incidence in the upper right-hand corner,
  • 48:11projected from March 1st onward.
  • 48:14Hospitalizations, with the dashed line,
  • 48:17showing expanded hospital capacity in Connecticut.
  • 48:27We see projections of deaths under this scenario,
  • 48:31cumulative incidence as a proportion of the population size
  • 48:34among people who are alive.
  • 48:35So this is what you would get if you conducted
  • 48:37a seroprevalence study in the future.
  • 48:39We hope this is useful for planning those types of studies,
  • 48:42and estimates of the affective reproduction number
  • 48:46in Connecticut over time.
  • 48:48There are two scenarios in particular that we want to show
  • 48:54policy makers that correspond to slow and fast reopening.
  • 48:57Really, this is not reopening scenarios.
  • 48:59I'm not sure what happened with this green annotation.
  • 49:03I don't know if you can see it.
  • 49:04If I did that or somebody else did, but just ignore that.
  • 49:09I'm not sure where it came from.
  • 49:11Under slow reopening, we imagine that people
  • 49:16release 10% of their latent suppressed contact
  • 49:19every month and under a scenario like this,
  • 49:22where everybody keeps distancing and everything goes
  • 49:24very well in the state, new infections continue their drop
  • 49:29and rise very slowly into the late summer and fall,
  • 49:32hospitalization stays low throughout the summer.
  • 49:39Deaths sort of begin to plateau and do not rise above
  • 49:4310,000 by the end of the summer.
  • 49:48Right, so this is the scenario that the state
  • 49:50is really hoping for.
  • 49:51It's a slow reopening that does not substantially increase
  • 49:55new infections with very slow rise
  • 49:59in new infections as the state reopens.
  • 50:02In contrast, a more pessimistic scenario,
  • 50:05which I think corresponds more to
  • 50:12a fast reopening, is one in which contact increases by 10%
  • 50:17or 10% of suppressed contact is released every two weeks.
  • 50:22This results in a very fast resurgence of new cases,
  • 50:26new hospitalizations and deaths by the end of the summer.
  • 50:29This is what the governor would like to avoid
  • 50:33when school children are scheduled
  • 50:34to go back to school in the fall.
  • 50:40There is a lot of interest right now in seroprevalence
  • 50:42because of competing claims about herd immunity
  • 50:45and how many people have been already infected
  • 50:47and have evidence of prior infection.
  • 50:50Under these scenarios, we can produce projections
  • 50:53of the proportion of people in a random sample
  • 50:56in the state, who might have evidence of prior infection.
  • 51:01So this is very important for designing seroprevalence
  • 51:04studies that we can use to further calibrate these models
  • 51:07and that can be used to guide policy.
  • 51:14I'm going to try to finish up very quickly here.
  • 51:17There are a couple of key messages from this work
  • 51:19that we tried to convey to policy makers.
  • 51:21The first is that the state is doing pretty well,
  • 51:23in terms of suppression of contact, closure of schools
  • 51:26and the stay at home order have effectively reduce
  • 51:28transmission and hospitalizations in Connecticut.
  • 51:32If contact increases quickly, the state's at serious risk
  • 51:37of big resurgence by later summer 2020.
  • 51:40Real time metrics that policy makers have access to
  • 51:42are really not going to serve as an early warning system
  • 51:46for that resurgence.
  • 51:49The state probably needs to be evaluating future projections
  • 51:53under realistic contact scenarios for the state.
  • 51:57We still have a lot of uncertainty that we tried to capture
  • 52:01in model projections about cumulative incidence,
  • 52:04asymptomatic fraction,
  • 52:06how things are going to go with children,
  • 52:10the effects of enhanced testing and contact tracing
  • 52:13and how contact patterns may change following reopening.
  • 52:18So we are issuing a series of reports, which you can read
  • 52:22online and we will be updating them in real time
  • 52:27as the summer goes on.
  • 52:28You can find them at this URL.
  • 52:29You can also email me and I'll point you to them.
  • 52:34These are sort of continuously updated research products
  • 52:36and I hope that they will represent
  • 52:38the latest information from Connecticut
  • 52:40and our latest predictions for the state as it reopens.
  • 52:45Also, there's a document here which summarizes
  • 52:47much more detail about the transmission model
  • 52:50that I have given here in this presentation.
  • 52:53I'm gonna skip over this stuff about our workflow.
  • 52:56We can talk about it later, if anybody is interested,
  • 52:59but this is just how we transition from regular research
  • 53:03to doing this type of very active software development.
  • 53:07I will end here.
  • 53:08I want to thank all of the people in the group
  • 53:11and beyond, who have been working on this tirelessly
  • 53:13over the last couple of months.
  • 53:15All of the products that I've told you about
  • 53:17are publicly available.
  • 53:18You can find the source code on Git
  • 53:21on our Git repositories and you can find the web application
  • 53:24and the reports online as well.
  • 53:27So I'd be happy to take any questions.
  • 53:29- Thanks, thanks Forrest for the last part.
  • 53:32I think some people have some questions using the chat box.
  • 53:39Ken asked, "Is the model used at currently proposing
  • 53:43"used at hospital or by your medical group?"
  • 53:50- The ICU planning app
  • 53:53has been used, we know, and possibly is being used
  • 53:56at Yale New Haven Hospital.
  • 53:58The projections for Connecticut are not intended for use
  • 54:01in any particular hospital systems,
  • 54:03though I think they will be of interest
  • 54:05to leaders of systems who are planning to accommodate
  • 54:10a potential second wave of infections
  • 54:12as it might occur later in the summer.
  • 54:15I hope that as we get farther in the summer,
  • 54:18if there is a second wave that appears to be coming,
  • 54:21that the projections will be useful in planning
  • 54:23capacity expansion efforts, possibly at or beyond levels
  • 54:27that we already saw in April.
  • 54:31So we will be generating any information
  • 54:33that decision makers at those hospital systems
  • 54:37think would be useful as they plan their response.
  • 54:40That's a great question.
  • 54:42- Thanks.
  • 54:43And...
  • 54:49Let me see and Sherry asked,
  • 54:55"In the first reopening model, what amount was the reopening
  • 54:58"assumed to start in?"
  • 55:01- Exactly on May 20th,
  • 55:03which is when the governor began the process of reopening.
  • 55:08It is also true that the governor has been
  • 55:13giving information about potential reopening plans
  • 55:15for a very long time
  • 55:17and that there is some change in contact as people
  • 55:19begin to anticipate those changes in policy.
  • 55:23I think that if you are looking at human mobility data
  • 55:28from cell phones and other sources,
  • 55:31you will see that people have been moving around
  • 55:33for a while, increasing their level of activity
  • 55:36outside of the home, even before May 20th in Connecticut.
  • 55:41Whether that has actually resulted
  • 55:46in a substantial increase in transmission remains to be seen
  • 55:50but I don't think we should assume that just because
  • 55:52people are moving around and possibly returning
  • 55:54to some types of work that there will be a corresponding
  • 55:58increase in transmission.
  • 56:02- Okay thanks.
  • 56:03Daniel asks, "Is the increase in incidence starting
  • 56:07"in September a cumulative effect of prolonged increase
  • 56:14"in contact."
  • 56:15- Can I just ask the question directly?
  • 56:16So I'm wondering, in the parts where you're showing
  • 56:18the two reopening models, it looked like the curve
  • 56:21starts to go back up around August,
  • 56:23September in the slow one.
  • 56:25I'm wondering if that's because you reach a threshold
  • 56:28above a certain percentage of contact
  • 56:30or if it's a cumulative effect?
  • 56:32Like, if we were to keep contact at .2 for example,
  • 56:36throughout all of this time and it weren't to increase
  • 56:39above a threshold, is there a situation which you don't see
  • 56:42that tail come up again?
  • 56:45- Yes, great question.
  • 56:47If you like to think in terms of the effective
  • 56:49reproduction number, this increase just corresponds
  • 56:52to a time about three weeks after
  • 56:55that number goes above one.
  • 56:58So there is a threshold effect and to answer your question,
  • 57:00if contact were to remain below a level
  • 57:06that would give that value of one,
  • 57:08then you would not see this type of resurgence.
  • 57:11I think as a practical matter, it is very unlikely
  • 57:14that the state can avoid a situation where the effective
  • 57:17reproduction number does above one.
  • 57:21I think this is not the stated strategy of anyone
  • 57:24and it's probably not, but I think it is the realistic
  • 57:29expectation about what will happen in reality.
  • 57:32The reality is that the state is going to try very hard
  • 57:35to increase a level of contact just about to that level,
  • 57:39where they would see some local outbreaks
  • 57:41that can be extinguished but they will try to maximize
  • 57:45the level of contact, meaning economic activity
  • 57:49and social mobility
  • 57:53that the state can achieve.
  • 57:54So they'll try to get as much economic productivity
  • 57:57and contact as they can without causing resurgence
  • 58:03or large outbreak or an overrun of hospital capacity.
  • 58:08- Thank you.
  • 58:09- Thanks.
  • 58:12- Akil here have two questions.
  • 58:15So the first one is are there any assumptions
  • 58:17of the proposed population who have Covid-19
  • 58:20but have not been tested?
  • 58:24- There are implicit and explicit assumptions
  • 58:26about that proportion.
  • 58:28I think we can produce predictions
  • 58:31for the current prevalence and also cumulative incidence
  • 58:37but those predictions depend quite a lot on our prior
  • 58:39assumptions about the asymptomatic faction.
  • 58:44We don't have very precise information about how many
  • 58:48or what proportion of infections are totally asymptomatic
  • 58:51and would go undetected by the healthcare system
  • 58:54because people don't seek testing or seek care of any kind
  • 58:58when they're not feeling sick.
  • 59:02So certainly, we can try to learn about those things.
  • 59:03There's some information in the available case counts
  • 59:07and in hospitalizations and deaths about that stuff,
  • 59:12but we still have a lot of uncertainty about current
  • 59:15cumulative incidence.
  • 59:17I think it's fair to say that currently prevalence
  • 59:18is quite low in Connecticut.
  • 59:22- Okay thanks.
  • 59:22I guess I saw something new saying they test the people
  • 59:26(unclear speaking).
  • 59:29Because they can test other people that have the ability
  • 59:32and then they have some estimate of the asymptomatic case,
  • 59:37the rate of them?
  • 59:38- Yes, that's true.
  • 59:41In some very specific settings, like institutional settings
  • 59:44like nursing homes and correctional institutions,
  • 59:47you can test everybody and then you can learn how many
  • 59:51infections are asymptomatic.
  • 59:53The question then becomes of how representative
  • 59:57those samples are compared to the rest of the state.
  • 01:00:01Is it safe to take situations where people
  • 01:00:05are living in very close proximity
  • 01:00:08and possibly poor health conditions and to generalize
  • 01:00:13all of that information to the state?
  • 01:00:15I think there is some very good anecdotal evidence
  • 01:00:17from prisons, from nursing homes
  • 01:00:18and also testing systematic testing of healthcare workers
  • 01:00:23that we can try to take into account,
  • 01:00:26but it remains unclear how generalizable
  • 01:00:29that information is.
  • 01:00:30For example, healthcare workers may be immunologically
  • 01:00:33somewhat unlike members of the general population
  • 01:00:36who are not continuously exposed to different types
  • 01:00:40of illness and to coronaviruses in particular.
  • 01:00:43So I would hesitate to take large screening studies
  • 01:00:45of nurses for example, and apply the asymptomatic fraction
  • 01:00:51or prevalence or incidence in that sample
  • 01:00:53to the general population.
  • 01:00:56- Thanks.
  • 01:00:58And the second question of Akil is can Covid-19 models
  • 01:01:03from different states learn from each other?
  • 01:01:06I have relay the question is because currently your model
  • 01:01:10is most about stating the data
  • 01:01:11and you can validate how good the model is.
  • 01:01:15Because states, maybe they have their reopening plan
  • 01:01:18at different times, can this provide useful information
  • 01:01:21about how good the model is by learning
  • 01:01:24from different states.
  • 01:01:25- Yes, great question.
  • 01:01:27It is always true that information from other contexts
  • 01:01:30can be very useful if you know what is different
  • 01:01:33in those other contexts.
  • 01:01:35I would love to be able to use more granular information
  • 01:01:37from neighboring states throughout the northeast
  • 01:01:40to inform projections from Connecticut,
  • 01:01:42'cause as we know, Connecticut is not an island
  • 01:01:45and as soon as New York opens up and people start working
  • 01:01:47in New York, then everything will change quite a lot,
  • 01:01:50quite quickly in Connecticut.
  • 01:01:53So I would like to share information.
  • 01:01:56We have focused on Connecticut here because we have very
  • 01:01:58detailed information about Connecticut but no special access
  • 01:02:02in Massachusetts, Rhode Island and New York.
  • 01:02:07So that's why we've done it, but I think it will become
  • 01:02:09very important and I always thought
  • 01:02:13it would be the job of the CDC and the US
  • 01:02:16to synthesize a national and local projections
  • 01:02:18and to gather all the granular local information
  • 01:02:21and to put it all together.
  • 01:02:22That has not happened in this particular pandemic.
  • 01:02:26So I think everyone else is trying to scramble
  • 01:02:31to aggregate information at the right levels
  • 01:02:33to produce predictions that are actionable locally.
  • 01:02:37But there's not coordination right now
  • 01:02:39between groups that are doing state-specific
  • 01:02:42reopening plans, unfortunately.
  • 01:02:44As for whether the differences or staggered reopening
  • 01:02:48can be used as a kind of instrument to identify
  • 01:02:50the causal effects of reopening, I assume that's the subtext
  • 01:02:54of the question, the answer is yes.
  • 01:02:57I think people are very interested in doing that.
  • 01:02:59The problem is that reopening is somewhat endogenous.
  • 01:03:03The states to reopening as a function of the conditions
  • 01:03:06currently in the states and also obviously,
  • 01:03:09as a function of the political considerations
  • 01:03:12of the leadership and of the population.
  • 01:03:16Right now I don't think it's safe to say that reopening
  • 01:03:19occurs randomly in some time interval
  • 01:03:22and that we can exploit that randomness in a simple way
  • 01:03:24to assess the effect of reopening.
  • 01:03:27Certainly, some of the states that we observe
  • 01:03:29reopening quickly, take Georgia for example.
  • 01:03:34Those states are likely to see at least local
  • 01:03:36and possibly very broad resurgences and outbreaks
  • 01:03:42that may result in reversion to more restrictive movement
  • 01:03:45conditions in those states.
  • 01:03:48So I think really, there's this going to be a long,
  • 01:03:50longitudinal sequence of treatments,
  • 01:03:53meaning changes in state regulations and then outcomes,
  • 01:03:57which the regulators will observe
  • 01:04:00and then this kink of cat and mouse game,
  • 01:04:01where decision makers try to tamp down on local outbreaks
  • 01:04:06and then respond to ones that occur in the future.
  • 01:04:11So we will try to learn about the effects of all those
  • 01:04:13interventions and changes in policies
  • 01:04:16but I think that there is cause for some skepticism
  • 01:04:20in really learning a generalizable causable effects
  • 01:04:24just from the time series.
  • 01:04:27- Thanks.
  • 01:04:29I guess one last very specific question about a talk.
  • 01:04:33So Paul asked, "Have you considered how real time data
  • 01:04:36"metrics, such as oxygen sensors from fitness trackers
  • 01:04:39"could effect your predictions?"
  • 01:04:41- Very interested in distributed measurements
  • 01:04:44at the population level that could be helpful
  • 01:04:46to inform some of these things.
  • 01:04:49I think that we have not yet seen widespread adoption
  • 01:04:53of mobile apps
  • 01:04:58for self monitoring for contact tracing.
  • 01:05:04There is some adoption of thermometers and oxygen sensors
  • 01:05:10but as far as I know, there are no data streams
  • 01:05:12that are publicly available.
  • 01:05:14- This is Paul Forcher, I asked the question.
  • 01:05:18There are some, there's--
  • 01:05:20I'm participating in two studies.
  • 01:05:22One that's run out of by Mike Snider,
  • 01:05:26who use to be at Yale who's head of Stanford Genomics.
  • 01:05:31The other one's institute
  • 01:05:34and any of you can sign up for these things
  • 01:05:37and if you have a fitness tracker that's tracking
  • 01:05:41oxygen levels, there's emerging evidence that changing
  • 01:05:46oxygen levels can be predictive of Covid infection
  • 01:05:49before the patients are symptomatic and there's some...
  • 01:05:54So I would, those are two studies that you could
  • 01:05:57connect with and I wouldn't be surprised at all
  • 01:05:59if they would share all of their realtime data
  • 01:06:01that they're collecting with you.
  • 01:06:02- Yeah, that is a great idea, thank you.
  • 01:06:04With these--
  • 01:06:05- Mike Snider's a former Yale person,
  • 01:06:08so you already have an inroad with that guy.
  • 01:06:12- Yeah, thank you, that's a great idea.
  • 01:06:15- Okay thanks, I guess that's all questions for the talk.
  • 01:06:20If you have any questions, I guess they can talk to you.
  • 01:06:24Like the audience can talk to Forrest offline.
  • 01:06:28- Please feel free to email me, anybody who has questions.
  • 01:06:30- Some people want to hear more about the talks,
  • 01:06:34like you didn't have time to cover,
  • 01:06:35that I guess the interest you can talk to Forrest offline.
  • 01:06:40Also, this talk will be recorded
  • 01:06:41and will be publicly available.
  • 01:06:44Also, on the previous talk are also recorded.
  • 01:06:47I'll also send out a link to everyone
  • 01:06:51in the School of Public Health,
  • 01:06:54so if you want you can access it.
  • 01:06:56Okay thank, thanks Forrest.
  • 01:06:58- Thanks everyone.
  • 01:06:59- And thanks for everyone. - Thanks everyone.