I’m not trying to do optimization. I’m trying to get an updated distribution of observed values of which I don’t have that many by incorporating historical values (that are a proxy for the measured data) which I do have a lot of.
How do you know that kappa and b follow HalfNormal with the parameters you specified? Also, the priors should be based on the data in df_for_prior, not the observed.