Biostar Beta. Not for public use.
Estimating LD score in PLINK
0
Entering edit mode
17 months ago
United Kingdom

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?

ADD COMMENTlink
0
Entering edit mode

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

ADD REPLYlink
0
Entering edit mode

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
0
Entering edit mode

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

ADD REPLYlink
0
Entering edit mode

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

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.3.1