Python SciPy TransformedDensityRejection.ppf_hat用法及代碼示例

本文簡要介紹 python 語言中scipy.stats.sampling.TransformedDensityRejection.ppf_hat的用法。

用法:

TransformedDensityRejection.ppf_hat(u)#

計算 u 處帽子分布的 CDF 的倒數。

參數 ：：

u： array_like

百分位數數組

返回 ：：

ppf_hat： array_like

對應於給定百分位數的分位數數組。

例子：

>>> from scipy.stats.sampling import TransformedDensityRejection
>>> from scipy.stats import norm
>>> import numpy as np
>>> from math import exp
>>>
>>> class MyDist:
...     def pdf(self, x):
...         return exp(-0.5 * x**2)
...     def dpdf(self, x):
...         return -x * exp(-0.5 * x**2)
...
>>> dist = MyDist()
>>> rng = TransformedDensityRejection(dist)
>>>
>>> rng.ppf_hat(0.5)
-0.00018050266342393984
>>> norm.ppf(0.5)
0.0
>>> u = np.linspace(0, 1, num=1000)
>>> ppf_hat = rng.ppf_hat(u)