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
omega~dt(0,1,1)T(0.001,) #scaled T-distr. (half), with phi=1 and df=1 (Johnson example script)

###draw cluster level values from mutivariate DP base distributions
for (k in 1:n){
delta[k,1:npar] ~ dmnorm.vcov(mu[1:npar],Omega.mat[1:npar,1:npar])
}

for (xx in 1:npar){
for (yy in 1:npar){
Omega.mat[xx,yy]<-omega^2 * R[xx,yy] ##R = (H'H)^-1; data; g-prior 
}}


##fixed parameters
beta[1:npar]~dmnorm.vcov(mu.b[1:npar],R.beta[1:npar,1:npar])


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])

}

}


















