Archive for model choice

The Search for Certainty

Posted in Books, Statistics with tags , , , , , , , , on January 25, 2010 by xi'an

“The expected value is hardly ever expected.” The Search for Certainty, page 207.

Following the advice of an Og’s reader, Erkan, I read Krzysztof Burdzy’s The Search for Certainty over the past weeks (and mostly in the métro). The author is primarily a probabilist, but he has set upon launching a radical criticism of the philosophical foundations of Statistics, the book subtitle being (somehow misleadingly) put as On the Clash of Science and Philosophy of Probability. Or so he thinks. (The introduction to the book is available on the author’s webpage and he has also started a blog related with the book. The cover is made of seven dices all showing sixes, with on top the probability p=0.0000036… of such an event, provided the dice is fair. And provided one does not take into account the suspicious fives all facing the reader!)

“Of the four well crystallized philosophies of probability, two chose the certainty as their intellectual holy grail. Those are the failed theories of von Mises and de Finetti.” The Search for Certainty, page 30.

Although I have read the book with a pen, I will not go at this stage into a detailed analysis of The Search for Certainty. Indeed, I found the book both annoying and unconvincing, for reasons not very different from the criticisms addressed at Taleb’s The Black Swan. The book aims at demonstrating that the philosophical arguments underlying both frequentist and Bayesian Statistics are wrong. Unsurprisingly (!), I find the book lacking in this demonstration and overall poor from a scholarly perspective. It compares with Taleb’s The Black Swan in that the attempts at philosophy are more related to everyday “common sense” than to deep (and scholarly) philosophy (and they also involve the apparently inevitable Karl Popper, “the champion of the propensity theory of probability”, p.43!). The main point made in The Search for Certainty is very narrow in that Burdzy concentrates on two very specific entries to frequentism and subjectivism, namely von Mises’ and de Finetti’s, respectively, while those are not your average statistician’s references. For instance, von Mises bases his definition of frequency properties on the notion of collectives, a notion I had not previously encountered. Similarly, de Finetti’s statement “Probability does not exist” cannot be seen as the core principle of many Bayesian statisticians and I certainly do not relate to his all-subjective perspective for conducting Bayesian inference.

“Mr. Winston is unique because we know something about him that we do not know about any other individual in the population.” The Search for Certainty, page 66.

The style of The Search for Certainty is fairly annoying, if of a different sort than Taleb‘s. It is predominantly non-technical, the worked-out examples being always of the coin and balls-in-an-urn type. Arguments are never more than one paragraph long and metaphors and weak analogies are more confusing than helpful. The book criticises a lot decision-theory and the related coherence of Bayesian procedures, ie against the Dutch Book argument, but it misses the connections between frequentist optimality(ies) and Bayesian procedures, like Wald’s complete class theorem, Welch and Peers (1964) matching priors, or the more recent Berger’s frequentist-Bayesian perspective. There are also (in my opinion) confusions, as when the basis for the frequency approach to probability is criticised for being connected with an unrealistic infinity of events, thus confusing concepts with experiments (indeed, the Large Law of Large Numbers cannot be proved by an experiment), the existence of a model versus its assessment (“Kolmogorov’s axioms say nothing on how to match the mathematical results with reality”, p.31), the use of a probability against its “truth”, the inclusion of time (and thus model shifts) into mathematical axioms, the assimilation of frequentist statistics to unbiased estimators, a somehow diffuse belief that some priors can be proven to be better than others, an argument that  they can be evaluated by their predictive performances… The insistence in adding new axioms to Kolmogorov’s is furthermore puzzling:

(L4) If there exists a symmetry on the space of possible outcomes which maps an event A onto an event B then the two events have equal probabilities, that is, P(A) = P(B).
(L5) An event has probability 0 if and only if it cannot occur. An event has probability 1 if and only if it must occur.

The axiom (L4) relates to both the Principle of Insufficient Reasons, whose limitations are the cause for much debate in the selection of prior distributions, and to invariance principles that lead to Haar measures as default noninformative priors. But I do not see the point in adding such an axiom into the tenets of probability. (And (L5) is more psychological than mathematical…)

“There is no justification for the use of the Bayes theorem in the subjective theory.” The Search for Certainty, page 144.

In conclusion, it is hopefully obvious that I did not overly enjoy this The Search for Certainty and that I do not consider it makes a significant contribution to the foundations of statistical inference and in particular to Bayesian analysis. Being examined by an outsider to our discipline certainly has a strong appeal, but only if done at a deep enough level.

Ps-Be warned that there is a homonymous book with full titlde The Search for Certainty: A Philosophical Account of Foundations of Mathematics on issues related with Gödel’s incompleteness theorem, by Marcus Gianquinto, book that I mistakenly bought for The Search for Certainty discussed here. Professor Gianquinto is actually a professor of philosophy at UCL, who specialises in epistemology so he could also comment on this book.

Model uncertainty in Warwick

Posted in Statistics, University life with tags , , , , , on November 19, 2009 by xi'an

Yet another workshop: Prior to the Valencia 9/ISBA 2010 meeting, there will be a workshop on model uncertainty organised by the Centre for Research in Statistical Methodology (CRiSM) at the University of Warwick, on May 30-June 1. The purpose of the meeting is to discuss recent developments and research on topics related to model uncertainty and model choice. The workshop largely follows the format of the 2006 CRiSM workshop on Bayesian Inference for Stochastic Systems. As in 2006, the timing of the conference has been chosen to ease the transition, for all interested participants, to the Valencia International Meetings on Bayesian Statistics, starting on June 3. The invited speakers are
Jim Berger (Duke)
Carlos Carvalho (Chicago)
David Dunson (Duke)
Jon Forster (Southampton)
Arnoldo Frigessi (Oslo)
Alan Gelfand (Duke)
Ed George (Pennsylvania)
Chris Holmes (Oxford)
Michael Jordan (Berkeley)
Robert Kohn (Sydney)
Athanasios Kottas (California)
Antonio Lijoi (Pavia)
David Madigan (Columbia)
Peter Müller (Texas)
Christian Robert (Paris)
David Spiegelhalter (Cambridge)
Yee Whye Teh (UCL)
Nanny Wermuth (Gothenburg)
Henry Wynn (LSE)
and I am looking forward this topical meeting because this is an area central to my research interests. Since most attendees will also be in Benidorm the following week, we can envision ten days of exciting debates about the nature of Bayesian model choice!

Yes III

Posted in Mountains, Statistics, University life with tags , , on October 4, 2009 by xi'an

No, I haven’t mixed my exclamation marks and my I’s on the keyboard: I am attending early this week the YES III meeting, where YES stands for Young European statisticians. The workshop is centred on Paradigms of Model Choice and is organised by Eurandom in Eindhoven, Holland. It is directed at young statisticians (mainly Ph.D. students and postdocs) who are interested in the problem of model choice and there are four short seminars of three 45 minute talks each on various aspects of model choice given by Laurie Davies, Duisburg-Essen, on approximate models and regularization, Peter Grünwald, Amsterdam, on the minimum description length principle,  Nils Hjort, Oslo on focused information criterion, and myself. I will cover my recent coverage of computational tool for Bayesian model choice and have the opportunity to expose the new results on the Dickey-Savage ratio for the first time! I am looking forward this meeting, having a fond memory of a previous Young Statisticians meeting in Ovronnaz, Switzerland a few years ago (including an epic hike with Jean-Michel Marin to the Muveran, where he took this great sunset picture of the Petit Muveran).

Ps- Another Yes!!! moment is that PNAS Editorial Board has reconsidered the decision to decline our letter, thus that the submission of our comments on Ratman et al. (June 30, 2009) has been allowed, with a further reduction of its length to 500 words and the standard provisions “Your Letter will be forwarded to the Editorial Board for consideration. The Board may decide to accept, reject, or ask for revisions to your Letter. Additionally, the Board may solicit a response from the authors. Accepted Letters will be published online only.”

JSM 2009 impressions [day 2]

Posted in Books, Statistics, University life with tags , , , , , , , , , , on August 4, 2009 by xi'an

Julien Cornebise wrote his impressions on yesterday [day 2] as comments to day 1 and he is welcome as a guest editor! I completely agree with his views on George Casella’s Medallion Lecture on design, which emphasized the need to reconsider this somehow neglected part of the Statistics curriculum. George’s lecture was both passionate and broad, which made it accessible to the large audience there. It was based on his Statistical Design book, on sale at the Springer Verlag booth in the Exhibit hall when you go to check the Enigma machine at the NSA booth. Along with a whole table of new books in the Use R! series, soon to be augmented by our book Introducing Monte Carlo Methods with R with George Casella, which is available in a draft version at the booth. (We actually signed the contract for Introducing Monte Carlo Methods with R with Springer yesterday afternoon.) The Springer editor, John Kimmel, is one of the ASA Fellows this year, in recognition of his support of the dissemination of new ideas in Statistics (my wording) and this is a great initiative from the ASA committee on Fellows as he unreservedly deserves it, if only for launching the Use R! series!

PICT6678

As mentioned by Julien, the session on the future of Statistics was reserved to the happy “fews” who managed to get a seat and others had to stay in the “present” thanks to this safety regulation that seems to be implemented on some talks/rooms and not others. I passed the first people being stopped by a fierce guard on my way to the “past”, ie to the cosmology and astrophysics session. There, I enjoyed very much Larry Wasserman’s talk on Nonparametric estimation of filaments for uncovering a challenging problem as well as for his elegant resolution of the problem. As well as the presentation by Laura Cayon of Detection of weak lensing, where I discovered that my old Purdue friend Anirban das Gupta was also involved in cosmology. I also went to the Monte Carlo and Sequential Analyses: Methods and Applications session, organised by Mike West, but the talks were too short to make much of an impact on me, even though I appreciated the talk by Minghui Shi on Particle stochastic search for high-dimensional variable selection that linked with Nicolas Chopin’s early work on exploring a large dataset and I was also intrigued by the talk of Ioanna Manolopoulou on Targeted sequential resampling from large Data sets in mixture modeling for using proxies to the real mixture model. The day ended up with a Board meeting for ISBA, that unfortunately took place outside in a hot humid weather… I now have to get ready for the Gertrude Cox Scholarship 5k race, since it starts at 6:15am (yes, am!).

JSM 2009 impressions [day 1]

Posted in Statistics, University life with tags , , , on August 3, 2009 by xi'an

Last afternoon, I attended the Bayesian model choice invited session, with talks by David Madigan, Hani Doss and Hugh Chipman. The room was packed and the talks were quite interesting. Madigan et al’s Sequential Bayesian Model Selection uses Mike West’s dynamic parameters to allow for more adaptive models, with the drawback of a lack of stabilisation in the estimation of the parameter which reflects the time-varying modelling. (I could not really understand why the model choice approach could not handle the hidden Markov chain of the model labels, though.)

The second talk, Doss’ Estimation of Large Families of Bayes Factors from Markov Chain Output, was more central to my interests (and to my own talk). and had several innovations worth mentioning. One was the use of the Bayes factor for prior comparison, which is something I never thought of before. In a strictly Bayesian perspective, this is not surprising as it means using the data to compare your priors! Completely un-orthodox! At a deeper level, I am still wondering at the validation of the Bayes factor in this setting… The second innovation was that, thanks to the same likelihood appearing in both numerator and denominator, Hani Doss was able to use a new type of bridge sampling estimator due to the identity

\mathbb{E}_{\pi_2(\cdot|x)} \left[ \dfrac{\pi_1(\theta)}{\pi_2(\theta)} \mid x \right] = \dfrac{m_1(x)}{m_2(x)}

and thus, using a posterior sample from one posterior is sufficient for approximating the Bayes factor in this case. A third innovation was using the control variates

\dfrac{\pi_i(\theta) - \pi_j(\theta)}{\pi_1(\theta) f(x|\theta)}

since they always are unbiased estimators of zero. A point I need to investigate further is how using improper priors (the example in Hani Doss’ talk was based on Zellner’s g-prior) is possible in this regard.

The talk by Hugh Chipman was about the BART (Bayesian additive regression tree) “machine learning” technique of Chpiman, George and McCulloch, which is always impressive as a completely non-parametric regression method. The innovation was in using many trees (50,200) simultaneously in the model, not as in model averaging, but as a basis for more complex dependences. The differences in the influence of the covariates seems to vane away as the number of trees goes up, but, as discussed during this session, this may be due to the frequency of uses of a covariate being too crude an indicator.

I then went to Aad van der Vart’s Le Cam lecture, on Some Frequentist Results on Posterior Distributions on Infinite-Dimensional Parameter Spaces, who, as usual, managed to give a very methodical and clear overview of the results on the asymptotics of Bayesian non-parametric estimators.

The meeting is, as forecasted, a monster of a meeting, but the conference center is so huge that it is not overwhelming. The JSM staff is quite efficient and they managed to solve the issue of the Series B editors meeting almost immediately.