Question: Ldheatmap Does Not Include All Snps In The Analysis (Tri-Allelic Snp Excluded)
gravatar for Agatha
7.0 years ago by
Agatha350 wrote:

I am trying to generate an LD plot using LDheatmap in R.

I have 4 SNPs geotyped for around 80 000 individuals :


   p1205 p0864 p3077 p1244
1     G/G     C/C     T/T     C/C
2     A/A     C/C     T/T     C/C
3     G/G     T/T     T/T     T/T
4     A/A     C/C     T/T     C/C
5     G/G     C/C     T/T     C/C
6     G/A     C/C     T/T     C/C

MyHeatmap <- LDheatmap(m.genos, genetic.distances = as.numeric(m.dist),distances="physical",LDmeasure="r",title="LD analysis",    

However, the LDheatmap function computes the ld matrix using ONLY 3 of them.

> MyHeatmap$LDmatrix
           p3077    p1205     p0864
p3077      NA 0.02379055 0.02184959
p1205      NA         NA 0.21808013
p0864      NA         NA         NA

I have managed to succesfully generate the LD matrix for other SNPs, so I was wondering what could me the problem in this case?

The excluded SNP is tri allelic : => possibly causing the LDheatmap function to fail in generating the complete LD matrix

  Number of samples typed: 87058 (99.8%)

  Allele Frequency: (3 alleles)
  Count Proportion
 C  110853       0.64
 T   54539       0.31
 A    8724       0.05
 NA    272         NA

 Genotype Frequency:
    Count Proportion
   C/C 35354       0.41
   C/T 34647       0.40
   C/A  5498       0.06
   T/T  8552       0.10
   T/A  2788       0.03
   A/A   219       0.00
    NA    136         NA

  Heterozygosity (Hu)  = 0.4940388
  Poly. Inf. Content   = 0.4119684
R snp • 2.4k views
ADD COMMENTlink modified 7.0 years ago by Fabio Marroni2.6k • written 7.0 years ago by Agatha350
gravatar for Fabio Marroni
7.0 years ago by
Fabio Marroni2.6k
Fabio Marroni2.6k wrote:

I think that several tools for calculating LD for SNPs will exclude triallelic SNPs (because I think they make computation somewhat more difficult). If you use the r package genetics, for example (which is the one used by LDheatmap for the calculation, if you let it do the calc) it will remove triallelic SNPs. An option might be to use tools that calculate LD for multiallelic markers. If you only have 4 SNPs it should not be difficult. For example, Arlequin can do that. I guess many other tools can do that.

ADD COMMENTlink written 7.0 years ago by Fabio Marroni2.6k

that's probably it...thought so as well. thanks for the alternative suggestions.

ADD REPLYlink written 6.9 years ago by Agatha350
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: 1275 users visited in the last hour