## 10w2170, Banff [2]

**O**ver the two days of the **Hierarchical Bayesian Methods in Ecology** workshop, we managed to cover normal models, testing, regression, Gibbs sampling, generalised linear models, Metropolis-Hastings algorithms and of course a fair dose of hierarchical modelling. At the end of the Saturday marathon session, we spent one and half discussing some models studied by the participants, which were obviously too complex to be solved on the spot but well-defined so that we could work on MCMC implementation and analysis. And on Sunday morning, a good example of Poisson regression proposed by Devin Goodman led to an exciting on-line programming of a random effect generalised model, with the lucky occurrence of detectable identifiability issues that we could play with… I am impressed at the resilience of the audience given the gruesome pace I pursued over those two days, covering the five first chapters of * Bayesian Core*, all the way to the mixtures! In retrospect, I think I need to improve my coverage of testing as the noninformative case presumably sounded messy. And unconvincing. I also fear the material on hierarchical models was not sufficiently developed. But, overall, the workshop provided a wonderful opportunity to exchange with bright PhD students from Ecology and Forestry about their models and (hierarchical) Bayesian modelling.

**A**s an aside, I also noticed that several of the participants from the University of Alberta were (considering) using data cloning, which is a MCMC method for computing maximum likelihood estimates that was locally developed by Subhash Lele (University of Alberta). The method is based on using MCMC data augmentation methods to optimise likelihood functions by considering a larger and larger number of replicates of the missing data. Data cloning is thus essentially a user-friendly simulated annealing technique for missing data models and one of the multiple replicas of what I called the prior feedback technique in 1991. The idea of prior feedback (also exposed in our 1996 JASA paper with Gene Hwang and in * Monte Carlo Statistical Methods*) is that when the likelihood is raised to a high power, any prior distribution with an unrestricted support leads to a “posterior” distribution concentrated around the MLE. (I actually got this idea from reading the 1991 Read Paper by Murray Aitkin on posterior Bayes factors. The following discussion was quite critical, including Dennis Lindley’s now famous “

*One hardly advances the respect with which statisticians are held in society by making such declarations”*!, but instead of solving the problem with improper priors in Bayesian testing it suggested to me a computational technique to compute MLE’s…) With Arnaud Doucet and Simon Godsill, we developed in 2002 an algorithm for missing data and state-space models called SAME (for state augmentation for marginal estimation) that got published in

*Statistics and Computing*but resurfaces periodically in the literature, the latest avatar being data cloning… Judging from a web search, the data cloning approach is getting increasingly popular in Ecology, including an R package named

**and developed by Peter Solymos.**

*dcr*
May 11, 2011 at 7:27 am

[…] June 11, is a tutorial session given by the organisers, as usual, and I will give a short talk on computational methods (Monte Carlo, MCMC and ABC) for Bayesian inference. The regular talks start on Sunday, June 12, and […]

September 25, 2010 at 12:10 am

[…] the third movie is not out yet in the US, since the second one was played in a movie theatre in Banff last week… It is nonetheless possible to see it in Swedish or in […]

September 22, 2010 at 12:16 am

[…] some discussions I had last week at Banff about data cloning, I re-read the 2007 “Data cloning” paper in Ecology Letters […]

September 16, 2010 at 9:58 am

[…] ridge too many… On Monday morning, Devin and I left Banff very early to try one more scramble before my plane (as the pouring rain of the previous night had […]