Вопросы по теме 'winbugs'
Логит-регрессия с WinBUGS/R2WinBUGS, не сходящейся к оценкам правдоподобия
Я пользователь Win-7 с R 2.15.2.
Может ли кто-нибудь помочь мне, почему следующая модель плохо сходится с оценками простой логит-модели?
Отредактировано
Mydata <- structure(list(gg = c(13.659955, 6.621436486, 3.017166776, 2.516795069,...
1108 просмотров
schedule
18.06.2022
rJAGS-версия write.model()?
В пакете R2WinBUGS есть функция write.model() . Пакет R rjags не имеет такой функции, о которой я знаю. write.model создает временный текстовый файл, который может быть прочитан WinBUGS как модель.
Я знаю, что могу ввести write.model в...
672 просмотров
schedule
27.06.2022
Ошибка сообщения WinBUGS: ожидаемая структура ключевого слова
Я новичок в WinBUGS, и мне не удалось заставить работать приведенный ниже код. Модель syntactically correct (это иерархическая логит-модель со случайными эффектами), но когда я нажимаю load data , появляется сообщение об ошибке expected key word...
3062 просмотров
schedule
27.03.2023
Несколько определений узлов winbugs
У меня есть проблема с этим кодом в winbugs. Модель синтаксически верна, и данные загружаются, но когда я компилирую, программный вывод представляет собой «несколько определений узла Z». Я не знаю, как решить проблему.
Это модель:
#BUGS...
1601 просмотров
schedule
19.11.2022
Будет ли JAGS оценивать все родительские узлы dcat или только нужный?
Скажем, у нас есть следующее утверждение:
for (i in 1:N) {
pi[i,1] <- ....
pi[i,2] <- ....
pi[i,3] <- ....
...
pi[i,100] <- ...
Y[i] ~ dcat(p[i,])
}
Допустим, что Y[1] = 5 . Будут ли jags оценивать все...
1251 просмотров
schedule
01.05.2023
Ошибка Openbugs слишком много констант
Я пытаюсь запустить приведенный ниже код, в котором я использовал нулевой трюк для своего выражения логарифмического правдоподобия (фи - это мое логарифмическое правдоподобие):
model{
for (l in 1:k) {
d.1[l] ~ dbern(p.1)
d.2[l] ~ dbern(p.2)...
130 просмотров
schedule
14.03.2023
Модель BUGS для (вложенных?) повторных измерений ANOVA
Мне было интересно, есть ли у кого-нибудь код для модели BUGS/JAGS для повторных измерений ANOVA? По сути, у меня есть ответ ( y ), который я хочу смоделировать с учетом времени суток , дня и лечения . >. Я также хотел бы включить два термина...
694 просмотров
schedule
08.02.2024
Явно укажите функцию логита в WinBUGS/OpenBUGS
Я новичок в OpenBUGS, и у меня возникла проблема с подгонкой модели под функцию logit().
Почитав вокруг, я обнаружил, что одним из возможных решений для этого было бы явное указание логит-функции без использования собственной логит-функции WinBUGS:...
1025 просмотров
schedule
02.04.2023
Ошибка WinBUGS: векторное отношение z должно включать последовательные элементы переменной
Я пытаюсь смоделировать многомерную модель Probit для своих двоичных данных. Я пробовал все, но WinBUGS, в свою очередь, выдает мне эту ошибку. Любые идеи или предложения тепло приветствуются.
model{ for (i in 1:ns){ ## цикл по исследованиям...
330 просмотров
schedule
06.02.2023
Ошибка выполнения в JAGS
Я пытаюсь сделать это в JAGS:
z[l] ~ dbeta(0.5,0.5)
y[i,l] ~ z[l]*dnorm(0,10000) + inprod(1-z[l],dnegbin(exp(eta_star[i,l]),alpha[l]))
( dnorm(0,10000) моделирует дельту Дирака в 0: см. здесь если вам интересна модель).
Но я получаю:...
1281 просмотров
schedule
25.05.2023
Предварительная предиктивная проверка в JAGS - ошибка несоответствия размеров
Мне интересно, может ли кто-нибудь помочь мне решить эту проблему. Я хотел бы сделать оценку соответствия для модели нормального логарифма Пуассона, представленной ниже (это всего лишь простая тестовая модель). Когда я закомментирую строки fit ‹-...
277 просмотров
schedule
12.11.2022
Трехуровневая вложенная модель смешанных эффектов
Я пытаюсь смоделировать трехуровневую вложенную линейную модель смешанных эффектов в rjags (трехуровневая: несколько наблюдений за несколькими людьми в нескольких группах). В группах есть уникальные наборы людей.
Эквивалентная модель в lme4...
1658 просмотров
schedule
21.11.2022
Репликация функции prod() R в WinBUGS
Как с помощью WinBUGS вычислить произведение всех значений в одном векторе?
Я попытался использовать цикл for по тому же вектору.
Например:
В R, если A ‹- [1,2,3,4], prod(A) = 24.
Однако,
в BUGS, если a ‹- 2 и for (i in 1:n){ a ‹- a *...
144 просмотров
schedule
05.01.2023
индекс массива больше, чем верхняя граница массива для d
Может ли кто-нибудь помочь мне с этой проблемой, пожалуйста? Я использовал другие наборы данных с другим количеством исследований (NS) и обработок (NT), и все работало нормально. Любая помощь будет высоко оценена.
Набор данных выглядит следующим...
239 просмотров
schedule
26.07.2023
R, интервал HPD (наивысшая апостериорная плотность) на основе выборок из апостериорных данных, WinBUGS
Как рассчитать интервал HPD (наивысшая апостериорная плотность) из апостериорных выборок? У меня есть четыре параметра, и я генерирую 1000 образцов из распределения апостериорных параметров. Теперь Как рассчитать HPD в программном обеспечении R. Я...
440 просмотров
schedule
20.05.2023
Есть ли способ получить и сохранить позиции матричного элемента в JAGS?
Я разрабатываю байесовскую иерархическую модель в R с кодом BUGS в JAGS.
В моей модели у меня есть две матрицы, которые содержат релевантную информацию друг о друге в одной и той же позиции матрицы. Моя информация структурирована по строкам. Я...
34 просмотров
schedule
14.07.2023