Nonparametric estimation of the division rate of a size-structured population
arXiv:1103.4226 · doi:10.1137/110828344
Abstract
We consider the problem of estimating the division rate of a size-structured population in a nonparametric setting. The size of the system evolves according to a transport-fragmentation equation: each individual grows with a given transport rate, and splits into two offsprings of the same size, following a binary fragmentation process with unknown division rate that depends on its size. In contrast to a deterministic inverse problem approach, as in (Perthame, Zubelli, 2007) and (Doumic, Perthame, Zubelli, 2009), we take in this paper the perspective of statistical inference: our data consists in a large sample of the size of individuals, when the evolution of the system is close to its time-asymptotic behavior, so that it can be related to the eigenproblem of the considered transport-fragmentation equation (see \cite{PR} for instance). By estimating statistically each term of the eigenvalue problem and by suitably inverting a certain linear operator (see previously quoted articles), we are able to construct a more realistic estimator of the division rate that achieves the same optimal error bound as in related deterministic inverse problems. Our procedure relies on kernel methods with automatic bandwidth selection. It is inspired by model selection and recent results of Goldenschluger and Lepski.
References in corpus (4)
Cited by in corpus (10)
- Estimator selection: a new method with applications to kernel density estimation
- Nonparametric estimation for interacting particle systems : McKean-Vlasov models
- A General Inverse Problem for the Growth-Fragmentation Equation
- Estimating the division rate and kernel in the fragmentation equation
- A purely mechanical model with asymmetric features for early morphogenesis of rod-shaped bacteria micro-colony
- Nonparametric estimation of the fragmentation kernel based on a PDE stationary distribution approximation
- Cell lineage statistics with incomplete population trees
- Minimal penalty for Goldenshluger-Lepski method
- Adaptive nonparametric estimation of a component density in a two-class mixture model
- A convergent numerical algorithm for the stochastic growth-fragmentation problem