SBMSplitMerge (version 1.1.1)

param_gamma: Gamma parameter model

Description

A parammod with gamma-distributed parameters

Usage

param_gamma(a0, a1, b0, b1)

Arguments

a0

theta_0 ~ Gamma(a0,a1)

a1

theta_0 ~ Gamma(a0,a1)

b0

theta_k ~ Gamma(b0,b1)

b1

theta_k ~ Gamma(b0,b1)

Value

a parammod

Details

This model represents a prior on theta with: $$theta_0 ~ Gamma(a0,a1)$$ $$theta_k ~ Gamma(b0,b1)$$ for k = 1 ... kappa

Examples

# NOT RUN {
## theta0 ~ Gamma(1,1); thetak ~ Gamma(5,5)
pg <- param_gamma(1,1,5,5)
pg$r(5) ## a draw with 5 within-block parameters
# }