Post Made Community Wiki by whuber
source | link

Introductory examples in computational statistics class

I'm looking for an example of Bayesian inference for a class with the properties:

  1. The problem is easy to state, and the model & prior are both pretty reasonable, and
  2. R can't really calculate the integral that shows up when you try to get the posterior distribution (either it gives an error, or it gives an answer that is very far off).

Bonus points if it is known that the integral in question really doesn't have a `nice formula' in some sense (e.g. in the formalism introduced by the people who do differential Galois theory).

In case a reader is feeling especially kind, I'll say why I care about this:

I'll be teaching a computational statistics class next fall, and am looking for some tidy examples to give in the first 1-2 lectures that convey something about where problems in computational statistics come from and that they are somewhat unavoidable. I plan to have a moderate number of these (perhaps 5), and so I don't want any of them to require a lot of setup. Normally I don't worry so much about getting a 'hook' like this in the first class, but the first few weeks will otherwise be taken up with quite a bit of introductory material (intro to programming in R, then a review of notation, optimization, notions of complexity & basic computer science theory) and I worry that we won't get to any of the `good stuff' before the shopping period is over.

The class is at the advanced undergrad/early grad level, and I'll be teaching out of Computational Statistics by Givens and Hoeting.