生成的值，其总和为1的向量但是元素应该是大约但不完全均匀的

Question

假设我想使长度为N的向量我的目标是确保所有的元素是大约1 / N（但不是所有的元素应该是正好是1 / N）和矢量的总和应等于一。有没有一种方法R中做到这一点？

Answer 1

您可以使用狄利克雷分布，大和等参数。

> library(bayesm)
> rdirichlet(c(100,100,100))
          [,1]
[1,] 0.3360327
[2,] 0.3280179
[3,] 0.3359495
> rdirichlet(c(10000,10000,10000))
          [,1]
[1,] 0.3311760
[2,] 0.3353209
[3,] 0.3335032
> rdirichlet(c(1000000,1000000,1000000))
          [,1]
[1,] 0.3330997
[2,] 0.3334884
[3,] 0.3334119