model {
  for(i in 1:length(w)){
     logloss[i]~dnorm(mu[i],1/sig2[i])
  }
#
# set up sig2
#
  for (i in 1:length(w)){
    sig2[i]<-sigd2[d[i]]
  }
  for (j in 1:10){
    sigd2[j]<-sum(a[j:10])
  }
  for (k in 1:10){
    a[k]~dunif(0,1)
  }
  
#
# set up log-means     
# 
  for (i in 1:length(w)){
    mu[i]<-level[w[i]] + dev[d[i]]
  }
#
# specify priors
#
  for (i in 1:numlev){
    level[i]~dunif(0,maxlev)
    }
  dev[1]<-0 
  for (i in 2:10){
    dev[i]~dunif(-5,5)
    }
}