15.12.2012 Views

scipy tutorial - Baustatik-Info-Server

scipy tutorial - Baustatik-Info-Server

scipy tutorial - Baustatik-Info-Server

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

SciPy Reference Guide, Release 0.8.dev<br />

Examples<br />

loc : array-like, optional<br />

location parameter (default=0)<br />

scale : array-like, optional<br />

scale parameter (default=1)<br />

size : int or tuple of ints, optional<br />

shape of random variates (default computed from input arguments )<br />

moments : string, optional<br />

composed of letters [’mvsk’] specifying which moments to compute where ‘m’<br />

= mean, ‘v’ = variance, ‘s’ = (Fisher’s) skew and ‘k’ = (Fisher’s) kurtosis. (default=’mv’)<br />

>>> import matplotlib.pyplot as plt<br />

>>> numargs = truncnorm.numargs<br />

>>> [ a,b ] = [0.9,]*numargs<br />

>>> rv = truncnorm(a,b)<br />

Display frozen pdf<br />

>>> x = np.linspace(0,np.minimum(rv.dist.b,3))<br />

>>> h=plt.plot(x,rv.pdf(x))<br />

Check accuracy of cdf and ppf<br />

>>> prb = truncnorm.cdf(x,a,b)<br />

>>> h=plt.semilogy(np.abs(x-truncnorm.ppf(prb,a,b))+1e-20)<br />

Random number generation<br />

>>> R = truncnorm.rvs(a,b,size=100)<br />

Truncated Normal distribution.<br />

The standard form of this distribution is a standard normal truncated to the range [a,b] — notice<br />

that a and b are defined over the domain of the standard normal. To convert clip values for a specific<br />

mean and standard deviation use a,b = (myclip_a-my_mean)/my_std, (myclip_b-my_mean)/my_std<br />

596 Chapter 3. Reference

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!