Posts Tagged ‘ books ’

ratio-of-uniforms [#4]

December 1, 2016
By
ratio-of-uniforms [#4]

Possibly the last post on random number generation by Kinderman and Monahan’s (1977) ratio-of-uniform method. After fiddling with the Gamma(a,1) distribution when a<1 for a while, I indeed figured out a way to produce a bounded set with this method: considering an arbitrary cdf Φ with corresponding pdf φ, the uniform distribution on the set […]

Read more »

sampling by exhaustion

November 24, 2016
By
sampling by exhaustion

The riddle set by The Riddler of last week sums up as follows: Within a population of size N, each individual in the population independently selects another individual. All individuals selected at least once are removed and the process iterates until one or zero individual is left. What is the probability that there is zero […]

Read more »

Monty Python generator

November 22, 2016
By
Monty Python generator

By some piece of luck I came across a paper by the late George Marsaglia, genial contributor to the field of simulation, and Wai Wan Tang, entitled The Monty Python method for generating random variables. As shown by the below illustration, the concept is to flip the piece H outside the rectangle back inside the […]

Read more »

simulation under zero measure constraints

November 16, 2016
By
simulation under zero measure constraints

A theme that comes up fairly regularly on X validated is the production of a sample with given moments, either for calibration motives or from a misunderstanding of the difference between a distribution mean and a sample average. Here are some entries on that topic: How to sample from a distribution so that mean of […]

Read more »

copy code at your own peril

November 13, 2016
By
copy code at your own peril

I have come several times upon cases of scientists [I mean, real, recognised, publishing, senior scientists!] from other fields blindly copying MCMC code from a paper or website, and expecting the program to operate on their own problem… One illustration is from last week, when I read a X Validated question [from 2013] about an […]

Read more »

Example 7.3: what a mess!

November 12, 2016
By
Example 7.3: what a mess!

A rather obscure question on Metropolis-Hastings algorithms on X Validated ended up being about our first illustration in Introducing Monte Carlo methods with R. And exposing some inconsistencies in the following example… Example 7.2 is based on a [toy] joint Beta x Binomial target, which leads to a basic Gibbs sampler. We thought this was […]

Read more »

variance of an exponential order statistics

November 9, 2016
By
variance of an exponential order statistics

This afternoon, one of my Monte Carlo students at ENSAE came to me with an exercise from Monte Carlo Statistical Methods that I did not remember having written. And I thus “charged” George Casella with authorship for that exercise! Exercise 3.3 starts with the usual question (a) about the (Binomial) precision of a tail probability […]

Read more »

SAS on Bayes

November 7, 2016
By
SAS on Bayes

Following a question on X Validated, I became aware of the following descriptions of the pros and cons of Bayesian analysis, as perceived by whoever (Tim Arnold?) wrote SAS/STAT(R) 9.2 User’s Guide, Second Edition. I replied more specifically on the point It [Bayesian inference] provides inferences that are conditional on the data and are exact, […]

Read more »

ratio-of-uniforms [#3]

November 3, 2016
By
ratio-of-uniforms [#3]

Being still puzzled (!) by the ratio-of-uniform approach, mostly failing to catch its relevance for either standard distributions in a era when computing a cosine or an exponential is negligible, or non-standard distributions for which computing bounds and boundaries is out-of-reach, I kept searching for solutions that would include unbounded densities and still produce compact […]

Read more »

ratio-of-uniforms [#2]

October 30, 2016
By
ratio-of-uniforms [#2]

Following my earlier post on Kinderman’s and Monahan’s (1977) ratio-of-uniform method, I must confess I remain quite puzzled by the approach. Or rather by its consequences. When looking at the set A of (u,v)’s in R⁺×X such that 0≤u²≤ƒ(v/u), as discussed in the previous post, it can be represented by its parameterised boundary u(x)=√ƒ(x),v(x)=x√ƒ(x)    x […]

Read more »


Subscribe

Email:

  Subscribe