# Bayesian Methods for the Physical Science. Learning from
# Examples in Astronomy and Physics. By S. Andreon and B. Weaver.
model {
intrscat <- 1/sqrt(prec.intrscat)
prec.intrscat ~ dgamma(1.0E-5,1.0E-5)
alpha ~ dnorm(0.0,1.0E-4)
beta ~ dt(0,1,1)
for (i in 1:length(obstot)) {

 obsbkg[i] ~ dpois(nbkg[i])
 obstot[i] ~ dpois(nbkg[i]/C[i]+10^lgn[i]) 
 nbkg[i] ~ dunif(0,3000)
 # modeling data selection
 lgM[i] ~ dnorm(14.5,pow(0.33,-2))
 # decomment for ignoring data selection 
 # lgM[i] ~dunif(13,16)

 obslgM[i] ~ dnorm(lgM[i],pow(errlgM[i],2))
 # modeling mass -n relation 
 lgnm[i] <- alpha+1.5 +beta*(lgM[i]-14.5)
 lgn[i] ~ dnorm(lgnm[i], pow(intrscat,-2))
}
}