# This WinBUGS code is a companion to the paper “Bayesian models for pooling microarray studies with
# multiple sources of replications” by Erin M. Conlon, Joon J. Song and Jun S. Liu, BMC Bioinformatics
# (2006), 7:247.  Questions regarding this code can be addressed to the corresponding author Erin M. Conlon. 
# Please see the accompanying README file named README.BayesPoolMicro.pdf for additional information.

# Example1 Model
# Note: The parameters in this model are specifically for the data of Example1

model
	{
		for ( i in 1:N ) {
		
		
		y[1,i] ~ dnorm(mu[1,i], tau1[i])
		y[2,i] ~ dnorm(mu[1,i], tau1[i])
		y[3,i] ~ dnorm(mu[1,i], tau1[i])
		mu[1,i] ~ dnorm(theta1[i],sigma1[i])
		
		y[4,i] ~ dnorm(mu[2,i], tau1[i])
		mu[2,i] ~ dnorm(theta1[i],sigma1[i])
		
		y[5,i] ~ dnorm(mu[3,i], tau1[i])
		mu[3,i] ~ dnorm(theta1[i],sigma1[i])
		
		y[6,i] ~ dnorm(mu[4,i], tau2[i])
		y[7,i] ~ dnorm(mu[4,i], tau2[i])
		mu[4,i] ~ dnorm(theta2[i],sigma2[i])
		
		y[8,i] ~ dnorm(mu[5,i], tau2[i])
		mu[5,i] ~ dnorm(theta2[i],sigma2[i])
		
		y[9,i] ~ dnorm(mu[6,i], tau2[i])
		mu[6,i] ~ dnorm(theta2[i],sigma2[i])		
		
    	        theta1[i] ~ dnorm(0,eta1[i])
	
			eta1[i]<-etas1[T1[i]]	
			
	        theta2[i] ~ dnorm(0,eta2[i])
	
			eta2[i]<-etas2[T1[i]]			
			
		T[i] ~ dbern(p)
		T1[i]<-T[i]+1
           				
		}
		 
		p ~ dbeta(1,1)
			
		for ( i in 1:N ){

		tau1[i] ~ dgamma(1.5,0.1101) 
		tau2[i] ~ dgamma(1.5,0.03) 

		}
		
		for ( i in 1:N ){

		sigma1[i] ~ dgamma(1.5,0.04185) 
		sigma2[i] ~ dgamma(1.5,0.03045) 

		}
		
		eta21 ~ dgamma(.1,.1)
		eta22 ~ dgamma(.1,.1)

	        etas1[1]<-eta21
	        etas1[2]<-eta21*c1
                c1 ~ dgamma(1,100)
	
	        etas2[1]<-eta22
                etas2[2]<-eta22*c2
                c2 ~ dgamma(1,100)

	}