R/pltgamma.R

In RCASPAR: A package for survival time prediction based on a piecewise baseline hazard Cox regression model.

#### Plot the gamma distribution with the chosen parameters
pltgamma <- function (a, b) {
    x <- seq(0,20,000.1)
    y <- x ^ (a - 1) * exp((-x/b)) / ((b ^ a) * gamma(a)) 
    plot (x, y, type="s")
    }