我从 Dirichlet 分布生成 100 个随机数,然后我需要使用输出从 Gamma 分布生成。这是代码:
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 的点积(这是 Dirichlet 的输出)。
我收到此错误消息:
Error in b %*% a : non-conformable arguments