--- a/classic-bugs/vol1/epil/epil3.bug
+++ b/classic-bugs/vol1/epil/epil3.bug
@@ -9,11 +9,11 @@
 model {
    for(j in 1:N) {
       for(k in 1:T) {
-         log(mu[j,k]) <- a0 + alpha.Base0 * (log.Base4[j] - log.Base4.bar)   
-                            + alpha.Trt0 * (Trt[j] - Trt.bar)  
-                            + alpha.BT  * BT.resid[j]
-                            + alpha.Age * (log.Age[j] - log.Age.bar)  
-                            + alpha.V4  * (V4[k] - V4.bar) 
+         log(mu[j,k]) <- alpha0 + alpha.Base * log.Base4[j]   
+                            + alpha.Trt * Trt[j] 
+                            + alpha.BT  * log.Base4[j] * Trt[j]
+                            + alpha.Age * log.Age[j]  
+                            + alpha.V4  * V4[k] 
                             + b1[j] + b[j,k];
          y[j,k] ~ dpois(mu[j,k]);
          b[j,k] ~ dnorm(0.0, tau.b);
@@ -22,15 +22,10 @@
       log.Base4[j] <- log(Base[j]/4);
       log.Age[j] <- log(Age[j]);
    }
-# covariate means:
-   log.Age.bar <- mean(log.Age[]);                
-   Trt.bar       <- mean(Trt[]);                   
-   log.Base4.bar <- mean(log.Base4[]);           
-   V4.bar        <- mean(V4[]);                  
 # priors:
-   a0          ~ dnorm(0.0,1.0E-4); 		           
-   alpha.Base0 ~ dnorm(0.0,1.0E-4);            
-   alpha.Trt0  ~ dnorm(0.0,1.0E-4);            
+   alpha0      ~ dnorm(0.0,1.0E-4); 		           
+   alpha.Base  ~ dnorm(0.0,1.0E-4);            
+   alpha.Trt   ~ dnorm(0.0,1.0E-4);            
    alpha.BT    ~ dnorm(0.0,1.0E-4);            
    alpha.Age   ~ dnorm(0.0,1.0E-4);            
    alpha.V4    ~ dnorm(0.0,1.0E-4);
@@ -38,11 +33,4 @@
    sigma.b1 <- 1.0/sqrt(tau.b1);
    tau.b       ~ dgamma(1.0E-3,1.0E-3); 
    sigma.b <- 1.0/sqrt(tau.b);
-# re-calculate coefficients on original scale: 
-   alpha0 <- a0 - alpha.Base0 * log.Base4.bar 
-                - alpha.Trt0 * Trt.bar 
-                - alpha.Age * log.Age.bar 
-                - alpha.V4 * V4.bar;
-   alpha.Trt <- alpha.Trt0 - beta.Trt * alpha.BT;
-   alpha.Base <- alpha.Base0 - beta.Base * alpha.BT;
 }