sample函数的一个bug

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

想请教，当查看sample函数的例子时，发现有一个bug，即生成的一个长度为10的样本，这是为什么呀？谢谢~

# sample()'s surprise：

x <- 1:10   
sample(x[x >  8]) # length 2   
sample(x[x >  9]) # oops -- length 10!   
sample(x[x > 10]) # length 0


## safer version:
resample <- function(x, ...) x[sample.int(length(x), ...)]
resample(x[x >  8]) # length 2
resample(x[x >  9]) # length 1
resample(x[x > 10]) # length 0

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：Sample AMPL amp bug SAM surprise version 样本

实际生成sample时是下面这样的

> x <- 1:10   
> sample(x[x >  8])
[1] 10  9
> sample(x[x >  9])
[1]  6  5  4  2  1  3  9  8 10  7

1. sample(x[x>9])

复制代码

实际上就是:

1. sample(10)

复制代码

生成10个随机数，没有问题的

jiangbeilu 发表于 2016-5-26 17:45
实际上就是:

生成10个随机数，没有问题的

谢谢！刚刚试了一下，确实。sample(c(10))是生成10个随机数，这种只有一个元素的情况应该就变成是sample(10)了。