Analysis of comparative data with hierarchical autocorrelation
arXiv:0804.3166 · doi:10.1214/08-AOAS173
Abstract
The asymptotic behavior of estimates and information criteria in linear models are studied in the context of hierarchically correlated sampling units. The work is motivated by biological data collected on species where autocorrelation is based on the species' genealogical tree. Hierarchical autocorrelation is also found in many other kinds of data, such as from microarray experiments or human languages. Similar correlation also arises in ANOVA models with nested effects. I show that the best linear unbiased estimators are almost surely convergent but may not be consistent for some parameters such as the intercept and lineage effects, in the context of Brownian motion evolution on the genealogical tree. For the purpose of model selection I show that the usual BIC does not provide an appropriate approximation to the posterior probability of a model. To correct for this, an effective sample size is introduced for parameters that are inconsistently estimated. For biological studies, this work implies that tree-aware sampling design is desirable; adding more sampling units may not help ancestral reconstruction and only strong lineage effects may be detected with high power.
Published in at http://dx.doi.org/10.1214/08-AOAS173 the Annals of Applied Statistics (http://www.imstat.org/aoas/) by the Institute of Mathematical Statistics (http://www.imstat.org)
Cited by in corpus (10)
- Detecting adaptive evolution in phylogenetic comparative analysis using the Ornstein-Uhlenbeck model
- Asymptotic theory with hierarchical autocorrelation: Ornstein-Uhlenbeck tree models
- Inferring evolutionary histories of pathway regulation from transcriptional profiling data
- Phylogenetic effective sample size
- Phylogenetic confidence intervals for the optimal trait value
- A Consistent Estimator of the Evolutionary Rate
- On the convergence of the maximum likelihood estimator for the transition rate under a 2-state symmetric model
- Stochastic models in phylogenetic comparative methods: analytical properties and parameter estimation
- Generalization and Robustness of Batched Weighted Average Algorithm with V-geometrically Ergodic Markov Data
- Learning From Non-iid Data: Fast Rates for the One-vs-All Multiclass Plug-in Classifiers