Inferring extinction date of a species using non-homogeneous Poisson processes with a change-point The code uses utility programs from a book which has the following citation: John K. Kruschke's Doing Bayesian Data Analysis: A Tutorial with R and BUGS (1e) / A Tutorial with R, JAGS, and Stan The homepage for the book is here .