Articles by xi'an

Random generators for parallel processing

October 28, 2010 | xi'an

Given the growing interest in parallel processing through GPUs or multiple processors, there is a clear need for a proper use of (uniform) random number generators in this environment. We were discussing the issue yesterday with Jean-Michel Marin and briefly looked at a few solutions: given p parallel streams/threads/...
[Read more...]

Le Monde puzzle [42]

October 24, 2010 | xi'an

An interesting suduko-like puzzle for this week puzzle in Le Monde thi A 10×10 grid is filled by a random permutation of {0,…,99}. The 4 largest figures in each row are coloured in yellow and the 4 largest values in each column are coloured in red. What is the range of the number of ... [Read more...]

abc

October 21, 2010 | xi'an

Michael Blum and Olivier François, along with Katalin Csillery, just released an R package entitled abc. (I am surprised the name was not already registered!) Its aim is obviously to implement ABC approximations for Bayesian inference: Description The ’abc’ package provides various functions for parameter estimation and model selection ... [Read more...]

Coincidence in lotteries

October 19, 2010 | xi'an

Last weekend, my friend and coauthor Jean-Michel Marin was interviewed (as Jean-Claude Marin, sic!) by a national radio about the probability of the replication of a draw on the Israeli Lottery. Twice the same series of numbers appeared within a month. This lotery operates on a principle of 6/37 + 1/8: 6 numbers are ... [Read more...]

Adap’skiii program

October 19, 2010 | xi'an

We have just posted the (mostly definitive) program for Adap’skii, January 3-4, The Canyons, Utah. This is taking place just before and as a satellite of the larger MCMSki III conference, January 4-7, same location. The registration for the conference and for lodging is available through the  MCMCSki III ... [Read more...]

Le Monde puzzle [41]

October 17, 2010 | xi'an

The current puzzle in Le Monde this week is again about prime numbers: The control key on a credit card is an integer η(a) associated with the card number a such that, if the card number is c=ab, its key η(c) satisfies η(c)=η(a)+η(b)-1. There is only ... [Read more...]

Postdoc position in computational Bayesian statistics

October 14, 2010 | xi'an

Here is an announcement I received that should interest potential postdocs (willing to come to Paris). The location is on the Orsay campus, south of Paris. In the framework of the ANR-funded Metacoli project which aims at identifying the metabolic underpinnings of the lifestyle diversity in the E. coli species, ... [Read more...]

Parallel processing of independent Metropolis-Hastings algorithms

October 11, 2010 | xi'an

With Pierre Jacob, my PhD student, and Murray Smith, from National Institute of Water and Atmospheric Research, Wellington, who actually started us on this project at the last and latest Valencia meeting, we have completed a paper on using parallel computing in independent Metropolis-Hastings algorithms. The paper is arXived and ...
[Read more...]

Le Monde puzzle [40]

October 10, 2010 | xi'an

The puzzle in Le Monde this week[end] is called the “square five” (sic!): Two players each have twenty-five cards with five times each of the digits 1,2,3,4,5. They alternate putting one card on top of the pile, except that they can instead take an arbitrary number of consecutive cards from ... [Read more...]

Pre-ordinary meeting

October 7, 2010 | xi'an

Those are the slides for the (basic) introduction of the paper by Mark Girolami and Ben Calderhead at the RSS next week. Not to be confused with my comments on the paper. Filed under: R, Statistics, Travel, University life Tagged: Hamiltonian, Langevi... [Read more...]

Typos…

October 5, 2010 | xi'an

Edward Kao just sent another typo found both in  Monte Carlo Statistical Methods (Problem 3.21) and in Introducing Monte Carlo Methods with R (Exercise 3.17), namely that should be I also got another email from Jerry Sin mentioning that matrix summation in the matrix commands of Figure 1.2 of Introducing Monte Carlo Methods ...
[Read more...]

Le Monde puzzle [34]

October 3, 2010 | xi'an

Since the puzzle in this week (-end) edition of Le Monde is not (easily) solvable via an R program, I chose to go back to an older puzzle that my students can solve. Eleven [distinguishable] token are [arbitrarily] distributed around a 200 meter perimeter-long ring. They all start moving at the ... [Read more...]

Typo in Example 5.18

October 2, 2010 | xi'an

Edward Kao pointed out several typos in Example 5.18 of Monte Carlo Statistical Methods. First, the customers in area i should be double-indexed, i.e. which implies in turn that . Then the summary T should be defined as and as given that the first m customers have the fifth plan missing. ...
[Read more...]

Le Monde puzzle [38]

September 29, 2010 | xi'an

Since I have resumed my R class, I will restart my resolution of Le Monde mathematical puzzles…as they make good exercises for the class. The puzzle this week is not that exciting: Find the four non-zero different digits a,b,c,d such that abcd is equal to the ... [Read more...]

Galton & simulation

September 27, 2010 | xi'an

Stephen Stigler has written a paper in the Journal of the Royal Statistical Society Series A on Francis Galton’s analysis of (his cousin) Charles Darwin’ Origin of Species, leading to nothing less than Bayesian analysis and accept-reject algorithms! “On September 10th, 1885, Francis Galton ushered in a new era of ... [Read more...]

Riemann, Langevin & Hamilton [reply]

September 27, 2010 | xi'an

Here is a (prompt!) reply from Mark Girolami corresponding to the earlier post: In preparation for the Read Paper session next month at the RSS, our research group at CREST has collectively read the Girolami and Calderhead paper on Riemann manifold Langevin and Hamiltonian Monte Carlo methods and I hope ... [Read more...]

Effective sample size

September 23, 2010 | xi'an

In the previous days I have received several emails asking for clarification of the effective sample size derivation in “Introducing Monte Carlo Methods with R” (Section 4.4, pp. 98-100). Formula (4.3) gives the Monte Carlo estimate of the variance of a self-normalised importance sampling estimator (note the change from the original version ... [Read more...]

Monte Carlo Statistical Methods third edition

September 23, 2010 | xi'an

Last week, George Casella and I worked around the clock on starting the third edition of Monte Carlo Statistical Methods by detailing the changes to make and designing the new table of contents. The new edition will not see a revolution in the presentation of the material but rather a ...
[Read more...]

R tee-shirt

September 21, 2010 | xi'an

I gave my introduction to the R course in a crammed amphitheatre of about 200 students today. Had to wear my collectoR teeshirt from Revolution Analytics, even though it only made the kids pay attention for about 30 seconds… The other few “lines” that worked were using the Proctor & Gamble “car 54″ poster ...
[Read more...]

Typo in Example 3.6

September 16, 2010 | xi'an

Edward Kao pointed out the following difficulty about Example 3.6 in Chapter 3 of “Introducing Monte Carlo Methods with R”: I have two questions that have puzzled me for a while. I hope you can shed some lights. They are all about Example 3.6 of your book. 1. On page 74, there is a term [...] [Read more...]
1 37 38 39 40 41 45

Never miss an update!
Subscribe to R-bloggers to receive
e-mails with the latest R posts.
(You will not see this message again.)

Click here to close (This popup will not appear again)