beta(a, b) lbeta(a, b) gamma(a) lgamma(a) digamma(a) trigamma(a) tetragamma(a) pentagamma(a) choose(n,k) lchoose(n,k)
beta and lbeta return the beta function
and the natural logarithm of the beta function.
The functions gamma and lgamma return the gamma function
and the natural logarithm of the gamma function.
The functions digamma, trigamma, tetragamma and
pentagamma
return the first, second, third and fourth derivatives of the logarithm
of the gamma function.
The functions choose and lchoose
return binomial coefficients and
their logarithms.
x <- seq(.1, 4, length=201); dx <- diff(x)[1] par(mfrow=c(2,3)) for(ch in c("", "l","di","tri","tetra","penta")) { deriv <- nchar(ch)>=2; if(deriv) dy <- diff(y) / dx nm <- paste(ch,"gamma",sep=''); y <- get(nm)(x) plot(x, y, type='l', main=nm, col=1) abline(h=0,col="lightgray") if(deriv) lines(x[-1], dy, col=2, lty=2) }