# Used by RcodeDiasorinExample.txt
 model
 {
 for(i in 1:n1){
   low[i] ~ dlnorm(muL1, tauL1)
   }
 for(i in 1:n2){
   normal[i] ~ dlnorm(muN1, tauN1)
   }
 muL1 ~ dnorm(4.87, 0.001)
 muN1 ~ dnorm(5.39, 0.001)
 tauL1 ~ dgamma(1.0376,0.001)
 tauN1 ~ dgamma(1.04653,0.001)
 medianL1 <- exp(muL1)
 medianN1 <- exp(muN1)
 relmedian1 <- medianN1 / medianL1
 test1 <- step(medianN1-medianL1)
 lowf1 ~ dlnorm(muL1, tauL1)
 normalf1 ~ dlnorm(muN1, tauN1)
 }