dma {SBMSplitMerge}

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

Details

This model posits:

kappa-1 ~ Pois(delta)

omega|kappa, gamma ~ Dirichlet(gamma)

Z_i|omega ~ Multinomial(omega) for i=1 .. n

Value

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

Examples

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