Hi,
I'm just wondering whether it's deliberate that there aren't any log
probability density functions in the GSL? I've had a quick search
through the archives and haven't seen anything obvious, but if you
could point me in the right direction that would be great!
As a bit of background, I'm using the GSL for doing Bayesian MCMC
inference and the starting point of a chain tends to be somewhere with
(very) low density. The lack of log density functions means that
densities will be truncated to zero before I can take the log of them.
For all of the distributions I've looked at in GSL it's been trivial
to rearrange the maths so it's working with log densities and I'd be
happy to tidy up my hacks into a patch if it would be useful to other
people. Which, as far as I can tell, seems to be anyone doing much
with likelihoods--the product of lots of very small numbers gets
truncated to zero very quickly and working with the log of these is an
easy way of getting sensible results!
Of course, any other suggestions would be appreciated!
Thanks in advance,
Sam