-
Notifications
You must be signed in to change notification settings - Fork 10
Add prefix argument variant_id to plink conversion #390
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
Add prefix argument variant_id to plink conversion #390
Conversation
plink1.9 --bfile example --recode vcf-iid --out tmp mv tmp example.vcf
This is ready for a look now and should get us to full initial plink support. I think the definition of correctness here being equivalent to running plink's own VCF conversion through vcf2zarr makes things much easier, and it fairly clear we're doing the correct thing now. |
7ceaf57
to
8a2192c
Compare
c4a86b6
to
8202bdc
Compare
Ok, I think that's plink support done now. The coverage misses are incidental I think and we can ignore them. |
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.
Very nice, good to nail down and use plink --recode
as the validation
# bed and bim files, but should share a .fam | ||
self.prefix = str(prefix) | ||
paths = PlinkPaths( | ||
self.prefix + ".bed", |
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.
Might be nicer to have the PlinkPaths __init__
just take the prefix?
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.
Pros and cons - I tend to try and keep these dataclasses as simple as possible for reasoning about serialising and so on.
("tests/data", 5045032), | ||
# NOTE: removed this after additional plink files broke it in #390. | ||
# It's not worth the trouble - see #330 | ||
# ("tests/data", 5053610), |
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.
Ha, glad you removed this. Was driving me potty.
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.
LGTM
Looks like CI is complaining about a missing docs section |
Fixup tests for plink Covert some corner cases in the plink outout
Disable the top-level test as it's not working well. Not clear we want this tbh
Add the docs files
8202bdc
to
b51c8d5
Compare
Ha, whoops, forgot to add the docs files! |
Docs look good. |
No description provided.