rbinom(1000, size=74, prob=0.89)
我对rbinom函数仍然有点生疏。在上面的代码中,我生成了对74个试验的1000个观察结果,每个试验的成功概率= 0.89。我希望我生成的偏差小于或等于74.如何限制rbinom函数来解决这个问题?虽然我很少,但是从上面的代码中偏离了> = 74,这导致我的程序出现问题。
答案 0 :(得分:2)
您偶尔会有一些等于74的偏差,因此>=
可能会返回TRUE
。他们从不大于74但是因为你只运行了74次试验。请改为any( x > 74 )
。
set.seed(5895)
x <- rbinom(1000, size=74, prob=0.89)
table( x >= 74 )
# FALSE TRUE
# 999 1
table( x > 74 )
# FALSE
# 1000
max(x)
# [1] 74
是的,我确实需要经过一些.Random.seed
才能找到一个可以得到这个结果! : - )