编辑代码

Ex <- 1.09
p <- 0.00109

r <- Ex*(1-p)/p
print(r) # r= 998.91

# P(X>=10) = 1-P(X=0)-P(X=1)-...-P(X=9)
p09 <- pnbinom(9,r,1-p) # r要与后面的概率值相对应，例如：r若表示失败次数，则prob应为失败概率

cat("calculate Ex by pnbinom:",1-p09,"\n")