SBMSplitMerge (version 1.1.1)

dma: Dirichlet Multinomial Allocation

Description

A blockmod for Dirichlet Multinomial Allocation (DMA)

Usage

dma(gamma, delta)

Arguments

gamma

parameter for Dirichlet component

delta

parameter for Poison component

Value

a block model representing a dma(gamma, delta) distribution

Details

This model posits: $$kappa-1 ~ Pois(delta)$$ $$omega|kappa, gamma ~ Dirichlet(gamma)$$ $$Z_i|omega ~ Multinomial(omega) for i=1 .. n$$

Examples

# NOT RUN {
## simulate from a DMA(2, 5) prior
## This models the `number of blocks-1` as Poisson(5)
## and block assignments as Dirichlet-Multinomial(2, 2, ...)
m <- dma(2, 5)
print(m)
m$r(10)
# }