在 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
最佳答案
大小参数的参数回收是主要原因。
因为 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 - 它可能看起来并不一开始也很明显。]
关于r - 从具有意外结果的二项式(K,p)中采样,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/14920734/