  Birth-death processes track the size of a univariate population, but many
biological systems involve interaction between populations, necessitating
models for two or more populations simultaneously. A lack of efficient methods
for evaluating finite-time transition probabilities of bivariate processes,
however, has restricted statistical inference in these models. Researchers rely
on computationally expensive methods such as matrix exponentiation or Monte
Carlo approximation, restricting likelihood-based inference to small systems,
or indirect methods such as approximate Bayesian computation. In this paper, we
introduce the birth(death)/birth-death process, a tractable bivariate extension
of the birth-death process. We develop an efficient and robust algorithm to
calculate the transition probabilities of birth(death)/birth-death processes
using a continued fraction representation of their Laplace transforms. Next, we
identify several exemplary models arising in molecular epidemiology,
macro-parasite evolution, and infectious disease modeling that fall within this
class, and demonstrate advantages of our proposed method over existing
approaches to inference in these models. Notably, the ubiquitous stochastic
susceptible-infectious-removed (SIR) model falls within this class, and we
emphasize that computable transition probabilities newly enable direct
inference of parameters in the SIR model. We also propose a very fast method
for approximating the transition probabilities under the SIR model via a novel
branching process simplification, and compare it to the continued fraction
representation method with application to the 17th century plague in Eyam.
Although the two methods produce similar maximum a posteriori estimates, the
branching process approximation fails to capture the correlation structure in
the joint posterior distribution.
