Publications by xi'an
Imperial postdoc in Bayesian nonparametrics
Here is another announcement for a post-doctoral position in London (UK) to work with Sarah Filippi. In the Department of Mathematics at Imperial College London. (More details on the site or in this document. Hopefully, the salary is sufficient for staying in London, if not in South Kensington!) The post holder will work on developing a novel Bay...
1583 sym 2 img
the riddle of the stands
The simple riddle of last week on The Riddler, about the minimum number of urinals needed for n men to pee if the occupation rule is to stay as far as possible from anyone there and never to stand next to another man, is quickly solved by an R code: ocupee=function(M){ ok=rep(0,M) ok[1]=ok[M]=1 ok[trunc((1+M/2))]=1 while (max(diff((1:M)...
1068 sym R (233 sym/1 pcs) 4 img
maximal spacing around order statistics
The riddle from the Riddler for the coming weeks is extremely simple to express in mathematical terms, as it summarises into characterising the distribution of when the n-sample is made of iid Normal variates. I however had a hard time finding a result connected with this quantity since most available characterisations are for either Uniform or ...
2274 sym 10 img
Le Monde puzzle [#1051]
A combinatoric Le Monde mathematical puzzle of limited size: When the only allowed move is to switch two balls from adjacent boxes, what is the minimal number of moves to return all balls in the above picture to their respective boxes? Same question with six boxes and 12 balls. The question is rather interesting to code as I decided to use recur...
1798 sym R (1155 sym/3 pcs) 4 img
maximal spacing around order statistics [#2]
The proposed solution of the riddle from the Riddler discussed here a few weeks ago is rather approximative, in that the distribution of when the n-sample is made of iid Normal variates is (a) replaced with the distribution of one arbitrary minimum and (b) the distribution of the minimum is based on an assumption of independence between the abso...
1049 sym 4 img
a chain of collapses
A quick riddler resolution during a committee meeting (!) of a short riddle: 36 houses stand in a row and collapse at times t=1,2,..,36. In addition, once a house collapses, the neighbours if still standing collapse at the next time unit. What are the shortest and longest lifespans of this row? Since a house with index i would collapse on its own...
1365 sym 2 img
Le Monde puzzle [#1053]
An easy arithmetic Le Monde mathematical puzzle again: If coins come in units of 1, x, and y, what is the optimal value of (x,y) that minimises the number of coins representing an arbitrary price between 1 and 149? If the number of units is now four, what is the optimal choice? The first question is fairly easy to code coinz <- function(x,y)...
1793 sym R (498 sym/3 pcs) 2 img
Bayesian GANs [#2]
As an illustration of the lack of convergence of the Gibbs sampler applied to the two “conditionals” defined in the Bayesian GANs paper discussed yesterday, I took the simplest possible example of a Normal mean generative model (one parameter) with a logistic discriminator (one parameter) and implemented the scheme (during an ISBA 2018 sessio...
1291 sym 6 img
hitting a wall
Once in a while, or a wee bit more frequently (!), it proves impossible to communicate with a contributor of a question on X validated. A recent instance was about simulating from a multivariate kernel density estimate where the kernel terms at x¹,x²,… are Gaussian kernels applied to the inverses of the norms |x-x¹|, |x-x²|,… rather than ...
1986 sym 2 img
seven digit addition
Another quick riddle from the riddler: solve the equation EXMREEK + EHKREKK = ?K?H?X?E which involves every digit between 0 and 9. While the puzzle can be unravelled by considering first E and K, which must be equal to 6 and 3, a simple R code also leads to the conclusion isok <- function(a,b){ s=as.numeric(unlist(strsplit(as.character(sum(10^(...
693 sym R (647 sym/1 pcs)