View source: R/invgamma_couplings.R
Sample from maximally coupled inverse gamma with given parameters alpha1, alpha2, beta1, beta2, where the parametrization is that the log-density of IG(alpha, beta) at x is
alpha * log(beta) - lgamma(alpha) - (alpha+1) * log(x) - beta / x
1 | rigamma_coupled(alpha1, alpha2, beta1, beta2)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.