ディリクレ分布から 100 個の乱数を生成していますが、出力を使用してガンマ分布から生成する必要があります。コードは次のとおりです。
a <- rdirichlet(100, c(1,1,1))
b <- c(3,3,3)
sapply(a, function(x) {rgamma(100, shape=2, rate =(b%*%a)) })
ここで、ガンマ分布の割合はベクトル b と a の内積であることに注意してください (これはディリクレからの出力です)。
次のエラー メッセージが表示されます。
Error in b %*% a : non-conformable arguments