Lấy mẫu từ nhị thức (K, p) với kết quả không mong muốn

Question

Trong tệp trợ giúp cho rbinom, đối số kích thước là một số thử nghiệm (bao gồm số không) nhưng không nói nếu điều này cũng có thể là vectơ.

Cách đúng của việc sử dụng chức năng này là

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

    0 1 
809 191 

Nhưng những gì đang xảy ra ở đây?

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

    0 1 
894 106 

Answer 1

Tái chế đối số kích thước là nguyên nhân chính.

Vì n là 1000, 0:1 được tái chế cho đến khi bạn nhận được 500 0 và 500 1 của (xen kẽ).

Tất cả những cái 0-kích thước cho 0:

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

Kết quả là 500 0 's + 500 thử nghiệm Bernoulli với p = 0,2, kết quả là khoảng 100 1' s ra 1000 yếu tố.

[Kết quả của bạn dường như không đáng ngạc nhiên đối với tôi, nhưng đối số tái chế có thể bị cắn nếu bạn không tìm kiếm, và - trong khi có lý do tại sao số lần thành công trong 0 thử nghiệm Bernoulli phải được định nghĩa là 0 - nó có vẻ không rõ ràng lúc đầu hoặc]

Answer 2

Tài liệu lỗi:.

Nếu 'kích thước' không phải là một số nguyên , 'NaN' được trả về. [nhấn mạnh của tôi]

Bạn đang cung cấp nhiều hơn một số nguyên, vì vậy tài liệu sẽ ngụ ý rằng bạn sẽ nhận được NaN.

Điều khó hiểu của nó bởi vì nó nêu rõ nơi các đối số khác có thể là vectơ nhưng không phải là size. Tôi sẽ nộp một lỗi tài liệu với người bảo trì, mà trong trường hợp này có thể có nghĩa là trình theo dõi lỗi R chính.