Schnelle Alternative für numpy.median.reduceat

Gibt es in Bezug auf diese Antwort eine schnelle Möglichkeit, Mediane über ein Array zu berechnen, das Gruppen mit einer ungleichen Anzahl von Elementen enthält? Z.B: data = [1.00, 1.05, 1.30, 1.20, 1.06, 1.54, 1.33, 1.87, 1.67, ... ] index = [0, 0, 1, 1, 1, 1, 2, 3, 3, ... ] Und dann möchte ich...