# Simulate Bivariate Beta Distribution , BIBETA(6, 20, 2) in R [duplicate]

I need to simulate bivariate beta distribution, $$BIBETA(6, 20, 2)$$ in r. I am looking for a package/ code that would generate bivariate beta distribution. I couldn't find the r function for this distribution.

The probability density function of the Bivariate beta distribution is given as below.

$$f(x,y) = \dfrac{1}{B(a,b,c)} \dfrac{x^{a-1}y^{b-1}(1-x)^{b+c-1}(1-y)^{a+c-1}}{(1-xy)^{a+b+c}}$$ Is anyone know which r function can be used for this simulation?

There is a package called MultiRNG that implements this sort of multivariate simulation for a wide class of multivariate distributions (in your particular case, you are interested in the draw.dirichlet function).
• Do they have dirichlet cdf in r? Commented Oct 4, 2018 at 13:17
• I am trying to use the above function. I have a <- c(6, 20,2). I generated x and y random variables using draw.dirichlet. Now what is t in the above function. is t 2 by 2 matrix? Commented Oct 4, 2018 at 13:36
• The draw.dirichlet implements a Dirichlet distribution, which is NOT the bivariate beta asked for. Commented Jul 26, 2021 at 1:07