Вопросы по теме 'jags'
dinterval () для данных с интервальной цензурой?
Я новичок в JAGS и пытаюсь понять, как dinterval () работает в JAGS для цензурированных данных. Я моделирую грубые данные, где у меня есть только верхняя и нижняя границы для каждой точки данных (а не истинное значение). Вот простой пример того,...
2804 просмотров
schedule
25.10.2023
Как написать цикл for, когда функция увеличивается с каждой итерацией?
Я пытаюсь оценить вероятность обнаружения животных из n.sites за несколько периодов наблюдения, когда животные удаляются, а обнаружение меняется во времени и пространстве. Это работает, если я делаю что-то подобное для 5 периодов наблюдения:...
479 просмотров
schedule
03.11.2022
Нижний и верхний пределы детерминированного узла в JAGS
Я пытаюсь усечь детерминированный узел theta.t , чтобы убедиться, что его значения находятся между 0 и 1. Я пытался решить проблему, используя T(0,1) в конце строки 6, но это не работает, и я Я не знаю, как использовать ~ dinterval() для решения...
259 просмотров
schedule
25.08.2022
Регрессия Пуассона с несколькими ковариатами с использованием JAGS: как упростить модель
У меня есть набор данных количества заявок с y в качестве количества требований, 16 ковариат, а именно от x1 до x16 (состоит из 0 и 1), которые я расположил в матрице дизайна под названием X и E в качестве воздействия (также называемого смещением). Я...
1308 просмотров
schedule
06.02.2022
Определение новых функций и дистрибутивов на языке BUGS/JAGS/STAN
Я новичок в мире статистического анализа и недавно заинтересовался языком моделирования BUGS/JAGS/STAN. Что меня действительно удивляет, так это то, что я не видел ни одного примера определения новых функций или дистрибутивов, чтобы избежать...
2034 просмотров
schedule
12.06.2023
Оценка неизвестной переменной отклика в JAGS — обучение без учителя
Я пытаюсь оценить значения отклика процентного покрытия ( COV ) по известным параметрам распределения. Я могу сделать это, указав данные ответа как NA в OpenBUGS (например, код ниже), но JAGS этого не допустит. Кто-нибудь знает, как я могу добиться...
185 просмотров
schedule
03.05.2023
Почему JAGS требует не менее двух цепочек для расчета DIC?
Я пытаюсь понять, как JAGS рассчитывает отклонение и информационный критерий отклонения (DIC).
Функция dic.samples в rjags в R выдает ошибку, если у вас есть только одна цепочка.
stop("2 or more parallel chains required")
Формула...
1071 просмотров
schedule
18.03.2023
set.seed в WinBUGS через R
Можно ли установить начальное значение в WinBUGS для воспроизведения оценок параметров, как это можно сделать с set.seed в R ?
503 просмотров
schedule
08.02.2023
Как прогнозировать значения, используя оценки из rjags/JAGS
После настройки модели и обучения ее с помощью Gibbs Sampling я получил результат всего предсказания скрытых значений с помощью:
jags <- jags.model('example.bug',
data = data,
n.chains = 4,
n.adapt =...
3212 просмотров
schedule
15.02.2024
Почему мне нужно определять параметр верхнего уровня в JAGS и как?
Согласно руководству пользователя r-jags (раздел Компиляция):
Любой узел, который используется в правой части отношения, но не определен в левой части какого-либо отношения, считается постоянным узлом. Его значение должно быть указано в файле...
91 просмотров
schedule
16.11.2023
R, JAGS, r2jags: доступ к последнему элементу в начале цикла for
Я работаю над проблемой дизайна эксперимента и пытаюсь подогнать модель JAGS через R и r2jags .
Чтобы измерить эффект переноса, я должен получить доступ к элементу i-1 в списке для одной из переменных. Когда i=1 , эта переменная должна...
58 просмотров
schedule
11.03.2023
Как структурировать модель и найти бета-версию байесовской модели в коде R?
Я пытаюсь структурировать модель Jags и найти бета-версию, а также априорные значения для байесовского анализа данных. В моей модели есть 3 предиктора x1, x2, x3 , а результатом является переменная распределения Бернулли Y . Как определить...
46 просмотров
schedule
26.11.2022
JAGS - применить функцию ко всем узлам параметров
Я новичок в JAGS, и я запускаю модель в R через пакет R2jags . Код модели основан на коде, взятом из Kéry & Schaub 2012 («Байесовский анализ населения с использованием WinBUGS»), стр. 399.
Вычисляется мера расхождения хи-квадрат.
model {...
192 просмотров
schedule
25.05.2023
Как я могу смоделировать переменную отклика из уже установленной модели?
Я уже подогнал регрессионную модель с помощью JAGS.
model{
for(i in 1:n) {
y[i] ~ dbeta(alpha[i], beta[i])
alpha[i] <- mu[i] * phi[i]
beta[i] <- (1 - mu[i]) * phi[i]
log(phi[i]) <- -inprod(X2[i, ],...
158 просмотров
schedule
01.05.2023
ошибка rjags Недопустимый аргумент вектора для ilogit
Я хотел бы сравнить регрессию betareg с той же регрессией с использованием rjags.
library(betareg)
d = data.frame(p= sample(c(.1,.2,.3,.4),100, replace= TRUE),
id = seq(1,100,1))
# I am looking to reproduce this regression with...
1284 просмотров
schedule
03.11.2022
Модели логистической регрессии JAGS с нецелыми весами
Для модели логистической регрессии ниже я хочу иметь возможность производить выборку из апостериорного, используя нецелочисленные значения для n (и y). Это может происходить в такой модели, когда доступны частичные данные или желательно снизить вес....
469 просмотров
schedule
17.05.2023
Ошибка: слайсер застрял на значении с бесконечной плотностью, запущенной биномиальной бета-моделью в JAGS
Я пытаюсь запустить биномиальную бета-модель в JAGS (см. Пример кода ниже). Я все время получаю сообщение об ошибке: Ошибка: при попытке запустить модель JAGS возникла следующая ошибка:
Error in node a0
Slicer stuck at value with infinite...
1734 просмотров
schedule
13.01.2024
Run Jags — извлечь несколько реализаций из объекта mcmc
У меня есть скрипт runjags, который генерирует предсказанную плотность нор для каждой клетки на острове. Я хочу получить несколько рисунков (около 100) от объекта mcmc для каждой ячейки. Мой научный руководитель считает, что я смогу сделать это с...
176 просмотров
schedule
08.11.2022
Извлечение коэффициентов из объекта mcmc
Извлечение данных из объектов всегда было для меня одним из самых запутанных аспектов R. Я установил модель байесовской линейной регрессии с использованием rjags и получил следующий объект mcmc:
summary(m_csim)
Iterations = 1:150000
Thinning...
493 просмотров
schedule
18.03.2023
Как я могу решить несоответствие размеров в модели Jags?
Я новичок в байесовском анализе и пытаюсь попрактиковаться на примере моделей Classic Capture-recapture: Mh2
это мой код
nind <- dim(venados)[1]
K <- 43
ntraps <- 13
M <- 150
nz <- M - nind
Yaug <- array(0, dim = c(M,...
74 просмотров
schedule
20.01.2023