Biostar Beta. Not for public use.
Question: "hrc-1000g-check-bim.py" Unable to open specified bim file
0
Entering edit mode

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
ADD COMMENTlink 9 months ago zhangcc89 • 0 • updated 9 months ago Tao • 300
Entering edit mode
1

Isn't there a suffix missing? iha_QCQC7

ADD REPLYlink 9 months ago
ATpoint
17k
Entering edit mode
0

Nope. I have double checked.

ADD REPLYlink 9 months ago
zhangcc89
• 0
Entering edit mode
0

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?

ADD REPLYlink 9 months ago
ATpoint
17k
Entering edit mode
0

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
ADD REPLYlink 9 months ago
zhangcc89
• 0
Entering edit mode
0

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.

ADD REPLYlink 9 months ago
WouterDeCoster
39k
Entering edit mode
0

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

ADD REPLYlink 9 months ago
WouterDeCoster
39k
2
Entering edit mode

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
ADD COMMENTlink 9 months ago Tao • 300
Entering edit mode
0

Cool! Thanks very much.

ADD REPLYlink 9 months ago
zhangcc89
• 0
Entering edit mode
1

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.

ADD REPLYlink 9 months ago
ATpoint
17k
Entering edit mode
0

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

ADD REPLYlink 9 months ago
WouterDeCoster
39k

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.0