xavierdidelot / TransPhylo

Reconstruction of transmission trees using genomic data
http://xavierdidelot.github.io/TransPhylo/
GNU General Public License v2.0
60 stars 22 forks source link

Precision fixes #6

Closed gtonkinhill closed 6 years ago

gtonkinhill commented 6 years ago

Updated probTTree.cpp to keep probabilities in log space. Adjusted time interval to 0.01, and added exceptions to bail out if the function calculating wbar0 looks like it is going to run forever (lines 129 & 211). Added in behaviour to match line 67 of the original function probTTree.R (lines 169-175).

xavierdidelot commented 6 years ago

Thanks Gerry!