Bonjour, je vous explique le contexte et mon problème.
En génomique évolutive, pour calculer le polymorphisme nucléotidique d'un échantillon, on peut utiliser l'approche de Watterson:
watterson 2.PNG
Cette approche prend donc en compte le S soit le nombre de site polymorphes / la taille en pb de la séquence étudiée.
Ensuite on calcul le théta qui est égale à S / alpha
Ma question est là, qu'est-ce que ce alpha? Je crois reconnaitre une somme harmonique.
Du coup cela reviendrait à faire : 1+ 1/2 + 1/3 +1/4 + 1/5 + 1/6 + 1/7 + 1/8 + 1/9 si on en crois la formule :
harmonic.PNG
Seulement cela n'est pas égale à 1.99 comme dans le papier.
Quelqu'un reconnait-il la formule et est capable de calculer ce alpha?
Merci et bonne soirée.
-----