Stochastic Models for the Inference of Life Evolution


SMILE is an interdisciplinary research group gathering probabilists, statisticians, bio-informaticians and biologists.
SMILE is affiliated to the Stochastics and Biology group of LPSM (Lab of Probability, Statistics and Modeling) at Sorbonne Université (ex Université Pierre et Marie Curie Paris 06).
SMILE is hosted within the CIRB (Center for Interdisciplinary Research in Biology) at Collège de France.
SMILE is supported by Collège de France and CNRS.
Visit also our homepage at CIRB.

Recent contributions of the SMILE group related to SARS-Cov2 and COVID-19.


SMILE is hosted at Collège de France in the Latin Quarter of Paris. To reach us, go to 11 place Marcelin Berthelot (stations Luxembourg or Saint-Michel on RER B).
Our working spaces are rooms 107, 121 and 122 on first floor of building B1 (ask us for the code). Building B1 is facing you upon exiting the traversing hall behind Champollion's statue.


You can reach us by email (amaury.lambert - at - or (smile - at -

Light on



The genealogical decomposition of a matrix population model with applications to the aggregation of stages

Matrix projection models are a central tool in many areas of population biology. In most applications, one starts from the projection matrix to quantify the asymptotic growth rate of the population (the dominant eigenvalue), the stable stage distribution, and the reproductive values (the dominant right and left eigenvectors, respectively). Any primitive projection matrix also has an associated ergodic Markov chain that contains information about the genealogy of the population. In this paper, we show that these facts can be used to specify any matrix population model as a triple consisting of the ergodic Markov matrix, the dominant eigenvalue and one of the corresponding eigenvectors. This decomposition of the projection matrix separates properties associated with lineages from those associated with individuals. It also clarifies the relationships between many quantities commonly used to describe such models, including the relationship between eigenvalue sensitivities and elasticities. We illustrate the utility of such a decomposition by introducing a new method for aggregating classes in a matrix population models to produce a simpler model with a smaller number of classes. Unlike the standard method, our method has the advantage of preserving reproductive values and elasticities. It also has conceptually satisfying properties such as commuting with changes of units.



Coagulation-transport equations and the nested coalescents

The nested Kingman coalescent describes the dynamics of particles (called genes) contained in larger components (called species), where pairs of species coalesce at constant rate and pairs of genes coalesce at constant rate provided they lie within the same species. We prove that starting from \$$rn\$$ species, the empirical distribution of species masses (numbers of genes\$$/n\$$) at time \$$t/n\$$ converges as \$$n\to\infty\$$ to a solution of the deterministic coagulation-transport equation $$ \partial_t d \ = \ \partial_x ( \psi d ) \ + \ a(t)\left(d\star d - d \right), $$ where \$$\psi(x) = cx^2\$$, \$$\star\$$ denotes convolution and \$$a(t)= 1/(t+\delta)\$$ with \$$\delta=2/r\$$. The most interesting case when \$$\delta =0\$$ corresponds to an infinite initial number of species. This equation describes the evolution of the distribution of species of mass \$$x\$$, where pairs of species can coalesce and each species' mass evolves like \$$\dot x = -\psi(x)\$$. We provide two natural probabilistic solutions of the latter IPDE and address in detail the case when \$$\delta=0\$$. The first solution is expressed in terms of a branching particle system where particles carry masses behaving as independent continuous-state branching processes. The second one is the law of the solution to the following McKean-Vlasov equation $$ dx_t \ = \ - \psi(x_t) \,dt \ + \ v_t\,\Delta J_t $$ where \$$J\$$ is an inhomogeneous Poisson process with rate \$$1/(t+\delta)\$$ and \$$(v_t; t\geq0)\$$ is a sequence of independent rvs such that \$$\mathcal L(v_t) = \mathcal L(x_t)\$$. We show that there is a unique solution to this equation and we construct this solution with the help of a marked Brownian coalescent point process. When \$$\psi(x)=x^\gamma\$$, we show the existence of a self-similar solution for the PDE which relates when \$$\gamma=2\$$ to the speed of coming down from infinity of the nested Kingman coalescent.

Upcoming seminars


Planning des salles du Collège de France.
Intranet du Collège de France.