我想计算多元dirichlet分布的密度，并根据这种分布生成随机实现。就像函数dmvnorm对多元正态分布所做的那样。我发现这是正态分布，我想知道是否有一个函数可以为Dirichlet和Gamma分布做到这一点：
g <- expand.grid(x = seq(-2,2,0.05), y = seq(-2,2,0.05)) ## x and y are the 2 normal distributions.
g$z <- dmvnorm(x=cbind(g$x,g$y),mean = c(0,0),sigma = diag(2),log = FALSE)