Welcome to BUGS on 5 th Jun 2005  at 18:51:9
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("molla.bug")
model fit.lineare;

const
   N  = 8, # numero di misure (da 3 a 10 dischetti) 
           # [per i primi due non c'e' allungamento]
   Nmin = 3,
   pi2 = 9.869604;   

var 
   l[N], t[N],        # dati
   m.disco, m.molla,  # dati 'costanti'
   c.l, m.l,          # parametri retta allungamento
   tau.l, sigma.l,
   c.t, m.t,          # parametri retta periodi
   tau.t, sigma.t,
   mu.l[N], mu.t[N],  # valori legati deterministicamente alle masse
   k, g;              # costante molla e accelerazione gravita'

data in "molla.dat.S";
inits in "molla.in";

{
   # l Vs m
   for (i in 1:N) {
      mu.l[i] <- c.l + m.l * (m.molla + (Nmin-1 + i) * m.disco);
      l[i] ~dnorm(mu.l[i], tau.l);
   }
   c.l ~ dnorm(0.0, 1.0E-4);
   m.l ~ dnorm(0.0, 1.0E-4);

   tau.l ~ dgamma(1.0E-3, 1.0E-6);
   sigma.l <- 1/sqrt(tau.l);

   # t vs sqrt(m)
   for (i in 1:N) {
      mu.t[i] <- c.t + m.t * sqrt(m.molla + (Nmin-1 + i) * m.disco);
      t[i] ~dnorm(mu.t[i], tau.t);
   }
   c.t ~ dnorm(0.0, 1.0E-4);
   m.t ~ dnorm(0.0, 1.0E-4);

   tau.t ~ dgamma(1.0E-3, 1.0E-5);
   sigma.t <- 1/sqrt(tau.t);

   # k e g
   k <- 4*pi2 / (m.t*m.t)
   g <- m.l * k 
}  

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(50000)      time for    50000  updates was  00:00:01 
Bugs>monitor(c.l)
Bugs>monitor(m.l)
Bugs>monitor(sigma.l)
Bugs>monitor(c.t)
Bugs>monitor(m.t)
Bugs>monitor(sigma.t)
Bugs>monitor(k)
Bugs>monitor(g)
Bugs>update(100000)      time for    100000 updates was  00:00:02 
Bugs>stats(c.l)
                  mean        sd       2.5% :  97.5%  CI    median      sample
               -5.275E-2   1.249E-3  -5.527E-2  -5.027E-2  -5.273E-2     100000
Bugs>stats(m.l)
                  mean        sd       2.5% :  97.5%  CI    median      sample
                2.229E-1   2.073E-3   2.188E-1   2.271E-1   2.229E-1     100000
Bugs>stats(sigma.l)
                  mean        sd       2.5% :  97.5%  CI    median      sample
                1.003E-3   3.623E-4   5.627E-4   1.925E-3   9.219E-4     100000
Bugs>stats(c.t)
                  mean        sd       2.5% :  97.5%  CI    median      sample
               -3.231E-2   1.317E-2  -5.872E-2  -5.996E-3  -3.223E-2     100000
Bugs>stats(m.t)
                  mean        sd       2.5% :  97.5%  CI    median      sample
                9.683E-1   1.736E-2   9.336E-1   1.003E+0   9.682E-1     100000
Bugs>stats(sigma.t)
                  mean        sd       2.5% :  97.5%  CI    median      sample
                5.507E-3   2.000E-3   3.072E-3   1.061E-2   5.062E-3     100000
Bugs>stats(k)
                  mean        sd       2.5% :  97.5%  CI    median      sample
                4.215E+1   1.517E+0   3.923E+1   4.529E+1   4.212E+1     100000
Bugs>stats(g)
                  mean        sd       2.5% :  97.5%  CI    median      sample
                9.396E+0   3.499E-1   8.731E+0   1.011E+1   9.388E+0     100000
Bugs>q()