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@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.
