r - Binomial(K, p) からのサンプリングと予期しない結果

Question

のヘルプファイルではrbinom、size 引数は試行回数 (0 を含む) ですが、これがベクトルになるかどうかはわかりません。

この関数の正しい使い方は

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

score 9 · Accepted Answer

サイズ引数の引数リサイクルが主な原因です。

nは 1000 なので、500代と 500代になるまで (交互に)0:1リサイクルされます。01

サイズが 0 のものはすべて次のようになります0。

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

0p=0.2 で 500 回 + 500 回のベルヌーイ試行が行わ1れ、1000 個の要素のうち約 100 個が得られます。

[あなたの結果は私には驚くべきものではありませんでしたが、探していなければ引数のリサイクルが噛み付く可能性があります。また、0 回のベルヌーイ試行の成功回数を 0 と定義する必要があるのには理由がありますが、そうは見えないかもしれません。最初は明らかです。]

score 4 · Accepted Answer

ドキュメントのバグ:

'size' が整数でない場合、'NaN' が返されます。【私のこだわり】

あなたはそれに複数の整数を与えているので、ドキュメントはあなたが得ることを暗示しますNaN.

他の引数をベクトルにすることはできますが、できないことを明示的に示しているため、混乱しますsize。メンテナにドキュメンテーションのバグを報告します。この場合、おそらくメインの R バグトラッカーを意味します。

2 に答える 2