小编Pau*_*ul 的帖子

我目前正在使用apTreeshape使用"Yule-Hardy"方法模拟系统发育树.我想要做的是为三个不同的分组(小型,中型和大型树)随机生成20到25个不同的数字,然后为从分组中选择的每个随机数生成大约40棵树.

我知道如何在Python的Python中做到这一点,但在R中,事情看起来有点不同.

我的想法是,如果我要生成一个充满随机数的向量(每个大小分组一个),然后使用它生成一个基本上包含每个随机数的所有重复值的向量.

这是我有的:

sm_leaves<-c(sample(3:50,25,replace=F));
s_leafy<-numeric();

for (i in 1:length(sm_leaves)) { 
    for (j in 1:10) {
        s_leafy[j+i-1]=sm_leaves[i];
    }
}

这给我的输出如下:

> s_leafy
[1]  5 38  6 22 29 20 19 46  9 18 39 50 34 11 43  7  8 32 10 42 14 37
[23] 23 13 28 28 28 28 28 28 28 28 28 28

但我想要的更像是:

> s_leafy
[1]  5  5  5  5  5  5  5  5  5  5 38 38 38 38 38 38 38 38 …