Spread of the plague in Venice, 1630-1631: epidemic entropy in a "natural experiment"

Precise modeling of epidemic spread is difficult. One explanation is that disease spread is inherently stochastic. This would suggest that the distribution of cases across geographic regions would progress towards that more favored by chance. If the epidemic proceeds long enough, the allocation of cases could approach that most expected, maximizing Boltzmann-Gibbs-Shannon entropy. Here, we tested these hypotheses on mortality data from the Venetian 1630-1631 plague epidemic. Entropy per case (intensive) of the quantile function (distribution of parishes ranked by case rates) increased from an effective number of 7.32 parishes (95% CI 3.32-12.55 parishes) to 47.9 parishes (47.5-48.9 parishes) out of 50 total, indicating that the quantile function approached a uniform maximum entropy distribution. Intensive entropy of the probability density function (parishes categorized by cumulative case rate) increased from 0.63 nats (0.32-0.93 nats) to 1.75 nats (1.53-1.87 nats). The PDF approached a Gaussian distribution. The Kullback-Leibler divergence decreased from 0.84 nats (0.71-1.42 nats) to 0.12 nats (0.083-0.35 nats). These findings quantify how disease spreads and demonstrate that observed heterogeneity in infections between regions may in some circumstances be explained by chance alone.