从n组数字生成x个随机数

问题描述 投票:2回答:1

假设我在R中有以下数字列表,如何通过指定每组需要的数量来选择6个随机数?

r<-c(1,2,3,4,5,6,7,8,9,10)
b<-c(11,12,13,14,15,16,17,18,19,20)
y<-c(21,22,23,24,25,26,27,28,29,30)
g<-c(31,32,33,34,35,36,37,38,39,40)

我怎样才能选择r中的2个随机数,b中的2个随机数,y中的1个随机数和g中的1个随机数?

我需要我的代码足够灵活,可以手动输入每个列表需要多少个数字。根据我的选择,我还希望最终输出只是6个随机数的列表。

r random sample
1个回答
3
投票

将向量放入一个列表中

foo <- list(
    c(1,2,3,4,5,6,7,8,9,10),
    c(11,12,13,14,15,16,17,18,19,20),
    c(21,22,23,24,25,26,27,28,29,30),
    c(31,32,33,34,35,36,37,38,39,40)
)

Sample x times

如果您不介意替换,则生成另一个绘制向量(指定从哪个列表条目进行采样)并迭代此向量。

bar <- c(1, 1, 2, 2, 3, 4)
sapply(bar, function(x) sample(foo[[x]], 1))

Sample n times

如果您不想更换,则必须从每个列表条目中仅采样一次。要执行此操作,您必须为列表条目命名您需要的绘制数量并迭代列表(对于每个条目提取名称并对此数字进行采样)。

names(foo) <- c(2, 2, 1, 1)
unlist(sapply(seq_along(foo), function(x) sample(foo[[x]], names(foo[x]))))

Sample 1 time

另一个解决方案是生成概率向量并将其作为prob中的sample参数传递。

sample(unlist(foo), length(bar),
       prob = rep(table(bar) / length(bar) / lengths(foo), lengths(foo)))

说明:

我们需要每个条目的样本数量:table(bar)

1 2 3 4 
2 2 1 1 

为每个条目分配概率:table(bar) / length(bar)

        1         2         3         4 
0.3333333 0.3333333 0.1666667 0.1666667 

给每个数字它的概率被采样:rep(table(bar) / length(bar) / lengths(foo), lengths(foo))

...

         2          2          2          2          2          2          2          2          2          2 
0.03333333 0.03333333 0.03333333 0.03333333 0.03333333 0.03333333 0.03333333 0.03333333 0.03333333 0.03333333 
         3          3          3          3          3          3          3          3          3          3 
0.01666667 0.01666667 0.01666667 0.01666667 0.01666667 0.01666667 0.01666667 0.01666667 0.01666667 0.01666667 

...
© www.soinside.com 2019 - 2024. All rights reserved.