2016-02-10 2 views
0

Я хотел бы использовать квантиль вектора в stan, но функция quantile не работает. См. ** ** в следующем примере.Quantile of a vector in stan

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

Почему вы хотите генерировать квантиль на каждом этапе процесса MCMC? Просто возьмите квантил из заднего распределения. –

ответ

1

quantile Функция не существует на языке СТ.

Язык Стэна фокусируется на выражении дифференцируемых статистических моделей. Функция квантиля не дифференцируема.