An efficient algorithm to compute marginal posterior genotype probabilities for every member of a pedigree with loops.

ABSTRACT:

Background

Marginal posterior genotype probabilities need to be computed for genetic analyses such as geneticcounseling in humans and selective breeding in animal and plant species.

Methods

In this paper, we describe a peeling based, deterministic, exact algorithm to compute efficiently genotype probabilities for every member of a pedigree with loops without recourse to junction-tree methods from graph theory. The efficiency in computing the likelihood by peeling comes from storing intermediate results in multidimensional tables called cutsets. Computing marginal genotype probabilities for individual i requires recomputing the likelihood for each of the possible genotypes of individual i. This can be done efficiently by storing intermediate results in two types of cutsets ca

SUBMITTER: Totir LR 

PROVIDER: S-EPMC2801663 | biostudies-literature | 2009 Dec

REPOSITORIES: biostudies-literature