(This article was first published on Freakonometrics - Tag - R-english, and kindly contributed to R-bloggers)
In the course, still introducing some concept of dependent distributions, we will talk about the Dirichlet distribution (which is a distribution over the simplex of
). Let
denote the Gamma distribution with density (on
)

Let
denote independent
random variables, with
. Then
where

has a Dirichlet distribution with parameter
. Note that
has a distribution in the simplex of
,

and has density

We will write
.
|
|
|
![]() |
![]() |
|
|
|
![]() |
![]() |
can be visualized below, e.g.
, with some kind of symmetry, 
or
and
, below
and finally, below,

Note that marginal distributions are also Dirichlet, in the sense that if


, and if
, then
's have Beta distributions, 
> library(MCMCpack) > alpha=c(2,2,5) > x=seq(0,1,by=.05) > vx=rep(x,length(x)) > vy=rep(x,each=length(x)) > vz=1-x-vy > V=cbind(vx,vy,vz) > D=ddirichlet(V, alpha) > persp(x,x,matrix(D,length(x),length(x))
To leave a comment for the author, please follow the link and comment on his blog: Freakonometrics - Tag - R-english.
R-bloggers.com offers daily e-mail updates about R news and tutorials on topics such as: visualization (ggplot2, Boxplots, maps, animation), programming (RStudio, Sweave, LaTeX, SQL, Eclipse, git, hadoop, Web Scraping) statistics (regression, PCA, time series, trading) and more...









Zero Inflated Models and Generalized Linear Mixed Models with R.
Zuur, Saveliev, Ieno (2012).