Archive for compatible conditional distributions

Back from Oxford

Posted in Books, pictures, Statistics, Travel, University life with tags , , , , , , , , , on February 25, 2012 by xi'an

Several interesting questions were raised during my seminar talk at Oxford. First, David Cox suggested I looked at the collection of the two p-values in the Poisson and geometric cases to check whether or not they could point out at a disagreement. (I am however unsure at how the p-values should be computed in this case, maybe as a likelihood-ratio test…) Chris Holmes asked about what happened to the ABC Bayes factor when both models were wrong. I had not thought of this earlier and will look into it: my first impression is that there is no reason for the same model to be chosen. It depends on the relative tail behaviour of the distribution of the summary statistics under both models… Stephen Lauritzen mentioned prior to the seminar a highly relevant book by a Copenhagen mathematician, whose definition of conditional densities was perfectly suited for constructing a convergence proof for ABC (to be incorporated in my Roma slides, if feasible). During the talk, he also pointed out at other (counter)examples of models where sufficient statistics remained sufficient across models: e.g., contingency tables with pairwise interactions. Arnaud Doucet got back to the Potts model (counter)example to stress that we needed perfect sampling to make it work and that our MCMC alternative could be adding another level of approximation to the process, which is quite right!

On the less academic side, I was in Oxford only for a short while, being due back in Paris for a presentation of our Statistics Master: I stil managed a short run in the morning in a nearby park where I saw a heron (blurred above!), as well as hints of the coming Spring (left) but I wish I had had more time (and indications!) to run along the rowers as I did in Cambridge. (I also wish I had had time to visit Tolkien’s favourite pub! Although I had a beer at the Lamb and Flag, which served as a meeting place for the later members of the Inklings…)

workshop in Columbia [day 3]

Posted in pictures, R, Running, Statistics, Travel, University life with tags , , , , , , , , , , on September 27, 2011 by xi'an

Although this was only a half-day of talks, the third day of the workshop was equally thought-challenging and diverse.  (I managed to miss the ten first minutes by taking a Line 3 train to 125th street, having overlooked the earlier split from Line 1… Crossing south Harlem on a Sunday morning is a fairly mild experience though.) Jean-Marc Azaïs gave a personal recollection on the work of Mario Wschebor, who passed away a week ago and should have attended the workshop. Nan Chen talked about the Monte Carlo approximation of quantities of the form

\mathbb{E}[f(\mathbb{E}[Y|X])]

which is a problem when f is non linear. This reminded me (and others) of the Bernoulli factory and of the similar trick we use in the vanilla Rao-Blackwellisation paper with Randal Douc. However, the approach was different in that the authors relied on a nested simulation algorithm that did not adapt against f. And did not account for the shape of f. Peter Glynn, while also involved in the above, delivered a talk on the initial transient that showed possibilities for MCMC convergence assessment (even though this is a much less active area than earlier). And, as a fitting conclusion, the conference organiser, Jingchen Liu gave a talk on non-compatible conditionals he and Andrew are using to handle massively-missing datasets. It reminded me of Hobert and Casella (1996, JASA) of course and also of discussions we had in Paris with Andrew and Nicolas. Looking forward to the paper (as I have missed some points about the difference between true and operational models)! Overall, this was thus a terrific workshop (I just wish I could have been able to sleep one hour more each night to be more alert during all talks!) and a fantastic if intense schedule fitting the start of the semester and of teaching (even though Robin had to teach my first R class in English on Friday). I also discovered that several of the participants were attending the Winter Simulation Conference later this year, hence another opportunity to discuss simulation strategies together.