Changes

From Genome Analysis Wiki
Jump to navigationJump to search
Line 7: Line 7:  
; lsr : phred-scaled likelihood of observed bases given heterozygous genotype
 
; lsr : phred-scaled likelihood of observed bases given heterozygous genotype
 
; lss : phred-scaled likelihood of observed bases given non-reference genotype
 
; lss : phred-scaled likelihood of observed bases given non-reference genotype
 +
 +
== Calculation ==
    
:<math>
 
:<math>
nref = {0.5 * depth * (1 + (lss - lrr)\over(lss + lrr - 2 * lrs + 6 * depth))} \  
+
nref_i = 0.5 * depth_i * (1 + {(lss_i - lrr_i)\over lss_i + lrr_i - 2 * lrs_i + 6 * depth_i}) \  
 
</math>
 
</math>
    
:<math>
 
:<math>
AB = (0.05 + {\sum_i P_i(het) * nref)\over(0.10 + \sum_i P_i(het)*depth)}
+
AB = {0.05 + \sum_i P_i(het) * nref_i\over 0.10 + \sum_i P_i(het)*depth_i}
 
</math>
 
</math>
 +
 +
== Points ==
 +
 +
Three points for correctly deciphering why this formula works.

Navigation menu