This should not happen. It is either a general bug (not clear whether in the R or Julia code) or a bug arising from different parameterization of the underlying distributions from which Gibbs is sampling (e.g., Gamma, InversePrior using shape and scale vs shape/rate).