model{
	for( i in 1:n ){
		y[i] ~ dbin( p[i], 1)
		logit( p[i] ) <- beta[1] + beta[2] * c.temp[i]
	}
	for( j in 1:2 ){
		ptilde[j] ~ dbeta( a[j], b[j] )
	}
	beta[1] <- 0.2717391 * logit( ptilde[1] ) + 0.7282609  * logit( ptilde[2] )
	beta[2] <- -0.5 * logit( ptilde[1] ) + 0.05 * logit( ptilde[2] )
}