The Human papillomaviruses (HPV) vaccine induces a herd immunity effect in genital warts when a large number of the population is vaccinated. This aspect should be taken into account when devising new vaccine strategies, like vaccination at older ages or male vaccination. Therefore, it is important to develop mathematical models with good predictive capacities. We devised a sexual contact network that was calibrated to simulate the Spanish epidemiology of different HPV genotypes. Through this model, we simulated the scenario that occurred in Australia in 2007, where 12–13 year-old girls were vaccinated with a three-dose schedule of a vaccine containing genotypes 6 and 11, which protect against genital warts, and also a catch-up program in women up to 26 years of age. Vaccine coverage were 73% in girls with three doses and with coverage rates decreasing with age until 52% for 20–26 year-olds. A fast 59% reduction in the genital warts diagnoses occurred in the model in the first years after the start of the program, similar to what was described in the literature.

a b s t r a c tSeasonal fluctuations in the incidence of several respiratory infections are a feature of epidemiological surveys all around the world. This phenomenon is characteristic of influenza and respiratory syncytial virus pandemics. However, the explanation of the seasonal outbreaks of these diseases remains poorly understood. Many statistical studies have been carried out in order to provide a correlation of the outbreaks with climatic or social factors without achieving a definitive conclusion. Here we show that, in a random social network, self-sustained seasonal epidemics emerge as a process modulated by the infection probability and the immunity period after recovering from the infection. This is a purely endogenous phenomenon that does not require any exogenous forcing. Assuming that this is the dominant mechanism for seasonal epidemics, many implications for public health policies for infectious respiratory diseases could be drawn.

Working in large networks applied to epidemiological-type models has led us to design a simple but effective computed distributed environment to perform a large amount of model simulations in a reasonable time in order to study the behavior of these models and to calibrate them. Finding the model parameters that best fit the available data in the designed distributed computing environment becomes a challenge and it is necessary to implement reliable algorithms for model calibration. In this article, we have adapted the random particle swarm optimization algorithm to our distributed computing environment to be applied to the calibration of a papillomavirus transmission dynamics model on a lifetime sexual partners network. And we have obtained a good fitting saving time and calculations compared with the exhaustive searching strategy we have been using so far.

No abstract

Recently, the transmission dynamics of the Human Papillomavirus (HPV) has been studied. In previous works, we have designed and implemented a computational model (agent-based simulation model) where the contagion of the HPV is described on a network of lifetime sexual partners. The run of a single simulation of this computational model, composed of a network with 500 000 nodes, takes about one hour and a half. In addition to set an adequate model, finding out the model parameters that best fit the proposed model to the available data of prevalence is a crucial goal. Taking into account that the necessary number of simulations to perform the calibration of the model may be very high, the aforementioned goal may become unaffordable. In this paper, we present a procedure to fit the proposed HPV model to the available data and the design of an asynchronous version of the Particle Swarm Optimization (PSO) algorithm adapted to the distributed computing environment. In the process, the number of particles used in PSO should be set carefully looking for a compromise between quality of the solutions and computation time. Another feature of the procedure presented here is that we want to capture the intrinsic uncertainty in the data (data come from a survey) when calibrating the model. To do so, we also propose the design of an algorithm to select the model parameter sets obtained during the calibration that best capture the data uncertainty.