R/invGauss.density.R

In MixedPoisson: Mixed Poisson Models

Defines functions invGauss.density

Documented in invGauss.density

invGauss.density <-
function(theta, delta) {
theta <- mpfr(theta, 99)
y=(delta*exp(delta^2)*theta^(-3/2)*exp(-((delta^2)/2)*(1/theta + theta)))/(sqrt(2*pi))
return(asNumeric(y))
}

Try the MixedPoisson package in your browser