mars.tensor.random.wald¶

mars.tensor.random.
wald
(mean, scale, size=None, chunk_size=None, gpu=None, dtype=None)[source]¶ Draw samples from a Wald, or inverse Gaussian, distribution.
As the scale approaches infinity, the distribution becomes more like a Gaussian. Some references claim that the Wald is an inverse Gaussian with mean equal to 1, but this is by no means universal.
The inverse Gaussian distribution was first studied in relationship to Brownian motion. In 1956 M.C.K. Tweedie used the name inverse Gaussian because there is an inverse relationship between the time to cover a unit distance and distance covered in unit time.
 Parameters
mean (float or array_like of floats) – Distribution mean, should be > 0.
scale (float or array_like of floats) – Scale parameter, should be >= 0.
size (int or tuple of ints, optional) – Output shape. If the given shape is, e.g.,
(m, n, k)
, thenm * n * k
samples are drawn. If size isNone
(default), a single value is returned ifmean
andscale
are both scalars. Otherwise,np.broadcast(mean, scale).size
samples are drawn.chunk_size (int or tuple of int or tuple of ints, optional) – Desired chunk size on each dimension
gpu (bool, optional) – Allocate the tensor on GPU if True, False as default
dtype (datatype, optional) – Datatype of the returned tensor.
 Returns
out – Drawn samples from the parameterized Wald distribution.
 Return type
Tensor or scalar
Notes
The probability density function for the Wald distribution is
\[P(x;mean,scale) = \sqrt{\frac{scale}{2\pi x^3}}e^ \frac{scale(xmean)^2}{2\cdotp mean^2x}\]As noted above the inverse Gaussian distribution first arise from attempts to model Brownian motion. It is also a competitor to the Weibull for use in reliability modeling and modeling stock returns and interest rate processes.
References
 1
Brighton Webs Ltd., Wald Distribution, http://www.brightonwebs.co.uk/distributions/wald.asp
 2
Chhikara, Raj S., and Folks, J. Leroy, “The Inverse Gaussian Distribution: Theory : Methodology, and Applications”, CRC Press, 1988.
 3
Wikipedia, “Wald distribution” http://en.wikipedia.org/wiki/Wald_distribution
Examples
Draw values from the distribution and plot the histogram:
>>> import matplotlib.pyplot as plt >>> import mars.tensor as mt >>> h = plt.hist(mt.random.wald(3, 2, 100000).execute(), bins=200, normed=True) >>> plt.show()