model {
	# this just has to be large enough to ensure all phi[i]'s > 0
	C <- 1000
	zeros <- 0
	nlogl <- 1/2*pow(1/delta*(log(p2)+log(1-p1)-log(1-p2)-log(p1)),2)-
		(alpha1-1)*log(p1)-(beta1-1)*log(1-p1)-
		(alpha2-1)*log(p2)-(beta2-1)*log(1-p2) + C
	zeros ~ dpois(nlogl)
	OR<-p2/(1-p2)*(1-p1)/p1
        DOR<-p1/(1-p1)*p2/(1-p2)
        RR <- p2/p1
        RD <- p2 - p1
	p1~dbeta(1,1)
	p2~dbeta(1,1)
}
