[R] Binary data on JAGS
Massimiliano Tripoli
m@@@|m|||@no@tr|po|| @end|ng |rom gm@||@com
Tue Jun 1 19:52:13 CEST 2021
Thanks, for suggestions (To define p before).
I thought it and I put this line of code inside the loop:
p[i] <- y[i] / n[i]
but I received a RUNTIME ERROR and:
(Attempt to redefine node p[1]).
Thanks
Massimiliano
Il giorno mar 1 giu 2021 alle ore 17:22 Bert Gunter <bgunter.4567 using gmail.com>
ha scritto:
> Where is p defined before it is used? (Is this part of what jags provides
> somehow?)
>
> Bert Gunter
>
> "The trouble with having an open mind is that people keep coming along and
> sticking things into it."
> -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip )
>
>
> On Tue, Jun 1, 2021 at 7:57 AM Massimiliano Tripoli <
> massimiliano.tripoli using gmail.com> wrote:
>
>> Dear R users,
>> I'm trying to reproduce the example 6.5.1 (Dobson (1983)) in BUGS book in
>> linux using JAGS.
>> Below the code as
>>
>>
>>
>> https://www.mrc-bsu.cam.ac.uk/software/bugs/the-bugs-project-the-bugs-book/bugs-book-examples/the-bugs-book-examples-chapter-6-6-5-1/
>>
>> # By R code:
>> library('rjags')
>>
>> jags <- jags.model(
>> file = "651.bug",
>> list (x = c(1.6907, 1.7242, 1.7552, 1.7842, 1.8113, 1.8369, 1.8610,
>> 1.8839),
>> n = c(59, 60, 62, 56, 63, 59, 62, 60),
>> y = c(6, 13, 18, 28, 52, 53, 61, 60)) ,
>> inits = list(beta = 0, alpha = 50),
>> n.chains = 3
>> )
>>
>> # By JAGS code: (651.bug)
>>
>> model {
>> for (i in 1:8) {
>> y[i] ~ dbin(p[i], n[i])
>> logit(p[i]) <- alpha + beta*(x[i] - mean(x[]))
>> phat[i] <- y[i]/n[i]
>> yhat[i] <- n[i]*p[i]
>> }
>> alpha ~ dnorm(0, 0.0001)
>> beta ~ dnorm(0, 0.0001)
>> }
>> # I received this error in R:
>>
>> # Resolving undeclared variables
>>
>> # Allocating nodes
>>
>> #Graph information:
>>
>> # Observed stochastic nodes: 8
>>
>> # Unobserved stochastic nodes: 2
>>
>> # Total graph size: 78
>>
>>
>>
>> #Initializing model
>>
>> #Deleting model
>>
>>
>>
>> #Error in jags.model(file = "651.bug", list(x = c(1.6907, 1.7242, 1.7552,
>> :
>> # Error in node y[1]
>>
>> #Node inconsistent with parents
>>
>> Thanks in advance
>> M. Tripoli
>>
>> [[alternative HTML version deleted]]
>>
>> ______________________________________________
>> R-help using r-project.org mailing list -- To UNSUBSCRIBE and more, see
>> https://stat.ethz.ch/mailman/listinfo/r-help
>> PLEASE do read the posting guide
>> http://www.R-project.org/posting-guide.html
>> and provide commented, minimal, self-contained, reproducible code.
>>
>
[[alternative HTML version deleted]]
More information about the R-help
mailing list