model
{
for(j in l:N){
z[j]~dbin(p[j],n[j])
logit(p[j])<-beta[1]+beta[2]*lst[j]+beta[3]*ndvi[j]+beta[4]*dist[j]+beta[5]*pipe[j]+u[loc[j]]+v[year[j]]
}
#Spatial random effects：
u[1:NL]~spatial.exp(mu[],x[],y[],tau.u,phi,1)
for(i in 1:NL){
mu<-0
}
#AR(1)prior distribution for temporal random effects：
v[1]~dnorm(0.0,tau2)
for(i in 2:NY){
vmean[i-1]<-ro*v[i-1];
v~dnorm(vmean[i-1],tau.v)
}
tau2<-1/((1-pow(ro,2))*pow(sigma.v,2))
ro~dunif(-1,1)
tau.v~dgamma(0.1,0.1)
sigma.v<-1/sqrt(tau.v)
#validation or prediction
for(j in l:m){
z.p[j]~dbin(p.p[j],n.p[j])
logit(p[j])<-beta[1]+beta[2]*lst.p[j]+beta[3]*ndvi.p[j]+beta[4]*dist.p[j]+
beta[5]*pipe.p[j]+ u.p[loc.p[j]]+v[year.p[j]]
}
for(j in 1:m1){
u.p[j]~spatial.unipred(mu.p[j],x.p[j],y.p[j],u[])
mu.p[j]<-0
}
#priors
Tau.u~dgamma(0.1,0.1)
#wior range for correlation at min distance(O．6kin)is O．01 to O．99
#prior range for correlation at max distance(49kin)is 0 to 0．2
for(i in 1:10){
beta~dnorm(0,0.1)
}
}

这个是二项分布的时空模型，想做poisson怎么修改呀，请各位大神赐教啊