Skip to main content

Table 4 Joint constrained sparse haplotype frequency estimation algorithm in the presence of noisy measurements

From: Maximum-parsimony haplotype frequencies inference based on a joint constrained sparse representation of pooled DNA

1 Set k = 0
2 Set ρ > 0
3 Set X0 = 0, Z0 = 0, U 1 0 = 0 , U 2 0 = 0
4 U4 = (IHT(I+H HT)−1H)
5 Repeat
6   k = k + 1
7    u 3 T = 1 T U 4 H H U 2 k Z 1 k + U 1 Z 2 k 1 T 1 T U 4 1
8    X k + 1 = U 4 U 1 k Z 1 k + H T ( U 2 k Z 2 k ) 1 u 3 T
9    Z 1 k + 1 =max Shrink X k + 1 + 1 ρ U 1 k , 1 ρ , 0
10    Z 2 k + 1 =proj(H, X k + 1 , U 2 k ,A,δ)
11    U 1 k + 1 = U 1 k + X k + 1 Z 1 k + 1
12    U 2 k + 1 = U 2 k +H X k + 1 Z 2 k + 1
13 until convergence