A new estimator of the discovery probability
AbstractSpecies sampling problems have a long history in ecological and biological studies and a number of issues, including the evaluation of species richness, the design of sampling experiments, the estimation of rare species variety, are to be addressed. Such inferential problems have recently emerged also in genomic applications, however exhibiting some peculiar features that make them more challenging: specifically, one has to deal with very large populations (genomic libraries) containing a huge number of distinct species (genes) and only a small portion of the library has been sampled (sequenced). These aspects motivate the Bayesian nonparametric approach we undertake, since it allows to achieve the degree of flexibility typically needed in this framework. Basing on an observed sample of size n, focus will be on prediction of a key aspect of the outcome from an additional sample of size m, namely the so–called discovery probability. In particular, conditionally on an observed basic sample of size n, we derive a novel estimator of the probability of detecting, at the (n + m + 1)–th observation, species that have been observed with any given frequency in the enlarged sample of size n + m. Such an estimator admits a closed form expression that can be exactly evaluated. The result we obtain allows us to quantify both the rate at which rare species are detected and the achieved sample coverage of abundant species, as m increases. Natural applications are represented by the estimation of the probability of discovering rare genes within genomic libraries and the results are illustrated by means of two Expressed Sequence Tags datasets.
Download InfoIf you experience problems downloading a file, check if you have the proper application to view it first. In case of further problems read the IDEAS help page. Note that these files are not on the IDEAS site. Please be patient as the files may be large.
Bibliographic InfoPaper provided by University of Pavia, Department of Economics and Management in its series DEM Working Papers Series with number 007.
Length: 28 pages
Date of creation: Oct 2012
Date of revision:
Bayesian nonparametrics; Gibbs–type priors; Rare species discovery; Species sampling models; Two–parameter Poisson–Dirichlet process.;
This paper has been announced in the following NEP Reports:
Please report citation or reference errors to , or , if you are the registered author of the cited work, log in to your RePEc Author Service profile, click on "citations" and make appropriate adjustments.:
- Chang Xuan Mao, 2002. "A Poisson model for the coverage problem with a genomic application," Biometrika, Biometrika Trust, vol. 89(3), pages 669-682, August.
- Antonio Lijoi & Ramsés H. Mena & Igor Prünster, 2007. "A Bayesian Nonparametric Method for Prediction in EST Analysis," ICER Working Papers - Applied Mathematics Series 16-2007, ICER - International Centre for Economic Research.
- Antonio Lijoi & Ramsés H. Mena & Igor Prünster, 2007. "Bayesian Nonparametric Estimation of the Probability of Discovering New Species," Biometrika, Biometrika Trust, vol. 94(4), pages 769-786.
- Kolossiatis, M. & Griffin, J.E. & Steel, M.F.J., 2011. "Modeling overdispersion with the normalized tempered stable distribution," Computational Statistics & Data Analysis, Elsevier, vol. 55(7), pages 2288-2301, July.
- Alejandro Jara & Timothy Hanson & Fernando A. Quintana & Peter Müller & Gary L. Rosner, . "DPpackage: Bayesian Semi- and Nonparametric Modeling in R," Journal of Statistical Software, American Statistical Association, vol. 40(i05).
- Antonio Lijoi & Igor Pruenster & Stephen G. Walker, 2008. "Bayesian nonparametric estimators derived from conditional Gibbs structures," ICER Working Papers - Applied Mathematics Series 06-2008, ICER - International Centre for Economic Research.
For technical questions regarding this item, or to correct its authors, title, abstract, bibliographic or download information, contact: (Alice Albonico).
If references are entirely missing, you can add them using this form.