model
{
	for(i in 1:37) 
	{
		resp1[i] ~ dnorm(mu1,prec)
	}
	for(j in 1:42) 
	{
		resp2[j] ~ dnorm(mu2,prec)
	}
	mu1 ~ dunif(50,150)
	mu2 ~ dunif(50,150)
        v ~ dunif(0.01,5000)
        prec <- 1/v
        std  <- sqrt(1/prec)

}
