R/dgamma.R

In Rlab: Functions and Datasets Required for ST370 Class

"dgamma" <-

function (x, shape = 1, rate = 1, scale = 1/rate, alpha = shape, beta = scale, log = FALSE) 



{



    if (any(alpha <= 0)) 



            stop("alpha (or shape) must be strictly positive")

            

    if (any(beta <= 0) || any(scale <= 0) || any(rate <= 0)) 



            stop("beta (or scale or rate) must be strictly positive")



     stats::dgamma(x,shape=alpha,scale=beta,log=log)



}