Issues with partitioning Heritability

269 views
Skip to first unread message

Bakhtawar Noor

unread,
Oct 7, 2017, 6:09:40 AM10/7/17
to ldsc_users
I tried partitioning heritability using the following command:

python ldsc.py 
	--h2 BMI.sumstats.gz\
	--ref-ld-chr baseline.\ 
	--w-ld-chr weights.\
	--overlap-annot\
	--frqfile-chr 1000G.mac5eur.\
	--out BMI_baseline

However, I get the following error:

IOError: Could not open baseline.1.l2.ldscore[./gz/bz2]


I've attached the image of the traceback error as well.


And just so you know I have a ldsc folder which I cloned from github. In that I have the folder baseline and in the baseline folder I have files like baseline.16.annot.gz, baseline.16.l2.ldscore.gz, and baseline.16.l2.M_5_50. I have attached the screenshot of baseline folder too. And I am using mac. I would really appreciate if someone can help me in this regard. Thank you in advance.

Error traceback.png
Screen Shot 2017-10-07 at 12.06.45 PM.png

Raymond Walters

unread,
Oct 9, 2017, 11:38:35 AM10/9/17
to Bakhtawar Noor, ldsc_users
Hi,
If I understand your folder structure correctly the problem is that you aren’t specifying the full folder + file stem of the reference files. You probably want something like:

python ldsc.py \
--h2 BMI.sumstats.gz \
--ref-ld-chr ./baseline/baseline. \
--w-ld-chr ./weights_hm3_no_hla/weights. \
--overlap-annot \
--frqfile-chr ./1000G_frq/1000G.acm5eur. \
--out baseline

Cheers,
Raymond


--
You received this message because you are subscribed to the Google Groups "ldsc_users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to ldsc_users+...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/ldsc_users/3f28846a-71cf-49b8-963d-5a2732647b2f%40googlegroups.com.
For more options, visit https://groups.google.com/d/optout.
<Error traceback.png><Screen Shot 2017-10-07 at 12.06.45 PM.png>

Reply all
Reply to author
Forward
0 new messages