的数学表达式西皮·统计·瓦尔德.pdf文件？

Question

我使用scipy.stats.wald.fit来拟合我的数据数组，在没有任何警告或错误的情况下，我得到了两个参数（-1.004123，0.11517）。我绘制了分布图和直方图图，发现它们匹配得很好。你知道吗

但是，我的数据中有负值，wald分布应该只适用于正值http://en.wikipedia.org/wiki/Inverse_Gaussian_distribution。这个scipy.stats.wald.fit是如何工作的？有人能给我scipy.stats.wald.pdf的数学表达式吗？你知道吗

实际上我想用Cuda-C编写wald-pdf函数，所以需要一个数学表达式。你知道吗

Answer 1

如果有疑问，请看the source！你知道吗

第5900..5906行显示

The probability density function for `wald` is::

    wald.pdf(x, a) = 1/sqrt(2*pi*x**3) * exp(-(x-1)**2/(2*x))

for ``x > 0``.

`wald` is a special case of `invgauss` with ``mu == 1``.

Answer 2

{a1}的位置参数。你知道吗

使用stats.wald.pdf(x, loc=-1.004123, scale=0.11517)或者直接将(x - loc)/scale插入@Hugh Bothwell为您找到的公式。你知道吗

>>> data = stats.wald.rvs(-1.004123, 0.11517, size=1000)
>>> res = stats.wald.fit(data)  # cheating, I know
>>> x = np.linspace(min(data), max(data), 100)
>>> ax.plot(x, stats.wald.pdf(x, *res), 'r-')
>>> ax.plot(x, stats.wald.pdf(x, loc=res[0], scale=res[1]), 'b ')
>>> ax.hist(data, normed=True, histtype='stepfilled', alpha=0.3)
>>> plt.show()