bbricks (version 0.1.1)

MAP.CatDP: Maximum A Posteriori(MAP) estimate of a "CatDP" object

Description

Generate the MAP estimate of "pi" in following model structure: pi|alpha ~ DP(alpha,U) x|pi ~ Categorical(pi) where DP(alpha,U) is a Dirichlet Process on positive integers, alpha is the "concentration parameter" of the Dirichlet Process, U is the "base measure" of this Dirichlet process, it is an uniform distribution on all positive integers. In the case of CatDP, x can only be positive integers. The model structure and prior parameters are stored in a "CatDP" object. MAP is pi_MAP = argmax_pi p(pi|alpha,x).

Usage

# S3 method for CatDP
MAP(obj, ...)

Arguments

obj

A "CatDP" object.

...

Additional arguments to be passed to other inherited types.

Value

numeric.

References

Teh, Yee W., et al. "Sharing clusters among related groups: Hierarchical Dirichlet processes." Advances in neural information processing systems. 2005.

See Also

CatDP

Examples

# NOT RUN {
x <- sample(1L:10L,size = 40,replace = TRUE)
obj <- CatDP()
posterior(obj = obj,ss = x)
MAP(obj)
# }