R语言负二项(nbinom)分布

假设有一组独立的伯努利数列，每次实验有两种结果“成功”和“失败”。每次实验的成功概率是p，失败的概率是1－p。我们得到一组数列，当预定的“非成功”次数达到r次，那么结果为“成功”的随机次数会服从负二项分布：X~NB(r;P)

帕斯卡分布，负二项分布的正整数形式，描述第n次成功发生在第x次的概率，是统计学上一种离散概率分布，常用于描述生物群聚性，医学上用来描述传染性或非独立性疾病的分布和致病生物的分布。

####负二项(nbinom)分布
# 1.负二项(nbinom)分布中抽样函数rnbinom
# r: size; p: prob
size = 10
prob <- 0.1
rnbinom(100, size=size, prob=prob)

# 2.负二项分布概率密度函数
x <- seq(1,100)
y <- dnbinom(x,size=size, prob=prob)
plot(x,y)

# 3.负二项分布累积概率
# P[X ≤ x]
pnbinom(60,size=size, prob=prob)
# P[X > x]
pnbinom(60,size=size, prob=prob,lower.tail = FALSE)

# probabilities p are given as log(p).
pnbinom(60,size=size, prob=prob,log.p = TRUE)

# 4.qnbinom函数(pnbinom的反函数)
# 累积概率为0.95时的x值
# x <- seq(1,200)
# plot(x,pnbinom(x,size=size, prob=prob))
qnbinom(0.95, size=size, prob=prob)
qnbinom(0.995, size=size, prob=prob)