*data from Cohlen, BJ, Velde, ERT, 
Looman, CWN, Eijckemans, R, Habbema,
JDF. Crossover or parallel design in infertility trials? The discussion continues, 
Fertility and Sterility 1998; 70: 40-45.

#Model involving period effects only, period as a factor

	model{
	   for( i in 1 : N1 ) {
	   	outcome1[i] ~ dbern(p1[i])
	       logit(p1[i]) <- intercept+beta2*period2[i]+beta3*period3[i]+beta4*period4[i]+beta5*period5[i]+beta6*period6[i]+b[i]
		b[i] ~ dnorm(0, tau)
	   }
	
	
	for( i in 1 : N2 ){
		outcome2[i] ~ dbern(p2[i])
		logit(p2[i]) <- intercept+beta2*period22[i]+beta3*period3[i]+beta4*period4[i]+beta5*period5[i]+beta6*period6[i]+b[i]
		
		}
		
		for( i in 1 : N3 ){
		outcome3[i] ~ dbern(p3[i])
		logit(p3[i]) <- intercept+beta2*period2[i]+beta3*period33[i]+beta4*period4[i]+beta5*period5[i]+beta6*period6[i]+b[i]
		
		 }
	
	for( i in 1 : N4 ){
		outcome4[i] ~ dbern(p4[i])
		logit(p4[i]) <- intercept+beta2*period2[i]+beta3*period3[i]+beta4*period44[i]+beta5*period5[i]+beta6*period6[i]+b[i]
		}
	
	for( i in 1 : N5 ){
		outcome5[i] ~ dbern(p5[i])
		logit(p5[i]) <- intercept+beta2*period2[i]+beta3*period3[i]+beta4*period4[i]+beta5*period55[i]+beta6*period6[i]+b[i]
		
		 }
	
	for( i in 1 : N6 ){
		outcome6[i] ~ dbern(p6[i])
		logit(p6[i]) <- intercept+beta2*period2[i]+beta3*period3[i]+beta4*period4[i]+beta5*period5[i]+beta6*period66[i]+b[i]
		}
	
	  sigma<- 1/sqrt(tau) 
	   intercept ~ dnorm(0,1.0E-4)
	  beta2 ~ dnorm(0,1.0E-4)
	  beta3 ~ dnorm(0,1.0E-4)
	  beta4 ~ dnorm(0,1.0E-4)
	  beta5 ~ dnorm(0,1.0E-4)
	  beta6 ~ dnorm(0,1.0E-4)
	   tau ~ dgamma(1,1.0E-3)
	  log.sigma<- log(sigma)
	   
	   
	}	
	
	list(outcome1=c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0, 0, 0, 0, 0,
 0, 0, 0 ,0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1) ,
		
		
		period2=c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0, 0, 0, 0, 0,
 0, 0, 0 ,0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0),
		
		period3=c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0, 0, 0, 0, 0,
 0, 0, 0 ,0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0),
		
		period4=c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0, 0, 0, 0, 0,
 0, 0, 0 ,0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0),
		
		period5=c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0, 0, 0, 0, 0,
 0, 0, 0 ,0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0),
		
		period6=c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0, 0, 0, 0, 0,
 0, 0, 0 ,0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0),
		
		outcome2 =c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0, 0, 0, 0, 0,
 0, 0, 0 ,0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,0)  ,
		
		period22=c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1),
		
		outcome3 =c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0, 0, 0, 0, 0,
 0, 0, 0 ,0, 0, 0, 0, 0, 0 ,0 ,0, 1, 1,1,1,1,1,1,1,1,1,1,0,0 )  , 


period33=c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1),

outcome4 =c( 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0, 0, 0, 0, 0,
 0, 0, 0 ,0, 0, 0, 0, 1, 1,0,0,0,0)  , 

period44=c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1),


outcome5 =c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0,0,0,0,1, 1, 1, 1,0,0)  , 


period55=c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1),

outcome6 =c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 1, 1, 1 )  , 


period66=c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1),



		N1 = 74, N2=67, N3=59,N4=46, N5=40, N6=34),
		
		
		
		list(intercept=-4,beta2=1,beta3=1,beta4=1,beta5=1,beta6=1,tau=15,b=c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0, 0, 0, 0, 0,0, 0, 0 ,0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0, 0, 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0)) 
		
		
		
		
	
list(intercept=-1,beta2=0.3,beta3=0.3,beta4=0.3,beta5=0.3,beta6=0.3,tau=10, b=c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1))		



*********************************************************************************
*********************************************************************************
#Model involving both treatment and period effects, period as a factor


	model{
	   for( i in 1 : N1 ) {
	   	outcome1[i] ~ dbern(p1[i])
	       logit(p1[i]) <- intercept+beta1*treatment1[i]+beta2*period2[i]+beta3*period3[i]+beta4*period4[i]+beta5*period5[i]+beta6*period6[i]+b[i]
		b[i] ~ dnorm(0, tau)
	   }
	
	
	for( i in 1 : N2 ){
		outcome2[i] ~ dbern(p2[i])
		logit(p2[i]) <- intercept+beta1*treatment2[i]+beta2*period22[i]+beta3*period3[i]+beta4*period4[i]+beta5*period5[i]+beta6*period6[i]+b[i]
		
		}
		
		for( i in 1 : N3 ){
		outcome3[i] ~ dbern(p3[i])
		logit(p3[i]) <- intercept+beta1*treatment3[i]+beta2*period2[i]+beta3*period33[i]+beta4*period4[i]+beta5*period5[i]+beta6*period6[i]+b[i]
		
		 }
	
	for( i in 1 : N4 ){
		outcome4[i] ~ dbern(p4[i])
		logit(p4[i]) <- intercept+beta1*treatment4[i]+beta2*period2[i]+beta3*period3[i]+beta4*period44[i]+beta5*period5[i]+beta6*period6[i]+b[i]
		}
	
	for( i in 1 : N5 ){
		outcome5[i] ~ dbern(p5[i])
		logit(p5[i]) <- intercept+beta1*treatment5[i]+beta2*period2[i]+beta3*period3[i]+beta4*period4[i]+beta5*period55[i]+beta6*period6[i]+b[i]
		
		 }
	
	for( i in 1 : N6 ){
		outcome6[i] ~ dbern(p6[i])
		logit(p6[i]) <- intercept+beta1*treatment6[i]+beta2*period2[i]+beta3*period3[i]+beta4*period4[i]+beta5*period5[i]+beta6*period66[i]+b[i]
		}
	
	  sigma<- 1/sqrt(tau) 
	   intercept ~ dnorm(0,1.0E-4)
	   beta1 ~ dnorm(0,1.0E-4)
	  beta2 ~ dnorm(0,1.0E-4)
	  beta3 ~ dnorm(0,1.0E-4)
	  beta4 ~ dnorm(0,1.0E-4)
	  beta5 ~ dnorm(0,1.0E-4)
	  beta6 ~ dnorm(0,1.0E-4)
	   tau ~ dgamma(1,1.0E-4)
	  log.sigma<- log(sigma)
	   
	   
	}	
	
	list(outcome1=c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0, 0, 0, 0, 0,
 0, 0, 0 ,0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1) ,
		
		treatment1=c(0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,1,1,1,0,1,0,0,0,0,1,1,0,0,0,0,0,0,1,1,1,1,1,0,0,0,0,0,0,0,1,1,1,0,0,0,0,1,1,1),
		
		period2=c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0, 0, 0, 0, 0,
 0, 0, 0 ,0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0),
		
		period3=c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0, 0, 0, 0, 0,
 0, 0, 0 ,0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0),
		
		period4=c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0, 0, 0, 0, 0,
 0, 0, 0 ,0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0),
		
		period5=c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0, 0, 0, 0, 0,
 0, 0, 0 ,0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0),
		
		period6=c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0, 0, 0, 0, 0,
 0, 0, 0 ,0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0),
		
		outcome2 =c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0, 0, 0, 0, 0,
 0, 0, 0 ,0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,0)  ,
		
		treatment2=c(1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,0,0,0,1,0,1,1,1,1,0,0,1,1,1,1,1,1,0,0,0,0,0,1,1,1,1,1,1,1,0,0,0),
		
		period22=c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1),
		
		outcome3 =c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0, 0, 0, 0, 0,
 0, 0, 0 ,0, 0, 0, 0, 0, 0 ,0 ,0, 1, 1,1,1,1,1,1,1,1,1,1,0,0 )  , 

treatment3=c(0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,1,1,1,0,1,0,0,0,0,1,1,0,0,0,0,0,0,1,1,1,1,1,0,0),

period33=c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1),

outcome4 =c( 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0, 0, 0, 0, 0,
 0, 0, 0 ,0, 0, 0, 0, 1, 1,0,0,0,0)  , 

treatment4=c(1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,0,0,0,1,0,1,1,1,1,0,0), 

period44=c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1),


outcome5 =c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0,0,0,0,1, 1, 1, 1,0,0)  , 

treatment5=c(0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,1,1,1,0,1),

period55=c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1),

outcome6 =c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 1, 1, 1 )  , 

treatment6=c(1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1),

period66=c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1),



		N1 = 74, N2=67, N3=59,N4=46, N5=40, N6=34),
		
		
		
		list(intercept=20,beta1=20,beta2=20,beta3=20,beta4=40,beta5=40,beta6=40,tau=3,b=c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0, 0, 0, 0, 0,0, 0, 0 ,0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0, 0, 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0)) 
		
		
		
		
	
list(intercept=-1,beta1=0.3,beta2=0.3,beta3=0.3,beta4=0.3,beta5=0.3,beta6=0.3,tau=1, b=c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1))		

******************************************************************************************************
******************************************************************************************************
#Model involving period effects only, period having linear effect


	model{
	   for( i in 1 : N1 ) {
	   	outcome1[i] ~ dbern(p1[i])
	       logit(p1[i]) <- intercept+beta2*period1[i]+b[i]
		b[i] ~ dnorm(0, tau)
	   }
	
	
	for( i in 1 : N2 ){
		outcome2[i] ~ dbern(p2[i])
		logit(p2[i]) <- intercept+beta2*period2[i]+b[i]
		
		}
		
		for( i in 1 : N3 ){
		outcome3[i] ~ dbern(p3[i])
		logit(p3[i]) <- intercept+beta2*period3[i]+b[i]
		
		 }
	
	for( i in 1 : N4 ){
		outcome4[i] ~ dbern(p4[i])
		logit(p4[i]) <- intercept+beta2*period4[i]+b[i]
		}
	
	for( i in 1 : N5 ){
		outcome5[i] ~ dbern(p5[i])
		logit(p5[i]) <- intercept+beta2*period5[i]+b[i]
		
		 }
	
	for( i in 1 : N6 ){
		outcome6[i] ~ dbern(p6[i])
		logit(p6[i]) <- intercept+beta2*period6[i]+b[i]
		}
	
	  sigma<- 1/sqrt(tau) 
	   intercept ~ dnorm(0,1.0E-4)
	  beta2 ~ dnorm(0,1.0E-4)
	   tau ~ dgamma(1,1.0E-3)
	  log.sigma<- log(sigma)
	   
	   
	}	
	
	
	
	list(outcome1=c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0, 0, 0, 0, 0,
 0, 0, 0 ,0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1) ,
		
		
		
		period1=c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1),

		
		
		
		outcome2 =c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0, 0, 0, 0, 0,
 0, 0, 0 ,0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,0)  ,
		
		
period2=c(2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2),		
		
		
outcome3 =c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0, 0, 0, 0, 0,
 0, 0, 0 ,0, 0, 0, 0, 0, 0 ,0 ,0, 1, 1,1,1,1,1,1,1,1,1,1,0,0 )  , 


period3=c(3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3),

outcome4 =c( 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0, 0, 0, 0, 0,
 0, 0, 0 ,0, 0, 0, 0, 1, 1,0,0,0,0)  , 


period4=c(4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4),


outcome5 =c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0,0,0,0,1, 1, 1, 1,0,0)  , 



period5=c(5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5),

outcome6 =c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 1, 1, 1 )  , 


period6=c(6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6),



		N1 = 74, N2=67, N3=59,N4=46, N5=40, N6=34),
		
		
		
		list(intercept=-4,beta2=1,tau=15,b=c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0, 0, 0, 0, 0,0, 0, 0 ,0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0, 0, 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0)) 
		
		
		
		
	
list(intercept=-1,beta2=0.3,tau=10, b=c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1))		

**************************************************************************
**************************************************************************
#Model involving both tretment and period effects, period having linear effect


	model{
	   for( i in 1 : N1 ) {
	   	outcome1[i] ~ dbern(p1[i])
	       logit(p1[i]) <- intercept+beta1*treatment1[i]+beta2*period1[i]+b[i]
		b[i] ~ dnorm(0, tau)
	   }
	
	
	for( i in 1 : N2 ){
		outcome2[i] ~ dbern(p2[i])
		logit(p2[i]) <- intercept+beta1*treatment2[i]+beta2*period2[i]+b[i]
		
		}
		
		for( i in 1 : N3 ){
		outcome3[i] ~ dbern(p3[i])
		logit(p3[i]) <- intercept+beta1*treatment3[i]+beta2*period3[i]+b[i]
		
		 }
	
	for( i in 1 : N4 ){
		outcome4[i] ~ dbern(p4[i])
		logit(p4[i]) <- intercept+beta1*treatment4[i]+beta2*period4[i]+b[i]
		}
	
	for( i in 1 : N5 ){
		outcome5[i] ~ dbern(p5[i])
		logit(p5[i]) <- intercept+beta1*treatment5[i]+beta2*period5[i]+b[i]
		
		 }
	
	for( i in 1 : N6 ){
		outcome6[i] ~ dbern(p6[i])
		logit(p6[i]) <- intercept+beta1*treatment6[i]+beta2*period6[i]+b[i]
		}
	
	  sigma<- 1/sqrt(tau) 
	   intercept ~ dnorm(0,1.0E-4)
	   beta1 ~ dnorm(0,1.0E-4)
	  beta2 ~ dnorm(0,1.0E-4)
	   tau ~ dgamma(1,1.0E-4)
	  log.sigma<- log(sigma)
	   
	   
	}
		
	list(outcome1=c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0, 0, 0, 0, 0,
 0, 0, 0 ,0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1) ,
		
		treatment1=c(0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,1,1,1,0,1,0,0,0,0,1,1,0,0,0,0,0,0,1,1,1,1,1,0,0,0,0,0,0,0,1,1,1,0,0,0,0,1,1,1),
		
		period1=c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1),

		
		
		
		outcome2 =c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0, 0, 0, 0, 0,
 0, 0, 0 ,0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,0)  ,
		
		treatment2=c(1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,0,0,0,1,0,1,1,1,1,0,0,1,1,1,1,1,1,0,0,0,0,0,1,1,1,1,1,1,1,0,0,0),

period2=c(2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2),		
		
		
outcome3 =c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0, 0, 0, 0, 0,
 0, 0, 0 ,0, 0, 0, 0, 0, 0 ,0 ,0, 1, 1,1,1,1,1,1,1,1,1,1,0,0 )  , 

treatment3=c(0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,1,1,1,0,1,0,0,0,0,1,1,0,0,0,0,0,0,1,1,1,1,1,0,0),

period3=c(3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3),

outcome4 =c( 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0, 0, 0, 0, 0,
 0, 0, 0 ,0, 0, 0, 0, 1, 1,0,0,0,0)  , 

treatment4=c(1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,0,0,0,1,0,1,1,1,1,0,0), 

period4=c(4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4),


outcome5 =c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0,0,0,0,1, 1, 1, 1,0,0)  , 

treatment5=c(0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,1,1,1,0,1),

period5=c(5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5),

outcome6 =c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 1, 1, 1 )  , 

treatment6=c(1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1),

period6=c(6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6),



		N1 = 74, N2=67, N3=59,N4=46, N5=40, N6=34),
		
		
		
		list(intercept=-4,beta1=1,beta2=1,tau=15,b=c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0, 0, 0, 0, 0,0, 0, 0 ,0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0, 0, 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0)) 
		
		
		
		
	
list(intercept=-1,beta1=0.3,beta2=0.3,tau=10, b=c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1))		

********************************************************************************************
********************************************************************************************
#The null Model


	model{
	   for( i in 1 : N1 ) {
	   	outcome1[i] ~ dbern(p1[i])
	       logit(p1[i]) <- intercept+b[i]
		b[i] ~ dnorm(0, tau)
	   }
	
	
	for( i in 1 : N2 ){
		outcome2[i] ~ dbern(p2[i])
		logit(p2[i]) <- intercept+b[i]
		
		}
		
		for( i in 1 : N3 ){
		outcome3[i] ~ dbern(p3[i])
		logit(p3[i]) <- intercept+b[i]
		
		 }
	
	for( i in 1 : N4 ){
		outcome4[i] ~ dbern(p4[i])
		logit(p4[i]) <- intercept+b[i]
		}
	
	for( i in 1 : N5 ){
		outcome5[i] ~ dbern(p5[i])
		logit(p5[i]) <- intercept+b[i]
		
		 }
	
	for( i in 1 : N6 ){
		outcome6[i] ~ dbern(p6[i])
		logit(p6[i]) <- intercept+b[i]
		}
	
	  sigma<- 1/sqrt(tau) 
	   intercept ~ dnorm(0,1.0E-3)
	   tau ~ dgamma(1.0E-4,1.0E-4)
	  log.sigma<- log(sigma)
	   
	   
	}	
	
	
	list(outcome1=c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0, 0, 0, 0, 0,
 0, 0, 0 ,0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1) ,
		
		
		
		

		
		
		
		outcome2 =c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0, 0, 0, 0, 0,
 0, 0, 0 ,0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,0)  ,
		
		
		
		
outcome3 =c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0, 0, 0, 0, 0,
 0, 0, 0 ,0, 0, 0, 0, 0, 0 ,0 ,0, 1, 1,1,1,1,1,1,1,1,1,1,0,0 )  , 




outcome4 =c( 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0, 0, 0, 0, 0,
 0, 0, 0 ,0, 0, 0, 0, 1, 1,0,0,0,0)  , 





outcome5 =c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0,0,0,0,1, 1, 1, 1,0,0)  , 





outcome6 =c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 1, 1, 1 )  , 





		N1 = 74, N2=67, N3=59,N4=46, N5=40, N6=34),
		
		
		
		list(intercept=-4,tau=15,b=c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0, 0, 0, 0, 0,0, 0, 0 ,0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0, 0, 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0)) 
		
		
		
		
	
list(intercept=-1,tau=10, b=c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1))		

****************************************************************************************
****************************************************************************************
#Model involving treatment effect only


	model{
	   for( i in 1 : N1 ) {
	   	outcome1[i] ~ dbern(p1[i])
	       logit(p1[i]) <- intercept+beta1*treatment1[i]+b[i]
		b[i] ~ dnorm(0, tau)
	   }
	
	
	for( i in 1 : N2 ){
		outcome2[i] ~ dbern(p2[i])
		logit(p2[i]) <- intercept+beta1*treatment2[i]+b[i]
		
		}
		
		for( i in 1 : N3 ){
		outcome3[i] ~ dbern(p3[i])
		logit(p3[i]) <- intercept+beta1*treatment3[i]+b[i]
		
		 }
	
	for( i in 1 : N4 ){
		outcome4[i] ~ dbern(p4[i])
		logit(p4[i]) <- intercept+beta1*treatment4[i]+b[i]
		}
	
	for( i in 1 : N5 ){
		outcome5[i] ~ dbern(p5[i])
		logit(p5[i]) <- intercept+beta1*treatment5[i]+b[i]
		
		 }
	
	for( i in 1 : N6 ){
		outcome6[i] ~ dbern(p6[i])
		logit(p6[i]) <- intercept+beta1*treatment6[i]+b[i]
		}
	
	  sigma<- 1/sqrt(tau) 
	   intercept ~ dnorm(0,1.0E-3)
	   beta1 ~ dnorm(0,1.0E-3)
	   tau ~ dgamma(1.0E-4,1.0E-4)
	  log.sigma<- log(sigma)
	   
	   
	}	
	
	
	
	list(outcome1=c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0, 0, 0, 0, 0,
 0, 0, 0 ,0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1) ,
		
		treatment1=c(0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,1,1,1,0,1,0,0,0,0,1,1,0,0,0,0,0,0,1,1,1,1,1,0,0,0,0,0,0,0,1,1,1,0,0,0,0,1,1,1),
		
		

		
		
		
		outcome2 =c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0, 0, 0, 0, 0,
 0, 0, 0 ,0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,0)  ,
		
		treatment2=c(1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,0,0,0,1,0,1,1,1,1,0,0,1,1,1,1,1,1,0,0,0,0,0,1,1,1,1,1,1,1,0,0,0),

	
		
		
outcome3 =c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0, 0, 0, 0, 0,
 0, 0, 0 ,0, 0, 0, 0, 0, 0 ,0 ,0, 1, 1,1,1,1,1,1,1,1,1,1,0,0 )  , 

treatment3=c(0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,1,1,1,0,1,0,0,0,0,1,1,0,0,0,0,0,0,1,1,1,1,1,0,0),



outcome4 =c( 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0, 0, 0, 0, 0,
 0, 0, 0 ,0, 0, 0, 0, 1, 1,0,0,0,0)  , 

treatment4=c(1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,0,0,0,1,0,1,1,1,1,0,0), 




outcome5 =c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0,0,0,0,1, 1, 1, 1,0,0)  , 

treatment5=c(0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,1,1,1,0,1),



outcome6 =c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 1, 1, 1 )  , 

treatment6=c(1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1),




		N1 = 74, N2=67, N3=59,N4=46, N5=40, N6=34),
		
		
		
		list(intercept=-4,beta1=1,tau=15,b=c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0 ,0 ,0, 0, 0 ,0, 0, 0, 0, 0,0, 0, 0 ,0, 0, 0, 0, 0, 0 ,0 ,0, 0, 0, 0, 0, 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0)) 
		
		
		
		
	
list(intercept=-1,beta1=0.3,tau=10, b=c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1))