[Read more...]

[Read more...]

TOMACS (ACM Transactions on Modeling and Computer Simulation) is launching a call for paper submission. The special topic is Monte Carlo Methods in Statistics and Arnaud Doucet and myself are the special issue editors. Here are the details.: Over the last two decades Monte Carlo methods have attracted much attention ... [Read more...]

A casualty of cut-and-paste in Chapter 3 of Introducing Monte Carlo Methods with R. Brad McNeney from Simon Fraser sent me a nice email about the end of Example 3.6 missing a marginal estimate. Indeed, it does. And it should have been obvious from the “estimates” we derived, 19 and 16, which are not ... [Read more...]

__ printSudoku(z) +-------+-------+-------+ |   9   |       | 7   5 | |     6 |       |   9   | | 4 5 3 | 1 7   | 2 8   | +-------+-------+-------+ |     5 |     7 |   6   | | 1   9 | 6 8   |       | |   8   |   3   |     1 | +-------+-------+-------+ | 7   2 | [...] [Read more...]

I just threw a post up on Revolutions, which got a lot longer than I planned. And got me thinking. And reading (see refs in previous post). And trying. Turns out that it was way easier than I thought! The problem:From the blog post: " OpenSSH is now available on ... [Read more...]

Statistics and Computing has received several papers on ABC and plans to make a special ABC issue out of these. All submissions prior to June 2011 that will be accepted will be published in this special issue. The special issue is identified as an article type on the on-line page. In ... [Read more...]

On Friday, a 2008 paper on Stochastic Approximation and Newton’s Estimate of a Mixing Distribution by Ryan Martin and J.K. Ghosh was posted on arXiv. (I do not really see why it took so long to post on arXiv a 2008 Statistical Science paper but given that it is not ... [Read more...]

Here is an email sent by one of my students a few days ago: Do you know how to integrate a function with an  “if”? For instance: __X=rnorm(100) __Femp=function(x){ +   return(sum(Xintegrate(Femp,0,1)$value does not work. My reply was that the fundamental reason it does not ... [Read more...]

This year useR! conference will take place in Warwick, on August 16-18.  It is being organised by the department of Statistics and funded by CRiSM and Revolution Analytics (providers of the R tee-shirt!). I wish I could attend but mid-August is usually associated with genuine (post-JSM) family vacations. Filed under: ... [Read more...]

François Perron is visiting me for two months from Montréal and, following a discussion about the parallel implementation of MCMC algorithms—to which he also contributed with Yves Atchadé in 2005—, he remarked that a deterministic choice of permutations with the maximal contrast should do better than random or ... [Read more...]

A simple challenge in Le Monde this week: find the group of four primes such that any sum of three terms in the group is prime and the overall sum is minimised. Here is a quick exploration by simulation, using the schoolmath package (with its imperfections): A=primes(start=1,end=53)[... [Read more...]

After the very exciting and I think quite successful ABC in Paris meeting two years ago, Michael Stumpf from Imperial College London suggested a second edition in London along the same lines. Michael kindly associated me with the planning of this meeting. It is (logically) called ABC in London (or ... [Read more...]

This is the 1000th post on the ‘Og! Here are the entries that have had above 1000 views (not viewers) so far: In{s}a(ne)!! 5,353 “simply start over and build something better” 4,345 Julien on R shortcomings 1,966 Sudoku via simulated annealing 1,762 Of black swans and bleak prospects 1,462 Do we need an ... [Read more...]

The paper “Exact sampling for intractable probability distributions via a Bernoulli factory” by James Flegal and Radu Herbei got posted on arXiv without me noticing, presumably because it came out just between Larry Brown’s conference in Philadelphia and my skiing vacations! I became aware of it only yesterday and ... [Read more...]

We have now received reports back from JCGS for our parallel MCMC paper and they all are very nice and supportive! The reviewers essentially all like the Rao-Blackwellisation concept we developed in the paper and ask for additions towards a more concrete feeling for the practical consequences of the method. ... [Read more...]

Another Sudoku-like puzzle from the weekend edition of Le Monde. The object it starts with is a 9×9 table where each entry is an integer and where neighbours take adjacent values. (Neighbours are defined as north, west, south and east of an entry.) The question is about whether or not it ... [Read more...]

An ‘Og reader. Emmanuel Charpentier, sent me the following email about model choice: I read with great interest your critique of Peter Congdon’s 2006 paper (CSDA, 50(2):346-357) proposing a method of estimation of posterior model probabilities based on improper distributions for parameters not present in the model inder examination, as ... [Read more...]

I am in Montpellier today for a general meeting of our EMILE ANR grant bi-yearly meeting, discussing ABC model choice and new applications of simulation methods in population genetics, including an upgrade of DIYABC into the 2.0 version. (The above is a picture comparing “true” Bayes factors with ABC versions, commented ... [Read more...]

The second day at the Latent Gaussian Models workshop in Zürich was equally interesting. Among the morning talks, let me mention Daniel Bové who gave a talk connected with the hyper-g prior paper he wrote with Leo Held (commented in an earlier post) and the duo of Janine Illian ... [Read more...]