cut&paste typo in R book

A casualty of cut-and-paste in Chapter 3 of Introducing Monte Carlo Methods with R. Brad McNeney from Simon Fraser sent me a nice email about the end of Example 3.6 missing a marginal estimate. Indeed, it does. And it should have been obvious from the “estimates” we derived, 19 and 16, which are not even on the support of the posterior distribution represented on Figure 3.5… The R code is given as

> mean(y[,1]*apply(y,1,f)/den)/mean(apply(y,1,h)/den)
[1] 19.33745
> mean(y[,2]*apply(y,1,f)/den)/mean(apply(y,1,h)/den)
[1] 16.54468

and should have been

> mean(y[,1]*f(y)/den)/mean(f(y)/den)
[1] 94.08314
> mean(y[,2]*f(y)/den)/mean(f(y)/den)
[1] 80.42832

As also suggested by Brad, a similar modification applies to the remark after eqn. (3.7):


should be


One Response to “cut&paste typo in R book”

Leave a Reply

Fill in your details below or click an icon to log in: Logo

You are commenting using your account. Log Out /  Change )

Google photo

You are commenting using your Google account. Log Out /  Change )

Twitter picture

You are commenting using your Twitter account. Log Out /  Change )

Facebook photo

You are commenting using your Facebook account. Log Out /  Change )

Connecting to %s

This site uses Akismet to reduce spam. Learn how your comment data is processed.

%d bloggers like this: