TYPICAL START

C:\E\BJ\TEX\classes\707\week15 > bugs
Welcome to BUGS on 1 st Dec 1998  at 11:24:20
BUGS : Copyright (c) 1992 .. 1995 MRC Biostatistics Unit.
All rights reserved.
Version 0.600 for 32 Bit PC.
For general release: please see documentation for disclaimer.
The support of the Economic and Social Research Council (UK)
is gratefully acknowledged.
Bugs>compile("outmix.bug")

[ program listing here ]

Parsing model declarations.
Loading data value file(s).
Warning -- expected data read before end of file
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:01

-------------------------------------------------------------------



RESULTS FOR NORMAL ERRORS MODEL:


#     Y[i] ~ dt(mu[i],ttau[i],df[i]); # uncomment this line for outlier detect
      Y[i] ~ dt(mu[i],tau,ndf);       # uncomment this line for std analysis


Bugs>update(500)
      time for    500    updates was  00:00:02
Bugs>monitor(b0)

Bugs>monitor(b)

Bugs>monitor(sigma)

Bugs>update(1000)
      time for    1000   updates was  00:00:03
Bugs>stats(b0)

                  mean        sd       2.5% :  97.5%  CI    median      sample
                1.945E+0   2.208E-1   1.460E+0   2.370E+0   1.953E+0     1000
Bugs>stats(b)

                  mean        sd       2.5% :  97.5%  CI    median      sample
[1]            -3.422E-1   3.252E-1  -9.829E-1   2.935E-1  -3.450E-1     1000
[2]            -5.818E-2   3.036E-1  -6.262E-1   5.661E-1  -6.208E-2     1000
[3]            -2.688E-2   3.122E-1  -6.321E-1   6.027E-1  -3.072E-2     1000
Bugs>stats(sigma)

                  mean        sd       2.5% :  97.5%  CI    median      sample
                7.928E-1   9.627E-2   6.275E-1   1.006E+0   7.825E-1     1000
Bugs>

-----------------------------------------------------------------------

RESULTS FOR OUTLIER DETECTION

      Y[i] ~ dt(mu[i],ttau[i],df[i]); # uncomment this line for outlier detect
#     Y[i] ~ dt(mu[i],tau,ndf);       # uncomment this line for std analysis


Bugs>update(500)
      time for    500    updates was  00:00:04
Bugs>monitor(b0)

Bugs>monitor(b)

Bugs>monitor(sigma)

Bugs>monitor(Z)

Bugs>monitor(p)

Bugs>update(1000)
      time for    1000   updates was  00:00:09

[ Z never seems to get updated ]

------------------------------------------------------------------------