This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Theta Watterson on SNP

Hi,

I would like to calculate Theta watterson on SNP data. Could you please suggest me any suitable software?

Thanks, YK

snp next-gen r
• 6,012 views • link
updated by Gabriel R. Gabriel R. • written by yp • yp

3 answers

One answer here, another option is dadi.

• 0 views • link
written by h.mon h.mon

Thanks for the response. I have an SNP file in the VCF format. I am looking for a software which can take SNP as an input. The software you suggested, it takes BAM file as an input. I would like to calculate Theta Watterson using all the SNPs.

Thanks, YK

• 0 views • link
written by yp • yp

I would also suggest ANGSD http://www.popgen.dk/angsd/index.php/Thetas,Tajima,Neutrality_tests if you have access to the original BAM files.

• 0 views • link
written by Gabriel R. Gabriel R.

Since you tagged R on your question I would suggest pegas R package. If you want to estimate theta based on the number of segregating sites e.g. SNPs you just have to use: theta.s(s, n) , where s is the number of SNPs and n is the number of sequences.

• 0 views • link
written by GabrielMontenegro GabrielMontenegro

Log in to answer this question.

More posts like this