Kalkatungu Country 05 by Cheryl Perez. Used with permission.
In this post, I’m going to share the first model I built. This one imagines the campus where I work. At the time, there was still a hope that in-person classes would resume and that everything would go back to normal. Colleges and universities wanted to get students back on campus for pedagogical as well as financial reasons.
All of this made me nervous. When people are desperate to open, it makes it difficult to assess interpret information well. We cease to think rationally so we’re not making high quality decisions. I was nervous that desperate schools and communities would open prematurely and that cases would spiral out of control.
So built this model to visualize what could happen. It assumes a lot of things:
A campus of 1200 students and 100 faculty (no administrators, part time instructors, or other staff)
Faculty teaching 3 classes that have between 5 and 40 students, and a mode of 15 students
Each day consists of 6 teaching periods
No lunch or extracurricular meetings
At the end of the day students go to single, double, and triple sized dorm suites
Obviously there are a lot of simplifications here. And I have no data to support these numbers, but they seemed okay-ish.
I assumed that students that weren’t in class would leave campus with probability 5% or that they would meet up with another random student on campus with probability 5%. No data for this.
In terms of infections, I assumed that the campus had no cases and that all cases started somewhere in the community. I assumed that the community had 223 cases per 100k. Believe it or not, that was about where it was in April, 2020. Once infected, I assumed that a person would infect a contact with probability 5%. Again, just made this up. I also assumed instantaneous incubation. At the time, there was a lot of uncertainty around how long it took for an infected person to become symptomatic and to infectious to others.
This little exercise was an opportunity to play around with some R packages that I hadn’t used in a long while. Specifically, the igraph and ndtv packages. The networkDynamic documentation was super helpful as was Katya Ognyanova’s tutorial.
I ended up simulating a period of 3 days, which turned out like this:
(Please be patient, it may take a minute to load)