R/gam1.2.R

gam1.2 <-
function (y, lambda) 
{
    gamma <- double(1)
    gamma <- exp((y + lambda - 1/2) * log(y + lambda) + (1/(12 * 
        (y + lambda))) - (y + 1/2) * log(y) - (lambda - 1/2) * 
        log(lambda) - (1/2) * log(2 * pi) - 1/(12 * lambda))
    gamma
}

Try the poistweedie package in your browser

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

poistweedie documentation built on May 2, 2019, 3:25 a.m.