Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Filter M2 calls that are near other filtered calls on the same haplotype #5092

Merged
merged 6 commits into from
Aug 8, 2018

Conversation

davidbenjamin
Copy link
Contributor

@LeeTL1220 This gets rid of a bunch of false positives with no effect on sensitivity. It also does some useful refactoring of the filtering engine in order to exploit Takuto's two-pass formalism for more filters than just OB.

Copy link
Contributor

@LeeTL1220 LeeTL1220 left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@davidbenjamin Minor stuff. Feel free to merge after you address the comments.

* will not exceed the requested false positive rate
*
* @param posteriors A list of posterior probabilities, which gets sorted
* @param requestedFPR We set the filtering threshold such that the FPR doesn't exceed this value
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Since this is a public static method, please document and have checks for requestedFPR. For example, can it be negative?

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

done

import java.io.IOException;
import java.util.*;

public class FilteringFirstPass {
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Could you add some class javadoc? Would like to know how this fits in the grand scheme of things...

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

done

@codecov-io
Copy link

Codecov Report

Merging #5092 into master will increase coverage by 0.006%.
The diff coverage is 96.33%.

@@               Coverage Diff               @@
##              master     #5092       +/-   ##
===============================================
+ Coverage     86.478%   86.485%   +0.006%     
- Complexity     29155     29168       +13     
===============================================
  Files           1813      1813               
  Lines         135228    135256       +28     
  Branches       15023     15027        +4     
===============================================
+ Hits          116943    116976       +33     
+ Misses         12825     12818        -7     
- Partials        5460      5462        +2
Impacted Files Coverage Δ Complexity Δ
...ute/hellbender/utils/variant/GATKVCFConstants.java 80% <ø> (ø) 4 <0> (ø) ⬇️
...r/tools/walkers/mutect/Mutect2FilteringEngine.java 84.211% <100%> (-1.155%) 50 <4> (-4)
...bender/tools/walkers/mutect/FilterMutectCalls.java 97.222% <100%> (ø) 12 <2> (ø) ⬇️
...ols/walkers/mutect/FilteringFirstPassUnitTest.java 100% <100%> (ø) 3 <2> (?)
...ls/walkers/mutect/M2FiltersArgumentCollection.java 100% <100%> (ø) 1 <0> (ø) ⬇️
...e/hellbender/utils/variant/GATKVCFHeaderLines.java 99.333% <100%> (+0.004%) 10 <0> (ø) ⬇️
...ender/tools/walkers/mutect/FilteringFirstPass.java 95.652% <95.652%> (ø) 21 <21> (?)
...ithwaterman/SmithWatermanIntelAlignerUnitTest.java 90% <0%> (+30%) 2% <0%> (ø) ⬇️
...utils/smithwaterman/SmithWatermanIntelAligner.java 90% <0%> (+40%) 3% <0%> (+2%) ⬆️

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

Successfully merging this pull request may close these issues.

3 participants