R语言 超几何分布

R语言中的超几何分布被定义为一种方法，用于计算无替换抽样时的概率，以获得密度值。

在R语言中，有4个内置函数来生成超几何分布。

• dhyper()

dhyper(x, m, n, k)

• phyper()

phyper(x, m, n, k)

• qhyper()

qhyper(x, m, n, k)

• rhyper()

rhyper(N, m, n, k)

其中

x： 代表数值的数据集

m： 群体的大小

n： 抽取的样本数

k： 群体中的项目数

N： 超几何分布的数值

生成超几何分布的函数

dhyper() 函数

它被定义为超几何密度分布，用于获得密度值。

语法

dhyper(x_dhyper, m, n, k)

例1 :

   
# Specify x-values for dhyper function
x_dhyper <- seq(0, 22, by = 1.2)
  
# Apply dhyper function
y_dhyper <- dhyper(x_dhyper, m = 45, n = 30, k = 20)   
   
# Plot dhyper values
plot(y_dhyper) 

输出：

phyper() 函数

超几何累积分布函数，用于估计在测试或调试过程开始时最初驻留在程序中的故障数量，基于超几何分布，使用相应的值计算x中的每个值。

语法

phyper(x, m, n, k)

例1 :

# Specify x-values for phyper function
x_phyper <- seq(0, 22, by = 1)    
    
# Apply phyper function
y_phyper <- phyper(x_phyper, m = 40, n = 20, k = 31)  
  
# Plot phyper values
plot(y_phyper)  

输出

qhyper() 函数

它基本上是 Hypergeometric Quantile Function，用于指定0到1之间的概率序列。

语法

qhyper(x, m, n, k)

例1 :

# Specify x-values for qhyper function
x_qhyper <- seq(0, 1, by = 0.02)        
  
# Apply qhyper function
y_qhyper <- qhyper(x_qhyper, m = 49, n = 18, k = 30)    
  
# Plot qhyper values
plot(y_qhyper)

输出：

rhyper()函数

它一般是指通过指定种子和样本量来生成随机数的函数。

语法

rhyper(x, m, n, k)

例1 :

# Set seed for reproducibility
# Specify sample size
set.seed(400)                                 
N <- 10000                                      
  
# Draw N hypergeometrically distributed values
y_rhyper <- rhyper(N, m = 50, n = 20, k = 30) 
y_rhyper         
  
# Plot of randomly drawn hyper density
hist(y_rhyper,                                          
     breaks = 50,
     main = "")

输出: