# Negative Binomial mean model

model{

	# Likelihood

	for(i in 1:n.cc){

	log(mu[i]) <- beta[1]*intercept.cc[i] + beta[2]*share[i] +
		beta[3]*owner[i] + beta[4]*accounts[i]

	p[i] <- r/(r + mu[i])
	
	reports[i] ~ dnegbin(p[i], r)
	}

	# prior distributions

	for(i in 1:4){beta[i] ~ dnorm(0, 0.001)}
	
	r ~ dgamma(1, 1)
}