R语言负二项分布函数实践

R语言负二项分布函数实践

在统计学中，负二项分布（Negative Binomial Distribution）是一种描述离散随机变量的概率分布。它常用于表示在一系列独立的伯努利试验中，成功次数固定而实验次数不确定的情况。本文将展示如何使用R语言中的负二项分布函数进行实践，并附带相应的源代码。让我们开始吧！

负二项分布的概率密度函数为：

P(X=k) = (k+r-1)C(r-1) * p^r * (1-p)^k

其中，X是成功次数，k是实验次数，r是成功的个数，p是单次实验成功的概率。

首先，我们需要安装并加载R语言中的extraDistr库，该库提供了负二项分布函数的相关功能。可以使用以下命令进行安装和加载：

install.packages("extraDistr")  # 安装extraDistr库
library(extraDistr)  # 加载extraDistr库

接下来，我们定义一个函数，用于生成负二项分布的随机样本。函数的输入参数包括实验次数n、成功的个数r和单次实验成功的概率p。函数使用rnbinom()函数从负二项分布中生成指定数量的随机样本。

generate_negative_binomial <- function(n, r, p) {
  samples <- rnbinom(n, r, p)
  return(sampl