# Poisson mean model

model{

	# Likelihood

	for(i in 1:n.rec){

	trips[i] ~ dpois(lambda[i])

	log(lambda[i]) <- beta[1]*intercept.rec[i] + beta[2]*ski[i] +
		beta[3]*fee[i] + beta[4]*income[i]
	
	}

	# prior distributions

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

}