heterogeneity_param_prior function

Determine the prior distribution for the heterogeneity parameter