乱数を使う my.qdirichlet <- function(v,a,n=10000){ r <- rdirichlet(n,a) d <- ddirichlet(r,a) d0 <- ddirichlet(v,a) length(which(d<=d0))/n } ベータ分布で試してみる n <- 100 a <- c(1,1) a2 <- c(3,8) R <- rdirichlet(n,a) q <- rep(0,n) for(i i…