Question: Estimating LD score in PLINK
0
gravatar for GabrielMontenegro
6 months ago by
United Kingdom
GabrielMontenegro540 wrote:

LD score is defined as the sum of LD r2 between a variant and all the variants in a region.

I want to estimate the LD score for an arbitrary numbers of SNPs in a PLINK format dataset that I have. I was simply planning to get the LD between all snps to get the LDscore. I am trying to do this using this command, but I get an empty output:

plink --bfile myfile --extract snps_to_calculate_ldscore.txt \
      --r2  --ld-window-kb 9999999 \
      --ld-window-r2 0 --out ldscore

Any ideas?

ld plink ldscore • 369 views
ADD COMMENTlink modified 6 months ago by RamRS25k • written 6 months ago by GabrielMontenegro540

At minimum, you need to add --ld-window.

ADD REPLYlink modified 6 months ago by h.mon28k • written 6 months ago by chrchang5236.2k

Thanks for the quick reply. What --ld-window parameter should I use in order to always get r2 values between the SNPs that I am extracting with --extract?

ADD REPLYlink written 6 months ago by GabrielMontenegro540

"9999999" should be fine; with more, you'd run out of disk space anyway.

ADD REPLYlink written 6 months ago by chrchang5236.2k

When I add the --ld-window parameter with the other commands, I get this message: Error: Missing --ld-window parameter. What parameter for --ld-window should I use here? Thanks.

ADD REPLYlink written 6 months ago by GabrielMontenegro540
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2351 users visited in the last hour