This is a test version of Biostars. For the public version, visit https://www.biostars.org.
"hrc-1000g-check-bim.py" Unable to open specified bim file

Hi, I have met some problems about running hrc-1000g-check-bim.py. Can any experts give me some advice? Thanks very much.

Reference Panel:             HRC
Bim filename:                /Volumes/Mac8/GWAS_all20190509/GWAS/iha_QCQC7
Reference filename:          /Volumes/GWAS_all2019/GWAS/pipe_QC_to_Phase_githb/_Chr/HRC.r1-1.GRCh37.wgs.mac5.sites.tab
Allele frequencies filename: iha_QCQC7
Chromosome flag set:         No
Allele frequency threshold:  0.2

***ERROR****: Unable to open specified bim file: /Volumes/Mac8/GWAS_all20190509/GWAS/iha_QCQC7

Please check path and filename and try again
gwas
• 3,711 views • link
updated by Tao Tao • written by zhangcc89 • zhangcc89

Isn't there a suffix missing? iha_QCQC7

• 0 views • link
written by ATpoint ATpoint

I mean did you check that the file you find via ls /Volumes/Mac8/GWAS_all20190509/GWAS/iha_QCQC7* is not actually called /Volumes/Mac8/GWAS_all20190509/GWAS/iha_QCQC7.bim?

• 0 views • link
written by ATpoint ATpoint

The same errors with adding specified path:

Reference Panel:             HRC
Bim filename:                /Volumes/Mac8/GWAS_all/pipe_QC_to_Phase_githb/GWAS_Chr/iha_filter_QCQC7
Reference filename:          /Volumes/Mac8/GWAS_all/pipe_QC_to_Phase_githb/GWAS_Chr/HRC.r1-1.GRCh37.wgs.mac5.sites.tab
Allele frequencies filename: /Volumes/Mac8/GWAS_all/pipe_QC_to_Phase_githb/GWAS_Chr/iha_filter_QCQC7
Chromosome flag set:         No
Allele frequency threshold:  0.2
• 0 views • link
written by zhangcc89 • zhangcc89

Please use ADD COMMENT or ADD REPLY to answer to previous reactions, as such this thread remains logically structured and easy to follow. I have now moved your reaction but as you can see it's not optimal. Adding an answer should only be used for providing a solution to the question asked.

If you add new information, please update your initial question.

• 0 views • link
written by WouterDeCoster WouterDeCoster

Can you show the output of file /Volumes/Mac8/GWAS_all20190509/GWAS/iha_QCQC7 and ls /Volumes/Mac8/GWAS_all20190509/GWAS/iha_QCQC7*?

• 0 views • link
written by WouterDeCoster WouterDeCoster

1 answer

Hi, I remember it needs to specify the suffix of bim file. Such as:

perl HRC-1000G-check-bim.pl -b data.bim -f data.frq -r HRC.r1-1.GRCh37.wgs.mac5.sites.tab -h
• 0 views • link
written by Tao Tao

Me: Isn't there a suffix missing? iha_QCQC7

You: Nope. I have double checked.

Didn't you say you checked if suffixing changes things? Please see my new comment above.

• 0 views • link
written by ATpoint ATpoint

If an answer was helpful you should upvote it, if the answer resolved your question you should mark it as accepted.
Upvote|Bookmark|Accept

• 0 views • link
written by WouterDeCoster WouterDeCoster

Log in to answer this question.

More posts like this