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 at 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 at 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]]