相當於Matlab的Python中

Question

gaminv的，我想在Python翻譯下面Matlab代碼：

tau=40 %scale parameter 
a=3 %shape parameter 
t = gaminv(0.90,a,tau); 

代碼返回T = 212.8928。我試過了：

import scipy.stats.distributions as dist 
tau=40 
a=3 
t = dist.invgamma.cdf(90, a, scale = tau) 

但python代碼沒有返回相同的t值。難道我做錯了什麼？

Answer 1

您希望ppf（「百分比點功能」，這是用於CDF或位數函數的逆一個稍差通用名）的scipy.stats.gamma方法：

In [27]: from scipy.stats import gamma 

In [28]: gamma.ppf(0.9, 3, scale=40) 
Out[28]: 212.8928135133684 

參見How to calculate the inverse of the normal cumulative distribution function in python?;在那裏的分佈是不同的，但它基本上是同一個問題。

invgamma是不同的分佈。