需要在R中编写一个新的随机生成器

Question

我需要在-1,1之间生成7个随机数,它们的总和等于1.我用这个代码来做.

 diff(c(0, sort(round(runif(7,-1,1),2)), 1))

但是我有一个很大的问题.

该代码的一个输出是-0.89,0.21,0.00,0.21,0.30,0.19,0.61,-0.63.

问题是它是统一的我猜这样每次都会在第一个和最后一个数字中生成大的随机数,这不是我想要的.我需要将它传播到所有数字.恩.0.22 -.21 .33 -.12 0.11 0.35 -0.08(总和不等于1只是一个例子)

你知道我可以写一个代码来获得这种随机数吗？

Answer 1

您的总体想法可能受到随机描述中链接的答案的启发.标准问题是如何在0和1之间生成7个加1的数字.答案是:

diff(c(0, sort(runif(6, 0, 1)), 1))
#> [1] 0.27960792 0.02035231 0.02638626 0.09945877 0.25134002 0.03379598 0.28905874

获得-1到1之间数字的必要修改非常简单; 只是忽略了sort:

diff(c(0, runif(6, 0, 1), 1))
#> [1]  0.9961661 -0.6528227  0.5298829 -0.2087127 -0.2298045  0.2017705  0.3635203

这是如何运作的？我们再次划分零和一之间的空间.但是b省略了这种排序,我们考虑到了倒退的可能性,即负数是可能的.这是1000代的直方图:

这种方法的一个弱点是第一个和最后一个数字必然是正数.如果这困扰你,你可以添加一个额外的sample,例如:

sample(diff(c(0, runif(6, 0, 1), 1)), 7)
#> [1] -0.004242793 -0.725348335  0.385971491  0.320525822  0.389915347
#> [6]  0.053195271  0.579983197