如何计算python中正常累积分布函数的倒数?

时间:2013-12-17 05:57:49

标签: python scipy normal-distribution

如何计算Python中正态分布的累积分布函数(CDF)的倒数?

我应该使用哪个库?可能是scipy?

3 个答案:

答案 0 :(得分:92)

NORMSINV(在评论中提到)是标准正态分布的CDF的倒数。使用scipy,您可以使用scipy.stats.norm对象的ppf方法进行计算。首字母缩写ppf代表percent point function,这是quantile function的另一个名称。

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,这是“标准”正态分布。您可以通过分别指定locscale参数来使用不同的均值和标准差。

In [35]: norm.ppf(0.95, loc=10, scale=2)
Out[35]: 13.289707253902945

如果查看scipy.stats.norm的源代码,您会发现ppf方法最终会调用scipy.special.ndtri。因此,要计算标准正态分布的CDF的倒数,可以直接使用该函数:

In [43]: from scipy.special import ndtri

In [44]: ndtri(0.95)
Out[44]: 1.6448536269514722

答案 1 :(得分:10)

# given random variable X (house price) with population muy = 60, sigma = 40
import scipy as sc
import scipy.stats as sct
sc.version.full_version # 0.15.1

#a. Find P(X<50)
sct.norm.cdf(x=50,loc=60,scale=40) # 0.4012936743170763

#b. Find P(X>=50)
sct.norm.sf(x=50,loc=60,scale=40) # 0.5987063256829237

#c. Find P(60<=X<=80)
sct.norm.cdf(x=80,loc=60,scale=40) - sct.norm.cdf(x=60,loc=60,scale=40)

#d. how much top most 5% expensive house cost at least? or find x where P(X>=x) = 0.05
sct.norm.isf(q=0.05,loc=60,scale=40)

#e. how much top most 5% cheapest house cost at least? or find x where P(X<=x) = 0.05
sct.norm.ppf(q=0.05,loc=60,scale=40)

答案 2 :(得分:1)

Python 3.8开始,标准库提供NormalDist对象作为statistics模块的一部分。

它可以用来获取 逆累积分布函数 inv_cdf -cdf的逆),也称为给定平均值mu)和分位数函数百分点函数标准偏差sigma):

from statistics import NormalDist

NormalDist(mu=10, sigma=2).inv_cdf(0.95)
# 13.289707253902943

可以将其简化为标准正态分布mu = 0sigma = 1):

NormalDist().inv_cdf(0.95)
# 1.6448536269514715