-
Notifications
You must be signed in to change notification settings - Fork 14
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
edits related to production pipeline #114
edits related to production pipeline #114
Conversation
if wildcards.annots == "all_sites": | ||
# Adding selection to the whole contig | ||
contig.add_dfe(intervals=np.array([[0, int(contig.length)]]), DFE=dfe) | ||
else: | ||
## Adding annotation only seletion on exon region | ||
annot = species.get_annotations(wildcards.annots) | ||
annot_intervals = annot.get_chromosome_annotations(wildcards.chrms) | ||
contig.add_dfe(intervals=annot_intervals, DFE=dfe) |
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 someone check this? I was getting a dfe
is called before it's defined error. This logic resolves that, but I want to make sure it's correct otherwise.
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.
this looks good to me
I think this PR is good to merge now. Still having some msmc2 hiccups with Vaquita, but I'll troubleshoot and submit some solution in a different PR (perhaps not running msmc2 on Vaquita if nothing else). |
closes #115 |
Silas is this ready to go? should I merge this? |
Yeah I think so! I don't anticipate adding more changes. |
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.
i went through this and it looks correct to me. haven't run it tho!
Needed to make some edits to get production pipeline going. Caught a bug or two in the process.
Let's leave this open for a little bit. May need to add more commits.