Question: How to calculate LD values of SNP based on chromosomal positions and their call rate?
gravatar for blacktomato27
21 months ago by
United States
blacktomato2740 wrote:

Hi to all Good Afternoon

I have file containing SNPs (affy 44k chip) with their call rate and chromosomal positions with me and these belongs to rice crop. Now how can i calculate LD values (r2) for these SNPs, my basic objective is to select tag SNPs based on LD for genotyping work. any help with any R packages or own functions in this regard will be highly appreciated Thanks in Advance


snp R • 826 views
ADD COMMENTlink written 21 months ago by blacktomato2740

To calculate LD you need a set comprising many genotypes (similar hapmap or 1kG, see also 1000 Genomes Ld Calculation ) from possibly many populations, do you have that? A quick google search did not reveal a genome wide public resource like that, but some older papers studying a selection of genotypes.

ADD REPLYlink modified 21 months ago • written 21 months ago by Michael Dondrup46k

You can try to convert the data to PLINK format and then run a LD analysis there. Affymetrix PowerTools has a function that converts genotype data to PLINK format:

ADD REPLYlink written 21 months ago by Kevin Blighe46k
Please log in to add an answer.


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