param_gamma {SBMSplitMerge}

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)

Details

This model represents a prior on theta with:

theta_0 ~ Gamma(a0,a1)

theta_k ~ Gamma(b0,b1)

for k = 1 ... kappa

Value

a parammod

Examples

## 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