def rarmax(vector): m = np.amax(vector) indices = np.nonzero(vector == m) [0] return pr.choice(indices)