как определить следующую плотность в WinBUGS:
f(x|u)=abs(x)/(cos(pi*u)^2)e^{x^2/(cos(pi*u)^2)}
это не дистрибутивы WinBUGS по умолчанию. Я не знаю, как это определить. есть идеи?
model{
for(i in 1:n){
X[i]~f(x|u) ##abs(X[i]/pow(cos(pi*U[i]),2))*exp(-pow((X[i]/cos(pi*U[i])),2))
U[i]~dunif(aalpha,5)
}
#priors
aalpha~dunif(0,1)
}