Цензурная модель JAGS подсчитывает данные
Я кодирую иерархическую модель Пуассона в JAGS+R для данных подсчета для цензурированных данных.
A - это матрица, строки - это разные места, а столбцы - это разные интервалы времени, в которые я считаю дождливые дни. В качестве ковариат у меня есть набор матриц X_k$. Для каждого места и временного интервала у меня есть ковариаты, такие как средняя температура (сохраненная в X_1), число ветреных дней (в X_2), средняя влажность (в X_3).
Должен ли я отделить цензурированные данные от нецензурных? Как мне сделать это с матрицами?
Спасибо за вашу помощь!
Обновить:
У меня есть это в цикле до последнего, но один зарегистрированный интервал времени (для каждого места может быть цензура из-за отказа оборудования)
mu [i, j]<- a [1] * x [i, 1] + a [2] * x [i, 2] + b [1, j] * varx [i, j] + b [2, J] * varx [I, J]
N [i, j] <- dpois (лямбда [i,j])
log (лямбда [i,j]) <- mu [i, j] + alpha [i]
альфа [i]<- G0 [скрытый [i]]
скрытый [i]~ dcat(prob[])
повторил это за последний промежуток времени. добавленной
censoor[I]<- шаг (-censored [I])
в последнем бите интервала цензура [i] является вектором, который указывает, есть ли сбой оборудования
Я новичок, и это не работает, любая помощь? Спасибо