等效於 Python 中的“qt” R function

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

From my research, the qt function from R is related with to the t -distribution: R - qt function syntax

我的問題是：

如何在 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