bayesian - スタンのベクトルの分位数

Question

スタンでベクトルの分位数を使用したいのですが、機能分位数が機能していないようです。次の例の ** ** を参照してください。

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);**
}