Re: [R] generate random numbers that sum up to 1

From: Alberto Monteiro <>
Date: Mon 16 Oct 2006 - 15:22:04 GMT

Grant Izmirlian wrote:
> So, Alberto, you didn't see my post?
I think I didn't - but you are demanding too much from my memory; I can hardly remember what I saw yesterday!

> If Y has d independent
> components that are gamma distributed with common rate and shapes
> A_1, A_2, ..., A_d, then X, given by the components of Y divided by
> their sum has distribution Dirichlet(A_1, A_2, ..., A_d). If you
> want Uniform on the d-simplex, then use A_1 = A_2 = ... = A_d = 1
> (just as Duncan said)
The problem is that I wasn't aware that this was the Dirichlet distribution [R does not have this distribution, AFAIK, but I should have consulted the Borg of All Wisdom, the Wikipedia].

Alberto Monteiro mailing list PLEASE do read the posting guide and provide commented, minimal, self-contained, reproducible code. Received on Tue Oct 17 01:54:56 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Mon 16 Oct 2006 - 20:30:28 GMT.

Mailing list information is available at Please read the posting guide before posting to the list.