In our model the data are drawn from LogN(mu_i + log(c_ij), tau_i) The prior for c_ij is a categorical prior with category probabilities pi1, ..., pik, and c_ij can take values 1, ..., k where k is the length of pi. This function samples from the full conditional posterior of all c_ijs, given vectors of equal length yijs, muis, tauis
postCij(yijs, pi, muis, tauis)
Integer vector
Numeric Vector, cycle lengths
Numeric vector, must sum to 1. Sampled probabilities for c_ijs
Numeric vector, log of sampled mean for all individuals yijs
Numeric vector > 0, sampled precision for all individuals yijs