Statistics 解释似然Wattersonθ

Statistics 解释似然Wattersonθ,statistics,bioinformatics,Statistics,Bioinformatics,我第一次计算Wattersonθ,现在正在使用。据我所知,该值可以是0到1。然而,在ANGSD上,使用似然法运行它,这是我在一个BAM文件上计算Chr22上的θ时得到的结果: (indexStart,indexStop)(firstPos_withData,lastPos_withData)(WinStart,WinStop) Chr WinCenter tW tP tF tH tL Tajima fuf fud

我第一次计算Wattersonθ,现在正在使用。据我所知,该值可以是0到1。然而,在ANGSD上,使用似然法运行它,这是我在一个BAM文件上计算Chr22上的θ时得到的结果:

(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
所以我感兴趣的是tW(theta Watterson)。。但这个值对我来说毫无意义。我想我必须把它除以NSite,它会给我tW=0.0009671336118。这样做有意义吗

我还有每个站点的结果,日志规模如下:

#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
使用这些值计算总Wtheta是否更好?如果是的话,我是否仍要计算总和并除以总站点数

对不起,这是新来的新手问题。 非常感谢