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

GC lambda #35

Open
Cui4399 opened this issue Mar 11, 2024 · 0 comments
Open

GC lambda #35

Cui4399 opened this issue Mar 11, 2024 · 0 comments

Comments

@Cui4399
Copy link

Cui4399 commented Mar 11, 2024

“We provide a simple script (misc/estimate_lambdagc.py) to infer the genomic control factor to reinflate the summary statistics if it is unknown. The following script should be executed after step 1 completes.” I wonder if this estimate of GC lambda is the same as the GC lambda from ldsc.py?
For example:
Call:
./ldsc.py
--no-check-alleles
--h2 ./Trait5.sumstats.gz
--ref-ld-chr ./eur_w_ld_chr/
--out ./Trait5_h2
--w-ld-chr ./eur_w_ld_chr/
Total Observed scale h2: 0.002 (0.0013)
Lambda GC: 1.0046
Mean Chi^2: 1.0109
Intercept: 0.9948 (0.0071)
Ratio < 0 (usually indicates GC correction).
Analysis finished at Mon Mar 4 15:11:53 2024
Total time elapsed: 30.09s
I'm asking this because the GC lambda estimate from estimate_lambdagc.py is different from the GC lambda from ldsc

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