Publications by xi'an
About Fig. 4 of Fagundes et al. (2007)
Yesterday, we had a meeting of our EMILE network on statistics for population genetics (in Montpellier) and we were discussing our respective recent advances in ABC model choice. One of our colleagues mentioned the constant request (from referees) to include the post-ABC processing devised by Fagundes et al. in their 2007 ABC paper. (This paper c...
4497 sym 18 img
Accepted lack of confidence
I just got the following email from PNAS about our Lack of confidence in ABC model choice. Editor's Remarks to Author: both referees now find the manuscript acceptable for publication as do I. Each suggests small changes which I encourage the authors to make prior to having the manuscript go into production. Congratulations on an interesting...
1177 sym R (279 sym/1 pcs) 18 img
The foundations of Statistics [reply]
Shravan Vasishth has written a response to my review both published on the Statistics Forum. His response is quite straightforward and honest. In particular, he acknowledges not being a statistician and that he “should spend more time studying statistics”. I also understand the authors’ frustration at trying “to recruit several statistici...
2957 sym 18 img
Le Monde puzzle [#28]
The puzzle of last weekend in Le Monde was about finding the absolute rank of x9 when given the relative ranks of x1,….,x8 and the possibility to ask for relative ranks of three numbers at a time. In R terms, this means being able to use > rank(x[-9]) [1] 1 7 4 6 8 3 2 5 > rank(x[1:3]) [1] 1 3 2 or yet being able to sort the first 8 components...
2244 sym R (109 sym/2 pcs) 16 img
Ternary sorting
The last Le Monde puzzle made me wonder about the ternary version of the sorting algorithms, which all seem to be binary (compare x and y, then…). The problem is, given (only) a blackbox procedure that returns the relative order of three arbitrary numbers, how many steps are necessary to sort a series of n nnumbers? The heapsort entry in Wikipe...
1288 sym 16 img
Bayesian Core and loose logs
Jean-Michel (aka Jean-Claude!) Marin came for a few days so that we could make late progress on the revision of our book Bayesian Core towards an Use R! version. In one of the R programs in the mixture chapter, we were getting improbable answers, until we found an R mistake in the shape of > sum(c(1,2,3,log=TRUE)) [1] 7 > sum(c(1,2,3),log=TRUE...
986 sym R (67 sym/1 pcs) 16 img
A slice of infinity
Peng Yu sent me an email about the conditions for convergence of a Gibbs sampler: The following statement mentions convergence. But I’m not familiar what the regularity condition is. “But it is necessary to have a finite probability of moving away from the current state at all times in order to satisfy the regularity conditions on which th...
2480 sym R (637 sym/2 pcs) 16 img
Core not in CiRM
Despite not enjoying this year the optimal environment of CiRM, we are still making good progress on the revision (or the R vision) of Bayesian Core. In the past two days, we went over Chapters 1 (Introduction), 2 (Normal Models), 5 (Capture-Recapture Experiments), and 6 (Mixture Models), with Chapters 3 (Regression), 4 (Generalised Linear Models...
2023 sym 18 img
Le Monde puzzle [#29]
This week, the puzzle from the weekend edition of Le Monde was easy to state: in the sequence (8+17n), is there a 6th power? a 7th? an 8th? If so, give the first occurrence. So I first wrote an R code for a function testing whether an integer is any power: ispower=function(x){ ispo=FALSE logx=log(x) i=trunc(logx/log(2)) while((i>1)&&(!ispo)){ j=t...
1885 sym R (720 sym/5 pcs) 18 img
JSM 2011 [3]
Monday August 01 was the first full day of JSM 2011 and full is the appropriate word to describe the day! It started for me at 7am with a round table run by Marc Suchard on parallel computing (or at 3am if I am considering the time I woke up!). I was rather out of my depth there, given that my link with parallel computing is rather formal, having...
5581 sym 20 img