In our model the data are drawn from LogN(mu_i + log(c_ij), tau_i). The prior for tau_i is given as Gamma(thetai*phi, phi). This function draws from the conditional posterior of tau_i. Note that we parameterize with RATE, not SCALE.
postTaui(yij, cij, mui, thetai, phi = 1)
Numeric vector, repeated sampled value of length(yij)
Numeric vector, cycle lengths for a single individual
Positive Integer vector, a sampled vector of length(yij) where the corresponding values in cij indicate a sampled number of TRUE cycles in each cycle length given by yij
Numeric, log of sampled mean of this individual's yijs
Numeric, mean of prior (gamma) distribution on taui
Numeric, rate for Taui prior
Additionally, note that in order to vectorize the remainder of the MCMC algorithm this function returns the sampled value repeated for length(yij)