Computational aspects of DNA mixture analysis
Abstract
Statistical analysis of DNA mixtures is known to pose computational challenges due to the enormous state space of possible DNA profiles. We propose a Bayesian network representation for genotypes, allowing computations to be performed locally involving only a few alleles at each step. In addition, we describe a general method for computing the expectation of a product of discrete random variables using auxiliary variables and probability propagation in a Bayesian network, which in combination with the genotype network allows efficient computation of the likelihood function and various other quantities relevant to the inference. Lastly, we introduce a set of diagnostic tools for assessing the adequacy of the model for describing a particular dataset.
- Publication:
-
arXiv e-prints
- Pub Date:
- July 2013
- DOI:
- 10.48550/arXiv.1307.4956
- arXiv:
- arXiv:1307.4956
- Bibcode:
- 2013arXiv1307.4956G
- Keywords:
-
- Statistics - Methodology;
- Quantitative Biology - Quantitative Methods;
- 62-04 (Primary) 62P10 (Secondary)
- E-Print:
- doi:10.1007/s11222-014-9451-7