sampler.dp {SBMSplitMerge}

Dirichlet process sampler

Description

Dirichlet process sampler

Usage

sampler.dp(currsbm, edges, sbmmod, sigma)

Arguments

currsbm

the current state of the sampler

edges

an edges object

sbmmod

an sbmmod model

sigma

random walk parameter for theta

Value

next state of currsbm object

See Also

For full algorithm details see http://doi.org/10.17635/lancaster/thesis/296

Examples

model <- sbmmod(crp(4), param_norm(0,0,1,1,3,3,1,1), edges_norm())
trueSBM <- model$r(100)
Edges <- redges(trueSBM, model$edge)
dp_out <- sampler(Edges, model, 25, "dp", sigma=0.1)