以给定的概率生成随机数

Question

我有一个问题，基本上是以下 matlab 问题的矢量化 R 解决方案：

Generate random number with given probability matlab

我能够根据随机统一数和每个单个事件的给定概率（总和为 100% - 仅可能发生一个事件）生成随机事件结果：

sum(runif(1,0,1) >= cumsum(wdOff))

但是该函数只接受一个随机统一数字，而我希望它接受一个随机统一数字的向量并为这些条目输出相应的事件。

所以基本上我正在寻找 Matlab 中 Oleg 向量化解决方案的 R 解决方案（从 cmets 到 matlab 解决方案）：

“向量化解：sum(bsxfun(@ge, r, cumsum([0, prob]),2) 其中r 是列向量，prob 是行向量。--Oleg”

如果您需要更多信息，请告诉我。

Answer 1

你可以只做一个加权随机样本，而不用担心你的 cumsum 方法：

sample(c(1, 2, 3), size = 100, replace = TRUE, prob = c(0.5, 0.1, 0.4))

如果你已经有了数字，你也可以这样做：

x <- runif(10, 0, 1)
as.numeric(cut(x, breaks = c(0, 0.5, 0.6, 1)))