gravatar for yolek64754

2 hours ago by

Hi everyone,

I have been computing Watterson theta for the first time and I am using ANGSD . For my understanding this value can be from 0 to 1. However on ANGSD use likelihood to run it and this is the results I have when I compute theta on Chr22 on one bamfile:

(indexStart,indexStop)(firstPos_withData,lastPos_withData)(WinStart,WinStop)   Chr     WinCenter       tW      tP      tF      tH      tL      Tajima  fuf     fud     fayh    zeng    nSites
(0,32860927)(16050032,51244307)(0,51244307)     22      25622153        31780.907015    31780.907015    31780.907015    31780.907015    31780.907015    -nan    -nan    -nan    -nan    -nan    32860927

So what I am interested in is tW (theta Watterson).. but this value means no sense to me. I think I do have to divide it by the nSites which would give me tW=0.0009671336118. Does it make sense to do that ?

I also have the results per sites on a log scale that looks like that:

#Chromo Pos Watterson   Pairwise    thetaSingleton  thetaH  thetaL
22  16050032    -7.634670   -7.634670   -7.634670   -7.634670   -7.634670
22  16050033    -7.634670   -7.634670   -7.634670   -7.634670   -7.634670
22  16050034    -7.634670   -7.634670   -7.634670   -7.634670   -7.634670
22  16050035    -7.634670   -7.634670   -7.634670   -7.634670   -7.634670
22  16050036    -7.634670   -7.634670   -7.634670   -7.634670   -7.634670

Is it better do compute the overall Wtheta using those values ? If yes, should I still do the sum and dividing them by the total number of sites ?

Sorry for the newbie questions, new on this field.
Thanks a lot

link

modified 1 hour ago

written
2 hours ago
by

yolek6475410



Source link