Предварительное условие при условии другого предшествующего

Question

Я пытаюсь выбрать два параметра (предшествующих) из категориального распределения от 1 до 5000, тета [1] и тета [2] с требованием, чтобы theta1

Я пробовал (среди прочего):

theta[1] ~ dcat(p1[])
p1[1:n] <- 1/n 


theta[2] ~ dcat(p2[])
pi2[1:theta[1]] <- 0 
pi2[sum(theta[1],1):n] <- 1/sum(n, -pi1) 

с n = 5000, так что theta2 выбирается из категориального распределения от theta1 до n.

Ошибка: неизвестная переменная тета [1].

Любая помощь будет оценена.

Answer 1

Если в этой категориальной переменной с n=5000, единственное требование - theta1<theta2, вы можете использовать order() функция:

theta.star[1] ~ dcat(p1[])
theta.star[2] ~ dcat(p1[])
theta <- order(theta.star)

В order() функция - это способ наложить ограничения порядка в JAGS.