extend_us: Extend auxiliary variables to include additional empty...

In pierrejacob/montecarlodsm: Gibbs sampler for Dempster-Shafer inference in Categorical distributions

Description

Given a list of auxiliary variables Us, such as generated by gibbs_sampler, add empty categories, i.e. zero counts.

Usage

extend_us(Us, whichbefore, whichnew)

Arguments

Us	a list of auxiliary variables.
whichbefore	a vector of indices indicating location of previous categories
whichnew	a vector of indices indicating location of new (empty) categories

Value

A list with the following entries: #'

• "etas": an array of dimension niterations x K x K. For iteration 'iteration', etasiteration,, is a K x K matrix representing a convex polytope. The feasible set at that iteration, is the set of all vectors theta such that, theta_l/theta_k < etasiteration,k,l for all k,l in 1,...,K.

• "Us": a list of K arrays. For iteration 'iteration', and category 'k' in 1,...,K, if N_k = 0 then Us[k] is NA. If N_k > 0 then Us[k]iteration,, is a matrix with N_k rows, and K columns, representing the auxiliary variables "U"'s generated at that iteration.

pierrejacob/montecarlodsm documentation built on June 16, 2021, 1:06 p.m.