#
# Code will run in posterior sampling mode.
#
# A scale parameter will be fit to the bkg model.
#
# The total number of Gibbs draws is 100, every 1th draws will be saved.
#
# The model will be fit using the Multi Scale Prior.
#
# The data matrix is 32 by 32.
#
# The data file should contain a  2^5 by 2^5 matrix of counts.
#
# Starting Values for the smoothing parameter (alpha):
# Aggregation Level:  0,   alpha: 1 (Full Data)
# Aggregation Level:  1,   alpha: 1
# Aggregation Level:  2,   alpha: 1
# Aggregation Level:  3,   alpha: 1
# Aggregation Level:  4,   alpha: 1  (In the 2x2 table)
#
# The prior distribution on the total count from the multiscale component is
# Gamma(1.000000, 0.050000).
#
# The hyper-prior smoothing parameter (kappa 2) is 1000.
#
iteration logPost stepSize cycleSpinRow cycleSpinCol smoothingParam0 smoothingParam1 smoothingParam2 smoothingParam3 smoothingParam4 expectedMSCounts bkgScale 
1        nan              0 24 7 0.041498 0.069668 0.009699 0.070901 0.006929 932.247 0.996418 
2    66.5734            nan 17 9 0.083929 0.094157 0.075079 0.052572 0.090303 960.488 0.986589 
3   -109.667   -176.2399665 15 19 0.044499 0.012064 0.009548 0.091328 0.053779 962.989 0.957962 
4   -9.45613    100.2104795 12 22 0.070877 0.098171 0.065221 0.091063 0.067813 995.527 0.946676 
5   -26.0333   -16.57716426 13 12 0.039021 0.084046 0.016827 0.039142 0.068502 1016.26 0.916823 
6    59.4649    85.49823946 24 25 0.083611 0.050053 0.016664 0.078747 0.050088 987.514 0.948896 
7   -48.2457   -107.7106895 31 13 0.038766 0.094246 0.122111 0.049400 0.077196 1002.45 0.983044 
8    2.04551    50.29124926 25 12 0.003680 0.023495 0.057253 0.078452 0.071850 966.696 0.973481 
9       -inf           -inf 20 5 0.036658 0.064900 0.044893 0.028359 0.096293 920.149 0.959964 
10   -7.68269            inf 25 19 0.042937 0.105968 0.049982 0.118455 0.092208 1001.81 0.937206 
11     9.0086    16.69129447 0 14 0.053994 0.124554 0.087217 0.116594 0.076534 1000.45 0.960948 
12   -2.48216   -11.49075765 25 27 0.015602 0.014107 0.030248 0.094982 0.035316 960.974 0.983065 
13    53.5858    56.06792912 31 10 0.069321 0.048132 0.090939 0.116249 0.140315 1006.89 0.932607 
14    33.3158   -20.26993464 9 18 0.043866 0.077591 0.014238 0.065697 0.047648 993.859 0.946764 
15   -6.36457   -39.68041217 10 2 0.014391 0.021455 0.056514 0.043493 0.080810 971.182 0.932841 
16    1.09834    7.462909722 15 21 0.018304 0.059591 0.078740 0.025571 0.066341 972 0.952708 
17    1.81875    0.720415407 31 10 0.055331 0.092180 0.057259 0.065255 0.029652 955.89 0.956464 
18   -44.7083   -46.52705932 19 9 0.111673 0.020771 0.052810 0.100201 0.020471 889.881 0.987486 
19    -17.837    26.87130963 10 21 0.030181 0.071208 0.032027 0.028269 0.048950 1055.95 0.928222 
20    88.8721     106.709109 21 14 0.072016 0.040783 0.014885 0.066219 0.051210 1019.18 0.947735 
21   -2.32097   -91.19308106 20 2 0.019769 0.031293 0.058257 0.034182 0.059763 948.486 0.93514 
22   -4.16208   -1.841110707 18 9 0.060167 0.078597 0.033559 0.052262 0.039356 1013.93 0.958712 
23    59.0434     63.2054367 2 22 0.103970 0.071964 0.038809 0.103745 0.054655 964.846 0.954279 
24   -71.7422   -130.7855404 15 28 0.029227 0.098999 0.043369 0.069111 0.070413 974.179 0.992278 
25     63.688    135.4301597 24 31 0.038753 0.048788 0.026274 0.064537 0.115890 1000.4 0.973245 
26   -1.74389   -65.43187072 17 10 0.049942 0.004005 0.040790 0.093723 0.062178 965.582 0.967314 
27       -inf           -inf 17 3 0.070879 0.048348 0.056620 0.045891 0.118931 952.448 0.947931 
28  0.0402851            inf 16 2 0.038032 0.021686 0.034831 0.062329 0.093312 980.247 0.943168 
29    -34.479   -34.51930691 22 17 0.032978 0.018720 0.009687 0.073888 0.042415 990.366 0.960284 
30       -inf           -inf 30 3 0.029127 0.050854 0.047296 0.043762 0.111438 968.087 0.943268 
31     44.511            inf 22 10 0.066899 0.048326 0.049195 0.096622 0.051108 957.043 0.930972 
32    59.9534    15.44245126 2 1 0.102280 0.034073 0.022458 0.059598 0.092226 946.732 0.985412 
33   -34.7163   -94.66978638 15 30 0.074130 0.125252 0.042295 0.080836 0.053831 947.747 0.991618 
34    -42.057   -7.340696567 21 13 0.036662 0.030230 0.029061 0.086430 0.040477 984.374 0.939785 
35    33.3824    75.43938805 17 24 0.055000 0.054530 0.024576 0.035630 0.063781 1019.56 0.975581 
36   -31.7246   -65.10696415 26 1 0.071353 0.087134 0.106237 0.037832 0.045031 926.317 0.966152 
37   -102.428   -70.70310992 24 31 0.048547 0.036788 0.047724 0.059655 0.105540 986.457 0.979482 
38     -5.384    97.04372398 3 15 0.013067 0.049262 0.130219 0.011615 0.018899 918.983 0.931688 
39    17.2415    22.62550037 26 4 0.058719 0.102688 0.011737 0.076419 0.039024 998.141 0.974298 
40    41.9007    24.65920473 15 4 0.059313 0.073311 0.080249 0.047044 0.043027 957.294 0.93726 
41   -37.3787   -79.27936266 29 4 0.075544 0.047612 0.093903 0.019202 0.074097 961.509 0.978679 
42    54.7418    92.12043367 27 27 0.028364 0.069834 0.093903 0.068339 0.056216 1004.71 0.923247 
43   -3.16807   -57.90985489 8 29 0.075589 0.024170 0.064346 0.112657 0.043962 982.391 0.903265 
44   -53.7358   -50.56772255 12 30 0.014025 0.015825 0.009955 0.072316 0.129384 988.297 0.939654 
45   -25.6712    28.06463514 6 8 0.089107 0.042308 0.125129 0.015308 0.118094 986.003 0.95583 
46    15.7266    41.39775632 9 18 0.112442 0.069154 0.108902 0.078669 0.060683 969.396 0.949068 
47    32.1378    16.41121623 11 7 0.015918 0.098402 0.013047 0.093585 0.099032 967.196 0.965912 
48   -4.57243   -36.71024176 25 19 0.032207 0.065291 0.057883 0.060554 0.098077 1035.71 0.975526 
49    5.20641    9.778845802 21 17 0.043717 0.015670 0.061033 0.024580 0.073144 979.142 0.933092 
50    -8.9451   -14.15150986 31 11 0.063190 0.041114 0.025728 0.090979 0.100246 1045.44 0.97738 
51    26.4935    35.43862462 30 20 0.050335 0.024695 0.015147 0.086452 0.102340 999.777 0.982031 
52    8.66576    -17.8277684 13 13 0.005106 0.059493 0.062756 0.100018 0.149734 923.188 0.969605 
53       -inf           -inf 4 27 0.080332 0.043514 0.141029 0.105726 0.111623 944.122 0.971125 
54    15.2711            inf 3 7 0.037365 0.109705 0.062020 0.007517 0.026015 961.221 0.953864 
55   -47.1369   -62.40800491 31 31 0.105943 0.117230 0.031641 0.017102 0.080161 1005.07 0.929026 
56   -19.2564     27.8805467 27 15 0.043890 0.041292 0.006571 0.049236 0.066801 1045.91 0.983791 
57       -inf           -inf 10 3 0.110790 0.020554 0.051462 0.079765 0.059837 1021.48 0.980751 
58   -56.8234            inf 4 25 0.015173 0.051632 0.064531 0.080988 0.113355 1002.82 0.93854 
59   -23.6999    33.12347692 30 3 0.028116 0.080733 0.036802 0.061548 0.086527 981.309 0.959664 
60    51.7681    75.46806047 25 19 0.069405 0.088716 0.047466 0.074885 0.127599 941.141 0.951248 
61    -8.3198    -60.0879407 23 4 0.058851 0.075906 0.060683 0.074067 0.065289 1024.75 0.927192 
62   -27.8804   -19.56064248 31 3 0.075360 0.128779 0.058492 0.077866 0.098901 965.712 0.949899 
63   -10.9818    16.89864133 24 3 0.095606 0.047073 0.056014 0.023994 0.045787 1023.24 0.945374 
64   -69.9563   -58.97452552 25 5 0.076531 0.050022 0.028489 0.077643 0.029463 986.948 0.944349 
65   -61.1477    8.808599264 24 12 0.062600 0.043753 0.034261 0.043098 0.025969 1019.72 0.923168 
66   -18.1887    42.95904637 19 20 0.071846 0.102589 0.070084 0.075786 0.027198 986.658 0.951586 
67    -16.802    1.386686873 28 12 0.032148 0.012716 0.050050 0.025153 0.047949 1037.58 0.973288 
68   -16.1652    0.636753893 7 8 0.101169 0.023670 0.058118 0.029912 0.115057 1009.97 0.964564 
69    -16.142  0.02322880318 5 20 0.100853 0.078456 0.047243 0.054582 0.055144 1013.58 0.927249 
70   -58.0725   -41.93051606 7 18 0.049434 0.012171 0.014513 0.068020 0.038519 969.269 0.951866 
71    16.4199    74.49245789 21 18 0.060174 0.023074 0.052633 0.071967 0.052989 925.738 0.978281 
72   -10.8593   -27.27926246 0 8 0.022532 0.105192 0.056539 0.006918 0.101429 970.319 0.976501 
73    -21.692   -10.83265925 25 0 0.133552 0.036650 0.046151 0.008672 0.035675 905.019 0.985115 
74   -73.3571   -51.66514496 21 13 0.064964 0.004109 0.032848 0.081184 0.038667 914.426 0.975492 
75       -inf           -inf 27 5 0.076689 0.043751 0.046078 0.069784 0.011242 939.984 0.968996 
76    45.7287            inf 14 27 0.085099 0.027592 0.061540 0.101239 0.044059 943.454 0.953823 
77   -5.81194   -51.54062591 1 27 0.056160 0.024733 0.075159 0.025582 0.047820 987.678 0.993284 
78    26.2321     32.0440892 5 15 0.086052 0.036700 0.069905 0.048597 0.015619 1044.64 0.967825 
79   -18.1453   -44.37748723 21 0 0.071661 0.015733 0.023555 0.093921 0.055577 968.177 0.98633 
80    -58.986   -40.84066612 21 12 0.016188 0.099511 0.048245 0.023491 0.083770 961.19 0.937597 
81   -21.0535    37.93253815 13 1 0.019964 0.047080 0.108149 0.003985 0.109831 1042.63 0.960442 
82       -inf           -inf 0 31 0.037748 0.025018 0.018343 0.035028 0.064130 1021.27 0.95377 
83    24.9375            inf 19 30 0.031644 0.072280 0.017790 0.068271 0.132742 1028.31 0.955647 
84    7.23514   -17.70237958 6 24 0.114319 0.075102 0.072267 0.077894 0.107576 1035.09 0.944133 
85   -66.6324    -73.8675238 7 28 0.043998 0.087407 0.061186 0.072040 0.086195 1054.5 0.949617 
86    42.0086    108.6409829 12 26 0.038871 0.074238 0.027425 0.012027 0.109680 1034.39 0.92673 
87    29.8297   -12.17887188 24 13 0.058403 0.043671 0.079129 0.043497 0.052428 1064.19 0.930944 
88    28.5024   -1.327313804 22 24 0.075527 0.061188 0.066290 0.069854 0.116405 971.983 0.956431 
89   -23.1101   -51.61252309 12 31 0.017201 0.062835 0.026844 0.049994 0.057349 963.566 0.944969 
90    23.5018    46.61194388 12 20 0.016251 0.028741 0.049749 0.011432 0.016684 1005.37 0.933102 
91    6.85149   -16.65034438 31 4 0.031705 0.076696 0.082427 0.079548 0.049537 974.184 0.936931 
92    37.5863    30.73484409 17 24 0.046159 0.059958 0.054248 0.089417 0.026696 992.314 0.978489 
93   -23.9461   -61.53245975 21 11 0.014585 0.010110 0.036925 0.059542 0.037600 1097.48 0.919809 
94    43.4705    67.41660681 15 20 0.040617 0.038755 0.065204 0.039003 0.058016 972.109 0.974026 
95    4.45489   -39.01558478 6 29 0.078493 0.111230 0.053631 0.070901 0.058016 978.423 0.947338 
96   -25.4947   -29.94962865 22 14 0.068788 0.050981 0.112212 0.042464 0.056454 974.019 0.983925 
97    10.9128    36.40755032 9 22 0.070902 0.029724 0.045187 0.026469 0.105670 992.887 0.954473 
98   -3.69783   -14.61064968 18 21 0.003471 0.052857 0.102139 0.058046 0.065163 984.879 0.996713 
99       -inf           -inf 17 28 0.091326 0.090062 0.087249 0.075993 0.086718 1004.61 0.930667 
100   -86.6297            inf 22 25 0.062062 0.096942 0.027133 0.086058 0.012960 1043.94 0.976471 