R/ll.invGauss.R

Defines functions ll.invGauss

Documented in ll.invGauss

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

Try the MixedPoisson package in your browser

Any scripts or data that you put into this service are public.

MixedPoisson documentation built on May 2, 2019, 12:40 p.m.