The annual changing patterns of viral population sizes, defined as the effective number of infections, were estimated by Bayesian skyride analysis (Minin et al., 2008), which uses a piecewise demographic model with a Gaussian Markov random field (GMRF) smoothing prior.