powered by
Truncated Gamma draws are needed for the standard deviation of the random effects Gibbs conditional.
rtgamma(n, shape, rate, a)
Returns n truncated Gamma, i.e., \(Gam(shape, rate)I(a, infinity)\).
n
Number of samples.
Sampling from a truncated Gamma where \(E[x]=shape/rate\).
This parameter is the inverse of the scale which is an alternative representation for the Gamma distribution.
The truncation point, i.e., \(a<x\).
Gentle J. (2013) Random number generation and Monte Carlo methods. Springer, New York, NY.
set.seed(12) rtgamma(1, 3, 1, 4) rtgamma(1, 3, 1, 4) a=rtgamma(10000, 10, 2, 1) mean(a) min(a)
Run the code above in your browser using DataLab