如何计算Python中正态分布的累积分布函数(CDF)的反函数?
我应该使用哪个库?可能是卑鄙的?
NORMSINV(在注释中提到)是标准正态分布的CDF的倒数。使用scipy,您可以使用对象的ppf方法进行计算scipy.stats.norm。首字母缩写词ppf代表 百分比点函数 ,它是 分位数函数的 另一个名称。
scipy
ppf
scipy.stats.norm
In [20]: from scipy.stats import norm In [21]: norm.ppf(0.95) Out[21]: 1.6448536269514722
检查它是否与CDF相反:
In [34]: norm.cdf(norm.ppf(0.95)) Out[34]: 0.94999999999999996
默认情况下,norm.ppf使用mean = 0和stddev = 1,这是“标准”正态分布。您可以通过分别指定loc和scale参数来使用不同的均值和标准差。
norm.ppf
loc
scale
In [35]: norm.ppf(0.95, loc=10, scale=2) Out[35]: 13.289707253902945
如果您看一下的源代码scipy.stats.norm,就会发现该ppf方法最终会调用scipy.special.ndtri。因此,要计算标准正态分布的CDF的倒数,可以直接使用该函数:
scipy.special.ndtri
In [43]: from scipy.special import ndtri In [44]: ndtri(0.95) Out[44]: 1.6448536269514722