Adding log pdf methods to the distribution objects

[rasmusab]

I would be cool if the distributions also had a method for calculating the log pdf, like in R you can set log = TRUE for all distributions. The use case for this would be when you want to do maximum likelihood (or Bayesian statistics) but want to use the log likelihood for numerical accuracy.