R/gam1.1.R

gam1.1 <-
function (y, lambda) 
{
    gamma <- double(1)
    if (y == 0) {
        gamma <- 1
    }
    else {
        gamma <- 1
        for (i in 0:(y - 1)) {
            gamma <- (gamma * (lambda + i))/(i + 1)
        }
    }
    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.