|analyzing pooled sequenced data with selection|| Main project page|
Previous entry Next entry
Intro to Wen & Stephens in 2D
Suppose we have only summary-level data for haplotypes h1,h2,...,h2n. Specifically let the summary-level data be denoted by . We assume in this two locus model, the first locus is typed and the second locus is untyped. We hope to predict what the allele frequency is at the untyped SNP using information from panel data (perhaps this can be interpreted as our prior). Formally, let y1 denote the allele frequency at the typed SNP and y2 the allele frequency at the untyped SNP. We assume that h1,h2,...,h2n are independent and identically distributed from P(M) (our prior).
(assign as the vector of untyped SNPs) is a function of both the panel data (μ2) and the typed SNPs (y1).
Li & Stephens in 2D
We describe the Li & Stephens haplotype copying model: Let h1,h2,...,hk denote the k sampled haplotypes at 2 loci. Thus there are 4 possible haplotypes. The first haplotype is randomly chosen with equal probability from the four possible haplotypes.
Consider now the conditional distribution of hk + 1 given h1,h2,...,hk. Recall the intuition is that hk + 1 is a mosaic of h1,h2,..,hk.
Let Xj denote which hapolotype hk + 1 copies at site j (so ).
We model Xj as a markov chain on 1,..,k with . The transition probabilities are:
if x' = x and
otherwise. ρj and dj denote recombination and physical distances, respectively.
Now in a hidden markov model, there is also the transmission process. To mimic the effects of mutation, the copying process may be imperfect. if hx,j = a and otherwise. , where the motivation is the more haplotypes the less frequent mutation occurs.