model{ for(j in 1:Ns){ for(i in 1:Nt){ Resp[j,i]~dbin(pA[j,i],Nk) Pred[j,i]~dbin(pA[j,i],Nk) pA[j,i]<-pnorm(percept[j,i],c[j],tau) percept[j,i]<-a[j]*st[i,1]+b[j]*st[i,2] } } # priors tau<-pow(sigma,-2) sigma~dunif(0,50) for(j in 1:Ns){ a[j]~dnorm(0,0.001) b[j]~dnorm(0,0.001) c[j]~dnorm(0,0.001) } }