A pathogen can infect multiple hosts. For example, zoonotic diseases like rabies often colonize both humans and animals. Meanwhile, a single host can sometimes be infected with many pathogens, such as malaria and meningitis. Therefore, we studied two susceptible classes S 1 ( t ) and S 2 ( t ) , each of which can be infected when interacting with two different infectious groups I 1 ( t ) and I 2 ( t ) . The stochastic models were formulated through the continuous time Markov chain (CTMC) along with their deterministic analogues. The statistics for the developed model were studied using the multi-type branching process. Since each epidemic class was assumed to transmit only its own type of pathogen, two reproduction numbers were obtained, in addition to the probability-generating functions of offspring. Thus, these, together with the mean number of infections, were used to estimate the probability of extinction. The initial population of infectious classes can influence their probability of extinction. Understanding the disease extinctions and outbreaks could result in rapid intervention by the management for effective control measures.
Keywords: Gillespie algorithm; basic reproduction number; branching process; continuous time Markov chain; epidemic extinction; stochastic differential equation.