Welcome to BUGS on 13 th Jun 2006  at 1:12:12
BUGS : Copyright (c) 1992 .. 1995 MRC Biostatistics Unit.
All rights reserved.
Version 0.600 for unix systems.
For general release: please see documentation for disclaimer.
The support of the Economic and Social Research Council (UK)
is gratefully acknowledged.
Bugs>compile("pfit.bug")
model fit.piano;

const
   N = 19;    # number of observations

var
   x[N], y[N], Z[N],
   mu[N],
   a, b, c, tau, sigma;

data x, y, Z in "dati.dat";

inits in "pfit.in";

{
   for (i in 1:N) {
      mu[i] <- a + b*x[i] + c*y[i];
      Z[i]  ~ dnorm(mu[i],tau);
   }
   a ~ dnorm(0.0, 1.0E-14);
   b ~ dnorm(0.0, 1.0E-10);
   c ~ dnorm(0.0, 1.0E-10);
   tau   ~ dgamma(1.0E-5, 1.0E-5);
   sigma <- 1.0/sqrt(tau);
}

Parsing model declarations.
Loading data value file(s).
Loading initial value file(s).
Parsing model specification.
Checking model graph for directed cycles.
Generating code.
Generating sampling distributions.
Generating initial values
Checking model specification.
Choosing update methods.
compilation took  00:00:00 
Bugs>update(100000)      time for    100000 updates was  00:00:03 
Bugs>monitor(a)
Bugs>monitor(b)
Bugs>monitor(c)
Bugs>monitor(sigma)
Bugs>update(100000)      time for    100000 updates was  00:00:03 
Bugs>stats(a)
                  mean        sd       2.5% :  97.5%  CI    median      sample
                4.795E+5   3.458E+2   4.788E+5   4.802E+5   4.795E+5     100000
Bugs>stats(b)
                  mean        sd       2.5% :  97.5%  CI    median      sample
               -5.635E+0   5.717E-3  -5.646E+0  -5.624E+0  -5.635E+0     100000
Bugs>stats(c)
                  mean        sd       2.5% :  97.5%  CI    median      sample
                9.996E+2   9.327E-1   9.978E+2   1.001E+3   9.996E+2     100000
Bugs>stats(sigma)
                  mean        sd       2.5% :  97.5%  CI    median      sample
                9.543E+2   1.835E+2   6.763E+2   1.390E+3   9.272E+2     100000
Bugs>q()