Я кодирую иерархическую модель Пуассона в 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(lambda[i,j])
log(lambda[i,j]) ‹- mu[i,j] ] + альфа[i]
альфа[i]‹- G0[латентный[i]]
латентный[i]~ dcat(вероятность[])
повторил это для последнего временного интервала. Добавлен
censored[i]‹-step(-censored[i])
в последний бит временного интервала, censored[i] — это vec, указывающий на сбой оборудования.
Я новичок, и это не работает, любая помощь? Спасибо