Maggie (and other Stats PDLers) - I ran into some numerical stability issues with the Poisson probability mass function. The Poisson distribution is mathematically well-behaved, but is not numerically stable. Fortunately, I was able to use Stirling's approximation to increase the range of accuracy and the allowed input range to effectively work for anything.
I've issued a pull request. You can check it out and comment here: https://github.com/maggiexyz/PDL-Stats/pull/1 David P. S. Maggie's recent switch to Github didn't spur my interest in PDL::Stats (blame my research for that), but it sure made it a heck of a lot easier for me to manage this proposed modification. :-) -- "Debugging is twice as hard as writing the code in the first place. Therefore, if you write the code as cleverly as possible, you are, by definition, not smart enough to debug it." -- Brian Kernighan
_______________________________________________ Perldl mailing list [email protected] http://mailman.jach.hawaii.edu/mailman/listinfo/perldl
