Question: plink calculating LD using R2 command
0
gravatar for janhuang.cn
9 months ago by
janhuang.cn80
janhuang.cn80 wrote:

I am calculating the R2 of SNPs from 1000G phase 3 vcf files.

I first converted vcf to bed, and used R2 command. But I found some discrepancy of my results from the results on LD link, an online tool to calculate LD based on 1000G phase 3. The reason I calculated myself is I need a full LD table for all available SNPs from 1000G.

Both rs138281 and rs138285 are in the bed file, but they are not in LD. No R2 was returned by R2 command. I did not set a R2 cutoff, so it should returl all R2.

plink --bfile chr22_1000Gphase3_EUR_snp_maf_rmvsnp --r2 --out chr22_1000Gphase3_EUR_ldtable

However, they are highly correlated according to LD link (R2=0.9413 among EUR population)

Does anyone have an idea what might go wrong?

r2 ld plink • 641 views
ADD COMMENTlink modified 9 months ago by pfs230 • written 9 months ago by janhuang.cn80
1
gravatar for pfs
9 months ago by
pfs230
USA/Boston
pfs230 wrote:

Just a guess here. From the Plink documentation on LD.https://www.cog-genomics.org/plink/1.9/ld

The default output is the limited table format. Since you did not request a different format this should be the one you are getting. Using the default the LD will only be calculated for variants in close proximity (10-1). I am guessing you have more than 9 other SNPs between these two SNPs so LD for these two SNPs are not being calculated.

Below is copy paste from Plink documentation

By default, --r calculates and reports raw inter-variant allele count correlations, while --r2 reports squared correlations. You can request values for all pairs in matrix format (if you specify 'bin' and/or one of the matrix shape modifiers), all pairs in table format (with 'inter-chr'), or a limited window in table format (this is the default). Results are saved to plink.ld{.gz/.bin}.

By default, when a limited window report is requested, every pair of variants with at least (10-1) variants between them, or more than 1000 kilobases apart, is ignored. You can change the first threshold with --ld-window, and the second threshold with --ld-window-kb.

ADD COMMENTlink modified 9 months ago • written 9 months ago by pfs230

Thank you. That is the reason.

ADD REPLYlink written 9 months ago by janhuang.cn80
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: 1560 users visited in the last hour