Publications by xi'an
astroABC: ABC SMC sampler for cosmological parameter estimation
“…the chosen statistic needs to be a so-called sufficient statistic in that any information about the parameter of interest which is contained in the data, is also contained in the summary statistic.” Elise Jenningsa and Maeve Madigan arXived a paper on a new Python code they developed for implementing ABC-SMC, towards astronomy or rather c...
2547 sym 4 img
an inverse permutation test
A straightforward but probabilistic riddle this week in the Riddler, which is to find the expected order of integer i when the sequence {1,2,…,n} is partitioned at random into two sets, A and B, each of which is then sorted before both sets are merged. For instance, if {1,2,3,4} is divided in A={1,4} and B={2,3}, the order of 2 in {1,4,2,3} is ...
1405 sym 10 img
approximate lasso
Here is a representation of the precision of a kernel density estimate (second axis) against the true value of the density (first axis), which looks like a lasso of sorts, hence the title. I am not sure this tells much, except that the estimated values are close to the true values and that a given value of f(x) is associated with two different es...
855 sym 6 img
Journal of Open Source Software
A week ago, I received a request for refereeing a paper for the Journal of Open Source Software, which I have never seen (or heard of) before. The concept is quite interesting with a scope much broader than statistical computing (as I do not know anyone in the board and no-one there seems affiliated with a Statistics department). Papers are very ...
1232 sym 6 img
importance sampling by kernel smoothing [experiment]
Following my earlier post on Delyon and Portier’s proposal to replacing the true importance distribution ƒ with a leave-one-out (!) kernel estimate in the importance sampling estimator, I ran a simple one-dimensional experiment to compare the performances of the traditional method with this alternative. The true distribution is a N(0,½) with...
2070 sym R (794 sym/1 pcs) 6 img
grim knight [a riddle]
The Riddler of this week had a riddle that is a variation of the knight tour problem, namely “…how long is the longest path a knight can travel on a standard 8-by-8 chessboard without letting the path intersect itself?” the riddle being then one of a self-avoiding random walk [kind]… As I could not get back to sleep last night, I spent a ...
2419 sym R (174 sym/1 pcs) 10 img
tractable Bayesian variable selection: beyond normality
David Rossell and Francisco Rubio (both from Warwick) arXived a month ago a paper on non-normal variable selection. They use two-piece error models that preserve manageable inference and allow for simple computational algorithms, but also characterise the behaviour of the resulting variable selection process under model misspecification. Interest...
3934 sym 6 img
an attempt at EP-ABC from scratch, nothing more… [except for a few bugs]
Following a request from one of the reviewers of our chapter Likelihood-free model choice, I tried to run EP-ABC on a toy problem and to compare it with the outcome of a random forest ABC. Literally starting from scratch, namely from the description found in Simon and Nicolas’ JASA paper. To run my test, I chose as my modelled data an exponen...
3083 sym R (1459 sym/7 pcs) 6 img
ratio-of-uniforms
One approach to random number generation that had always intrigued me is Kinderman and Monahan’s (1977) ratio-of-uniform method. The method is based on the result that the uniform distribution on the set A of (u,v)’s in R⁺xX such that 0≤u²≤ƒ(v/u) induces the distribution with density proportional to ƒ on V/U. Hence the name. The proo...
2223 sym 10 img
a grim knight [cont’d]
As discussed in the previous entry, there are two interpretations to this question from The Riddler: “…how long is the longest path a knight can travel on a standard 8-by-8 board without letting the path intersect itself?” as to what constitutes a path. As a (terrible) chess player, I would opt for the version on the previous post, the kn...
1839 sym 14 img