总和为 1 的概率权重的有效随机数组

Question

我需要一些非常简单的东西，在我看来我的解决方案有太多步骤（考虑到 R 是为统计而完成的）：

我只需要一个包含 N 个概率权重（介于 0 和 1 之间）的数组，总和为 1。什么是简单和/或有效的解决方案？

我的解决方案：

N <- 40
weights <- runif(40)
f <- 1/sum(weights) 
weights <- f*weights
print(sum(weights))

Answer 1

解决方案：

你可以使用prop.table函数和runif函数：

weights <- prop.table(runif(N))

和跑步一样

w <- runif(N); weights <- w / sum(w)

说明：

1. runif 函数生成 N 0 到 1 范围内的随机数
2. prop.table 重新调整它们，使 sum 等于 1（与 w / sum(w) 相同）