-
Notifications
You must be signed in to change notification settings - Fork 0
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
AG try debugging w extend #1
Conversation
…tAussie/GERMLINE into ag-try-debugging-w-extend
… when answers are inconsistent
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Can you highlight what you changed in the logic of the program besides debug prints/variable/names/style? It's difficult to see how the behavior has changed.
@@ -90,7 +95,7 @@ int main(int argc, char* argv[]) | |||
|
|||
GERMLINE germline; | |||
germline.mine( params ); | |||
return 1; | |||
return 0; |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
There's a spot in the pipeline where I remember we were handling this bug, so we'll need to fix that.
{ | ||
if ( node[0]->getChromosome( a )->getMarkerSet()->getMarkerBits()[cur_marker] != node[1]->getChromosome( b )->getMarkerSet()->getMarkerBits()[cur_marker] ) | ||
err = true; | ||
found_mismatch = true; // err may actually mean end-of-homozygous-tract? |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
It looks that way
* `germline -haploid -w_extend -bits 41` (or any bits value > 1) now produces correct match tracts | ||
* exit codes follow shell-script convention: 1 for error, 0 for success | ||
* `germline -version` now reports the version number | ||
* added tests |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Nice
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Looks good to me, admittedly, I'm not fluent in C++ so don't understand all of this.
@svohr the one change was passing the marker_set_number parameter to the getMarkerSet function call in lines 85 and 86 of Match.cpp |
needs testing of haploid and homoz-only -w_extend with multiple dogs.
Edit: have tested this in bulk on dobermans and ref panel v7, looks good (eCOI is consistent with measured COI)