########################################################################

# Tobit models using 0's trick

#

# Mroz, T. A. 1987.

#  "The sensitivity of an empirical model of married women's hours to

#     work economic and statistical assumptions," Econometrica 55: 765–799.

#

# dependent variable: working hours (censored at 0)

# independent variables:

#             child, age, edu, hinc

#

########################################################################

model

{

              for(i in 1:n)

              {

ind[i] <- 1 - equals(hour[i], 0) # 0 if hour=0, 1 if hour >0                        

 

z[i]  <- 0

z[i] ~ dpois(psi[i])

 

l.lik[i] <- (1-ind[i])*log( phi( -mu[i]/sigma ) )

                                        +ind[i] *( -lsigma -0.5*pow( (hour[i]-mu[i])/sigma, 2) )

                                                    

psi[i] <- -l.lik[i] + 10000

mu[i]  <- b[1]+b[2]*(child[i]-mean(child[]))/sd(child[])

                                            +b[3]*(age[i]-mean(age[]))/sd(age[])

                                            +b[4]*(edu[i]-mean(edu[]))/sd(edu[])

                                            +b[5]*(hinc[i]-mean(hinc[]))/sd(hinc[])

              } 

 

              b[1]  ~ dnorm( 0, 1.0E-6)

              b[2]  ~ dnorm( 0, 1.0E-6)

              b[3]  ~ dnorm( 0, 1.0E-6)

              b[4]  ~ dnorm( 0, 1.0E-6)

              b[5]  ~ dnorm( 0, 1.0E-6)

                                         

              lsigma  ~ dnorm( 0, 1.0E-6)

              sigma  <- exp( lsigma )

 

              b.cons  <- b[1]-b[2]*mean(child[])/sd(child[])-b[3]*mean(age[])/sd(age[])

                           -b[4]*mean(edu[])/sd(edu[])-b[5]*mean(hinc[])/sd(hinc[])

              b.child <- b[2]/sd(child[])

              b.age   <- b[3]/sd(age[])

              b.edu   <- b[4]/sd(edu[])

              b.hinc  <- b[5]/sd(hinc[])

              }                                                                                                

# Initial values

list( b = c(0, 0, 0, 0, 0), lsigma = 0)

# Data 1 uncensored 428. censored 325.

list(n = 753)

# Data 2

hour[] child[]   age[]   edu[]  hinc[] 

  1610       1      32      12  4.0288 

  1656       0      30      12  8.4416 

  1980       1      35      12  3.5807 

   456       0      34      12  3.5417 

  1568       1      31      14      10 

  2032       0      54      12  6.7106 

  1440       0      37      16  3.4277 

  1020       0      54      12  2.5485 

  1458       0      48      12  4.2206 

  1600       0      39      12  5.7143 

  1969       0      33      12  9.7959 

  1960       0      42      11       8 

   240       1      30      12  5.3004 

   997       0      43      12  4.3413 

  1848       0      43      10   10.87 

  1224       0      35      11  9.1499 

  1400       0      43      12  6.1224 

   640       0      39      12  6.1498 

  2000       0      45      12   6.917 

  1324       0      35      12  4.7103 

  2215       0      42      16   3.131 

  1680       0      30      12       4 

  1600       0      48      13  7.2227 

   800       0      45      12  7.9652 

  1955       1      31      12  4.0884 

   660       0      43      17  14.181 

   525       0      59      12  6.5359 

  1904       0      32      12     8.5 

  1516       1      31      17  6.2762 

   346       0      42      12  5.2083 

  1040       0      50      11  2.7821 

   732       0      59      16  4.9107 

  1880       0      36      13  5.8669 

  1680       0      51      12    7.52 

  2081       0      45      16  7.5449 

   690       0      42      11   5.538 

  4210       0      46      12  6.9231 

  2205       0      46      10       5 

  1952       0      51      14  7.3064 

  1302       0      30      17  11.218 

   112       1      30      12  3.8462 

   893       0      57      12  5.8366 

   583       1      31      16  13.725 

   480       0      48      12  6.3493 

  1900       0      30      12  5.2528 

   576       0      34      12  1.3075 

  2056       0      48      16  2.7998 

  1984       0      45      12  2.6961 

  2640       0      51      12  7.5688 

   240       0      30      12  3.4077 

  1173       0      46      12  6.5401 

  3640       0      58      12  7.2139 

   340       0      37      12  6.2732 

   500       0      52       8  5.8824 

  1599       0      52      10  3.8095 

  1830       0      31      16  6.3776 

  1920       0      55      14  6.0453 

  2052       0      34      17  8.8119 

  2312       0      55      14  8.8765 

   196       0      39      12   3.342 

  2500       0      40      14  3.1836 

  1980       0      43      12  6.9149 

  1840       0      48       8  5.5147 

   320       0      47      12  5.2808 

   419       0      41      12     3.2 

  1880       0      36       8  5.8791 

    72       0      46      17    6.25 

   120       0      34      12  7.4879 

  1885       0      41      12  6.9767 

   240       0      51      12       8 

  1729       0      33      12  4.1319 

  1850       0      52      12  14.476 

  2033       0      58       9  5.7343 

   608       2      34      10  3.3742 

  1153       0      31      12  1.8333 

  2208       0      48      12  5.6025 

   252       0      32      12  4.2975 

   337       0      49      17  3.9161 

    90       2      32      15  4.8787 

  1174       0      58      12  9.2402 

   372       0      50       6  6.6519 

    30       0      60      14  4.6217 

  1800       0      50      12  9.9741 

   282       0      56      14  20.918 

   720       0      51       9    1.94 

  1440       0      54      17  8.8351 

  2100       0      59      13  4.6667 

  1000       0      46       9  5.1894 

   952       0      46      15  18.724 

  1413       1      39      12  10.417 

  2100       0      44      12  8.8458 

   120       2      33      12   7.569 

  3000       1      33      12  4.4507 

  1000       0      48      12  9.3204 

   336       0      31      12  9.1687 

  1216       0      45      12  6.4834  

   988       0      45      12  7.8125 

  2581       0      32      13  11.404 

  2030       0      47      12   5.087 

   413       0      34      13  6.9758 

   782       0      37      12  5.2863 

  1388       0      36      12  11.562 

  1450       1      47      12  8.6061 

  1720       0      48      16  12.805 

   800       0      42      12  6.6946 

   360       0      33      13  8.3333 

  2000       0      46      11  4.1667 

  1176       0      47      12  5.4759 

  2058       0      44      12   5.144 

   900       0      36      12  11.667 

   215       2      31      17  7.2917 

  2000       0      55      14  4.8611 

   757       0      45      16  12.164 

  1264       0      47      17  10.823 

  2064       0      46      12  12.478 

  1280       0      49      11   6.162 

  1715       0      49      12  7.1846 

  2000       0      45      12   7.093 

    12       1      38      17  17.826 

  1528       0      47      10   6.393 

   561       0      54      13  11.719 

  2058       0      41      11  4.7883 

  1823       0      43      12  4.2328 

  1680       1      31      16  7.1066 

  1440       0      47      17  8.3333 

  4950       0      35      12  6.9888 

  1176       0      45      16  9.5913 

  1100       1      33      12  6.0204 

  1516       0      54      16  2.9686 

   900       0      35       8       7 

  1080       1      31      12  7.1627 

   480       0      55      12  4.7364 

   288       0      34      12   8.781 

  1875       0      38      13  9.8315 

   630       0      45      11  6.6496 

   234       0      47      12  9.8485 

  1600       0      39      12  9.8039 

   960       1      36      14  8.2604 

   120       1      33      12  9.0676 

  2025       0      50      12  8.5082 

  1809       0      58      12  7.1565 

  3087       0      49      17  7.0423 

   910       0      41      14  6.5934 

  1840       0      51      12  4.3457 

   784       0      53       9  7.2362 

   400       1      36      12    1.74 

  1000       0      46      12  6.8966 

  1904       0      36      12  5.0455 

  1771       0      53      14   17.25 

  1486       0      40      16  8.3333 

   740       0      42      17   9.375 

  1820       1      33      15   8.279  

  1275       0      43      12     4.5 

   450       1      31      16  4.1995 

  1221       0      47      17  16.064 

  1550       0      54      17  11.277 

   180       1      33      12  0.5494 

  2090       0      43      16   9.548 

  1960       0      46      13  5.3333 

  1440       0      35      12  5.4545 

   794       0      37      11  5.9026 

   993       0      37      16  9.8576 

   160       0      34      14  11.685 

   105       1      43      16  6.2278 

  1200       0      46      12  9.1748 

   450       0      35       9  6.2972 

   996       0      46      17  10.471 

  1052       0      46      14  14.583 

  1128       0      43      12  6.2143 

  1840       0      30      12  6.6288 

  1910       0      41      11  3.8246 

   980       0      54      12  7.8125 

  2317       0      31      12  3.4806 

  1089       0      44      10    5.55 

   800       0      32      12  11.594 

  1230       0      47       5  9.5238 

  1158       0      46      17  7.2805 

  2272       0      37      11       3 

   528       0      51      12    10.4 

  1000       0      49      12   10.05 

   520       0      36      14  8.7866 

   760       0      39      11   10.14 

  1920       0      48      12       5 

  1220       0      38      14  11.888 

   200       0      40      12  13.333 

  2480       1      39      10  8.6806 

  2750       0      37      16  7.2917 

  2014       0      49      13  3.5242 

  1355       0      33      12  6.1818 

    80       0      30      12  2.9755 

  1670       0      54      12  4.5769 

   520       0      39      11  3.1056 

   288       0      43      12  5.8228 

  2014       0      31       9     7.5 

   800       0      33      13  8.3333 

  1984       0      40      12  9.2347 

  1823       0      36      12     5.5 

  1500       0      51      12       6 

  2261       0      44      13  6.6667 

  1728       0      42      16  13.061 

  1960       0      40      12  12.794 

  1578       1      34      16      12 

  1316       0      30      17    6.73 

  1530       0      54      12   6.679 

  2220       0      51      12  4.3945 

  1336       0      44       9  4.6875 

  1008       0      43      12  7.6989 

  1944       0      34      12  4.3287  

  2000       0      45      13  6.7416 

   258       0      39      12  9.3525 

  1785       0      50      12  9.1429 

   480       0      52      12  2.6418 

   772       0      41      12   7.365 

   900       0      59      10  1.4066 

  1428       0      52      12  1.5192 

   210       0      46      16  3.3439 

   239       1      41      12  2.9106 

  1878       0      33      11  4.3834 

   215       0      45      12  4.2098 

  2340       1      36      10  3.7129 

  1960       0      48      12  10.606 

   532       0      47      12  8.6382 

   394       0      45      12  3.8367 

   675       0      37      12  7.1671 

  1515       0      46      16  8.8634 

  1030       0      43      17  7.7369 

  1250       0      42      12  9.1702 

  1158       1      34      17   11.29 

   112       0      52      12  6.3269 

   336       0      37      12  4.7468 

  1984       0      37      12  4.8885 

   716       0      52       8  2.4842 

  1410       1      30      12   8.547 

  1300       0      31      13  6.6379 

  1640       0      38      12  3.1863 

  1202       0      43      12  4.5466 

   489       0      49       8  12.789 

  2076       0      55      12   6.528 

   526       0      38      17  16.403 

  1721       0      52      17       7 

  1327       0      48      12  13.542 

   584       0      32      13   6.713 

  1376       0      32      12  5.4487 

  1040       0      38      12  6.6872 

   548       0      46      12  9.0943 

    15       0      40      12  6.5337 

  1980       0      31       9   4.686 

  1520       0      43      10       8 

  1880       0      51      12   4.053 

  1260       1      30      16  8.1633 

  1092       0      52      13    16.5 

  1587       1      30       8  3.4615 

   156       0      51      16  11.361 

  1939       0      31      13  3.8911 

  1250       0      34      12  9.6059 

   610       0      49      11  2.6687 

   270       1      35      13  8.8865 

   660       1      53      12  6.3636 

  1000       0      32      12  17.728 

  1920       0      38      10  5.2923 

   200       0      54      12  15.686 

  1500       0      47      17  5.2423 

   868       0      45      15  3.0395 

  2318       0      47      16  16.745  

  2000       0      59      10    3.75 

    60       0      32      11   5.814 

  1715       0      45      12  10.833 

   550       0      40      12  9.0029 

  1960       0      47      14    6.54 

    44       1      36      16  5.9783 

  1920       0      56      14  8.6248 

  2540       0      41       8  4.8329 

   156       0      48       7  10.909 

   780       1      36      12  3.8864 

  3120       0      41      12  8.4574 

  2040       0      41      14       5 

  1610       0      36      12  7.4212 

   215       0      37      12  4.4938 

  1120       0      38      12  7.0252 

   846       0      43      14  3.2788 

  3225       0      54      16  1.2861 

  1376       0      38      12  5.5556 

   980       1      30      12  7.4951 

  1838       0      49      12  22.109 

  1494       0      45      13   9.058 

   450       0      51      13  11.538 

  1976       0      34      10  6.0468 

  2012       0      34      12  1.3653 

   561       0      41      12   16.55 

  1715       0      49      12  8.7719 

  1912       0      32      12  9.6354 

  3686       0      32      14  6.6788 

  1080       0      32      17  7.5727 

  1799       0      47      10  7.8125 

  1984       0      39       9  1.9835 

  1839       0      49      12  9.5238 

  1579       0      37      12   6.425 

    96       0      59      16  10.984 

  1920       0      50      12  12.766 

  1688       0      32      17  9.5545 

  1589       0      46      12  5.4967 

   345       0      43      17  9.7618 

  1521       0      37      11       5 

  1490       0      32      16  5.9995 

   989       0      39      11  6.9149 

   600       0      34      13  4.9669 

  2646       0      39      11  2.3896 

  2149       0      45       8  4.5122 

   320       0      50      11  9.2308 

  1600       0      40      12  11.458 

  2419       0      30      10  7.6993 

  2005       0      57      17   4.183 

  1960       0      39      12  8.3333 

  2147       0      53      12  5.1042 

  1207       0      48      17  8.2271 

  2000       0      46      14  5.6634 

  1260       0      47      12   5.222 

    90       0      43      12    7.19 

  1800       0      47      12  13.907 

   573       0      47      12  6.0163  

  1825       0      47      12  7.6782 

    75       0      46      12  5.3112 

  1348       0      34       9  4.1228 

  1880       0      48      10  7.8125 

  1240       0      30      12  4.6053 

   848       0      51      12  3.9746 

   150       0      52      12  15.063 

  2000       0      37      12  4.5064 

  1952       0      32      12  5.8824 

  1456       0      36      17  10.228 

  1740       0      35      12  6.5789 

  1400       0      45      17  8.8145 

  2000       0      56      12  2.9528 

  1750       0      40      10  4.5173 

  1101       1      45      12  5.1282 

  2000       0      32      12  4.1451 

  1877       0      45      12  8.7912 

   160       0      40      12  7.3061 

  1886       0      38      12  6.0069 

  1446       0      49      12  8.3333 

  1500       0      47      16  11.529 

   860       0      52      13  6.6667 

  1848       0      34      13  8.6806 

  1678       0      44      12   3.468 

   160       0      36      16  3.0321 

   108       0      50      17  13.807 

  1738       0      45      12  6.2174 

  1170       0      44      14  4.2308 

    15       0      57      12   11.47 

  2088       0      35      17  12.245 

  2490       0      46      12  9.2105 

   135       2      30      14  6.3309 

  1944       0      42      12  7.8039 

   690       0      34      12  7.6754 

   608       0      45      17  15.582 

    63       1      35      16  9.6436 

   154       0      40      16  12.886 

   420       0      32      12  6.5511 

   651       0      54       9    8.25 

   675       0      38      12  1.3792 

  1663       0      43      12  5.2632 

  1680       0      54      16  10.526 

   180       0      39      14      10 

  1581       0      37      12  12.548 

  1200       0      46      12  4.6207 

   450       0      56      11   1.854 

   547       0      41      12  7.4266 

   300       0      45      16  13.298 

   975       0      44      17  26.578 

  1621       0      50      17  13.959 

   300       0      37      14  13.126 

  1868       0      44      12  7.8883 

  1803       0      32      14   8.656 

  2143       1      34      12  3.4109 

  1080       0      32      10  4.6256 

  1352       0      37      12  4.2105  

   537       0      44      13  4.2708 

   352       0      34      16  3.3424 

   200       1      33      12  8.6359 

  2045       0      43       7  6.7581 

  1253       0      35      16  0.5843 

  1960       0      43      14   5.335 

  2000       0      34      12  4.7101 

  1960       0      36      10  4.9797 

  2000       0      41      12  2.9762 

  1568       0      41      16  12.744 

  1225       0      35      10  3.8462 

   780       1      32      12  7.2233 

   480       0      30      14  8.8077 

  1923       0      43      12  5.3864 

  2000       0      54       6  3.6846 

  2110       0      35      15  1.8437 

  1664       0      50      12  1.3013 

    48       1      34      17  7.2222 

  1791       0      52      14  7.3529 

  1404       0      35      13  11.864 

  1920       0      55       6  5.3327 

  2141       0      35      16     5.5 

  1720       0      49      14  7.2308 

  3533       2      38      15       7 

  2000       0      42      14  7.2137 

   800       0      48       8     4.5 

  3000       0      51      14  8.8632 

   293       0      43      12  9.3333 

  1872       0      43      12  3.8263 

  2058       0      38      12  5.4448 

  1832       0      44      12  8.7634 

   120       1      36      12     8.4 

  1632       0      38      12   5.102 

   778       0      47       8  7.4101 

  1984       0      34      12  5.3527 

   225       1      40      17  19.562 

  1960       0      31      12  5.0783 

   444       0      46      12  0.5128 

   384       0      36      14   5.402 

  1170       1      39      13  7.8646 

  1330       0      36      17  6.0531 

  1350       0      37       8       5 

   480       0      39      12  7.6829 

  1984       1      36      11  6.0386 

  1944       0      49      12   5.851 

    50       1      45      12  6.7416 

   460       2      32      17  8.1081 

   680       0      36      10  5.3061 

  2450       0      40      12  7.2709 

  2144       0      43      13  8.1776 

  1760       0      33      12  7.1006 

   490       0      30      12  6.5844 

     0       0      49      12  7.8529 

     0       2      30      16  11.929 

     0       1      30      12      18 

     0       0      41      12  10.022  

     0       0      45      12  9.3333 

     0       0      43      12   6.085 

     0       0      42      13  5.7054 

     0       0      60      12  9.1185 

     0       0      57      12  7.2072 

     0       0      38      10  5.3146 

     0       0      56      12  8.2828 

     0       0      32      12  7.0578 

     0       0      49       7  6.1659 

     0       0      55      12     2.7 

     0       1      36       9  7.6628 

     0       0      44      12  1.8491 

     0       0      44      10  3.8462 

     0       1      35      14  9.0543 

     0       2      44      14   9.375 

     0       0      45      12  8.3333 

     0       1      34      12  9.3733 

     0       2      30      17  4.4818 

     0       0      39       8  5.3476 

     0       0      36      12  3.2027 

     0       0      38      17  12.801 

     0       0      53      12  6.9103 

     0       0      36      12  4.2293 

     0       1      32      12  4.4479 

     0       0      51       9  8.0523 

     0       0      38      11  7.5057 

     0       2      33      12  7.4747 

     0       0      54      12  6.3218 

     0       0      38       9  2.1354 

     0       2      30      11  4.2092 

     0       2      34      12  0.7595 

     0       0      34       9  4.8077 

     0       0      50      12  12.788 

     0       2      30      17  7.1429 

     0       0      38      12  15.152 

     0       0      54      14  7.9082 

     0       1      30      12  6.9728 

     0       0      55      12  4.9181 

     0       0      51      10  8.3112 

     0       0      44      12  7.1429 

     0       0      53      12  5.6889 

     0       0      42      10  7.0833 

     0       0      38      12  5.6922 

     0       1      38      13    11.4 

     0       1      41      12  1.4152 

     0       0      35       8  2.7507 

     0       1      33      12  6.6182 

     0       0      48      13  9.2725 

     0       0      47      12  12.287 

     0       0      34      12  2.8549 

     0       2      33      13  7.9709 

     0       3      31      13  19.444 

     0       0      58       8  1.4806 

     0       0      49      12  1.8918 

     0       0      55       8  6.9215 

     0       0      44      14  10.692  

     0       0      44       9  9.9474 

     0       0      36      16      19 

     0       0      38      12  7.4882 

     0       0      37      16   9.259 

     0       0      47      12  11.322 

     0       0      47      12   6.447 

     0       1      32      12  6.5778 

     0       1      43      12     6.1 

     0       1      42      11  4.2308 

     0       0      56      12   2.217 

     0       0      38      13  6.5189 

     0       0      52      12  6.4923 

     0       0      50      12  11.364 

     0       0      33      16  9.8765 

     0       0      44      16  7.6576 

     0       0      41      12  8.3017 

     0       0      45      12  14.114 

     0       0      53      14      16 

     0       0      53      14  11.291 

     0       0      42      12  2.3395 

     0       2      32      13  7.7994 

     0       0      56      12       6 

     0       1      37      11  6.6893 

     0       1      40      12  8.5366 

     0       0      54      15   13.83 

     0       0      53       7    1.48 

     0       0      48      12  8.5714 

     0       1      36      12  8.2299 

     0       0      57      12  8.3877 

     0       0      51      12  7.3651 

     0       0      33      13  3.4415 

     0       0      52      12  3.3537 

     0       0      56      10  2.9762 

     0       1      36      12   5.181 

     0       1      36      14  10.293 

     0       0      46      12  13.243 

     0       0      31      10  3.6458 

     0       0      52      11  10.407 

     0       0      46      12  8.2117 

     0       2      35      12  8.6735 

     0       0      59      12  8.3208 

     0       0      36       8   4.287 

     0       1      51       7    10.5 

     0       1      31      16  7.9787 

     0       0      31      14  11.814 

     0       1      32      12  6.1111 

     0       1      35      16  12.578 

     0       0      40      12  11.244 

     0       1      33      10  2.7408 

     0       0      54       7  4.4444 

     0       1      36      12  5.8267 

     0       0      50      10  2.3478 

     0       0      54       8  4.4755 

     0       0      48      11  8.4986 

     0       0      41      15  7.1825 

     0       0      50      12   10.87  

     0       0      46      12  7.1222 

     0       0      42      13 10.6592 

     0       1      31       9  3.2762 

     0       0      53      12  7.7521 

     0       0      51      12  11.945 

     0       0      47      12  4.6154 

     0       0      50      12  13.021 

     0       0      37       6  2.6667 

     0       2      30      12  8.8333 

     0       0      49      12  2.3669 

     0       0      52      12       5 

     0       0      47      12  7.4597 

     0       0      49      12  7.7519 

     0       0      44      12  5.4054 

     0       0      53       8  6.1239 

     0       1      30      12  8.3369 

     0       0      54      12  7.1277 

     0       1      47       7  6.1204 

     0       0      56      15  3.1111 

     0       0      49      12  12.161 

     0       0      48       6  5.3846 

     0       0      49      12  3.9838 

     0       0      56      12  7.4697 

     0       0      46      12  11.411 

     0       0      45      12  12.936 

     0       0      32      12  5.6122 

     0       1      43      12  6.8099 

     0       1      34      12  8.5066 

     0       1      30      12   5.508 

     0       2      38      17  7.2285 

     0       1      33      16  9.4545 

     0       0      52      12  8.9835 

     0       0      43      11  5.9846 

     0       1      33      12  3.9629 

     0       0      45      10  7.0575 

     0       2      36      10  7.1565 

     0       1      34      12  4.0145 

     0       0      37      14  3.8829 

     0       0      46      10  7.9787 

     0       0      47      12  5.0366 

     0       2      31      16  7.0769 

     0       0      57       5  3.6184 

     0       1      30      12  13.466 

     0       0      30      12  14.531 

     0       0      44      12     6.8 

     0       0      53      13  6.6667 

     0       0      51       8  5.9172 

     0       1      39      12   10.24 

     0       0      52       8  8.1566 

     0       0      46       8  6.6775 

     0       0      47      12  0.9615 

     0       0      52       8  16.907 

     0       0      45      12   13.57 

     0       0      60      11  8.9329 

     0       0      41      13  4.7372 

     0       0      39       8  8.5106  

     0       0      49      12  5.2511 

     0       1      32      15  21.786 

     0       1      33      12  7.9125 

     0       0      36      10  7.1083 

     0       3      37      13  5.0434 

     0       1      30      12  2.2442 

     0       1      44      11  8.5069 

     0       0      48      12  7.5758 

     0       0      40      11  5.8824 

     0       0      47      13  40.441 

     0       0      36      12  6.2339 

     0       0      40      11  9.2971 

     0       0      46      12  0.7038 

     0       0      52      12  2.5425 

     0       0      44      12  5.0833 

     0       0      45      10  7.6843 

     0       2      30       7  4.0625 

     0       1      40      12  6.5106 

     0       0      43      12   2.951 

     0       0      49      12  5.6738 

     0       1      46      12  12.766 

     0       0      52      11       6 

     0       1      31      12  9.2955 

     0       1      42      10  5.1324 

     0       0      33      11  9.0659 

     0       0      57      16  11.931 

     0       0      49      10  6.7164 

     0       0      45      14  4.5851 

     0       0      56      11  9.2567 

     0       1      41      12  8.5656 

     0       0      56       5  6.7407 

     0       0      48      10  9.8788 

     0       0      52      16   9.375 

     0       0      51      12   6.277 

     0       0      35      11  6.4241 

     0       0      45      12  8.3718 

     0       0      54      12   12.22 

     0       0      54      12  3.5612 

     0       1      31      12   8.867 

     0       0      53       6  1.8515 

     0       2      35      14  8.0769 

     0       1      36      12  8.9796 

     0       0      59      12  6.6618 

     0       0      54      16  11.107 

     0       1      37      12  14.077 

     0       0      44      12  0.4121 

     0       1      34      17   23.81 

     0       0      49      12  9.5745 

     0       0      49      12  1.9286 

     0       0      60       9  2.5783 

     0       0      51      12  7.2274 

     0       1      30      12     8.8 

     0       0      47      12    7.75 

     0       0      36      12  14.688 

     0       1      35      12  10.487 

     0       0      58      12  4.0253  

     0       1      41      14  11.765 

     0       0      51      10  9.5238 

     0       0      47      12  6.9079 

     0       1      45       9   5.151 

     0       0      60      14  20.496 

     0       1      30      16  13.555 

     0       0      55      11  14.796 

     0       1      32      12  1.6905 

     0       0      36      12  7.2296 

     0       0      55      12  4.1152 

     0       0      47      12   2.999 

     0       0      47      12  8.0492 

     0       0      37      11   8.577 

     0       0      50      12  9.3235 

     0       0      30      12  4.7059 

     0       0      48      17  9.1346 

     0       0      43      10   7.767 

     0       1      48      11  4.1269 

     0       1      41      14      10 

     0       0      50      12  4.0727 

     0       0      58       8  7.8529 

     0       0      38      13  7.8168 

     0       0      37      12  9.1667 

     0       0      50      16  29.644 

     0       0      42       8  4.0761 

     0       1      37       9   8.362 

     0       0      41      16  12.062 

     0       0      31      12  2.3251 

     0       0      51      12   10.01 

     0       1      36      12  5.2734 

     0       0      54      15  10.417 

     0       0      49      12     7.5 

     0       1      48       9  6.3521 

     0       0      42       9  3.8071 

     0       1      41      12  9.9099 

     0       0      55      16  2.6701 

     0       0      42       9  4.3956 

     0       0      32      15   7.971 

     0       0      43      12  15.429 

     0       1      33      12  1.9854 

     0       0      48      15  7.2562 

     0       0      43      12  7.0352 

     0       1      47      17  8.1347 

     0       0      54      12  11.981 

     0       0      51      12  5.9922 

     0       0      51      10  4.1458 

     0       1      43      13  8.5102 

     0       0      53      12  2.9631 

     0       1      34      11  8.8149 

     0       1      31       8  2.4955 

     0       0      56      12  4.8843 

     0       0      42      16    10.2 

     0       0      32      12      10 

     0       1      35      12  5.3571 

     0       1      30      12  9.7959 

     0       0      51      10  1.3741  

     0       0      47      12   3.463 

     0       0      54      12  5.0408 

     0       3      31      15  16.667 

     0       0      47      10  4.3294 

     0       0      47      14  10.193 

     0       0      40      12  7.5833 

     0       0      48       8  4.1204 

     0       0      34       8  4.8077 

     0       0      38      12  3.1531 

     0       1      32      12      19 

     0       0      48      16  6.0535 

     0       0      41      12  3.7875 

     0       0      49       5  2.9287 

     0       0      59       8  9.2972 

     0       0      58      13  9.3827 

     0       0      41      12  26.531 

     0       0      45      12  40.509 

     0       1      30      14   15.16 

     0       0      41      12  18.964 

     0       2      30      12  5.4152 

     0       0      53      12  1.9988 

     0       0      31      12       8 

     0       0      43      14  11.153 

     0       1      31      12  5.9524 

     0       0      51      12  2.1194 

     0       0      43       9  3.8384 

     0       1      31      14  17.361 

     0       0      48      11  8.0092 

     0       1      31      12       3 

     0       0      44      12  10.589 

     0       0      48      11  10.913 

     0       0      53      12  5.6033 

     0       0      42      10  7.9365 

     0       2      39      12  2.9476 

     0       1      32      10  2.0562 

     0       0      36      12  1.3013 

     0       0      40      13  9.2715 

     0       2      31      12  4.8638 

     0       0      43      12  1.0898 

     0       0      60      12   12.44 

     0       0      39       9  6.0897

END