Question: Tajima'D and low-density SNP data
gravatar for chrys.delord
3.0 years ago by
chrys.delord0 wrote:

Dear all,

I am quite new to analysis of SNP data :-) My case is a bit particular because I only have genotypes from a few number of SNPs (~130) and more importantly, I do NOT have the associated DNA sequences because genotyping was performed through mass-spectrometry (and not sequencing!)

I was wondering if I still could use information from allele frequency spectrum and, for example, compute Tajima's D on my dataset?

Thank you very much for your enlightment!

Have a good day,


snp • 1.3k views
ADD COMMENTlink modified 3.0 years ago by Kevin Blighe69k • written 3.0 years ago by chrys.delord0
gravatar for Kevin Blighe
3.0 years ago by
Kevin Blighe69k
Republic of Ireland
Kevin Blighe69k wrote:

You could first extract sequence surround your SNPs, and there are various ways of doing this. There is a previous Biostars thread: extract sequence around snp

With the sequences, you then can create a DNAbin object using as.DNAbin() (from ape package) in R and run tajima.test() (from pegas package) on this. This will calculate Tajima's D.

ADD COMMENTlink written 3.0 years ago by Kevin Blighe69k
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: 963 users visited in the last hour