Question: Is it possible to specify a reference allele while computing LD in Plink?
2
gravatar for User 7754
3.5 years ago by
User 7754230
United Kingdom
User 7754230 wrote:

Hi,

 

I am using methods of analysis where it is necessary that the LD and effect sizes are computed with respect to the same reference allele, while in my case they are not since I am using an available panel to compute LD with Plink and my own data to compute effect size. To align the two measurements, I would have to compare each allele with those of the reference and flip the effect sizes in cases where the alleles are flipped, but things get messy when including both SNPs and INDELs in the comparisons. To avoid this step, I was hoping it could be possible to specify a reference allele directly when computing LD in Plink like this:

    plink --bfile 1000Genomes --r2  --extract list.snps --reference-allele list.reference.allele

If this is not, are there any tools out there to easily compare alleles from two datasets and adjust the affect size of the study if the alleles are flipped compared to a reference file such as 1000 Genomes? I have previously used 'QCGWAS' - CRAN which was really nice to do this, but the files I have now are all very large BED files that take too long to load into R, that is why if Plink could to this it would be much faster...  

Many thanks for any suggestions and ideas!

Fra

plink allele • 1.4k views
ADD COMMENTlink modified 2.7 years ago by Biostar ♦♦ 20 • written 3.5 years ago by User 7754230
2
gravatar for User 7754
3.5 years ago by
User 7754230
United Kingdom
User 7754230 wrote:

From Christopher Chang

Yes, --a2-allele lets you choose the reference alleles for the current run.

Based on  Christopher Chang's answer, to compute LD in Plink using a specific allele as a reference we can use:

plink --bfile 1000Genomes --r2  --extract list.snps ---a2-allele list.reference.allele 

 

ADD COMMENTlink written 3.5 years ago by User 7754230
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: 1151 users visited in the last hour