Parametric Bayesian inference for Y-linked two-sex branching models

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.

[1]  Miguel González,et al.  Parametric inference for Y-linked gene branching models: Expectation-maximization method , 2010 .

[2]  Andrew Gelman,et al.  General methods for monitoring convergence of iterative simulations , 1998 .

[3]  Rodrigo Martínez,et al.  Expectation-Maximization Algorithm for Determining Natural Selection of Y-Linked Genes Through Two-Sex Branching Processes , 2012, J. Comput. Biol..

[4]  M González,et al.  Bisexual branching processes in a genetic context: rates of growth for Y-linked genes. , 2008, Mathematical biosciences.

[5]  M González,et al.  Bisexual branching processes in a genetic context: the extinction problem for Y-linked genes. , 2006, Mathematical biosciences.

[6]  M. Jobling,et al.  Y-chromosomal diversity in the population of Guinea-Bissau: a multiethnic perspective , 2007, BMC Evolutionary Biology.

[7]  M. Hurles,et al.  European Y-chromosomal lineages in Polynesians: a contrast to the population structure revealed by mtDNA. , 1998, American journal of human genetics.

[8]  James O. Berger,et al.  Ordered group reference priors with application to the multinomial problem , 1992 .

[9]  D. Bates,et al.  Output Analysis and Diagnostics for MCMC , 2015 .

[10]  M. Hammer,et al.  Y-chromosome lineages trace diffusion of people and languages in southwestern Asia. , 2001, American journal of human genetics.

[11]  M. Molina,et al.  Bayesian estimation in the class of bisexual branching processes with population-size dependent mating , 2008 .

[12]  M. Molina,et al.  Two-sex Branching Models with Random Control on the Number of Progenitor Couples , 2012 .

[13]  M. González,et al.  Non-parametric Bayesian estimation for multitype branching processes through simulation-based methods , 2008, Comput. Stat. Data Anal..

[14]  G. Alsmeyer,et al.  Limiting genotype frequencies of Y-linked genes through bisexual branching processes with blind choice. , 2011, Journal of theoretical biology.

[15]  Peter Guttorp,et al.  Statistical inference for branching processes , 1991 .

[16]  Kateryna D. Makova,et al.  Evolution and Survival on Eutherian Sex Chromosomes , 2009, PLoS genetics.

[17]  C P Farrington,et al.  Branching process models for surveillance of infectious diseases controlled by mass vaccination. , 2003, Biostatistics.

[18]  Manuel Mota,et al.  Bisexual branching processes to model extinction conditions for Y-linked genes. , 2009, Journal of theoretical biology.

[19]  D. Rubin,et al.  Inference from Iterative Simulation Using Multiple Sequences , 1992 .