等效于 Python 中的“qt”R 函数

Question

我在 R 中有以下代码：

n <- 112 # Observations 
p <- 4 # Variables
alpha <- 0.05 # Alpha is alpha 
quant = qt(1-alpha/2, n-p-1) # which is 1.982383

根据我的研究，qtR 中的函数与t 分布有关：R - qt 函数语法

我的问题是：

如何在 Python 中获得等价物？

Answer 1

这是scipy.stats.t.ppf：

In [26]: import scipy.stats                                                                              
In [27]: n = 112                                                                                         
In [28]: p = 4                                                                                           
In [29]: alpha = 0.05                                                                                    
In [40]: scipy.stats.t.ppf(1 - alpha / 2, n - p - 1)                                                              
Out[40]: 1.9823833701230174

或者，如果您不想手动指定尾部是您感兴趣的内容，则有scipy.stats.t.isf：

In [26]: import scipy.stats                                                                              
In [27]: n = 112                                                                                         
In [28]: p = 4                                                                                           
In [29]: alpha = 0.05                                                                                    
In [40]: scipy.stats.t.ppf(1 - alpha / 2, n - p - 1)                                                              
Out[40]: 1.9823833701230174