Я пытаюсь использовать WinBUGS из R через BRugs и R2WinBUGS , код следующий: require(R2WinBUGS) require(BRugs) model<-function(){ for(i in 1:N){ y[i] <- x[i] + w[i] w[i] ~ dnorm(0, sigma.y) x[i] <- a - b*5 + v[i]...

Я пытаюсь оценить вероятность обнаружения животных из n.sites за несколько периодов наблюдения, когда животные удаляются, а обнаружение меняется во времени и пространстве. Это работает, если я делаю что-то подобное для 5 периодов наблюдения:...

Я новичок в мире статистического анализа и недавно заинтересовался языком моделирования BUGS/JAGS/STAN. Что меня действительно удивляет, так это то, что я не видел ни одного примера определения новых функций или дистрибутивов, чтобы избежать...

Я пытаюсь оценить значения отклика процентного покрытия ( COV ) по известным параметрам распределения. Я могу сделать это, указав данные ответа как NA в OpenBUGS (например, код ниже), но JAGS этого не допустит. Кто-нибудь знает, как я могу добиться...

Я пытаюсь использовать Rstan, чтобы соответствовать модели из примеров Байесовских идей и анализа данных Кристенсена, Джонсона, Бранскума и Хэнсона: Введение для ученых и статистиков. Авторы используют WinBUGS, поэтому необходима некоторая адаптация....

Можно ли установить начальное значение в WinBUGS для воспроизведения оценок параметров, как это можно сделать с set.seed в R ?

Я использую WinBUGS/R для разработки байесовской иерархической модели с тремя уровнями, но я изо всех сил пытаюсь найти достойные примеры кода. Кто-нибудь может предложить некоторые, пожалуйста? Я новичок в WinBUGS, но не в многоуровневом...

Я относительно новичок в R и OpenBugs и потратил много времени на устранение неполадок этой модели. Я смог выяснить изрядное количество из них самостоятельно с помощью онлайн-ресурсов, но я застрял на этой ошибке. В нем говорится, что это «несколько...

Я новичок в JAGS, и я запускаю модель в R через пакет R2jags . Код модели основан на коде, взятом из Kéry & Schaub 2012 («Байесовский анализ населения с использованием WinBUGS»), стр. 399. Вычисляется мера расхождения хи-квадрат. model {...

В настоящее время просматриваю примеры тома 1 и наткнулся на ошибку в примере с красителями. Когда я пытаюсь загрузить inits из примера, он возвращает «эта цепочка содержит неинициализированные переменные. Я не уверен, какая часть ее неверна, так...

Я хочу запустить иерархическую модель Пуассона для данных о хоккейных воротах. вот модель, настроенная в ошибках: modelString <- "model { for(i in 1:n_obs){ hockey_goals[i] ~ dpois(p[i]) log(p[i]) <- p_inter + p_age * age[i] +...

Как бороться с ошибкой компиляции, которая гласит: «среднее значение многомерного нормального Y [1,1] должно иметь то же количество компонентов, что и Y [1,1]»? Вот мой код и данные в WinBugs: model{ for (i in 1:10){ Y[i, 1:8] ~...