Two multi-sigmoidal diffusion models for the study of the evolution of the covid-19 pandemic

A proposal is made to employ stochastic models, based on diffusion processes, to represent the evolution of the SARS-CoV-2 virus pandemic. Specifically, two diffusion processes are proposed whose mean functions obey multi-sigmoidal Gompertz and Weibull-type patterns. Both are con-structed by introducing polynomial functions in the ordinary differential equations that originate the classical Gompertz and Weibull curves. The estimation of the parameters is approached by maximum likelihood. Various associated problems are analyzed, such as the determination of initial solutions for the necessary numerical methods in practical cases, as well as Bayesian methods to determine the degree of the polynomial. Additionally, strategies are suggested to determine the best model to fit specific data. A practical case is developed from data originating from several Spanish regions during the first two waves of the COVID-19 pandemic. The determination of the inflection time instants, which correspond to the peaks of infection and deaths, is given special attention. To deal with this particular issue, point estimation as well as first-passage times have been considered.