Stochastic models of kleptoparasitism

In this paper, we consider a model of kleptoparasitism amongst a small group of individuals, where the state of the population is described by the distribution of its individuals over three specific types of behaviour (handling, searching for or fighting over, food). The model used is based upon earlier work which considered an equivalent deterministic model relating to large, effectively infinite, populations. We find explicit equations for the probability of the population being in each state. For any reasonably sized population, the number of possible states, and hence the number of equations, is large. These equations are used to find a set of equations for the means, variances, covariances and higher moments for the number of individuals performing each type of behaviour. Given the fixed population size, there are five moments of order one or two (two means, two variances and a covariance). A normal approximation is used to find a set of equations for these five principal moments. The results of our model are then analysed numerically, with the exact solutions, the normal approximation and the deterministic infinite population model compared. It is found that the original deterministic models approximate the stochastic model well in most situations, but that the normal approximations are better, proving to be good approximations to the exact distribution, which can greatly reduce computing time.