data{
for (xx in 1:npar){
mu[xx]<-0  ##mean for G prior
mu.b[xx]<-0
}
}

model{

##prior for detection
p~dbeta(1,1)

##prior, g-prior parameter
tau~dgamma(0.1, 0.1)

###draw cluster level values from mutivariate DP base distributions
for (k in 1:n){
delta[k,1] ~ dnorm(mu[1],tau)
}


##fixed parameters
beta[1]~dnorm(mu.b[1],0.1)


for (i in 1:n){

	for(j in 1:J){
	lpsi[i,j]<-inprod(delta[i,],X[j,]) + inprod(beta,X[j,]) #frst term: species level effects, second term: fixed (avg) effects  
	psi[i,j]<-pnorm(lpsi[i,j],0,1)  #probit link
	z[i,j]~dbern(psi[i,j])
	p.eff[i,j]<-p*z[i,j]
	y[i,j]~dbin(p.eff[i,j],K)

	}

Nocc[i]<-sum(z[i,1:J])

}

}


















