On the exact and population bi-dimensional reproduction numbers in a stochastic SVIR model with imperfect vaccine

We aim to quantify the spread of a direct contact infectious disease that confers permanent immunity after recovery, within a non-isolated finite and homogeneous population. Prior to the onset of the infection and to prevent the spread of this disease, a proportion of individuals was vaccinated. But the administered vaccine is imperfect and can fail, which implies that some vaccinated individuals get the infection when being in contact with infectious individuals. We study the evolution of the epidemic process over time in terms of a continuous-time Markov chain, which represents a general SIR model with an additional compartment for vaccinated individuals. In our stochastic framework, we study two bi-dimensional variables recording infection events, produced by a single infectious individual or by the whole infected group, taking into account if the newly infected individual was previously vaccinated or not. Theoretical schemes and recursive algorithms are derived in order to compute joint probability mass functions and factorial moments for these random variables. We illustrate the applicability of our techniques by means of a set of numerical experiments.