dpoilog function

Returns the density for the Poisson lognormal distribution with parameters mu and sig