Skip to main content

Table 5 Joint constrained sparse haplotype frequency estimation algorithm with missing data

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   For i = 1,,O
8     u 3 , i = 1 T U 4 H i T u 2 , i k a ~ i + u i 1 z i k 1 1 T U 4 1
9     x i k + 1 = U 4 u 1 , i k z i k + H i T ( u 2 , i k a i ~ ) u 3 , i 1
10   end for;
11     Z k + 1 =max Shrink X k + 1 + 1 ρ U 1 k , 1 ρ , 0
12     U 1 k + 1 = U 1 k + X k + 1 Z k + 1
13   For i = 1,,O
14     u 2 , i k + 1 = u 2 , i k + H i x i k + 1 a ~ i
15   end for
16 until convergence