Репликация функции prod () R в WinBUGS

Привет и добро пожаловать на сайт!

Помните, что BUGS - это декларативный синтаксис, а не язык программирования, поэтому вы не можете перезаписывать значения переменных, как вы ожидаете, на таком языке, как R. Итак, вам нужно создать несколько промежуточных узлов, чтобы делать то, что вы вычисляете.

Если у вас есть следующие данные:

A <- [1,2,3,4]
nA <- 4

Затем вы можете включить в свою модель:

sumlogA[1] <- 0
for(i in 1:nA){
    sumlogA[i+1] <- sumlogA[i] + log(A[i])
}
prodA <- exp(sumlogA[nA+1])

Обратите внимание, что я работаю с логарифмической шкалой, а затем беру экспоненту суммы - это математически эквивалентно произведению, но является гораздо более стабильным в вычислительном отношении вычислением.

Надеюсь, это поможет,

Matt