Skip to main content
`Ctrl`+`K`
scipy.stats.mstats.

# plotting_positions

scipy.stats.mstats.plotting_positions(data, alpha=0.4, beta=0.4)

  Returns plotting positions (or empirical percentile points) for the data.

  Plotting positions are defined as `(i-alpha)/(n+1-alpha-beta)`, where:

  - i is the rank order statistics

  - n is the number of unmasked values along the given axis

  - alpha and beta are two parameters.

  Typical values for alpha and beta are:

  - (0,1) : `p(k) = k/n`, linear interpolation of cdf (R, type 4)

  - (.5,.5) : `p(k) = (k-1/2.)/n`, piecewise linear function (R, type 5)

  - (0,0) : `p(k) = k/(n+1)`, Weibull (R type 6)

  - (1,1) : `p(k) = (k-1)/(n-1)`, in this case, `p(k) = mode[F(x[k])]`. That’s R default (R type 7)

  - (1/3,1/3): `p(k) = (k-1/3)/(n+1/3)`, then `p(k) ~ median[F(x[k])]`. The resulting quantile estimates are approximately median-unbiased regardless of the distribution of x. (R type 8)

  - (3/8,3/8): `p(k) = (k-3/8)/(n+1/4)`, Blom. The resulting quantile estimates are approximately unbiased if x is normally distributed (R type 9)

  - (.4,.4) : approximately quantile unbiased (Cunnane)

  - (.35,.35): APL, used with PWM

  - (.3175, .3175): used in scipy.stats.probplot

  Parameters:
    data array_like

      Input data, as a sequence or array of dimension at most 2.

    alpha float, optional

      Plotting positions parameter. Default is 0.4.

    beta float, optional

      Plotting positions parameter. Default is 0.4.

  Returns:
    positions MaskedArray

      The calculated plotting positions.