Bugs>update(10000)
Bugs>monitor(mu1)
Bugs>monitor(mu2)
Bugs>monitor(sigma1)
Bugs>monitor(sigma2)
Bugs>monitor(y1)
Bugs>monitor(y2)
update(100000)
Bugs>stats(mu1)
                  mean        sd       2.5% :  97.5%  CI    median      sample
                1.161E+1   4.780E-1   1.067E+1   1.255E+1   1.161E+1     100000
Bugs>stats(mu2)
                  mean        sd       2.5% :  97.5%  CI    median      sample
                1.726E+1   4.593E-1   1.636E+1   1.816E+1   1.726E+1     100000
Bugs>stats(sigma1)
                  mean        sd       2.5% :  97.5%  CI    median      sample
                1.643E+0   2.840E-1   1.199E+0   2.305E+0   1.606E+0     100000
Bugs>stats(sigma2)
                  mean        sd       2.5% :  97.5%  CI    median      sample
                2.452E+0   2.537E-1   2.016E+0   3.011E+0   2.431E+0     100000
Bugs>stats(y1)
                  mean        sd       2.5% :  97.5%  CI    median      sample
                1.161E+1   1.802E+0   8.059E+0   1.518E+1   1.162E+1     100000
Bugs>stats(y2)
                  mean        sd       2.5% :  97.5%  CI    median      sample
                1.725E+1   2.710E+0   1.193E+1   2.262E+1   1.726E+1     100000


analisi con syst.R

> mysum(mu1)
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max.
  9.152  11.290  11.610  11.610  11.930  13.970
[1] media = 11.611828;   stddev = 0.477986
> mysum(mu2)
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max.
  15.24   16.95   17.26   17.26   17.56   19.54
[1] media = 17.257451;   stddev = 0.459352

> corr(mu1,mu2)
[1] cov = 0.089807;  rho = 0.409026

   => Nota sqrt(cov) = 0.3 , come deve:  -> sigma_z 


> mysum(sigma1)
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max.
 0.9055  1.4430  1.6060  1.6430  1.8020  4.7250
[1] media = 1.643035;   stddev = 0.283982
 > mysum(sigma2)
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max.
  1.621   2.274   2.431   2.452   2.606   4.079
[1] media = 2.452140;   stddev = 0.253751
> corr(sigma1,sigma2)
[1] cov = 0.000204;  rho = 0.002827

   => no corr, come deve


> mysum(y1)
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max.
 0.6481 10.4200 11.6200 11.6100 12.8000 21.7900
[1] media = 11.613516;   stddev = 1.802435
> mysum(y2)
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max.
  5.066  15.440  17.260  17.250  19.050  31.670
[1] media = 17.251885;   stddev = 2.710037

> corr(y1,y2)
[1] cov = 0.095038;  rho = 0.019456

   =>attese su devst di y1 e y2, indicando con s le stdev:
   
     s^2(y_i) =  s^2(mu_i) + sigma1^2 + sigma1z^2
1:      3.25 <-> 0.228 + 2.70 + 0.25 = 3.19  => ok
2:      7.34 <-> 0.228 + 6.01 + 1.00 = 7.24  => ok