我试图使用R从超几何分布中生成一些随机数。然而,当我有非常少量的白球和大量的黑球时,rhyper()表现得非常奇怪。这是我在电脑里得到的:
> sum(rhyper(100,1000,1e9-1000,1e6))
[1] 91
> sum(rhyper(100,2000,1e9-2000,1e6))
[1] 204
> sum(rhyper(100,10000,1e9-10000,1e6))
[1] 1016
> sum(rhyper(100,20000,1e9-20000,1e6))
[1] 1909
> sum(rhyper(100,50000,1e9-50000,1e6))
[1] 4968
> sum(rhyper(100,5000,1e9-5000,1e6))
[1] 60
> sum(rhyper(100,6000,1e9-6000,1e6))
[1] 164
> sum(rhyper(100,8000,1e9-8000,1e6))
[1] 0
> sum(rhyper(100,9000,1e9-9000,1e6))
[1] 45
前5个工作正常,但是对于第6个,我预计得到的数字大约为500,但不是60,而是7,7,8,9。
rhyper()函数或我的电脑有问题吗?