sufficientStatistics_Weighted.HDP: Weighted sufficient statistics of a "HDP" object

Description

For following model structure: $G | g a m m a \sim D P (g a m m a, U)$ $p i_{j} | G, a l p h a \sim D P (a l p h a, G), j = 1 : J$ $z | p i_{j} \sim C a t e g o r i c a l (p i_{j})$ $k | z, G \sim C a t e g o r i c a l (G), if z is a sample from the base measure G$ $t h e t a_{k} | p s i \sim H 0 (p s i)$ $x | t h e t a_{k}, k \sim F (t h e t a_{k})$ where DP(gamma,U) is a Dirichlet Process on positive integers, gamma is the "concentration parameter", U is the "base measure" of this Dirichlet process, U is an uniform distribution on all positive integers. DP(alpha,G) is a Dirichlet Process on integers with concentration parameter alpha and base measure G. The choice of F() and H0() can be described by an arbitrary "BasicBayesian" object such as "GaussianGaussian","GaussianInvWishart","GaussianNIW", "GaussianNIG", "CatDirichlet", and "CatDP". See ?BasicBayesian for definition of "BasicBayesian" objects, and see for example ?GaussianGaussian for specific "BasicBayesian" instances. As a summary, An "HDP" object is simply a combination of a "CatHDP" object (see ?CatHDP) and an object of any "BasicBayesian" type. In the case of HDP, z and k can only be positive integers. The sufficient statistics of a set of samples x in a "HDP" object is the same sufficient statistics of the "BasicBayesian" inside the "HDP", see examples.

Usage

# S3 method for HDP
sufficientStatistics_Weighted(obj, x, w, ...)

Arguments

obj

A "HDP" object.

Random samples of the "BasicBayesian" object.

numeric, sample weights.

...

further arguments passed to the corresponding sufficientStatistics method of the "BasicBayesian" object.

Value

Return the sufficient statistics of the corresponding BasicBayesian type, see examples.

References

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

Examples

Run this code

# NOT RUN {
## a HDP with Gaussian NIW observations
obj1 <- HDP(gamma=list(gamma=1,alpha=1,j=2,
                       H0aF="GaussianNIW",
                       parH0=list(m=0,k=1,v=2,S=1)))
## a HDP with Categorical-Dirichlet observations
obj2 <- HDP(gamma=list(gamma=1,alpha=1,j=2,
            H0aF="CatDirichlet",
            parH0=list(alpha=1,uniqueLabels=letters[1:3])))
x1 <- rnorm(100)
x2 <- sample(letters[1:3],100,replace = TRUE)
w <- runif(100)
sufficientStatistics_Weighted(obj = obj1,x=x1,w=w,foreach = FALSE)
sufficientStatistics_Weighted(obj = obj2,x=x2,w=w,foreach = FALSE)
sufficientStatistics_Weighted(obj = obj1,x=x1,w=w,foreach = TRUE)
# }

Run the code above in your browser using DataLab

Get 50% off unlimited learning