
model {
#eps is random effect mean zero variance sigma squared
#w is for spatial effects
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]+w[i] }
	
#spatial effects
w[1:N]~spatial.exp(mu1[],longitude[],latitude[],tau.sp,phi,1)

for (i in 1:N) {
	mu1[i]<-0}
tau.sp~dgamma(2.01,1.01)
sigma.sp<-1/tau.sp

phi~dunif(0.2,150)
range<-3/phi

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

}




