Phase transition on the convergence rate of parameter estimation under an Ornstein-Uhlenbeck diffusion on a tree

DOI10.1007/S00285-016-1029-XMaRDI QIDQ504087zbMATH OpenOpenAlexWikidataFDO

Authors Sebastien Roch, Cécile Ané, Lam Si Tung Ho

Publication date 25 January 2017

Published in Journal of Mathematical Biology (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1406.1568

zbMATH Keywords

consistencyphase transitionmaximum likelihood estimatorevolutionOrnstein-Uhlenbeckphylogenetic

Mathematics Subject Classification ID

Asymptotic properties of parametric estimators (62F12) Applications of statistics to biology and medical sciences; meta analysis (62P10) Non-Markovian processes: estimation (62M09) Problems related to evolution (92D15) Diffusion processes (60J60) Taxonomy, cladistics, statistics in mathematical biology (92B10)

Abstract: Diffusion processes on trees are commonly used in evolutionary biology to model the joint distribution of continuous traits, such as body mass, across species. Estimating the parameters of such processes from tip values presents challenges because of the intrinsic correlation between the observations produced by the shared evolutionary history, thus violating the standard independence assumption of large-sample theory. For instance Ho and An'e cite{HoAne13} recently proved that the mean (also known in this context as selection optimum) of an Ornstein-Uhlenbeck process on a tree cannot be estimated consistently from an increasing number of tip observations if the tree height is bounded. Here, using a fruitful connection to the so-called reconstruction problem in probability theory, we study the convergence rate of parameter estimation in the unbounded height case. For the mean of the process, we provide a necessary and sufficient condition for the consistency of the maximum likelihood estimator (MLE) and establish a phase transition on its convergence rate in terms of the growth of the tree. In particular we show that a loss of $sqrtn$-consistency (i.e., the variance of the MLE becomes $Omega(n^{-1})$, where $n$ is the number of tips) occurs when the tree growth is larger than a threshold related to the phase transition of the reconstruction problem. For the covariance parameters, we give a novel, efficient estimation method which achieves $sqrtn$-consistency under natural assumptions on the tree.

Recommendations

• Asymptotic theory with hierarchical autocorrelation: Ornstein-Uhlenbeck tree models
• On the convergence of the maximum likelihood estimator for the transition rate under a 2-state symmetric model
• Consistency of a phylogenetic tree maximum likelihood estimator
• Consistency and convergence rate of phylogenetic inference via regularization
• Phylogenetic confidence intervals for the optimal trait value

Cites work

• scientific article; zbMATH DE number 3643026 (Why is no real title available?)
• scientific article; zbMATH DE number 4062374 (Why is no real title available?)
• scientific article; zbMATH DE number 1865935 (Why is no real title available?)
• scientific article; zbMATH DE number 1418385 (Why is no real title available?)
• A phylogenetic comparative method for studying multivariate adaptation
• Asymptotic theory with hierarchical autocorrelation: Ornstein-Uhlenbeck tree models
• Broadcasting on trees and the Ising model.
• CLT for Ornstein-Uhlenbeck branching particle system
• Functionals of critical multitype branching processes
• Majority rule has transition ratio 4 on Yule trees under a 2-state symmetric model
• Mathematical Statistics
• Phylogenetic confidence intervals for the optimal trait value
• Robust Estimation of Latent Tree Graphical Models: Inferring Hidden States With Inexact Parameters
• \(U\)-statistics of Ornstein-Uhlenbeck branching particle system

Cited in

(11)

• Asymptotic theory with hierarchical autocorrelation: Ornstein-Uhlenbeck tree models
• Phylogenetic confidence intervals for the optimal trait value
• When can we reconstruct the ancestral state? Beyond Brownian motion
• When can we reconstruct the ancestral state? A unified theory
• A central limit theorem for punctuated equilibrium
• On the convergence of the maximum likelihood estimator for the transition rate under a 2-state symmetric model
• Normal approximation for mixtures of normal distributions and the evolution of phenotypic traits
• On the convergence order of a binary tree approximation of symmetrized diffusion processes
• A parametric interpretation of Bayesian nonparametric inference from gene genealogies: linking ecological, population genetics and evolutionary processes
• Species tree estimation under joint modeling of coalescence and duplication: sample complexity of quartet methods
• Using the Ornstein-Uhlenbeck process to model the evolution of interacting populations

Describes a project that uses

Uses Software

• ape
• GEIGER

This page was built for publication: Phase transition on the convergence rate of parameter estimation under an Ornstein-Uhlenbeck diffusion on a tree

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q504087)