Blog Archives

Research in pair next summer

April 29, 2010
By
Research in pair next summer

Today I received the very good news that our proposal with Jean-Michel Marin to undertake “research in pair” in CIRM, Luminy, a fortnight next summer was accepted! This research centre in Mathematics is a southern and French version of the renowned German centre of Oberwolfach and, while I would have prefered the cool Black Forest

Read more »

The Bernoulli factory

April 22, 2010
By
The Bernoulli factory

A few months ago, Latuszyński, Kosmidis, Papaspiliopoulos and Roberts arXived a paper I should have noticed earlier as its topic is very much related to our paper with Randal Douc on the vanilla Rao-Blackwellisation scheme. It is motivated by the Bernoulli factory problem, which aims at (unbiasedly) estimating f(p) from an iid sequence of Bernoulli

Read more »

Little R == r

April 21, 2010
By

There's big R, the R that I use to do most my work, the environment that makes pretty graphics, et. al. It's like matlab, only cooler. Or more cool. Or less uncool. You can see my prejudices here. Today i discovered little R. It's like big R, only little. Holy shit. Dirk gives...

Read more »

Sudokus more random than random!

April 18, 2010
By
Sudokus more random than random!

Darren Wraith pointed out this column about sudokus to me. It analyses the paper by Newton and De Salvo published in the Proceedings of the Royal Academy of Sciences A that I cannot access from home. The discussion contains this absurd sentence “Sudoku matrices are actually more random than randomly-generated matrices” which shows how mistreated

Read more »

An obscure integral

April 7, 2010
By
An obscure integral

Here is an email from Thomas I received yesterday about a computation in our book Introducing Monte Carlo Methods with R: I’m currently reading your book “Introduction to Monte Carlo Methods with R” and I quite highly appreciate your work. I’m not able to see how the integral on page 74, that describes the marginal

Read more »

Le Monde rank test (corr’d)

April 6, 2010
By
Le Monde rank test (corr’d)

Since my first representation of the rank statistic as paired was incorrect, here is the histogram produced by the simulation perm=sample(1:20) saple=sum(abs(sort(perm)-sort(perm))) when . It is obviously much closer to zero than previously. An interesting change is that the regression of the log-mean on produces > lm(log(memean)~log(enn)) Call: lm(formula = log(memean) ~ log(enn)) Coefficients: (Intercept)    

Read more »

Le Monde rank test (cont’d)

April 5, 2010
By
Le Monde rank test (cont’d)

Following a comment from efrique pointing out that this statistic is called Spearman footrule, I want to clarify the notation in namely (a) that the ranks of and are considered for the whole sample, i.e. instead of being computed separately for the ‘s and the ‘s, and then (b) that the ranks are reordered for

Read more »

Le Monde rank test

April 4, 2010
By
Le Monde rank test

In the puzzle found in Le Monde of this weekend, the mathematical object behind the silly story is defined as a pseudo-Spearman rank correlation test statistic, where the difference between the ranks of the paired random variables and is in absolute value instead of being squared as in the Spearman rank test statistic. I don’t

Read more »

Savage-Dickey [talk]

March 19, 2010
By
Savage-Dickey [talk]

Here are the slides for the Savage-Dickey paradox paper that I gave in San Antonio this morning: (Any suspected coincidence of the first part with earlier talks is for real!) I have tried to spell out as clearly as possible in the second part the issues of version choices that are at the core of

Read more »

Course in San Antonio, Texas

March 18, 2010
By
Course in San Antonio, Texas

Yesterday, I gave my short (3 hours) introduction to computational Bayesian statistics to a group of 25-30 highly motivated students. I managed to cover “only” the first three chapters, as I included some material on Bayes factor approximation and only barely reached Metropolis-Hastings. Here are the slides, modified from the original Bayesian Core slides: (It

Read more »