Квантиль вектора в стане

Я хотел бы использовать квантиль вектора в стане, но функция квантиля, похоже, не работает. См. ** ** в следующем примере.

data{
    vector[10] y;
    vector[10] x;
}

parameters{
    real a;
    real b;
}


model{
    vector[10] mu;
    real Q;
    mu = a*x+b;
    **Q = quantile(y-mu, 0.66);**
}