I would like to use the quantile of a vector in stan but the function quantile doesn't seem to work. See the ** ** in the following example.
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);**
}