model;
{
treat~dnorm(0,0.000001)
tau1~dgamma(0.001,0.001)
tau2~dgamma(0.001,0.001)
mean~dnorm(12,0.000001)
sd1<-sqrt(1/tau1)
sd2<-sqrt(1/tau2)
for (i in 1:20){
mu[i]<-treat*tr[i]+mucage[c[i]]
Y[i]~dnorm(mu[i],tau1)
}
for (i in 1:6){
mucage[i]~dnorm(mean,tau2)
}
}

c[] tr[]      Y[]	
1	0	12.99	
1	0	12.83	
1	1	14.10	
2	1	14.18	
2	1	14.37	
2	0	13.26	
2	0	13.94	
3	1	12.98	
3	1	13.54	
3	0	12.25	
3	0	12.73	
4	1	14.69	
4	1	14.27	
4	0	13.04	
4	0	13.60	
5	1	14.15	
5	1	15.01	
5	0	14.02	
6	1	13.22	
6	0	11.51	
END

list(tau1=1,tau2=1)