Question: Estimating LD score in PLINK
0
gravatar for GabrielMontenegro
18 months ago by
United Kingdom
GabrielMontenegro560 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 • 992 views
ADD COMMENTlink modified 18 months ago by _r_am31k • written 18 months ago by GabrielMontenegro560

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

ADD REPLYlink modified 18 months ago by h.mon31k • written 18 months ago by chrchang5237.4k

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 18 months ago by GabrielMontenegro560

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

ADD REPLYlink written 18 months ago by chrchang5237.4k

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 18 months ago by GabrielMontenegro560
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: 1919 users visited in the last hour