--- a/classic-bugs/vol1/epil/epil2.bug
+++ b/classic-bugs/vol1/epil/epil2.bug
@@ -6,11 +6,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  * Trt[j] * log.Base4[j]
+                            + alpha.Age * log.Age[j]  
+                            + alpha.V4  * V4[k]
                             + b1[j];
          y[j,k] ~ dpois(mu[j,k]);
       }
@@ -18,25 +18,13 @@
       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);
    tau.b1      ~ dgamma(1.0E-3,1.0E-3); 
    sigma.b1 <- 1.0/sqrt(tau.b1);
-# 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;
 }