update_s function

Update class weight vector

This function updates the class weight vector by drawing from its posterior distribution.

update_s(delta, m)

Arguments

• delta: A numeric for the concentration parameter vector rep(delta,C) of the Dirichlet prior for s. Per default, delta = 1. In case of Dirichlet process-based updates of the latent classes, delta = 0.1 per default.
• m: The vector of current class frequencies.

Returns

A vector, a draw from the Dirichlet posterior distribution for s.

Details

Let $m=(m_1,\dots,m_C)$ be the frequencies of $C$ classes. Given the class weight (probability) vector $s=(s_1,\dots,s_C)$, the distribution of $m$ is multinomial and its likelihood is

The conjugate prior $p(s)$ for $s$ is a Dirichlet distribution, which has a density function proportional to

where $\delta = (\delta_1,\dots,\delta_C)$

is the concentration parameter vector. Note that in {RprobitB}, $\delta_1=\dots=\delta_C$. This restriction is necessary because the class number $C$ can change. The posterior distribution of $s$ is proportional to

which in turn is proportional to a Dirichlet distribution with parameters $\delta+m$.

Examples

### number of classes
C <- 4
### current class sizes
m <- sample.int(C)
### concentration parameter for Dirichlet prior (single-valued)
delta <- 1
### updated class weight vector
update_s(delta = 1, m = m)