Estoy tratando de hacer unGráfico de estimación de densidad de kernel (KDE) con seaborn y ubicar la mediana. El código se parece a esto:

import seaborn as sns
import numpy as np
import matplotlib.pyplot as plt

sns.set_palette("hls", 1)
data = np.random.randn(30)
sns.kdeplot(data, shade=True)

# x_median, y_median = magic_function()
# plt.vlines(x_median, 0, y_median)

plt.show()

Como puedes ver necesito unmagic_function() para obtener los valores medios x e y de lakdeplot. Entonces me gustaría trazarlos con, p.vlines. Sin embargo, no puedo entender cómo hacer eso. El resultado debería verse más o menos así (obviamente, la barra mediana negra está mal aquí):

Supongo que mi pregunta no está estrictamente relacionada con los nacimientos marinos y también se aplica a otros tipos de parcelas de matplotlib. Cualquier idea es muy apreciada.