
model {
for (i in 1:N) {
	positives[i]~dbin(p[i], examined[i])
	logit(p[i])<-b[1]+b[2]*alt[i]+b[3]*temp[i]+b[4]*rain[i]+b[5]*ndvi[i] }


for (i in 1:5) {
	b[i]~dnorm(0.0,0.01) 
	or[i]<-exp(b[i]) }




#Predict
for (i in 1:P) {
	logit(p.pred[i])<-b[1]+b[2]*alt.pred[i]+b[3]*temp.pred[i]+b[4]*rain.pred[i]+b[5]*ndvi.pred[i] }

}