# model from Oring_model_a.txt
model{
  for( i in 1:23){
  y[i]  ~ dbin( p[i], 1)
  logit( p[i] ) <- beta[1] + beta[2] * temp[i]
  }
  for( j in 1:2 ){
  ptilde[j]  ~ dbeta( a[j], b[j] )
  }
  beta[1] <- (75/20) * logit( ptilde[1] )
                 - (55/20) * logit( ptilde[2] )
  beta[2] <- (-1/20) * logit( ptilde[1] )
                 + (1/20) * logit( ptilde[2] )}

# data redefined
list(y=c(0,1,0,0,0,0,0,0,1,1,1,0,0,1,0,0,0,0,
0,0,1,0,1),
temp=c(66,70,69,68,67,72,73,70,57,63,70,78,67,
53,67,75,70,81,76,79,75,76,58),
a=c(1.000,0.577), b=c(0.577,1.000))

# initial values (inits) redefined
list( ptilde = c( 0.5, 0.5 ) )