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

AG try debugging w extend #1

Merged
merged 28 commits into from
May 29, 2018
Merged

AG try debugging w extend #1

merged 28 commits into from
May 29, 2018

Conversation

adonomay
Copy link

@adonomay adonomay commented May 13, 2018

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)

@adonomay adonomay changed the base branch from master to ag-version-1-5-1-from-columbia-germline May 13, 2018 03:34
@adonomay adonomay changed the base branch from ag-version-1-5-1-from-columbia-germline to master May 15, 2018 01:45
Copy link

@svohr svohr left a 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;
Copy link

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?
Copy link

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
Copy link

Choose a reason for hiding this comment

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

Nice

Copy link

@ajsams ajsams left a 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.

@adonomay
Copy link
Author

@svohr the one change was passing the marker_set_number parameter to the getMarkerSet function call in lines 85 and 86 of Match.cpp

@adonomay adonomay merged commit a8fb1cc into master May 29, 2018
@adonomay adonomay deleted the ag-try-debugging-w-extend branch May 29, 2018 16:44
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Development

Successfully merging this pull request may close these issues.

3 participants