Simulation of stochastic models of structured population in population genetics under neutrality

Abstract

This paper describes some population genetic models under neutrality, involving genetic drift and mutations. Starting with Kingman’s coalescent we show how structured populations can be modeled. We detail these models by showing how simulation algorithms can be written. In particular we highlight the latent processes than rule out the explicit computation of the likelihood function on a dataset.