Skip to main content

Table 3 Joint constrained sparse haplotype frequency estimation algorithm

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 A + U 1 k Z k 1 T 1 T U 4 1

8

   X k + 1 = U 4 U 1 k Z k + H T ( U 2 k A ) 1 u 3 T

9

   Z k + 1 =max Shrink X k + 1 + 1 ρ U 1 k , 1 ρ , 0

10

   U 1 k + 1 = U 1 k + X k + 1 Z k + 1

11

   U 2 k + 1 = U 2 k +H X k + 1 A

12

until convergence