Takes a vector of collection indices to which individuals (vector elements) were assigned, and returns a Dirichlet random variable generated by adding the prior to the sum of each collection's occurrences, and simulating an alpha from a gamma distribution with this shape parameter.
dirch_from_allocations(C, lambda)
a vector giving different categories of individual (not counts of categories - untabulated)
priors for the categories
The categories are labeled in C from 1 up to n. n is the length of lambda
,
which is a vector of priors. Note that all elements of lambda
must be strictly greater than 0.