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(thetaiphi, phi). This function uses a MH step to draw a new sample of phi. Proposal distribution is Gamma(currentPhirhoPhi, rhoPhi). Note that we parameterize with RATE, not SCALE.
postPhi(taui, thetai, currentPhi, rhoPhi = 1000)
Numeric, new draw of phi
Numeric vector, individuals precisions.
Numeric vector. individuals precisions means (estimate)
Previous draw of phi
Proposal rate for gamma distribution that draws proposal for phi, default is 1000.