Difference between revisions of "Genotype Likelihood Based Allele Balance"
From Genome Analysis Wiki
Jump to navigationJump to searchLine 11: | Line 11: | ||
:<math> | :<math> | ||
− | + | nref_i = 0.5 * depth * (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) * | + | AB = {0.05 + \sum_i P_i(het) * nref_i\over 0.10 + \sum_i P_i(het)*depth_i} |
</math> | </math> | ||
Revision as of 15:21, 22 October 2010
A statistic developed by Dr. Tom Blackwell.
Inputs
- depth
- at each site
- lrr
- phred-scaled likelihood of observed bases given reference genotype
- lsr
- phred-scaled likelihood of observed bases given heterozygous genotype
- lss
- phred-scaled likelihood of observed bases given non-reference genotype
Calculation
Points
Three points for correctly deciphering why this formula works.