R/examples/gamma_params.R

prm <- gamma.params(mu=15.3, sigma=9.3)

prm[1]/prm[2] # mean
sqrt(prm[1]/(prm[2]^2)) # sd

plot(density(rgamma(n=1000, shape=prm[1], rate=prm[2])), main='')
abline(v=15.3, lty=2, col='red')
gilesjohnr/genpatch documentation built on May 12, 2019, 10:50 a.m.