x is age grouped in five years, it ranges from 0 to 110 (eg. 0,5,10,...110). B ranges from 0 to 1. The exponent is (x+B) not (x-B) so it is always defined. It rarely happens in the models that p is greater than 1 but anyway I put a condition that if it is greater than 1 then substitute with 0.999. I am trying to approximate a posterior distribution of the parameters (A to H) through slice sampling or Metropolis Hastings method. but the problem here is that the deaths is binomial and the prior distributions are not for the probability of success of the binomial, they are for the parameters that are in Heligman Pollard model (the p equation). Sorry for replying late and thank you for your time.
... View more