bayesian

0зной

1ответ

Что такое `trait: units` и` variance.function (trait): units` в `MCMCglmm`?

Я пытаюсь выполнить MCMCglmm подобного df=data.frame(y=rep(c(0:2),each=12),x=rnorm(3*12)) MCMCglmm(data=df,fixed=y~x, family="categorical") и я получаю сообщение об ошибке For error structures invo

r statistics glm bayesian mcmc 2016-10-23

---

5зной

1ответ

Стэн. Использование target + = синтаксис

Я начинаю изучать Стэна. Может ли кто-нибудь объяснить, когда и как использовать синтаксис, такой как ...? target += вместо того, чтобы просто: y ~ normal(mu, sigma) Например, в Стан руководстве в

bayesian stan 2016-10-27

---

1зной

1ответ

Вероятностная инструкция ifelse в JAGS

Я хочу, чтобы переменной внутри модели JAGS присваивалось значение с вероятностью p [i] и другое значение с вероятностью 1-p [i]. p [i] - переменная в модели. Я пробовал следующее, но я не уверен, что

if-statement bayesian jags 2016-10-27

---

1зной

1ответ

95% Регион HPD с использованием R TeachingDemos Пакет

Мне нужно найти 95% -ный регион HPD, используя пакет TeachingDemos в R. У меня есть последующее распределение, следующее за гамма-распределением. После установки пакета и ввод в библиотеке (TeachingDe

r bayesian 2016-11-01

---

0зной

1ответ

Избегание ловушек оптимизации при моделировании порядковой предсказанной переменной в PyMC3

Я пытаюсь моделировать порядковую предсказанную переменную с использованием PyMC3 на основе подхода в главе 23 Doing Bayesian Data Analysis. Я хотел бы определить хорошее начальное значение, используя

bayesian pymc pymc3 2016-11-04

---

1зной

2ответ

Ошибка выполнения в 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: см here, если вы заи

r bayesian jags winbugs 2016-11-10

---