-
Notifications
You must be signed in to change notification settings - Fork 1
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
Fix SNP writers related to pathlib. Refactor and document Ancestry Objects, SNPObject, and dim reduction classes #5
Conversation
- Rename average_strands to is_phased. - Add common_variants_intersection to subset_to_common_variants. - Add common_markers_intersection to subset_to_common_markers. - Rename filename to file. - Rename correct_snp_variants to correct_flipped_variants. - Refine docstrings for improved clarity and consistency.
…ean_ attributes. Remove repeated properties and setters. Solve bug in SNP_PCA.ipynb
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 tried doing:
import snputils as su
input_path = "/private/groups/ioannidislab/salcc/snputils_data/bed/"
filename = "ALL.chr22.shapeit2_integrated_v1a.GRCh38.20181129.phased"
output_path = "/private/home/mcomajoa/snputils/"
snpobj = su.read_bed(input_path + filename)
and then either of these options:
snpobj.save_vcf(output_path + filename)
snpobj.save_bed(output_path + filename)
snpobj.save_pgen(output_path + filename)
and all of them result in a different error. I think only the save_bed error is related to the changes in this PR, and I see that before the PR it also resulted in another error, but it would still be nice to have all the writers working. If you think it's too much work let me know, and I can approve this PR for the doc changes and we create another one later to fix these issues.
- Fixed issue with when handling filenames not ending in . - Resolved unintended modification of within the BEDWriter class by using a deep copy.
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.
Now the path handling seems to work correctly. Thanks!
No description provided.