r - 从具有意外结果的二项式(K，p)中采样

Question

在 rbinom 的帮助文件中, size 参数是多次试验(包括零)，但它没有说明这是否也可以是向量。

使用此功能的正确方法是

table(rbinom(n = 1000, size = 1, prob = 0.2))

  0   1 
809 191

但是这里发生了什么？

table(rbinom(n = 1000, size = 0:1, prob = 0.2))

  0   1 
894 106 

Answer 1

大小参数的参数回收是主要原因。

因为 n是 1000，0:1被回收直到你得到 500 0和 500 1的(交替)。

所有 0 尺寸的都给出 0 :

   > rbinom(10,size=0,prob=0.2)  
     [1] 0 0 0 0 0 0 0 0 0 0

结果是 500 0 's + 500 次伯努利试验，p=0.2，结果约为 100 1超过 1000 个元素。

[你的结果对我来说似乎并不令人惊讶，但如果你不寻找它，参数回收可能会咬人，而且 - 虽然有理由将 0 次伯努利试验的成功次数定义为 0 - 它可能看起来并不一开始也很明显。]