我想计算多元狄利克雷分布的密度，并从这种分布中生成随机实现。就像函数 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)