15.12.2012 Views

scipy tutorial - Baustatik-Info-Server

scipy tutorial - Baustatik-Info-Server

scipy tutorial - Baustatik-Info-Server

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

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

Parameters<br />

a : array<br />

moment : int<br />

axis : int or None<br />

Returns<br />

The appropriate moment along the given axis or over all values if axis is :<br />

None. :<br />

mquantiles(a, prob=, [0.25, 0.5, 0.75], alphap=0.40000000000000002, betap=0.40000000000000002,<br />

axis=None, limit=())<br />

Computes empirical quantiles for a data array.<br />

Samples quantile are defined by Q(p) = (1 − g).x[i] + g.x[i + 1], where x[j] is the *j*th order statistic, and i<br />

= (floor(n*p+m)), m=alpha+p*(1-alpha-beta) and g = n*p + m - i.<br />

Typical values of (alpha,beta) are:<br />

• (0,1) : p(k) = k/n : linear interpolation of cdf (R, type 4)<br />

• (.5,.5) : p(k) = (k+1/2.)/n : piecewise linear function (R, type 5)<br />

• (0,0) : p(k) = k/(n+1) : (R type 6)<br />

• (1,1) : p(k) = (k-1)/(n-1). In this case, p(k) = mode[F(x[k])]. That’s R default (R type 7)<br />

• (1/3,1/3): p(k) = (k-1/3)/(n+1/3). Then p(k) ~ median[F(x[k])]. The resulting quantile estimates are<br />

approximately median-unbiased regardless of the distribution of x. (R type 8)<br />

• (3/8,3/8): p(k) = (k-3/8)/(n+1/4). Blom. The resulting quantile estimates are approximately unbiased<br />

if x is normally distributed (R type 9)<br />

• (.4,.4) : approximately quantile unbiased (Cunnane)<br />

• (.35,.35): APL, used with PWM<br />

Parameters<br />

a : array-like<br />

Input data, as a sequence or array of dimension at most 2.<br />

prob : array-like, optional<br />

List of quantiles to compute.<br />

alpha : float, optional<br />

Plotting positions parameter, default is 0.4.<br />

beta : float, optional<br />

Plotting positions parameter, default is 0.4.<br />

axis : int, optional<br />

Axis along which to perform the trimming. If None (default), the input array is first<br />

flattened.<br />

limit : tuple<br />

Tuple of (lower, upper) values. Values of a outside this closed interval are ignored.<br />

Returns<br />

quants : MaskedArray<br />

680 Chapter 3. Reference

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

Saved successfully!

Ooh no, something went wrong!