My first package, zipcode, is now available on CRAN. It contains the CivicSpace database of 43,191 U.S. zip codes.

Another exciting day at Adap’skiii!!! Yves Atchadé presented a very recent work on the fundamental issue of estimating the asymptotic variance estimation for adaptive MCMC algorithms, with an intriguing experimental observation that a non-converging bandwidth with rate 1/n was providing better coverage than the converging rate. (I always found the issue of estimating the asymptotic

While the above heading sounds like a title in reverse, its words are in the “correct” order in that our paper with George Casella, A Short History of Markov Chain Monte Carlo, has been accepted for publication by Statistical Science. This publication may sound weird when considering that the paper is also scheduled to appear

The last puzzle of the year in Le Monde reads as follows (as far as I understand its wording!): Iter(n,x,y) is the function Iter=function(n,x,y){ if (n==1){ output=trunc(y/10)+x*(y%%10) }else{ output=Iter(n-1,x,Iter(1,x,y))} return output } Find the seven-digit number z such that Iter(6,1,z)=12, Iter(6,2,z)=19, Iter(6,3,z)=29, and Iter(6,-1,z)=Iter(6,-2,z)=Iter(6,-3,z)=0. Obviously, the brute-force solution of listing all 90 million seven digit

The blog year started in August and consists of 30-something posts. Here is a summary. Quant concepts backtesting: Backtesting — almost wordless cointegration: American TV does cointegration efficient frontier: Anomalies meet volatility implied alpha: Implied alpha — almost wordless portfolio theory: Ancient portfolio theory random walk: The tightrope of the random walk returns: A tale … Continue reading...

Following Ashley’s latest comments on Chapter 5 of Introducing Monte Carlo Methods with R, I realised Example 5.5 was totally off-the-mark! Not only the representation of the likelihood should have used prod instead of mean, not only the constant should call the val argument of integrate, not only integrate uses lower and upper rather than

Ashley put the following comment on Chapter 5 of Introducing Monte Carlo Methods with R”: I am reading chapter 5. I try to reproduced the result on page 128. The R codes don’t work on my laptop. When I try to run the following codes on page 128 > for (i in 1:(nlm(like,sta)$it)){ + mmu=rbind(mmu,nlm(like,sta,iter=i)$est)}