等价于 "qt" R 函数 Python

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

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

我的问题是：

How can I get the equivalent in 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 [40]: scipy.stats.t.isf(alpha / 2, n - p - 1)                                                              
Out[40]: 1.9823833701230174