model
{

	for (i in 1:n) {
		Y[i] ~ dnorm(mu[i], sigmasq.inv)
		mu[i] <- beta[i]*gamma[i]
		beta[i] ~ dnorm(0, tausq.inv)
		gamma[i] ~ dbern(q0)
	}

sigmasq.inv ~ dgamma(0.01, 0.01)
tausq.inv ~ dgamma(0.01, 0.01)

sigmasq <- 1/sigmasq.inv
tausq <- 1/tausq.inv

q0 <- 1-p0
p0 ~ dbeta(1,1)

}

### Inits and Data for n=25 ###

##Inits
list(
beta = 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),  
gamma = 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),  
p0 = 0.5, sigmasq.inv=1.0, tausq.inv=1.0)

##Data


### End for n=25 ###

### Inits and Data for n=100 ###

## Inits


## Data


### End for n=100###

### Inits and data for n=500 ###

##Inits


##Data


### End for n=500 ###

### Inits and Data for n=5000 ###

##Inits


##Data


### End for n=5000 ####