model {
	for (i in 1:N) {   # loop over cell-lines
			
		for (j in 1:t) {  # loop over concentrations
				### Level 1 Model: within cell line model, four parameter logistic
			y[j,i]~dnorm(eta[j,i], invsigma2)
			eta[j,i] <- beta[i,1] + (beta[i,2])/(1+exp(beta[i,4]*(log(x[j])-beta[i,3])))  
				}
									
			# Level 2 Model: between cell line model, genotype fixed effect for all parameters
			
		beta[i, 1] <- exp(alpha[i, 1])
		beta[i, 2] <- exp(alpha[i, 2])
		beta[i, 3] <- alpha[i, 3]+alpha1[1]*gene.exp[i]
		beta[i, 4]<- alpha[i,4]+alpha1[2]*gene.exp[i]
			
		alpha[i,1:4]~ dmnorm(mu[1:4], invV[1:4, 1:4])
			
		}
		

	mu[1:4] ~ dmnorm(mean[1:4], invSigma[1:4, 1:4])
	invV[1:4, 1:4] ~ dwish(R[1:4, 1:4], 4)
	V[1:4, 1:4] <- inverse(invV[1:4, 1:4]) 
			
	alpha1[1]~dnorm(0, 0.00001)
	alpha1[2]~dnorm(0, 0.00001)
	invsigma2 ~ dgamma(0.001, 0.001)		
}