# Used by RcodeDiasorinExample.txt
 model
 {
 for(i in 1:n1){
   low[i] ~ dlnorm(muL, tauL)
   }
 for(i in 1:n2){
   normal[i] ~ dlnorm(muN, tauN)
   }
 muL ~ dnorm(4.87, 347.22)
 muN ~ dnorm(5.39, 357.14)
 tauL ~ dgamma(1.0376,0.001)
 tauN ~ dgamma(1.04653,0.001)
 medianL <- exp(muL)
 medianN <- exp(muN)
 relmedian <- medianN / medianL
 relprec <- tauN/tauL
 test <- step(medianN-medianL)
 lowf ~ dlnorm(muL, tauL)
 normalf ~ dlnorm(muN, tauN)
 }