Sampling von einem Binomial (K, p) mit unerwartetem Ergebnis

Question

In den Hilfedateien für rbinom, Größe Argument ist eine Reihe von Studien (inkl. Eine Null), aber es sagt nicht, ob dies auch ein Vektor sein kann.

Der richtige Weg, um diese Funktion zu verwenden, ist

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

    0 1 
809 191 

Aber was hier geschieht?

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

    0 1 
894 106 

Answer 1

Argument Wiederverwertung der Größe Argument ist die Hauptursache.

Da n 1000, 0:1 zurückgeführt wird, bis Sie 500 0 's und 500 1' get s (abwechselnd).

die alle 0-size diejenigen geben 0:

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

Resultierende in 500 0 's + 500 Bernoulli-Versuchen mit p = 0,2, in etwa 100 1 resultierenden' von 1000 s-Elemente.

[Ihre Ergebnisse schienen mir nicht überraschend, aber Argument Recycling kann beißen, wenn Sie nicht danach suchen, und - obwohl es Gründe gibt, warum die Anzahl der Erfolge in 0 Bernoulli-Studien als 0 definiert werden sollte - es kann auch nicht auf den ersten offensichtlich scheinen]

Answer 2

Dokumentation Fehler.

Wenn ‚Größe‘ ist nicht eine ganze Zahl, ‚NaN‘ zurückgegeben. [mein Schwerpunkt]

Sie geben es mehr als eine ganze Zahl, so würde die Dokumentation bedeuten, dass Sie NaN erhalten würden.

Es ist verwirrend, weil es explizit angibt, wo andere Argumente Vektoren sein können, aber nicht size. Ich würde einen Dokumentationsfehler mit dem Betreuer einreichen, der in diesem Fall wahrscheinlich den Haupt-R-Bug-Tracker bedeutet.