calculate percentiles

[mathause]

How would I calculate percentiles over a certain dimension?

ds.groupby('time.month').apply(np.percentile, q=5)

Does not work.

[clarkfitzg]

ds = xray.tutorial.load_dataset('air_temperature')
percents = ds.groupby('time.month').reduce(np.percentile, q=5)

Producing:

In [80]: percents
Out[80]:
<xray.Dataset>
Dimensions:  (month: 12)
Coordinates:
  * month    (month) int64 1 2 3 4 5 6 7 8 9 10 11 12
Data variables:
    air      (month) float64 239.7 239.3 244.3 251.7 261.4 271.5 274.0 272.5 ...

[shoyer]

Adding a dedicated method for calculating percentiles/quantiles would also be a good idea: #303

[mathause]

Thanks