Python-Matplotlib boxplot.Как показать процентили 0,10,25,50,75,90 и 100? - PullRequest
Новая
       17

Python-Matplotlib boxplot.Как показать процентили 0,10,25,50,75,90 и 100?

7 голосов
/ 30 декабря 2011

Я хотел бы построить EPSgram (см. Ниже), используя Python и Matplotlib.

Функция boxplot отображает только квартили (0, 25, 50, 75, 100).Итак, как я могу добавить еще два поля?

EPSGram boxplot

1 Ответ

2 голосов
/ 09 февраля 2012

Я собрал образец, если вам все еще интересно. Он использует scipy.stats.scoreatpercentile, но вы можете получать эти цифры из других мест: 

from random import random
import numpy as np
import matplotlib.pyplot as plt
from scipy.stats import scoreatpercentile

x = np.array([random() for x in xrange(100)])

# percentiles of interest
perc = [min(x), scoreatpercentile(x,10), scoreatpercentile(x,25),
               scoreatpercentile(x,50), scoreatpercentile(x,75),
               scoreatpercentile(x,90), max(x)]
midpoint = 0 # time-series time

fig = plt.figure()
ax = fig.add_subplot(111)
# min/max
ax.broken_barh([(midpoint-.01,.02)], (perc[0], perc[1]-perc[0]))
ax.broken_barh([(midpoint-.01,.02)], (perc[5], perc[6]-perc[5]))
# 10/90
ax.broken_barh([(midpoint-.1,.2)], (perc[1], perc[2]-perc[1]))
ax.broken_barh([(midpoint-.1,.2)], (perc[4], perc[5]-perc[4]))
# 25/75
ax.broken_barh([(midpoint-.4,.8)], (perc[2], perc[3]-perc[2]))
ax.broken_barh([(midpoint-.4,.8)], (perc[3], perc[4]-perc[3]))

ax.set_ylim(-0.5,1.5)
ax.set_xlim(-10,10)
ax.set_yticks([0,0.5,1])
ax.grid(True)
plt.show()

Output of the code above

...