(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 .

or and , below

and finally, below,

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

> 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...