Abstract
A Y-linked two-sex branching process with blind choice is a suitable model for analyzing the evolution of the number of carriers of two alleles of a Y-linked gene in a two-sex monogamous population where each female chooses her partner from among the male population without caring about his type (i.e., the allele he carries). This work focuses on the development of Bayesian inference for this model, considering a parametric framework with the reproduction laws belonging to the power series family of distributions. A sample is considered given by the observation of the total number of females and males (regardless of their types) up to some generation as well as the number of each type of male in the last generation. Using a simulation method based on the Gibbs sampler, we approximate the posterior distributions of the main parameters of this model. The accuracy of the procedure based on this sample is illustrated by way of a simulated example.
This is a preview of subscription content, log in via an institution to check access.
Similar content being viewed by others
References
Alsmeyer, G., Gutiérrez, C., Martínez, R.: Limiting genotype frequencies of Y-linked genes through bisexual branching processes with blind choice. J. Theor. Biol. 275, 42–51 (2011)
Berger, J., Bernardo, J.: Ordered group reference prior with application to a multinomial problem. Biometrika 79, 25–37 (1992)
Brooks, S., Gelman, A.: General methods for monitoring convergence of iterative simulations. J. Comput. Graph. Stat. 7, 434–455 (1998)
Farrington, C.P., Kanaan, M.N., Gay, N.J.: Branching process models for surveillance of infectious diseases controlled by mass vaccination. Biostatistics 4(2), 279–295 (2003)
Gelman, A., Rubin, D.: Inference from iterative simulation using multiple sequences. Stat. Sci. 7(4), 457–511 (1992)
González, M., Hull, D.M., Martínez, R., Mota, M.: Bisexual branching processes in a genetic context: The extinction problem for Y-linked genes. Math. Biosci. 202, 227–247 (2006)
González, M., Martín, J., Martínez, R., Mota, M.: Non-parametric Bayesian estimation for multitype branching processes through simulation-based methods. Comput. Stat. Data Anal. 52, 1281–1291 (2008a)
González, M., Martínez, R., Mota, M.: Bisexual branching processes in a genetic context: Rates of growth for Y-linked genes. Math. Biosci. 215, 167–176 (2008b)
González, M., Martínez, R., Mota, M.: Bisexual branching processes to model extinction conditions for Y-linked genes. J. Theor. Biol. 258, 478–488 (2009)
González, M., Gutiérrez, C., Martínez, R.: Parametric inference for Y-linked gene branching models: expectation-maximization method. In: González, M., del Puerto, I., Martínez, R., Molina, M., Mota, M., Ramos, A. (eds.) Proceedings of Workshop on Branching Processes and Their Applications. Lecture Notes in Statistics, vol. 197, pp. 191–204. Springer, Berlin (2010)
González, M., Gutiérrez, C., Martínez, R.: Expectation-maximization algorithm for determining natural selection of Y-linked genes through two-sex branching processes. J. Comput. Biol. (2012, accepted)
Guttorp, P.: Statistical Inference for Branching Processes. Wiley, New York (1991)
Hurles, M., Irven, C., Nicholson, J., Taylor, P., Santos, F., Loughlin, J., Jobling, M., Sykes, B.: European Y-chromosomal lineages in Polynesians: a contrast to the population structure revealed by mtDNA. Am. J. Hum. Genet. 63(6), 1793–1806 (1998)
Lucy, D., Aykroyd, R.: GenKern: Functions for generating and manipulating binned kernel density estimates (2010). http://CRAN.R-project.org/package=GenKern. R package version 1.1-10
Molina, M., Mota, M., Ramos, A.: Bayesian estimation in the class of bisexual branching processes with population-size dependent mating. Test 17, 179–196 (2008)
Molina, M., Mota, M., Ramos, A.: Two-sex branching models with random control on the number of progenitor couples. Methodol. Comput. Appl. Probab. 14(1), 35–48 (2012)
Nature (2003). www.nature.com/nature/focus/ychromosome/#evolution
Plummer, M., Best, N., Cowles, K., Vines, K.: coda: output analysis and diagnostics for MCMC (2010). http://CRAN.R-project.org/package=coda. R package version 0.14-2
Quintana-Murci, L., Krausz, C., Zerjal, T., Sayar, S.H., Hammer, M.F., Mehdi, S., Ayub, Q., Qamar, R., Mohyuddin, A., Radhakrishna, U., Jobling, M., Tyler-Smith, C., McElreavey, K.: Y-chromosome lineages trace diffusion of people and languages in southwestern Asia. Am. J. Hum. Genet. 68(2), 537–542 (2001)
R Development Core Team: R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna (2011). http://www.R-project.org/. ISBN 3-900051-07-0
Rosa, A., Ornelas, C., Jobling, M., Brehm, A., Villems, R.: Y-chromosomal diversity in the population of Guinea-Bissau: a multiethnic perspective. BMC Evol. Biol. 27, 107–124 (2007)
Wilson, M., Makova, K.: Evolution and survival on eutherian sex chromosomes. PLoS Genet. 5(7), e1000568 (2009). doi: 10.1371/journal.pgen.1000568
Y-Chromosome Consortium (2008). ycc.biosci.arizona.edu
Acknowledgements
This research was supported by the Ministerio de Ciencia e Innovación and the FEDER, grant MTM2009-13248.
The authors would like to thank the associate editor and the reviewers for providing valuable comments and suggestions which have significantly improved this paper.
Moreover, they are also grateful to Professors Horacio González-Velasco and Carlos García-Orellana for providing them with the computational support.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
González, M., Gutiérrez, C. & Martínez, R. Parametric Bayesian inference for Y-linked two-sex branching models. Stat Comput 23, 727–741 (2013). https://doi.org/10.1007/s11222-012-9342-8
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11222-012-9342-8