CODEWITHSUNDEEP
pythonnumpyscipyquantile
dsaxton
dsaxton

Reputation: 1015

Quantile functions in Python

I'm having trouble finding quantile functions for well-known probability distributions in Python, do they exist? In particular, is there an inverse normal distribution function? I couldn't find anything in either Numpy or Scipy.

Upvotes: 9

Views: 13603

Answers (2)

MaxD
MaxD

Reputation: 377

It seems new but I've found this about numpy and quantile. Maybe you can have a look (not tested)

Upvotes: 0

dbouz
dbouz

Reputation: 919

Check the .ppf() method of any distribution class in scipy.stats. This is the equivalent of a quantile function (otherwise named as percent point function or inverse CDF)

An example with the exponential distribution from scipy.stats:

# analysis libs
import scipy
import numpy as np
# plotting libs
import matplotlib as mpl
import matplotlib.pyplot as plt

# Example with the exponential distribution
c = 0
lamb = 2

# Create a frozen exponential distribution instance with specified parameters
exp_obj = scipy.stats.expon(c,1/float(lamb))

x_in = np.linspace(0,1,200) # 200 numbers in [0,1], input for ppf()
y_out = exp_obj.ppf(x_in)
plt.plot(x_in,y_out) # graphically check the results of the inverse CDF

Upvotes: 7

Related Questions