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

Can seqminer reads and merges these large GWAS files? #23

Open
jielab opened this issue Sep 26, 2021 · 0 comments
Open

Can seqminer reads and merges these large GWAS files? #23

jielab opened this issue Sep 26, 2021 · 0 comments

Comments

@jielab
Copy link

jielab commented Sep 26, 2021

Hi, there:

A recent Nature Genetics paper posted 35 large GWAS files at https://doi.org/10.35092/yhjc.12355382. Compared with the regular .tsv format, I assume that these BED format files should be much faster to be queries and managed.

Now I want to combine the BETA and SE of all these 35 GWAS into a single file. Reading each of them into R using read.table() and then merge() would take impossible computational time and that did not utilize the BED format. Can you please let me know if seqminer or something else could do this merging?

Similarly, these GWAS files don’t have rsIDs, and I need to add them. Can you please also suggest a good approach to do this. If these were VCF files, I could use bcftools annotate to add rsID from a dbSNP reference file.

Thank you & Best regards,
Jie

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant