Purpose
To use EM in the context of Exercise 7.28

> x <- c(3540, 3560, 3739, 2784, 2571, 2729, 3952, 993, 1908, 948,
+     1172, 1047, 3138, 5485, 5554, 2943, 4969, 4828)
> y <- c(3, 4, 1, 1, 3, 1, 2, 0, 2, 0, 1, 3, 5, 4, 6, 2, 5, 4)
> beta.complete <- sum(y)/sum(x)
> tj.complete <- (x + y)/(beta.complete + 1)
> beta.incomplete <- sum(y[-1])/sum(x[-1])
> tj.incomplete <- c(y[1]/beta.incomplete, ((x + y)/(beta.incomplete +
+     1))[-1])

Parameters under complete data

> print(beta.complete)
[1] 0.0008413892
> print(tj.complete)
 [1] 3540.0215 3561.0038 3736.8558 2782.6587 2571.8361 2727.7049 3950.6759
 [8]  992.1652 1908.3943  947.2030 1172.0139 1049.1173 3140.3577 5484.3855
[15] 5555.3258 2942.5242 4969.8184 4827.9378

Parameters under incomplete data

> print(beta.incomplete)
[1] 0.0008409786
> print(tj.incomplete)
 [1] 3567.2727 3561.0053 3736.8574 2782.6598 2571.8371 2727.7061 3950.6776
 [8]  992.1656 1908.3951  947.2034 1172.0144 1049.1177 3140.3590 5484.3877
[15] 5555.3281 2942.5254 4969.8205 4827.9398

MLE logic

> x <- c(3540, 3560, 3739, 2784, 2571, 2729, 3952, 993, 1908, 948,
+     1172, 1047, 3138, 5485, 5554, 2943, 4969, 4828)
> y <- c(3, 4, 1, 1, 3, 1, 2, 0, 2, 0, 1, 3, 5, 4, 6, 2, 5, 4)
> x.incomplete <- x[-1]
> beta.0 <- 0.001
> tau.0 <- rep(2000, 18)
> param <- list(beta = beta.0, tau = tau.0)
> em <- function(y, x, param) {
+     temp.beta <- param$beta
+     temp.tau <- param$tau
+     new.beta <- sum(y)/(temp.tau[1] + sum(x))
+     new.tau1 <- (temp.tau[1] + y[1])/(new.beta + 1)
+     new.tau2 <- (y[-1] + x)/(new.beta + 1)
+     new.tau <- c(new.tau1, new.tau2)
+     param.new <- list(beta = new.beta, tau = new.tau)
+     return(param.new)
+ }
> param.start <- param
> counter <- 1
> param.temp <- em(y, x.incomplete, param)
> while (abs(param.temp$beta - param.start$beta) > 1e-05) {
+     param.start <- param.temp
+     param.temp <- em(y, x.incomplete, param.temp)
+     counter <- counter + 1
+     print(param.temp)
+ }
$beta
[1] 0.0008652228


$tau
 [1] 2002.5358 3560.9190 3736.7669 2782.5924 2571.7748 2727.6400 3950.5819
 [8]  992.1416 1908.3489  947.1805 1171.9860 1049.0923 3140.2830 5484.2549
[15] 5555.1935 2942.4541 4969.7001 4827.8229
> k <- 1
> for (k in 1:18) {
+     while (abs(param.temp$tau[k] - param.start$tau[k]) > 1) {
+         param.start <- param.temp
+         param.temp <- em(y, x.incomplete, param.temp)
+         print(param.temp)
+     }
+ }
$beta
[1] 0.0008652026


$tau
 [1] 2003.8021 3560.9191 3736.7669 2782.5925 2571.7749 2727.6400 3950.5819
 [8]  992.1416 1908.3489  947.1805 1171.9860 1049.0923 3140.2830 5484.2550
[15] 5555.1936 2942.4542 4969.7002 4827.8230


$beta
[1] 0.0008651824


$tau
 [1] 2005.0673 3560.9192 3736.7670 2782.5925 2571.7749 2727.6401 3950.5820
 [8]  992.1416 1908.3489  947.1805 1171.9860 1049.0923 3140.2831 5484.2551
[15] 5555.1937 2942.4542 4969.7003 4827.8231


$beta
[1] 0.0008651623


$tau
 [1] 2006.3315 3560.9192 3736.7671 2782.5926 2571.7750 2727.6401 3950.5821
 [8]  992.1416 1908.3490  947.1805 1171.9860 1049.0924 3140.2831 5484.2552
[15] 5555.1939 2942.4543 4969.7004 4827.8231


$beta
[1] 0.0008651422


$tau
 [1] 2007.5947 3560.9193 3736.7672 2782.5927 2571.7750 2727.6402 3950.5822
 [8]  992.1417 1908.3490  947.1806 1171.9861 1049.0924 3140.2832 5484.2553
[15] 5555.1940 2942.4544 4969.7005 4827.8232


$beta
[1] 0.000865122


$tau
 [1] 2008.8568 3560.9194 3736.7672 2782.5927 2571.7751 2727.6403 3950.5823
 [8]  992.1417 1908.3490  947.1806 1171.9861 1049.0924 3140.2833 5484.2554
[15] 5555.1941 2942.4544 4969.7006 4827.8233


$beta
[1] 0.000865102


$tau
 [1] 2010.1178 3560.9194 3736.7673 2782.5928 2571.7752 2727.6403 3950.5823
 [8]  992.1417 1908.3491  947.1806 1171.9861 1049.0924 3140.2833 5484.2556
[15] 5555.1942 2942.4545 4969.7007 4827.8234


$beta
[1] 0.0008650819


$tau
 [1] 2011.3778 3560.9195 3736.7674 2782.5928 2571.7752 2727.6404 3950.5824
 [8]  992.1417 1908.3491  947.1806 1171.9861 1049.0924 3140.2834 5484.2557
[15] 5555.1943 2942.4545 4969.7008 4827.8235


$beta
[1] 0.0008650618


$tau
 [1] 2012.6368 3560.9196 3736.7675 2782.5929 2571.7753 2727.6404 3950.5825
 [8]  992.1417 1908.3492  947.1806 1171.9862 1049.0925 3140.2835 5484.2558
[15] 5555.1944 2942.4546 4969.7009 4827.8236


$beta
[1] 0.0008650418


$tau
 [1] 2013.8946 3560.9197 3736.7675 2782.5929 2571.7753 2727.6405 3950.5826
 [8]  992.1418 1908.3492  947.1806 1171.9862 1049.0925 3140.2835 5484.2559
[15] 5555.1945 2942.4547 4969.7010 4827.8237


$beta
[1] 0.0008650217


$tau
 [1] 2015.1515 3560.9197 3736.7676 2782.5930 2571.7754 2727.6405 3950.5827
 [8]  992.1418 1908.3492  947.1807 1171.9862 1049.0925 3140.2836 5484.2560
[15] 5555.1946 2942.4547 4969.7011 4827.8238


$beta
[1] 0.0008650017


$tau
 [1] 2016.4073 3560.9198 3736.7677 2782.5931 2571.7754 2727.6406 3950.5827
 [8]  992.1418 1908.3493  947.1807 1171.9862 1049.0925 3140.2836 5484.2561
[15] 5555.1947 2942.4548 4969.7012 4827.8239


$beta
[1] 0.0008649817


$tau
 [1] 2017.6621 3560.9199 3736.7678 2782.5931 2571.7755 2727.6406 3950.5828
 [8]  992.1418 1908.3493  947.1807 1171.9863 1049.0926 3140.2837 5484.2562
[15] 5555.1949 2942.4548 4969.7013 4827.8240


$beta
[1] 0.0008649618


$tau
 [1] 2018.9158 3560.9199 3736.7678 2782.5932 2571.7755 2727.6407 3950.5829
 [8]  992.1418 1908.3494  947.1807 1171.9863 1049.0926 3140.2838 5484.2563
[15] 5555.1950 2942.4549 4969.7014 4827.8241


$beta
[1] 0.0008649418


$tau
 [1] 2020.1684 3560.9200 3736.7679 2782.5932 2571.7756 2727.6407 3950.5830
 [8]  992.1419 1908.3494  947.1807 1171.9863 1049.0926 3140.2838 5484.2564
[15] 5555.1951 2942.4549 4969.7015 4827.8242


$beta
[1] 0.0008649219


$tau
 [1] 2021.4201 3560.9201 3736.7680 2782.5933 2571.7756 2727.6408 3950.5831
 [8]  992.1419 1908.3494  947.1808 1171.9863 1049.0926 3140.2839 5484.2565
[15] 5555.1952 2942.4550 4969.7016 4827.8243


$beta
[1] 0.0008649019


$tau
 [1] 2022.6707 3560.9202 3736.7681 2782.5933 2571.7757 2727.6409 3950.5831
 [8]  992.1419 1908.3495  947.1808 1171.9863 1049.0926 3140.2840 5484.2567
[15] 5555.1953 2942.4551 4969.7017 4827.8244


$beta
[1] 0.000864882


$tau
 [1] 2023.9202 3560.9202 3736.7681 2782.5934 2571.7757 2727.6409 3950.5832
 [8]  992.1419 1908.3495  947.1808 1171.9864 1049.0927 3140.2840 5484.2568
[15] 5555.1954 2942.4551 4969.7018 4827.8245


$beta
[1] 0.0008648622


$tau
 [1] 2025.1687 3560.9203 3736.7682 2782.5934 2571.7758 2727.6410 3950.5833
 [8]  992.1419 1908.3495  947.1808 1171.9864 1049.0927 3140.2841 5484.2569
[15] 5555.1955 2942.4552 4969.7019 4827.8246


$beta
[1] 0.0008648423


$tau
 [1] 2026.4162 3560.9204 3736.7683 2782.5935 2571.7758 2727.6410 3950.5834
 [8]  992.1420 1908.3496  947.1808 1171.9864 1049.0927 3140.2841 5484.2570
[15] 5555.1956 2942.4552 4969.7020 4827.8247


$beta
[1] 0.0008648224


$tau
 [1] 2027.6626 3560.9204 3736.7684 2782.5936 2571.7759 2727.6411 3950.5834
 [8]  992.1420 1908.3496  947.1809 1171.9864 1049.0927 3140.2842 5484.2571
[15] 5555.1957 2942.4553 4969.7021 4827.8248


$beta
[1] 0.0008648026


$tau
 [1] 2028.9080 3560.9205 3736.7684 2782.5936 2571.7759 2727.6411 3950.5835
 [8]  992.1420 1908.3497  947.1809 1171.9865 1049.0927 3140.2843 5484.2572
[15] 5555.1959 2942.4554 4969.7022 4827.8249


$beta
[1] 0.0008647828


$tau
 [1] 2030.1524 3560.9206 3736.7685 2782.5937 2571.7760 2727.6412 3950.5836
 [8]  992.1420 1908.3497  947.1809 1171.9865 1049.0928 3140.2843 5484.2573
[15] 5555.1960 2942.4554 4969.7023 4827.8250


$beta
[1] 0.000864763


$tau
 [1] 2031.3957 3560.9206 3736.7686 2782.5937 2571.7760 2727.6412 3950.5837
 [8]  992.1420 1908.3497  947.1809 1171.9865 1049.0928 3140.2844 5484.2574
[15] 5555.1961 2942.4555 4969.7024 4827.8251


$beta
[1] 0.0008647432


$tau
 [1] 2032.6380 3560.9207 3736.7687 2782.5938 2571.7761 2727.6413 3950.5838
 [8]  992.1421 1908.3498  947.1809 1171.9865 1049.0928 3140.2845 5484.2575
[15] 5555.1962 2942.4555 4969.7025 4827.8252


$beta
[1] 0.0008647234


$tau
 [1] 2033.8792 3560.9208 3736.7687 2782.5938 2571.7761 2727.6413 3950.5838
 [8]  992.1421 1908.3498  947.1810 1171.9866 1049.0928 3140.2845 5484.2576
[15] 5555.1963 2942.4556 4969.7026 4827.8253


$beta
[1] 0.0008647037


$tau
 [1] 2035.1195 3560.9209 3736.7688 2782.5939 2571.7762 2727.6414 3950.5839
 [8]  992.1421 1908.3498  947.1810 1171.9866 1049.0928 3140.2846 5484.2577
[15] 5555.1964 2942.4556 4969.7027 4827.8254


$beta
[1] 0.000864684


$tau
 [1] 2036.3587 3560.9209 3736.7689 2782.5939 2571.7762 2727.6415 3950.5840
 [8]  992.1421 1908.3499  947.1810 1171.9866 1049.0929 3140.2846 5484.2579
[15] 5555.1965 2942.4557 4969.7028 4827.8255


$beta
[1] 0.0008646642


$tau
 [1] 2037.5968 3560.9210 3736.7689 2782.5940 2571.7763 2727.6415 3950.5841
 [8]  992.1421 1908.3499  947.1810 1171.9866 1049.0929 3140.2847 5484.2580
[15] 5555.1966 2942.4558 4969.7029 4827.8256


$beta
[1] 0.0008646446


$tau
 [1] 2038.8340 3560.9211 3736.7690 2782.5940 2571.7763 2727.6416 3950.5841
 [8]  992.1421 1908.3500  947.1810 1171.9866 1049.0929 3140.2848 5484.2581
[15] 5555.1967 2942.4558 4969.7030 4827.8256


$beta
[1] 0.0008646249


$tau
 [1] 2040.0701 3560.9211 3736.7691 2782.5941 2571.7764 2727.6416 3950.5842
 [8]  992.1422 1908.3500  947.1810 1171.9867 1049.0929 3140.2848 5484.2582
[15] 5555.1968 2942.4559 4969.7031 4827.8257


$beta
[1] 0.0008646052


$tau
 [1] 2041.3051 3560.9212 3736.7692 2782.5942 2571.7764 2727.6417 3950.5843
 [8]  992.1422 1908.3500  947.1811 1171.9867 1049.0929 3140.2849 5484.2583
[15] 5555.1969 2942.4559 4969.7032 4827.8258


$beta
[1] 0.0008645856


$tau
 [1] 2042.5392 3560.9213 3736.7692 2782.5942 2571.7765 2727.6417 3950.5844
 [8]  992.1422 1908.3501  947.1811 1171.9867 1049.0930 3140.2850 5484.2584
[15] 5555.1971 2942.4560 4969.7033 4827.8259


$beta
[1] 0.0008645659


$tau
 [1] 2043.7722 3560.9213 3736.7693 2782.5943 2571.7765 2727.6418 3950.5845
 [8]  992.1422 1908.3501  947.1811 1171.9867 1049.0930 3140.2850 5484.2585
[15] 5555.1972 2942.4561 4969.7034 4827.8260


$beta
[1] 0.0008645463


$tau
 [1] 2045.0042 3560.9214 3736.7694 2782.5943 2571.7766 2727.6418 3950.5845
 [8]  992.1422 1908.3501  947.1811 1171.9868 1049.0930 3140.2851 5484.2586
[15] 5555.1973 2942.4561 4969.7035 4827.8261


$beta
[1] 0.0008645267


$tau
 [1] 2046.2352 3560.9215 3736.7695 2782.5944 2571.7766 2727.6419 3950.5846
 [8]  992.1423 1908.3502  947.1811 1171.9868 1049.0930 3140.2851 5484.2587
[15] 5555.1974 2942.4562 4969.7036 4827.8262


$beta
[1] 0.0008645072


$tau
 [1] 2047.4651 3560.9216 3736.7695 2782.5944 2571.7767 2727.6419 3950.5847
 [8]  992.1423 1908.3502  947.1812 1171.9868 1049.0931 3140.2852 5484.2588
[15] 5555.1975 2942.4562 4969.7037 4827.8263


$beta
[1] 0.0008644876


$tau
 [1] 2048.6941 3560.9216 3736.7696 2782.5945 2571.7767 2727.6420 3950.5848
 [8]  992.1423 1908.3503  947.1812 1171.9868 1049.0931 3140.2853 5484.2589
[15] 5555.1976 2942.4563 4969.7038 4827.8264


$beta
[1] 0.0008644681


$tau
 [1] 2049.9220 3560.9217 3736.7697 2782.5945 2571.7768 2727.6420 3950.5848
 [8]  992.1423 1908.3503  947.1812 1171.9869 1049.0931 3140.2853 5484.2590
[15] 5555.1977 2942.4563 4969.7038 4827.8265


$beta
[1] 0.0008644485


$tau
 [1] 2051.1489 3560.9218 3736.7698 2782.5946 2571.7768 2727.6421 3950.5849
 [8]  992.1423 1908.3503  947.1812 1171.9869 1049.0931 3140.2854 5484.2591
[15] 5555.1978 2942.4564 4969.7039 4827.8266


$beta
[1] 0.000864429


$tau
 [1] 2052.3747 3560.9218 3736.7698 2782.5946 2571.7769 2727.6421 3950.5850
 [8]  992.1424 1908.3504  947.1812 1171.9869 1049.0931 3140.2854 5484.2592
[15] 5555.1979 2942.4565 4969.7040 4827.8267


$beta
[1] 0.0008644096


$tau
 [1] 2053.5996 3560.9219 3736.7699 2782.5947 2571.7769 2727.6422 3950.5851
 [8]  992.1424 1908.3504  947.1812 1171.9869 1049.0932 3140.2855 5484.2594
[15] 5555.1980 2942.4565 4969.7041 4827.8268


$beta
[1] 0.0008643901


$tau
 [1] 2054.8234 3560.9220 3736.7700 2782.5948 2571.7770 2727.6423 3950.5852
 [8]  992.1424 1908.3504  947.1813 1171.9869 1049.0932 3140.2856 5484.2595
[15] 5555.1981 2942.4566 4969.7042 4827.8269


$beta
[1] 0.0008643706


$tau
 [1] 2056.0462 3560.9220 3736.7700 2782.5948 2571.7770 2727.6423 3950.5852
 [8]  992.1424 1908.3505  947.1813 1171.9870 1049.0932 3140.2856 5484.2596
[15] 5555.1982 2942.4566 4969.7043 4827.8270


$beta
[1] 0.0008643512


$tau
 [1] 2057.2680 3560.9221 3736.7701 2782.5949 2571.7771 2727.6424 3950.5853
 [8]  992.1424 1908.3505  947.1813 1171.9870 1049.0932 3140.2857 5484.2597
[15] 5555.1984 2942.4567 4969.7044 4827.8271


$beta
[1] 0.0008643318


$tau
 [1] 2058.4888 3560.9222 3736.7702 2782.5949 2571.7771 2727.6424 3950.5854
 [8]  992.1425 1908.3506  947.1813 1171.9870 1049.0932 3140.2858 5484.2598
[15] 5555.1985 2942.4567 4969.7045 4827.8272


$beta
[1] 0.0008643124


$tau
 [1] 2059.7086 3560.9223 3736.7703 2782.5950 2571.7772 2727.6425 3950.5855
 [8]  992.1425 1908.3506  947.1813 1171.9870 1049.0933 3140.2858 5484.2599
[15] 5555.1986 2942.4568 4969.7046 4827.8272


$beta
[1] 0.000864293


$tau
 [1] 2060.9273 3560.9223 3736.7703 2782.5950 2571.7772 2727.6425 3950.5855
 [8]  992.1425 1908.3506  947.1814 1171.9871 1049.0933 3140.2859 5484.2600
[15] 5555.1987 2942.4569 4969.7047 4827.8273


$beta
[1] 0.0008642736


$tau
 [1] 2062.1451 3560.9224 3736.7704 2782.5951 2571.7773 2727.6426 3950.5856
 [8]  992.1425 1908.3507  947.1814 1171.9871 1049.0933 3140.2859 5484.2601
[15] 5555.1988 2942.4569 4969.7048 4827.8274


$beta
[1] 0.0008642542


$tau
 [1] 2063.3618 3560.9225 3736.7705 2782.5951 2571.7773 2727.6426 3950.5857
 [8]  992.1425 1908.3507  947.1814 1171.9871 1049.0933 3140.2860 5484.2602
[15] 5555.1989 2942.4570 4969.7049 4827.8275


$beta
[1] 0.0008642349


$tau
 [1] 2064.5775 3560.9225 3736.7706 2782.5952 2571.7774 2727.6427 3950.5858
 [8]  992.1426 1908.3507  947.1814 1171.9871 1049.0933 3140.2861 5484.2603
[15] 5555.1990 2942.4570 4969.7050 4827.8276


$beta
[1] 0.0008642156


$tau
 [1] 2065.7922 3560.9226 3736.7706 2782.5952 2571.7774 2727.6427 3950.5858
 [8]  992.1426 1908.3508  947.1814 1171.9872 1049.0934 3140.2861 5484.2604
[15] 5555.1991 2942.4571 4969.7051 4827.8277


$beta
[1] 0.0008641963


$tau
 [1] 2067.0059 3560.9227 3736.7707 2782.5953 2571.7775 2727.6428 3950.5859
 [8]  992.1426 1908.3508  947.1814 1171.9872 1049.0934 3140.2862 5484.2605
[15] 5555.1992 2942.4571 4969.7052 4827.8278


$beta
[1] 0.000864177


$tau
 [1] 2068.2186 3560.9227 3736.7708 2782.5953 2571.7775 2727.6428 3950.5860
 [8]  992.1426 1908.3508  947.1815 1171.9872 1049.0934 3140.2862 5484.2606
[15] 5555.1993 2942.4572 4969.7053 4827.8279


$beta
[1] 0.0008641577


$tau
 [1] 2069.4303 3560.9228 3736.7708 2782.5954 2571.7776 2727.6429 3950.5861
 [8]  992.1426 1908.3509  947.1815 1171.9872 1049.0934 3140.2863 5484.2607
[15] 5555.1994 2942.4573 4969.7054 4827.8280


$beta
[1] 0.0008641385


$tau
 [1] 2070.6410 3560.9229 3736.7709 2782.5955 2571.7776 2727.6429 3950.5861
 [8]  992.1427 1908.3509  947.1815 1171.9872 1049.0934 3140.2864 5484.2608
[15] 5555.1995 2942.4573 4969.7055 4827.8281


$beta
[1] 0.0008641193


$tau
 [1] 2071.8507 3560.9229 3736.7710 2782.5955 2571.7777 2727.6430 3950.5862
 [8]  992.1427 1908.3510  947.1815 1171.9873 1049.0935 3140.2864 5484.2609
[15] 5555.1996 2942.4574 4969.7056 4827.8282


$beta
[1] 0.0008641


$tau
 [1] 2073.0593 3560.9230 3736.7711 2782.5956 2571.7777 2727.6430 3950.5863
 [8]  992.1427 1908.3510  947.1815 1171.9873 1049.0935 3140.2865 5484.2610
[15] 5555.1998 2942.4574 4969.7057 4827.8283


$beta
[1] 0.0008640808


$tau
 [1] 2074.2670 3560.9231 3736.7711 2782.5956 2571.7778 2727.6431 3950.5864
 [8]  992.1427 1908.3510  947.1816 1171.9873 1049.0935 3140.2865 5484.2612
[15] 5555.1999 2942.4575 4969.7058 4827.8284


$beta
[1] 0.0008640616


$tau
 [1] 2075.4737 3560.9231 3736.7712 2782.5957 2571.7778 2727.6431 3950.5864
 [8]  992.1427 1908.3511  947.1816 1171.9873 1049.0935 3140.2866 5484.2613
[15] 5555.2000 2942.4575 4969.7059 4827.8285


$beta
[1] 0.0008640425


$tau
 [1] 2076.6793 3560.9232 3736.7713 2782.5957 2571.7779 2727.6432 3950.5865
 [8]  992.1427 1908.3511  947.1816 1171.9874 1049.0935 3140.2867 5484.2614
[15] 5555.2001 2942.4576 4969.7060 4827.8286


$beta
[1] 0.0008640233


$tau
 [1] 2077.8840 3560.9233 3736.7713 2782.5958 2571.7779 2727.6433 3950.5866
 [8]  992.1428 1908.3511  947.1816 1171.9874 1049.0936 3140.2867 5484.2615
[15] 5555.2002 2942.4576 4969.7061 4827.8286


$beta
[1] 0.0008640042


$tau
 [1] 2079.0876 3560.9233 3736.7714 2782.5958 2571.7780 2727.6433 3950.5867
 [8]  992.1428 1908.3512  947.1816 1171.9874 1049.0936 3140.2868 5484.2616
[15] 5555.2003 2942.4577 4969.7062 4827.8287


$beta
[1] 0.0008639851


$tau
 [1] 2080.2903 3560.9234 3736.7715 2782.5959 2571.7780 2727.6434 3950.5868
 [8]  992.1428 1908.3512  947.1816 1171.9874 1049.0936 3140.2868 5484.2617
[15] 5555.2004 2942.4578 4969.7062 4827.8288


$beta
[1] 0.000863966


$tau
 [1] 2081.4920 3560.9235 3736.7716 2782.5959 2571.7781 2727.6434 3950.5868
 [8]  992.1428 1908.3512  947.1817 1171.9874 1049.0936 3140.2869 5484.2618
[15] 5555.2005 2942.4578 4969.7063 4827.8289


$beta
[1] 0.0008639469


$tau
 [1] 2082.6926 3560.9236 3736.7716 2782.5960 2571.7781 2727.6435 3950.5869
 [8]  992.1428 1908.3513  947.1817 1171.9875 1049.0936 3140.2870 5484.2619
[15] 5555.2006 2942.4579 4969.7064 4827.8290


$beta
[1] 0.0008639278


$tau
 [1] 2083.8923 3560.9236 3736.7717 2782.5960 2571.7782 2727.6435 3950.5870
 [8]  992.1429 1908.3513  947.1817 1171.9875 1049.0937 3140.2870 5484.2620
[15] 5555.2007 2942.4579 4969.7065 4827.8291


$beta
[1] 0.0008639088


$tau
 [1] 2085.0910 3560.9237 3736.7718 2782.5961 2571.7782 2727.6436 3950.5871
 [8]  992.1429 1908.3514  947.1817 1171.9875 1049.0937 3140.2871 5484.2621
[15] 5555.2008 2942.4580 4969.7066 4827.8292


$beta
[1] 0.0008638897


$tau
 [1] 2086.2887 3560.9238 3736.7718 2782.5961 2571.7783 2727.6436 3950.5871
 [8]  992.1429 1908.3514  947.1817 1171.9875 1049.0937 3140.2871 5484.2622
[15] 5555.2009 2942.4580 4969.7067 4827.8293


$beta
[1] 0.0008638707


$tau
 [1] 2087.4853 3560.9238 3736.7719 2782.5962 2571.7783 2727.6437 3950.5872
 [8]  992.1429 1908.3514  947.1818 1171.9876 1049.0937 3140.2872 5484.2623
[15] 5555.2010 2942.4581 4969.7068 4827.8294


$beta
[1] 0.0008638517


$tau
 [1] 2088.6810 3560.9239 3736.7720 2782.5962 2571.7784 2727.6437 3950.5873
 [8]  992.1429 1908.3515  947.1818 1171.9876 1049.0937 3140.2873 5484.2624
[15] 5555.2011 2942.4582 4969.7069 4827.8295


$beta
[1] 0.0008638327


$tau
 [1] 2089.8757 3560.9240 3736.7721 2782.5963 2571.7784 2727.6438 3950.5874
 [8]  992.1430 1908.3515  947.1818 1171.9876 1049.0938 3140.2873 5484.2625
[15] 5555.2012 2942.4582 4969.7070 4827.8296


$beta
[1] 0.0008638138


$tau
 [1] 2091.0694 3560.9240 3736.7721 2782.5964 2571.7785 2727.6438 3950.5874
 [8]  992.1430 1908.3515  947.1818 1171.9876 1049.0938 3140.2874 5484.2626
[15] 5555.2013 2942.4583 4969.7071 4827.8297


$beta
[1] 0.0008637948


$tau
 [1] 2092.2621 3560.9241 3736.7722 2782.5964 2571.7785 2727.6439 3950.5875
 [8]  992.1430 1908.3516  947.1818 1171.9876 1049.0938 3140.2874 5484.2627
[15] 5555.2014 2942.4583 4969.7072 4827.8297


$beta
[1] 0.0008637759


$tau
 [1] 2093.4539 3560.9242 3736.7723 2782.5965 2571.7786 2727.6439 3950.5876
 [8]  992.1430 1908.3516  947.1818 1171.9877 1049.0938 3140.2875 5484.2628
[15] 5555.2016 2942.4584 4969.7073 4827.8298


$beta
[1] 0.000863757


$tau
 [1] 2094.6446 3560.9242 3736.7723 2782.5965 2571.7786 2727.6440 3950.5877
 [8]  992.1430 1908.3516  947.1819 1171.9877 1049.0938 3140.2876 5484.2629
[15] 5555.2017 2942.4584 4969.7074 4827.8299


$beta
[1] 0.0008637381


$tau
 [1] 2095.8343 3560.9243 3736.7724 2782.5966 2571.7787 2727.6440 3950.5877
 [8]  992.1430 1908.3517  947.1819 1171.9877 1049.0939 3140.2876 5484.2630
[15] 5555.2018 2942.4585 4969.7075 4827.8300


$beta
[1] 0.0008637192


$tau
 [1] 2097.0231 3560.9244 3736.7725 2782.5966 2571.7787 2727.6441 3950.5878
 [8]  992.1431 1908.3517  947.1819 1171.9877 1049.0939 3140.2877 5484.2631
[15] 5555.2019 2942.4585 4969.7076 4827.8301


$beta
[1] 0.0008637003


$tau
 [1] 2098.2109 3560.9244 3736.7725 2782.5967 2571.7788 2727.6441 3950.5879
 [8]  992.1431 1908.3518  947.1819 1171.9878 1049.0939 3140.2877 5484.2632
[15] 5555.2020 2942.4586 4969.7077 4827.8302


$beta
[1] 0.0008636815


$tau
 [1] 2099.3977 3560.9245 3736.7726 2782.5967 2571.7788 2727.6442 3950.5880
 [8]  992.1431 1908.3518  947.1819 1171.9878 1049.0939 3140.2878 5484.2633
[15] 5555.2021 2942.4587 4969.7078 4827.8303


$beta
[1] 0.0008636626


$tau
 [1] 2100.5835 3560.9246 3736.7727 2782.5968 2571.7789 2727.6442 3950.5880
 [8]  992.1431 1908.3518  947.1820 1171.9878 1049.0939 3140.2879 5484.2634
[15] 5555.2022 2942.4587 4969.7078 4827.8304


$beta
[1] 0.0008636438


$tau
 [1] 2101.7683 3560.9246 3736.7728 2782.5968 2571.7789 2727.6443 3950.5881
 [8]  992.1431 1908.3519  947.1820 1171.9878 1049.0940 3140.2879 5484.2635
[15] 5555.2023 2942.4588 4969.7079 4827.8305


$beta
[1] 0.000863625


$tau
 [1] 2102.9521 3560.9247 3736.7728 2782.5969 2571.7789 2727.6443 3950.5882
 [8]  992.1432 1908.3519  947.1820 1171.9878 1049.0940 3140.2880 5484.2637
[15] 5555.2024 2942.4588 4969.7080 4827.8306


$beta
[1] 0.0008636062


$tau
 [1] 2104.1350 3560.9248 3736.7729 2782.5969 2571.7790 2727.6444 3950.5882
 [8]  992.1432 1908.3519  947.1820 1171.9879 1049.0940 3140.2880 5484.2638
[15] 5555.2025 2942.4589 4969.7081 4827.8307


$beta
[1] 0.0008635875


$tau
 [1] 2105.3169 3560.9248 3736.7730 2782.5970 2571.7790 2727.6444 3950.5883
 [8]  992.1432 1908.3520  947.1820 1171.9879 1049.0940 3140.2881 5484.2639
[15] 5555.2026 2942.4589 4969.7082 4827.8307


$beta
[1] 0.0008635687


$tau
 [1] 2106.4978 3560.9249 3736.7730 2782.5970 2571.7791 2727.6445 3950.5884
 [8]  992.1432 1908.3520  947.1820 1171.9879 1049.0940 3140.2881 5484.2640
[15] 5555.2027 2942.4590 4969.7083 4827.8308


$beta
[1] 0.00086355


$tau
 [1] 2107.6777 3560.9250 3736.7731 2782.5971 2571.7791 2727.6445 3950.5885
 [8]  992.1432 1908.3520  947.1821 1171.9879 1049.0941 3140.2882 5484.2641
[15] 5555.2028 2942.4590 4969.7084 4827.8309


$beta
[1] 0.0008635312


$tau
 [1] 2108.8566 3560.9250 3736.7732 2782.5971 2571.7792 2727.6446 3950.5885
 [8]  992.1433 1908.3521  947.1821 1171.9880 1049.0941 3140.2883 5484.2642
[15] 5555.2029 2942.4591 4969.7085 4827.8310


$beta
[1] 0.0008635125


$tau
 [1] 2110.0346 3560.9251 3736.7732 2782.5972 2571.7792 2727.6446 3950.5886
 [8]  992.1433 1908.3521  947.1821 1171.9880 1049.0941 3140.2883 5484.2643
[15] 5555.2030 2942.4591 4969.7086 4827.8311


$beta
[1] 0.0008634938


$tau
 [1] 2111.2116 3560.9252 3736.7733 2782.5972 2571.7793 2727.6447 3950.5887
 [8]  992.1433 1908.3521  947.1821 1171.9880 1049.0941 3140.2884 5484.2644
[15] 5555.2031 2942.4592 4969.7087 4827.8312


$beta
[1] 0.0008634752


$tau
 [1] 2112.3876 3560.9252 3736.7734 2782.5973 2571.7793 2727.6447 3950.5888
 [8]  992.1433 1908.3522  947.1821 1171.9880 1049.0941 3140.2884 5484.2645
[15] 5555.2032 2942.4593 4969.7088 4827.8313


$beta
[1] 0.0008634565


$tau
 [1] 2113.5626 3560.9253 3736.7735 2782.5973 2571.7794 2727.6448 3950.5888
 [8]  992.1433 1908.3522  947.1821 1171.9880 1049.0942 3140.2885 5484.2646
[15] 5555.2033 2942.4593 4969.7089 4827.8314


$beta
[1] 0.0008634379


$tau
 [1] 2114.7366 3560.9254 3736.7735 2782.5974 2571.7794 2727.6448 3950.5889
 [8]  992.1433 1908.3523  947.1822 1171.9881 1049.0942 3140.2886 5484.2647
[15] 5555.2034 2942.4594 4969.7090 4827.8315


$beta
[1] 0.0008634193


$tau
 [1] 2115.9097 3560.9254 3736.7736 2782.5975 2571.7795 2727.6449 3950.5890
 [8]  992.1434 1908.3523  947.1822 1171.9881 1049.0942 3140.2886 5484.2648
[15] 5555.2035 2942.4594 4969.7091 4827.8316


$beta
[1] 0.0008634007


$tau
 [1] 2117.0818 3560.9255 3736.7737 2782.5975 2571.7795 2727.6449 3950.5891
 [8]  992.1434 1908.3523  947.1822 1171.9881 1049.0942 3140.2887 5484.2649
[15] 5555.2036 2942.4595 4969.7091 4827.8316


$beta
[1] 0.0008633821


$tau
 [1] 2118.2530 3560.9256 3736.7737 2782.5976 2571.7796 2727.6450 3950.5891
 [8]  992.1434 1908.3524  947.1822 1171.9881 1049.0942 3140.2887 5484.2650
[15] 5555.2037 2942.4595 4969.7092 4827.8317


$beta
[1] 0.0008633635


$tau
 [1] 2119.4231 3560.9256 3736.7738 2782.5976 2571.7796 2727.6451 3950.5892
 [8]  992.1434 1908.3524  947.1822 1171.9881 1049.0943 3140.2888 5484.2651
[15] 5555.2038 2942.4596 4969.7093 4827.8318


$beta
[1] 0.0008633449


$tau
 [1] 2120.5923 3560.9257 3736.7739 2782.5977 2571.7797 2727.6451 3950.5893
 [8]  992.1434 1908.3524  947.1823 1171.9882 1049.0943 3140.2888 5484.2652
[15] 5555.2039 2942.4596 4969.7094 4827.8319


$beta
[1] 0.0008633264


$tau
 [1] 2121.7606 3560.9258 3736.7739 2782.5977 2571.7797 2727.6452 3950.5894
 [8]  992.1435 1908.3525  947.1823 1171.9882 1049.0943 3140.2889 5484.2653
[15] 5555.2040 2942.4597 4969.7095 4827.8320


$beta
[1] 0.0008633079


$tau
 [1] 2122.9278 3560.9258 3736.7740 2782.5978 2571.7798 2727.6452 3950.5894
 [8]  992.1435 1908.3525  947.1823 1171.9882 1049.0943 3140.2890 5484.2654
[15] 5555.2041 2942.4598 4969.7096 4827.8321


$beta
[1] 0.0008632894


$tau
 [1] 2124.0941 3560.9259 3736.7741 2782.5978 2571.7798 2727.6453 3950.5895
 [8]  992.1435 1908.3525  947.1823 1171.9882 1049.0943 3140.2890 5484.2655
[15] 5555.2043 2942.4598 4969.7097 4827.8322


$beta
[1] 0.0008632709


$tau
 [1] 2125.2594 3560.9260 3736.7742 2782.5979 2571.7799 2727.6453 3950.5896
 [8]  992.1435 1908.3526  947.1823 1171.9883 1049.0943 3140.2891 5484.2656
[15] 5555.2044 2942.4599 4969.7098 4827.8323


$beta
[1] 0.0008632524


$tau
 [1] 2126.4238 3560.9260 3736.7742 2782.5979 2571.7799 2727.6454 3950.5896
 [8]  992.1435 1908.3526  947.1823 1171.9883 1049.0944 3140.2891 5484.2657
[15] 5555.2045 2942.4599 4969.7099 4827.8324


$beta
[1] 0.0008632339


$tau
 [1] 2127.5872 3560.9261 3736.7743 2782.5980 2571.7800 2727.6454 3950.5897
 [8]  992.1435 1908.3526  947.1824 1171.9883 1049.0944 3140.2892 5484.2658
[15] 5555.2046 2942.4600 4969.7100 4827.8325


$beta
[1] 0.0008632155


$tau
 [1] 2128.7496 3560.9262 3736.7744 2782.5980 2571.7800 2727.6455 3950.5898
 [8]  992.1436 1908.3527  947.1824 1171.9883 1049.0944 3140.2893 5484.2659
[15] 5555.2047 2942.4600 4969.7101 4827.8325


$beta
[1] 0.000863197


$tau
 [1] 2129.9111 3560.9262 3736.7744 2782.5981 2571.7800 2727.6455 3950.5899
 [8]  992.1436 1908.3527  947.1824 1171.9883 1049.0944 3140.2893 5484.2660
[15] 5555.2048 2942.4601 4969.7102 4827.8326


$beta
[1] 0.0008631786


$tau
 [1] 2131.0716 3560.9263 3736.7745 2782.5981 2571.7801 2727.6456 3950.5899
 [8]  992.1436 1908.3528  947.1824 1171.9884 1049.0944 3140.2894 5484.2661
[15] 5555.2049 2942.4601 4969.7103 4827.8327


$beta
[1] 0.0008631602


$tau
 [1] 2132.2311 3560.9263 3736.7746 2782.5982 2571.7801 2727.6456 3950.5900
 [8]  992.1436 1908.3528  947.1824 1171.9884 1049.0945 3140.2894 5484.2662
[15] 5555.2050 2942.4602 4969.7103 4827.8328


$beta
[1] 0.0008631419


$tau
 [1] 2133.3897 3560.9264 3736.7746 2782.5982 2571.7802 2727.6457 3950.5901
 [8]  992.1436 1908.3528  947.1824 1171.9884 1049.0945 3140.2895 5484.2663
[15] 5555.2051 2942.4602 4969.7104 4827.8329


$beta
[1] 0.0008631235


$tau
 [1] 2134.5474 3560.9265 3736.7747 2782.5983 2571.7802 2727.6457 3950.5902
 [8]  992.1437 1908.3529  947.1825 1171.9884 1049.0945 3140.2895 5484.2664
[15] 5555.2052 2942.4603 4969.7105 4827.8330


$beta
[1] 0.0008631051


$tau
 [1] 2135.7040 3560.9265 3736.7748 2782.5983 2571.7803 2727.6458 3950.5902
 [8]  992.1437 1908.3529  947.1825 1171.9885 1049.0945 3140.2896 5484.2665
[15] 5555.2053 2942.4603 4969.7106 4827.8331


$beta
[1] 0.0008630868


$tau
 [1] 2136.8597 3560.9266 3736.7748 2782.5984 2571.7803 2727.6458 3950.5903
 [8]  992.1437 1908.3529  947.1825 1171.9885 1049.0945 3140.2897 5484.2666
[15] 5555.2054 2942.4604 4969.7107 4827.8332


$beta
[1] 0.0008630685


$tau
 [1] 2138.0145 3560.9267 3736.7749 2782.5984 2571.7804 2727.6459 3950.5904
 [8]  992.1437 1908.3530  947.1825 1171.9885 1049.0946 3140.2897 5484.2667
[15] 5555.2055 2942.4605 4969.7108 4827.8332


$beta
[1] 0.0008630502


$tau
 [1] 2139.1683 3560.9267 3736.7750 2782.5985 2571.7804 2727.6459 3950.5904
 [8]  992.1437 1908.3530  947.1825 1171.9885 1049.0946 3140.2898 5484.2668
[15] 5555.2056 2942.4605 4969.7109 4827.8333


$beta
[1] 0.0008630319


$tau
 [1] 2140.3211 3560.9268 3736.7750 2782.5985 2571.7805 2727.6460 3950.5905
 [8]  992.1437 1908.3530  947.1826 1171.9885 1049.0946 3140.2898 5484.2669
[15] 5555.2057 2942.4606 4969.7110 4827.8334


$beta
[1] 0.0008630136


$tau
 [1] 2141.4730 3560.9269 3736.7751 2782.5986 2571.7805 2727.6460 3950.5906
 [8]  992.1438 1908.3531  947.1826 1171.9886 1049.0946 3140.2899 5484.2670
[15] 5555.2058 2942.4606 4969.7111 4827.8335


$beta
[1] 0.0008629954


$tau
 [1] 2142.6239 3560.9269 3736.7752 2782.5986 2571.7806 2727.6461 3950.5907
 [8]  992.1438 1908.3531  947.1826 1171.9886 1049.0946 3140.2899 5484.2671
[15] 5555.2059 2942.4607 4969.7112 4827.8336


$beta
[1] 0.0008629772


$tau
 [1] 2143.7739 3560.9270 3736.7752 2782.5987 2571.7806 2727.6461 3950.5907
 [8]  992.1438 1908.3531  947.1826 1171.9886 1049.0947 3140.2900 5484.2672
[15] 5555.2060 2942.4607 4969.7113 4827.8337


$beta
[1] 0.0008629589


$tau
 [1] 2144.9229 3560.9271 3736.7753 2782.5987 2571.7807 2727.6462 3950.5908
 [8]  992.1438 1908.3532  947.1826 1171.9886 1049.0947 3140.2901 5484.2673
[15] 5555.2061 2942.4608 4969.7113 4827.8338


$beta
[1] 0.0008629407


$tau
 [1] 2146.0710 3560.9271 3736.7754 2782.5988 2571.7807 2727.6462 3950.5909
 [8]  992.1438 1908.3532  947.1826 1171.9886 1049.0947 3140.2901 5484.2674
[15] 5555.2062 2942.4608 4969.7114 4827.8339


$beta
[1] 0.0008629225


$tau
 [1] 2147.2181 3560.9272 3736.7755 2782.5988 2571.7808 2727.6463 3950.5909
 [8]  992.1439 1908.3532  947.1827 1171.9887 1049.0947 3140.2902 5484.2675
[15] 5555.2063 2942.4609 4969.7115 4827.8340


$beta
[1] 0.0008629044


$tau
 [1] 2148.3642 3560.9273 3736.7755 2782.5989 2571.7808 2727.6463 3950.5910
 [8]  992.1439 1908.3533  947.1827 1171.9887 1049.0947 3140.2902 5484.2676
[15] 5555.2064 2942.4609 4969.7116 4827.8340


$beta
[1] 0.0008628862


$tau
 [1] 2149.5095 3560.9273 3736.7756 2782.5989 2571.7808 2727.6464 3950.5911
 [8]  992.1439 1908.3533  947.1827 1171.9887 1049.0948 3140.2903 5484.2677
[15] 5555.2065 2942.4610 4969.7117 4827.8341


$beta
[1] 0.0008628681


$tau
 [1] 2150.6537 3560.9274 3736.7757 2782.5990 2571.7809 2727.6464 3950.5912
 [8]  992.1439 1908.3533  947.1827 1171.9887 1049.0948 3140.2903 5484.2678
[15] 5555.2066 2942.4610 4969.7118 4827.8342


$beta
[1] 0.0008628499


$tau
 [1] 2151.7971 3560.9275 3736.7757 2782.5990 2571.7809 2727.6465 3950.5912
 [8]  992.1439 1908.3534  947.1827 1171.9887 1049.0948 3140.2904 5484.2679
[15] 5555.2067 2942.4611 4969.7119 4827.8343


$beta
[1] 0.0008628318


$tau
 [1] 2152.9394 3560.9275 3736.7758 2782.5991 2571.7810 2727.6465 3950.5913
 [8]  992.1439 1908.3534  947.1827 1171.9888 1049.0948 3140.2905 5484.2680
[15] 5555.2068 2942.4612 4969.7120 4827.8344


$beta
[1] 0.0008628137


$tau
 [1] 2154.0809 3560.9276 3736.7759 2782.5991 2571.7810 2727.6465 3950.5914
 [8]  992.1440 1908.3534  947.1828 1171.9888 1049.0948 3140.2905 5484.2681
[15] 5555.2069 2942.4612 4969.7121 4827.8345


$beta
[1] 0.0008627957


$tau
 [1] 2155.2213 3560.9276 3736.7759 2782.5992 2571.7811 2727.6466 3950.5914
 [8]  992.1440 1908.3535  947.1828 1171.9888 1049.0948 3140.2906 5484.2682
[15] 5555.2070 2942.4613 4969.7122 4827.8346


$beta
[1] 0.0008627776


$tau
 [1] 2156.3609 3560.9277 3736.7760 2782.5992 2571.7811 2727.6466 3950.5915
 [8]  992.1440 1908.3535  947.1828 1171.9888 1049.0949 3140.2906 5484.2683
[15] 5555.2071 2942.4613 4969.7122 4827.8347


$beta
[1] 0.0008627595


$tau
 [1] 2157.4995 3560.9278 3736.7761 2782.5993 2571.7812 2727.6467 3950.5916
 [8]  992.1440 1908.3535  947.1828 1171.9889 1049.0949 3140.2907 5484.2684
[15] 5555.2072 2942.4614 4969.7123 4827.8347


$beta
[1] 0.0008627415


$tau
 [1] 2158.6371 3560.9278 3736.7761 2782.5993 2571.7812 2727.6467 3950.5917
 [8]  992.1440 1908.3536  947.1828 1171.9889 1049.0949 3140.2907 5484.2685
[15] 5555.2073 2942.4614 4969.7124 4827.8348


$beta
[1] 0.0008627235


$tau
 [1] 2159.7738 3560.9279 3736.7762 2782.5994 2571.7813 2727.6468 3950.5917
 [8]  992.1441 1908.3536  947.1828 1171.9889 1049.0949 3140.2908 5484.2686
[15] 5555.2074 2942.4615 4969.7125 4827.8349


$beta
[1] 0.0008627055


$tau
 [1] 2160.9096 3560.9280 3736.7763 2782.5994 2571.7813 2727.6468 3950.5918
 [8]  992.1441 1908.3537  947.1829 1171.9889 1049.0949 3140.2909 5484.2687
[15] 5555.2075 2942.4615 4969.7126 4827.8350


$beta
[1] 0.0008626875


$tau
 [1] 2162.0444 3560.9280 3736.7763 2782.5995 2571.7814 2727.6469 3950.5919
 [8]  992.1441 1908.3537  947.1829 1171.9889 1049.0950 3140.2909 5484.2688
[15] 5555.2076 2942.4616 4969.7127 4827.8351


$beta
[1] 0.0008626695


$tau
 [1] 2163.1783 3560.9281 3736.7764 2782.5995 2571.7814 2727.6469 3950.5919
 [8]  992.1441 1908.3537  947.1829 1171.9890 1049.0950 3140.2910 5484.2689
[15] 5555.2077 2942.4616 4969.7128 4827.8352


$beta
[1] 0.0008626516


$tau
 [1] 2164.3113 3560.9282 3736.7765 2782.5996 2571.7814 2727.6470 3950.5920
 [8]  992.1441 1908.3538  947.1829 1171.9890 1049.0950 3140.2910 5484.2690
[15] 5555.2078 2942.4617 4969.7129 4827.8353


$beta
[1] 0.0008626336


$tau
 [1] 2165.4433 3560.9282 3736.7765 2782.5996 2571.7815 2727.6470 3950.5921
 [8]  992.1441 1908.3538  947.1829 1171.9890 1049.0950 3140.2911 5484.2691
[15] 5555.2079 2942.4617 4969.7130 4827.8353


$beta
[1] 0.0008626157


$tau
 [1] 2166.5744 3560.9283 3736.7766 2782.5997 2571.7815 2727.6471 3950.5922
 [8]  992.1442 1908.3538  947.1829 1171.9890 1049.0950 3140.2911 5484.2692
[15] 5555.2080 2942.4618 4969.7130 4827.8354


$beta
[1] 0.0008625978


$tau
 [1] 2167.7045 3560.9284 3736.7767 2782.5997 2571.7816 2727.6471 3950.5922
 [8]  992.1442 1908.3539  947.1830 1171.9890 1049.0951 3140.2912 5484.2693
[15] 5555.2081 2942.4618 4969.7131 4827.8355


$beta
[1] 0.0008625799


$tau
 [1] 2168.8337 3560.9284 3736.7767 2782.5998 2571.7816 2727.6472 3950.5923
 [8]  992.1442 1908.3539  947.1830 1171.9891 1049.0951 3140.2912 5484.2694
[15] 5555.2082 2942.4619 4969.7132 4827.8356


$beta
[1] 0.000862562


$tau
 [1] 2169.9620 3560.9285 3736.7768 2782.5998 2571.7817 2727.6472 3950.5924
 [8]  992.1442 1908.3539  947.1830 1171.9891 1049.0951 3140.2913 5484.2695
[15] 5555.2083 2942.4619 4969.7133 4827.8357


$beta
[1] 0.0008625442


$tau
 [1] 2171.0893 3560.9285 3736.7769 2782.5999 2571.7817 2727.6473 3950.5924
 [8]  992.1442 1908.3540  947.1830 1171.9891 1049.0951 3140.2914 5484.2696
[15] 5555.2084 2942.4620 4969.7134 4827.8358


$beta
[1] 0.0008625263


$tau
 [1] 2172.2158 3560.9286 3736.7769 2782.5999 2571.7818 2727.6473 3950.5925
 [8]  992.1442 1908.3540  947.1830 1171.9891 1049.0951 3140.2914 5484.2697
[15] 5555.2085 2942.4620 4969.7135 4827.8359


$beta
[1] 0.0008625085


$tau
 [1] 2173.3412 3560.9287 3736.7770 2782.6000 2571.7818 2727.6474 3950.5926
 [8]  992.1443 1908.3540  947.1830 1171.9891 1049.0951 3140.2915 5484.2698
[15] 5555.2086 2942.4621 4969.7136 4827.8360


$beta
[1] 0.0008624907


$tau
 [1] 2174.4658 3560.9287 3736.7771 2782.6000 2571.7819 2727.6474 3950.5927
 [8]  992.1443 1908.3541  947.1831 1171.9892 1049.0952 3140.2915 5484.2699
[15] 5555.2087 2942.4622 4969.7137 4827.8360


$beta
[1] 0.0008624729


$tau
 [1] 2175.5894 3560.9288 3736.7771 2782.6001 2571.7819 2727.6475 3950.5927
 [8]  992.1443 1908.3541  947.1831 1171.9892 1049.0952 3140.2916 5484.2700
[15] 5555.2088 2942.4622 4969.7138 4827.8361


$beta
[1] 0.0008624551


$tau
 [1] 2176.7121 3560.9289 3736.7772 2782.6001 2571.7820 2727.6475 3950.5928
 [8]  992.1443 1908.3541  947.1831 1171.9892 1049.0952 3140.2916 5484.2701
[15] 5555.2089 2942.4623 4969.7138 4827.8362


$beta
[1] 0.0008624374


$tau
 [1] 2177.8338 3560.9289 3736.7773 2782.6002 2571.7820 2727.6476 3950.5929
 [8]  992.1443 1908.3542  947.1831 1171.9892 1049.0952 3140.2917 5484.2702
[15] 5555.2090 2942.4623 4969.7139 4827.8363


$beta
[1] 0.0008624196


$tau
 [1] 2178.9547 3560.9290 3736.7773 2782.6002 2571.7820 2727.6476 3950.5929
 [8]  992.1444 1908.3542  947.1831 1171.9893 1049.0952 3140.2918 5484.2703
[15] 5555.2091 2942.4624 4969.7140 4827.8364


$beta
[1] 0.0008624019


$tau
 [1] 2180.0746 3560.9290 3736.7774 2782.6003 2571.7821 2727.6477 3950.5930
 [8]  992.1444 1908.3542  947.1831 1171.9893 1049.0953 3140.2918 5484.2704
[15] 5555.2092 2942.4624 4969.7141 4827.8365


$beta
[1] 0.0008623841


$tau
 [1] 2181.1935 3560.9291 3736.7775 2782.6003 2571.7821 2727.6477 3950.5931
 [8]  992.1444 1908.3543  947.1832 1171.9893 1049.0953 3140.2919 5484.2705
[15] 5555.2093 2942.4625 4969.7142 4827.8366


$beta
[1] 0.0008623664


$tau
 [1] 2182.3116 3560.9292 3736.7775 2782.6004 2571.7822 2727.6478 3950.5931
 [8]  992.1444 1908.3543  947.1832 1171.9893 1049.0953 3140.2919 5484.2705
[15] 5555.2094 2942.4625 4969.7143 4827.8366


$beta
[1] 0.0008623487


$tau
 [1] 2183.4287 3560.9292 3736.7776 2782.6004 2571.7822 2727.6478 3950.5932
 [8]  992.1444 1908.3543  947.1832 1171.9893 1049.0953 3140.2920 5484.2706
[15] 5555.2095 2942.4626 4969.7144 4827.8367


$beta
[1] 0.0008623311


$tau
 [1] 2184.5449 3560.9293 3736.7777 2782.6005 2571.7823 2727.6479 3950.5933
 [8]  992.1444 1908.3544  947.1832 1171.9894 1049.0953 3140.2920 5484.2707
[15] 5555.2096 2942.4626 4969.7145 4827.8368


$beta
[1] 0.0008623134


$tau
 [1] 2185.6602 3560.9294 3736.7777 2782.6005 2571.7823 2727.6479 3950.5934
 [8]  992.1445 1908.3544  947.1832 1171.9894 1049.0954 3140.2921 5484.2708
[15] 5555.2097 2942.4627 4969.7145 4827.8369


$beta
[1] 0.0008622958


$tau
 [1] 2186.7745 3560.9294 3736.7778 2782.6006 2571.7824 2727.6480 3950.5934
 [8]  992.1445 1908.3544  947.1832 1171.9894 1049.0954 3140.2921 5484.2709
[15] 5555.2098 2942.4627 4969.7146 4827.8370


$beta
[1] 0.0008622781


$tau
 [1] 2187.8880 3560.9295 3736.7779 2782.6006 2571.7824 2727.6480 3950.5935
 [8]  992.1445 1908.3545  947.1833 1171.9894 1049.0954 3140.2922 5484.2710
[15] 5555.2099 2942.4628 4969.7147 4827.8371


$beta
[1] 0.0008622605


$tau
 [1] 2189.0005 3560.9296 3736.7779 2782.6007 2571.7825 2727.6481 3950.5936
 [8]  992.1445 1908.3545  947.1833 1171.9894 1049.0954 3140.2922 5484.2711
[15] 5555.2100 2942.4628 4969.7148 4827.8371


$beta
[1] 0.0008622429


$tau
 [1] 2190.1121 3560.9296 3736.7780 2782.6007 2571.7825 2727.6481 3950.5936
 [8]  992.1445 1908.3545  947.1833 1171.9895 1049.0954 3140.2923 5484.2712
[15] 5555.2101 2942.4629 4969.7149 4827.8372


$beta
[1] 0.0008622253


$tau
 [1] 2191.2227 3560.9297 3736.7781 2782.6008 2571.7825 2727.6482 3950.5937
 [8]  992.1445 1908.3546  947.1833 1171.9895 1049.0954 3140.2924 5484.2713
[15] 5555.2102 2942.4629 4969.7150 4827.8373


$beta
[1] 0.0008622078


$tau
 [1] 2192.3325 3560.9297 3736.7781 2782.6008 2571.7826 2727.6482 3950.5938
 [8]  992.1446 1908.3546  947.1833 1171.9895 1049.0955 3140.2924 5484.2714
[15] 5555.2103 2942.4630 4969.7151 4827.8374


$beta
[1] 0.0008621902


$tau
 [1] 2193.4413 3560.9298 3736.7782 2782.6009 2571.7826 2727.6482 3950.5938
 [8]  992.1446 1908.3546  947.1833 1171.9895 1049.0955 3140.2925 5484.2715
[15] 5555.2104 2942.4630 4969.7152 4827.8375


$beta
[1] 0.0008621727


$tau
 [1] 2194.5492 3560.9299 3736.7783 2782.6009 2571.7827 2727.6483 3950.5939
 [8]  992.1446 1908.3547  947.1834 1171.9895 1049.0955 3140.2925 5484.2716
[15] 5555.2104 2942.4631 4969.7152 4827.8376


$beta
[1] 0.0008621552


$tau
 [1] 2195.6562 3560.9299 3736.7783 2782.6010 2571.7827 2727.6483 3950.5940
 [8]  992.1446 1908.3547  947.1834 1171.9896 1049.0955 3140.2926 5484.2717
[15] 5555.2105 2942.4631 4969.7153 4827.8377


$beta
[1] 0.0008621377


$tau
 [1] 2196.7623 3560.9300 3736.7784 2782.6010 2571.7828 2727.6484 3950.5940
 [8]  992.1446 1908.3547  947.1834 1171.9896 1049.0955 3140.2926 5484.2718
[15] 5555.2106 2942.4632 4969.7154 4827.8377


$beta
[1] 0.0008621202


$tau
 [1] 2197.8675 3560.9301 3736.7784 2782.6011 2571.7828 2727.6484 3950.5941
 [8]  992.1447 1908.3548  947.1834 1171.9896 1049.0956 3140.2927 5484.2719
[15] 5555.2107 2942.4632 4969.7155 4827.8378


$beta
[1] 0.0008621027


$tau
 [1] 2198.9718 3560.9301 3736.7785 2782.6011 2571.7829 2727.6485 3950.5942
 [8]  992.1447 1908.3548  947.1834 1171.9896 1049.0956 3140.2927 5484.2720
[15] 5555.2108 2942.4633 4969.7156 4827.8379


$beta
[1] 0.0008620852


$tau
 [1] 2200.0751 3560.9302 3736.7786 2782.6012 2571.7829 2727.6485 3950.5943
 [8]  992.1447 1908.3548  947.1834 1171.9896 1049.0956 3140.2928 5484.2721
[15] 5555.2109 2942.4633 4969.7157 4827.8380


$beta
[1] 0.0008620678


$tau
 [1] 2201.1776 3560.9302 3736.7786 2782.6012 2571.7829 2727.6486 3950.5943
 [8]  992.1447 1908.3549  947.1835 1171.9897 1049.0956 3140.2929 5484.2722
[15] 5555.2110 2942.4634 4969.7158 4827.8381


$beta
[1] 0.0008620503


$tau
 [1] 2202.2791 3560.9303 3736.7787 2782.6013 2571.7830 2727.6486 3950.5944
 [8]  992.1447 1908.3549  947.1835 1171.9897 1049.0956 3140.2929 5484.2723
[15] 5555.2111 2942.4634 4969.7159 4827.8382


$beta
[1] 0.0008620329


$tau
 [1] 2203.3797 3560.9304 3736.7788 2782.6013 2571.7830 2727.6487 3950.5945
 [8]  992.1447 1908.3549  947.1835 1171.9897 1049.0956 3140.2930 5484.2724
[15] 5555.2112 2942.4635 4969.7159 4827.8382


$beta
[1] 0.0008620155


$tau
 [1] 2204.4794 3560.9304 3736.7788 2782.6014 2571.7831 2727.6487 3950.5945
 [8]  992.1448 1908.3550  947.1835 1171.9897 1049.0957 3140.2930 5484.2725
[15] 5555.2113 2942.4636 4969.7160 4827.8383


$beta
[1] 0.0008619981


$tau
 [1] 2205.5782 3560.9305 3736.7789 2782.6014 2571.7831 2727.6488 3950.5946
 [8]  992.1448 1908.3550  947.1835 1171.9897 1049.0957 3140.2931 5484.2726
[15] 5555.2114 2942.4636 4969.7161 4827.8384


$beta
[1] 0.0008619808


$tau
 [1] 2206.6761 3560.9305 3736.7790 2782.6015 2571.7832 2727.6488 3950.5947
 [8]  992.1448 1908.3550  947.1835 1171.9898 1049.0957 3140.2931 5484.2727
[15] 5555.2115 2942.4637 4969.7162 4827.8385


$beta
[1] 0.0008619634


$tau
 [1] 2207.7731 3560.9306 3736.7790 2782.6015 2571.7832 2727.6489 3950.5947
 [8]  992.1448 1908.3551  947.1836 1171.9898 1049.0957 3140.2932 5484.2728
[15] 5555.2116 2942.4637 4969.7163 4827.8386


$beta
[1] 0.0008619461


$tau
 [1] 2208.8691 3560.9307 3736.7791 2782.6015 2571.7833 2727.6489 3950.5948
 [8]  992.1448 1908.3551  947.1836 1171.9898 1049.0957 3140.2932 5484.2729
[15] 5555.2117 2942.4638 4969.7164 4827.8387


$beta
[1] 0.0008619287


$tau
 [1] 2209.9643 3560.9307 3736.7792 2782.6016 2571.7833 2727.6490 3950.5949
 [8]  992.1448 1908.3551  947.1836 1171.9898 1049.0958 3140.2933 5484.2729
[15] 5555.2118 2942.4638 4969.7165 4827.8387


$beta
[1] 0.0008619114


$tau
 [1] 2211.0586 3560.9308 3736.7792 2782.6016 2571.7834 2727.6490 3950.5949
 [8]  992.1449 1908.3552  947.1836 1171.9898 1049.0958 3140.2933 5484.2730
[15] 5555.2119 2942.4639 4969.7165 4827.8388


$beta
[1] 0.0008618941


$tau
 [1] 2212.1519 3560.9309 3736.7793 2782.6017 2571.7834 2727.6491 3950.5950
 [8]  992.1449 1908.3552  947.1836 1171.9899 1049.0958 3140.2934 5484.2731
[15] 5555.2120 2942.4639 4969.7166 4827.8389


$beta
[1] 0.0008618769


$tau
 [1] 2213.2444 3560.9309 3736.7794 2782.6017 2571.7834 2727.6491 3950.5951
 [8]  992.1449 1908.3552  947.1836 1171.9899 1049.0958 3140.2935 5484.2732
[15] 5555.2121 2942.4640 4969.7167 4827.8390


$beta
[1] 0.0008618596


$tau
 [1] 2214.3359 3560.9310 3736.7794 2782.6018 2571.7835 2727.6491 3950.5951
 [8]  992.1449 1908.3553  947.1837 1171.9899 1049.0958 3140.2935 5484.2733
[15] 5555.2122 2942.4640 4969.7168 4827.8391


$beta
[1] 0.0008618423


$tau
 [1] 2215.4266 3560.9310 3736.7795 2782.6018 2571.7835 2727.6492 3950.5952
 [8]  992.1449 1908.3553  947.1837 1171.9899 1049.0958 3140.2936 5484.2734
[15] 5555.2123 2942.4641 4969.7169 4827.8392


$beta
[1] 0.0008618251


$tau
 [1] 2216.5163 3560.9311 3736.7795 2782.6019 2571.7836 2727.6492 3950.5953
 [8]  992.1449 1908.3553  947.1837 1171.9899 1049.0959 3140.2936 5484.2735
[15] 5555.2124 2942.4641 4969.7170 4827.8392


$beta
[1] 0.0008618079


$tau
 [1] 2217.6052 3560.9312 3736.7796 2782.6019 2571.7836 2727.6493 3950.5953
 [8]  992.1450 1908.3554  947.1837 1171.9900 1049.0959 3140.2937 5484.2736
[15] 5555.2125 2942.4642 4969.7171 4827.8393


$beta
[1] 0.0008617907


$tau
 [1] 2218.6931 3560.9312 3736.7797 2782.6020 2571.7837 2727.6493 3950.5954
 [8]  992.1450 1908.3554  947.1837 1171.9900 1049.0959 3140.2937 5484.2737
[15] 5555.2126 2942.4642 4969.7171 4827.8394


$beta
[1] 0.0008617735


$tau
 [1] 2219.7802 3560.9313 3736.7797 2782.6020 2571.7837 2727.6494 3950.5955
 [8]  992.1450 1908.3554  947.1837 1171.9900 1049.0959 3140.2938 5484.2738
[15] 5555.2127 2942.4643 4969.7172 4827.8395


$beta
[1] 0.0008617563


$tau
 [1] 2220.8663 3560.9313 3736.7798 2782.6021 2571.7837 2727.6494 3950.5955
 [8]  992.1450 1908.3555  947.1838 1171.9900 1049.0959 3140.2938 5484.2739
[15] 5555.2128 2942.4643 4969.7173 4827.8396


$beta
[1] 0.0008617392


$tau
 [1] 2221.9516 3560.9314 3736.7799 2782.6021 2571.7838 2727.6495 3950.5956
 [8]  992.1450 1908.3555  947.1838 1171.9901 1049.0960 3140.2939 5484.2740
[15] 5555.2129 2942.4644 4969.7174 4827.8397


$beta
[1] 0.000861722


$tau
 [1] 2223.0360 3560.9315 3736.7799 2782.6022 2571.7838 2727.6495 3950.5957
 [8]  992.1450 1908.3555  947.1838 1171.9901 1049.0960 3140.2939 5484.2741
[15] 5555.2130 2942.4644 4969.7175 4827.8397


$beta
[1] 0.0008617049


$tau
 [1] 2224.1194 3560.9315 3736.7800 2782.6022 2571.7839 2727.6496 3950.5958
 [8]  992.1451 1908.3556  947.1838 1171.9901 1049.0960 3140.2940 5484.2742
[15] 5555.2130 2942.4645 4969.7176 4827.8398


$beta
[1] 0.0008616878


$tau
 [1] 2225.2020 3560.9316 3736.7801 2782.6023 2571.7839 2727.6496 3950.5958
 [8]  992.1451 1908.3556  947.1838 1171.9901 1049.0960 3140.2940 5484.2743
[15] 5555.2131 2942.4645 4969.7177 4827.8399


$beta
[1] 0.0008616707


$tau
 [1] 2226.2837 3560.9316 3736.7801 2782.6023 2571.7840 2727.6497 3950.5959
 [8]  992.1451 1908.3556  947.1838 1171.9901 1049.0960 3140.2941 5484.2744
[15] 5555.2132 2942.4646 4969.7177 4827.8400


$beta
[1] 0.0008616536


$tau
 [1] 2227.3645 3560.9317 3736.7802 2782.6024 2571.7840 2727.6497 3950.5960
 [8]  992.1451 1908.3557  947.1839 1171.9902 1049.0960 3140.2942 5484.2745
[15] 5555.2133 2942.4646 4969.7178 4827.8401


$beta
[1] 0.0008616365


$tau
 [1] 2228.4443 3560.9318 3736.7803 2782.6024 2571.7841 2727.6498 3950.5960
 [8]  992.1451 1908.3557  947.1839 1171.9902 1049.0961 3140.2942 5484.2745
[15] 5555.2134 2942.4647 4969.7179 4827.8402


$beta
[1] 0.0008616194


$tau
 [1] 2229.5233 3560.9318 3736.7803 2782.6025 2571.7841 2727.6498 3950.5961
 [8]  992.1451 1908.3557  947.1839 1171.9902 1049.0961 3140.2943 5484.2746
[15] 5555.2135 2942.4647 4969.7180 4827.8402


$beta
[1] 0.0008616024


$tau
 [1] 2230.6015 3560.9319 3736.7804 2782.6025 2571.7841 2727.6499 3950.5962
 [8]  992.1452 1908.3558  947.1839 1171.9902 1049.0961 3140.2943 5484.2747
[15] 5555.2136 2942.4648 4969.7181 4827.8403


$beta
[1] 0.0008615854


$tau
 [1] 2231.6787 3560.9320 3736.7804 2782.6026 2571.7842 2727.6499 3950.5962
 [8]  992.1452 1908.3558  947.1839 1171.9902 1049.0961 3140.2944 5484.2748
[15] 5555.2137 2942.4648 4969.7182 4827.8404


$beta
[1] 0.0008615684


$tau
 [1] 2232.7550 3560.9320 3736.7805 2782.6026 2571.7842 2727.6499 3950.5963
 [8]  992.1452 1908.3558  947.1839 1171.9903 1049.0961 3140.2944 5484.2749
[15] 5555.2138 2942.4649 4969.7182 4827.8405


$beta
[1] 0.0008615514


$tau
 [1] 2233.8304 3560.9321 3736.7806 2782.6026 2571.7843 2727.6500 3950.5964
 [8]  992.1452 1908.3559  947.1840 1171.9903 1049.0961 3140.2945 5484.2750
[15] 5555.2139 2942.4649 4969.7183 4827.8406


$beta
[1] 0.0008615344


$tau
 [1] 2234.9050 3560.9321 3736.7806 2782.6027 2571.7843 2727.6500 3950.5964
 [8]  992.1452 1908.3559  947.1840 1171.9903 1049.0962 3140.2945 5484.2751
[15] 5555.2140 2942.4650 4969.7184 4827.8406


$beta
[1] 0.0008615174


$tau
 [1] 2235.9787 3560.9322 3736.7807 2782.6027 2571.7844 2727.6501 3950.5965
 [8]  992.1452 1908.3559  947.1840 1171.9903 1049.0962 3140.2946 5484.2752
[15] 5555.2141 2942.4650 4969.7185 4827.8407


$beta
[1] 0.0008615004


$tau
 [1] 2237.0514 3560.9323 3736.7808 2782.6028 2571.7844 2727.6501 3950.5966
 [8]  992.1453 1908.3560  947.1840 1171.9903 1049.0962 3140.2946 5484.2753
[15] 5555.2142 2942.4651 4969.7186 4827.8408


$beta
[1] 0.0008614835


$tau
 [1] 2238.1233 3560.9323 3736.7808 2782.6028 2571.7845 2727.6502 3950.5966
 [8]  992.1453 1908.3560  947.1840 1171.9903 1049.0962 3140.2947 5484.2754
[15] 5555.2143 2942.4651 4969.7187 4827.8409


$beta
[1] 0.0008614666


$tau
 [1] 2239.1943 3560.9324 3736.7809 2782.6029 2571.7845 2727.6502 3950.5967
 [8]  992.1453 1908.3560  947.1840 1171.9904 1049.0962 3140.2947 5484.2755
[15] 5555.2144 2942.4652 4969.7188 4827.8410


$beta
[1] 0.0008614497


$tau
 [1] 2240.2645 3560.9324 3736.7810 2782.6029 2571.7845 2727.6503 3950.5968
 [8]  992.1453 1908.3560  947.1840 1171.9904 1049.0963 3140.2948 5484.2756
[15] 5555.2145 2942.4652 4969.7188 4827.8411


$beta
[1] 0.0008614328


$tau
 [1] 2241.3337 3560.9325 3736.7810 2782.6030 2571.7846 2727.6503 3950.5968
 [8]  992.1453 1908.3561  947.1841 1171.9904 1049.0963 3140.2948 5484.2757
[15] 5555.2146 2942.4653 4969.7189 4827.8411


$beta
[1] 0.0008614159


$tau
 [1] 2242.4021 3560.9326 3736.7811 2782.6030 2571.7846 2727.6504 3950.5969
 [8]  992.1454 1908.3561  947.1841 1171.9904 1049.0963 3140.2949 5484.2758
[15] 5555.2146 2942.4653 4969.7190 4827.8412


$beta
[1] 0.000861399


$tau
 [1] 2243.4695 3560.9326 3736.7811 2782.6031 2571.7847 2727.6504 3950.5970
 [8]  992.1454 1908.3561  947.1841 1171.9904 1049.0963 3140.2950 5484.2759
[15] 5555.2147 2942.4654 4969.7191 4827.8413


$beta
[1] 0.0008613822


$tau
 [1] 2244.5361 3560.9327 3736.7812 2782.6031 2571.7847 2727.6505 3950.5970
 [8]  992.1454 1908.3562  947.1841 1171.9905 1049.0963 3140.2950 5484.2759
[15] 5555.2148 2942.4654 4969.7192 4827.8414


$beta
[1] 0.0008613653


$tau
 [1] 2245.6019 3560.9327 3736.7813 2782.6032 2571.7848 2727.6505 3950.5971
 [8]  992.1454 1908.3562  947.1841 1171.9905 1049.0963 3140.2951 5484.2760
[15] 5555.2149 2942.4655 4969.7193 4827.8415


$beta
[1] 0.0008613485


$tau
 [1] 2246.6667 3560.9328 3736.7813 2782.6032 2571.7848 2727.6505 3950.5972
 [8]  992.1454 1908.3562  947.1841 1171.9905 1049.0964 3140.2951 5484.2761
[15] 5555.2150 2942.4655 4969.7193 4827.8415


$beta
[1] 0.0008613317


$tau
 [1] 2247.7307 3560.9329 3736.7814 2782.6033 2571.7848 2727.6506 3950.5972
 [8]  992.1454 1908.3563  947.1842 1171.9905 1049.0964 3140.2952 5484.2762
[15] 5555.2151 2942.4656 4969.7194 4827.8416


$beta
[1] 0.0008613149


$tau
 [1] 2248.7937 3560.9329 3736.7815 2782.6033 2571.7849 2727.6506 3950.5973
 [8]  992.1455 1908.3563  947.1842 1171.9905 1049.0964 3140.2952 5484.2763
[15] 5555.2152 2942.4656 4969.7195 4827.8417


$beta
[1] 0.0008612981


$tau
 [1] 2249.8559 3560.9330 3736.7815 2782.6033 2571.7849 2727.6507 3950.5974
 [8]  992.1455 1908.3563  947.1842 1171.9906 1049.0964 3140.2953 5484.2764
[15] 5555.2153 2942.4657 4969.7196 4827.8418


$beta
[1] 0.0008612814


$tau
 [1] 2250.9173 3560.9330 3736.7816 2782.6034 2571.7850 2727.6507 3950.5974
 [8]  992.1455 1908.3564  947.1842 1171.9906 1049.0964 3140.2953 5484.2765
[15] 5555.2154 2942.4657 4969.7197 4827.8419


$beta
[1] 0.0008612646


$tau
 [1] 2251.9777 3560.9331 3736.7816 2782.6034 2571.7850 2727.6508 3950.5975
 [8]  992.1455 1908.3564  947.1842 1171.9906 1049.0965 3140.2954 5484.2766
[15] 5555.2155 2942.4658 4969.7198 4827.8420


$beta
[1] 0.0008612479


$tau
 [1] 2253.0373 3560.9332 3736.7817 2782.6035 2571.7851 2727.6508 3950.5976
 [8]  992.1455 1908.3564  947.1842 1171.9906 1049.0965 3140.2954 5484.2767
[15] 5555.2156 2942.4658 4969.7198 4827.8420


$beta
[1] 0.0008612312


$tau
 [1] 2254.0960 3560.9332 3736.7818 2782.6035 2571.7851 2727.6509 3950.5976
 [8]  992.1455 1908.3565  947.1843 1171.9906 1049.0965 3140.2955 5484.2768
[15] 5555.2157 2942.4659 4969.7199 4827.8421


$beta
[1] 0.0008612144


$tau
 [1] 2255.1538 3560.9333 3736.7818 2782.6036 2571.7851 2727.6509 3950.5977
 [8]  992.1455 1908.3565  947.1843 1171.9907 1049.0965 3140.2955 5484.2769
[15] 5555.2158 2942.4659 4969.7200 4827.8422


$beta
[1] 0.0008611978


$tau
 [1] 2256.2108 3560.9333 3736.7819 2782.6036 2571.7852 2727.6510 3950.5978
 [8]  992.1456 1908.3565  947.1843 1171.9907 1049.0965 3140.2956 5484.2770
[15] 5555.2159 2942.4660 4969.7201 4827.8423


$beta
[1] 0.0008611811


$tau
 [1] 2257.2669 3560.9334 3736.7820 2782.6037 2571.7852 2727.6510 3950.5978
 [8]  992.1456 1908.3566  947.1843 1171.9907 1049.0965 3140.2956 5484.2770
[15] 5555.2160 2942.4660 4969.7202 4827.8424


$beta
[1] 0.0008611644


$tau
 [1] 2258.3221 3560.9335 3736.7820 2782.6037 2571.7853 2727.6510 3950.5979
 [8]  992.1456 1908.3566  947.1843 1171.9907 1049.0966 3140.2957 5484.2771
[15] 5555.2160 2942.4661 4969.7203 4827.8424


$beta
[1] 0.0008611478


$tau
 [1] 2259.3764 3560.9335 3736.7821 2782.6038 2571.7853 2727.6511 3950.5980
 [8]  992.1456 1908.3566  947.1843 1171.9907 1049.0966 3140.2957 5484.2772
[15] 5555.2161 2942.4661 4969.7203 4827.8425


$beta
[1] 0.0008611311


$tau
 [1] 2260.4299 3560.9336 3736.7821 2782.6038 2571.7854 2727.6511 3950.5980
 [8]  992.1456 1908.3567  947.1844 1171.9908 1049.0966 3140.2958 5484.2773
[15] 5555.2162 2942.4662 4969.7204 4827.8426


$beta
[1] 0.0008611145


$tau
 [1] 2261.4825 3560.9336 3736.7822 2782.6039 2571.7854 2727.6512 3950.5981
 [8]  992.1456 1908.3567  947.1844 1171.9908 1049.0966 3140.2958 5484.2774
[15] 5555.2163 2942.4662 4969.7205 4827.8427


$beta
[1] 0.0008610979


$tau
 [1] 2262.5342 3560.9337 3736.7823 2782.6039 2571.7854 2727.6512 3950.5981
 [8]  992.1457 1908.3567  947.1844 1171.9908 1049.0966 3140.2959 5484.2775
[15] 5555.2164 2942.4662 4969.7206 4827.8428


$beta
[1] 0.0008610813


$tau
 [1] 2263.5851 3560.9337 3736.7823 2782.6040 2571.7855 2727.6513 3950.5982
 [8]  992.1457 1908.3567  947.1844 1171.9908 1049.0966 3140.2959 5484.2776
[15] 5555.2165 2942.4663 4969.7207 4827.8428


$beta
[1] 0.0008610647


$tau
 [1] 2264.6351 3560.9338 3736.7824 2782.6040 2571.7855 2727.6513 3950.5983
 [8]  992.1457 1908.3568  947.1844 1171.9908 1049.0967 3140.2960 5484.2777
[15] 5555.2166 2942.4663 4969.7207 4827.8429


$beta
[1] 0.0008610482


$tau
 [1] 2265.6842 3560.9339 3736.7825 2782.6040 2571.7856 2727.6514 3950.5983
 [8]  992.1457 1908.3568  947.1844 1171.9909 1049.0967 3140.2961 5484.2778
[15] 5555.2167 2942.4664 4969.7208 4827.8430


$beta
[1] 0.0008610316


$tau
 [1] 2266.7325 3560.9339 3736.7825 2782.6041 2571.7856 2727.6514 3950.5984
 [8]  992.1457 1908.3568  947.1844 1171.9909 1049.0967 3140.2961 5484.2779
[15] 5555.2168 2942.4664 4969.7209 4827.8431


$beta
[1] 0.0008610151


$tau
 [1] 2267.7799 3560.9340 3736.7826 2782.6041 2571.7857 2727.6515 3950.5985
 [8]  992.1457 1908.3569  947.1845 1171.9909 1049.0967 3140.2962 5484.2780
[15] 5555.2169 2942.4665 4969.7210 4827.8432


$beta
[1] 0.0008609986


$tau
 [1] 2268.8265 3560.9340 3736.7826 2782.6042 2571.7857 2727.6515 3950.5985
 [8]  992.1458 1908.3569  947.1845 1171.9909 1049.0967 3140.2962 5484.2780
[15] 5555.2170 2942.4665 4969.7211 4827.8432


$beta
[1] 0.0008609821


$tau
 [1] 2269.8721 3560.9341 3736.7827 2782.6042 2571.7857 2727.6515 3950.5986
 [8]  992.1458 1908.3569  947.1845 1171.9909 1049.0967 3140.2963 5484.2781
[15] 5555.2171 2942.4666 4969.7212 4827.8433


$beta
[1] 0.0008609656


$tau
 [1] 2270.9170 3560.9342 3736.7828 2782.6043 2571.7858 2727.6516 3950.5987
 [8]  992.1458 1908.3570  947.1845 1171.9910 1049.0968 3140.2963 5484.2782
[15] 5555.2171 2942.4666 4969.7212 4827.8434


$beta
[1] 0.0008609491


$tau
 [1] 2271.9609 3560.9342 3736.7828 2782.6043 2571.7858 2727.6516 3950.5987
 [8]  992.1458 1908.3570  947.1845 1171.9910 1049.0968 3140.2964 5484.2783
[15] 5555.2172 2942.4667 4969.7213 4827.8435


$beta
[1] 0.0008609326


$tau
 [1] 2273.0040 3560.9343 3736.7829 2782.6044 2571.7859 2727.6517 3950.5988
 [8]  992.1458 1908.3570  947.1845 1171.9910 1049.0968 3140.2964 5484.2784
[15] 5555.2173 2942.4667 4969.7214 4827.8436


$beta
[1] 0.0008609162


$tau
 [1] 2274.0463 3560.9343 3736.7829 2782.6044 2571.7859 2727.6517 3950.5989
 [8]  992.1458 1908.3571  947.1846 1171.9910 1049.0968 3140.2965 5484.2785
[15] 5555.2174 2942.4668 4969.7215 4827.8436


$beta
[1] 0.0008608997


$tau
 [1] 2275.0876 3560.9344 3736.7830 2782.6045 2571.7860 2727.6518 3950.5989
 [8]  992.1459 1908.3571  947.1846 1171.9910 1049.0968 3140.2965 5484.2786
[15] 5555.2175 2942.4668 4969.7216 4827.8437


$beta
[1] 0.0008608833


$tau
 [1] 2276.1282 3560.9345 3736.7831 2782.6045 2571.7860 2727.6518 3950.5990
 [8]  992.1459 1908.3571  947.1846 1171.9911 1049.0969 3140.2966 5484.2787
[15] 5555.2176 2942.4669 4969.7216 4827.8438


$beta
[1] 0.0008608669


$tau
 [1] 2277.1678 3560.9345 3736.7831 2782.6045 2571.7860 2727.6519 3950.5991
 [8]  992.1459 1908.3572  947.1846 1171.9911 1049.0969 3140.2966 5484.2788
[15] 5555.2177 2942.4669 4969.7217 4827.8439


$beta
[1] 0.0008608505


$tau
 [1] 2278.2066 3560.9346 3736.7832 2782.6046 2571.7861 2727.6519 3950.5991
 [8]  992.1459 1908.3572  947.1846 1171.9911 1049.0969 3140.2967 5484.2789
[15] 5555.2178 2942.4670 4969.7218 4827.8439


$beta
[1] 0.0008608341


$tau
 [1] 2279.2446 3560.9346 3736.7832 2782.6046 2571.7861 2727.6519 3950.5992
 [8]  992.1459 1908.3572  947.1846 1171.9911 1049.0969 3140.2967 5484.2789
[15] 5555.2179 2942.4670 4969.7219 4827.8440


$beta
[1] 0.0008608178


$tau
 [1] 2280.2817 3560.9347 3736.7833 2782.6047 2571.7862 2727.6520 3950.5993
 [8]  992.1459 1908.3573  947.1846 1171.9911 1049.0969 3140.2968 5484.2790
[15] 5555.2180 2942.4671 4969.7220 4827.8441


$beta
[1] 0.0008608014


$tau
 [1] 2281.3179 3560.9347 3736.7834 2782.6047 2571.7862 2727.6520 3950.5993
 [8]  992.1460 1908.3573  947.1847 1171.9911 1049.0969 3140.2968 5484.2791
[15] 5555.2181 2942.4671 4969.7221 4827.8442


$beta
[1] 0.0008607851


$tau
 [1] 2282.3533 3560.9348 3736.7834 2782.6048 2571.7862 2727.6521 3950.5994
 [8]  992.1460 1908.3573  947.1847 1171.9912 1049.0970 3140.2969 5484.2792
[15] 5555.2182 2942.4672 4969.7221 4827.8443


$beta
[1] 0.0008607688


$tau
 [1] 2283.3878 3560.9349 3736.7835 2782.6048 2571.7863 2727.6521 3950.5994
 [8]  992.1460 1908.3573  947.1847 1171.9912 1049.0970 3140.2969 5484.2793
[15] 5555.2182 2942.4672 4969.7222 4827.8443


$beta
[1] 0.0008607525


$tau
 [1] 2284.4215 3560.9349 3736.7836 2782.6049 2571.7863 2727.6522 3950.5995
 [8]  992.1460 1908.3574  947.1847 1171.9912 1049.0970 3140.2970 5484.2794
[15] 5555.2183 2942.4673 4969.7223 4827.8444


$beta
[1] 0.0008607362


$tau
 [1] 2285.4543 3560.9350 3736.7836 2782.6049 2571.7864 2727.6522 3950.5996
 [8]  992.1460 1908.3574  947.1847 1171.9912 1049.0970 3140.2970 5484.2795
[15] 5555.2184 2942.4673 4969.7224 4827.8445


$beta
[1] 0.0008607199


$tau
 [1] 2286.4863 3560.9350 3736.7837 2782.6050 2571.7864 2727.6523 3950.5996
 [8]  992.1460 1908.3574  947.1847 1171.9912 1049.0970 3140.2971 5484.2796
[15] 5555.2185 2942.4674 4969.7225 4827.8446


$beta
[1] 0.0008607036


$tau
 [1] 2287.5174 3560.9351 3736.7837 2782.6050 2571.7865 2727.6523 3950.5997
 [8]  992.1461 1908.3575  947.1848 1171.9913 1049.0970 3140.2971 5484.2797
[15] 5555.2186 2942.4674 4969.7225 4827.8447


$beta
[1] 0.0008606874


$tau
 [1] 2288.5477 3560.9351 3736.7838 2782.6050 2571.7865 2727.6523 3950.5998
 [8]  992.1461 1908.3575  947.1848 1171.9913 1049.0971 3140.2972 5484.2797
[15] 5555.2187 2942.4675 4969.7226 4827.8447


$beta
[1] 0.0008606711


$tau
 [1] 2289.5771 3560.9352 3736.7839 2782.6051 2571.7865 2727.6524 3950.5998
 [8]  992.1461 1908.3575  947.1848 1171.9913 1049.0971 3140.2972 5484.2798
[15] 5555.2188 2942.4675 4969.7227 4827.8448


$beta
[1] 0.0008606549


$tau
 [1] 2290.6057 3560.9353 3736.7839 2782.6051 2571.7866 2727.6524 3950.5999
 [8]  992.1461 1908.3576  947.1848 1171.9913 1049.0971 3140.2973 5484.2799
[15] 5555.2189 2942.4676 4969.7228 4827.8449


$beta
[1] 0.0008606387


$tau
 [1] 2291.6334 3560.9353 3736.7840 2782.6052 2571.7866 2727.6525 3950.6000
 [8]  992.1461 1908.3576  947.1848 1171.9913 1049.0971 3140.2973 5484.2800
[15] 5555.2190 2942.4676 4969.7229 4827.8450


$beta
[1] 0.0008606225


$tau
 [1] 2292.6603 3560.9354 3736.7840 2782.6052 2571.7867 2727.6525 3950.6000
 [8]  992.1461 1908.3576  947.1848 1171.9914 1049.0971 3140.2974 5484.2801
[15] 5555.2191 2942.4676 4969.7229 4827.8450


$beta
[1] 0.0008606063


$tau
 [1] 2293.6864 3560.9354 3736.7841 2782.6053 2571.7867 2727.6526 3950.6001
 [8]  992.1462 1908.3577  947.1848 1171.9914 1049.0971 3140.2974 5484.2802
[15] 5555.2191 2942.4677 4969.7230 4827.8451


$beta
[1] 0.0008605901


$tau
 [1] 2294.7116 3560.9355 3736.7842 2782.6053 2571.7867 2727.6526 3950.6002
 [8]  992.1462 1908.3577  947.1849 1171.9914 1049.0972 3140.2975 5484.2803
[15] 5555.2192 2942.4677 4969.7231 4827.8452


$beta
[1] 0.000860574


$tau
 [1] 2295.7359 3560.9356 3736.7842 2782.6054 2571.7868 2727.6527 3950.6002
 [8]  992.1462 1908.3577  947.1849 1171.9914 1049.0972 3140.2975 5484.2804
[15] 5555.2193 2942.4678 4969.7232 4827.8453


$beta
[1] 0.0008605578


$tau
 [1] 2296.7594 3560.9356 3736.7843 2782.6054 2571.7868 2727.6527 3950.6003
 [8]  992.1462 1908.3577  947.1849 1171.9914 1049.0972 3140.2976 5484.2805
[15] 5555.2194 2942.4678 4969.7233 4827.8454


$beta
[1] 0.0008605417


$tau
 [1] 2297.7821 3560.9357 3736.7843 2782.6055 2571.7869 2727.6527 3950.6003
 [8]  992.1462 1908.3578  947.1849 1171.9915 1049.0972 3140.2976 5484.2805
[15] 5555.2195 2942.4679 4969.7233 4827.8454


$beta
[1] 0.0008605256


$tau
 [1] 2298.8039 3560.9357 3736.7844 2782.6055 2571.7869 2727.6528 3950.6004
 [8]  992.1462 1908.3578  947.1849 1171.9915 1049.0972 3140.2977 5484.2806
[15] 5555.2196 2942.4679 4969.7234 4827.8455


$beta
[1] 0.0008605095


$tau
 [1] 2299.8249 3560.9358 3736.7845 2782.6055 2571.7870 2727.6528 3950.6005
 [8]  992.1462 1908.3578  947.1849 1171.9915 1049.0972 3140.2977 5484.2807
[15] 5555.2197 2942.4680 4969.7235 4827.8456


$beta
[1] 0.0008604934


$tau
 [1] 2300.8450 3560.9358 3736.7845 2782.6056 2571.7870 2727.6529 3950.6005
 [8]  992.1463 1908.3579  947.1850 1171.9915 1049.0973 3140.2978 5484.2808
[15] 5555.2198 2942.4680 4969.7236 4827.8457


$beta
[1] 0.0008604774


$tau
 [1] 2301.8643 3560.9359 3736.7846 2782.6056 2571.7870 2727.6529 3950.6006
 [8]  992.1463 1908.3579  947.1850 1171.9915 1049.0973 3140.2978 5484.2809
[15] 5555.2199 2942.4681 4969.7237 4827.8457


$beta
[1] 0.0008604613


$tau
 [1] 2302.8828 3560.9360 3736.7846 2782.6057 2571.7871 2727.6530 3950.6007
 [8]  992.1463 1908.3579  947.1850 1171.9915 1049.0973 3140.2979 5484.2810
[15] 5555.2199 2942.4681 4969.7237 4827.8458


$beta
[1] 0.0008604452


$tau
 [1] 2303.9004 3560.9360 3736.7847 2782.6057 2571.7871 2727.6530 3950.6007
 [8]  992.1463 1908.3580  947.1850 1171.9916 1049.0973 3140.2979 5484.2811
[15] 5555.2200 2942.4682 4969.7238 4827.8459


$beta
[1] 0.0008604292


$tau
 [1] 2304.9172 3560.9361 3736.7848 2782.6058 2571.7872 2727.6530 3950.6008
 [8]  992.1463 1908.3580  947.1850 1171.9916 1049.0973 3140.2980 5484.2812
[15] 5555.2201 2942.4682 4969.7239 4827.8460


$beta
[1] 0.0008604132


$tau
 [1] 2305.9331 3560.9361 3736.7848 2782.6058 2571.7872 2727.6531 3950.6009
 [8]  992.1463 1908.3580  947.1850 1171.9916 1049.0973 3140.2980 5484.2813
[15] 5555.2202 2942.4683 4969.7240 4827.8461


$beta
[1] 0.0008603972


$tau
 [1] 2306.9482 3560.9362 3736.7849 2782.6059 2571.7872 2727.6531 3950.6009
 [8]  992.1464 1908.3581  947.1850 1171.9916 1049.0974 3140.2981 5484.2813
[15] 5555.2203 2942.4683 4969.7241 4827.8461


$beta
[1] 0.0008603812


$tau
 [1] 2307.9625 3560.9362 3736.7849 2782.6059 2571.7873 2727.6532 3950.6010
 [8]  992.1464 1908.3581  947.1851 1171.9916 1049.0974 3140.2981 5484.2814
[15] 5555.2204 2942.4684 4969.7241 4827.8462


$beta
[1] 0.0008603652


$tau
 [1] 2308.9759 3560.9363 3736.7850 2782.6059 2571.7873 2727.6532 3950.6010
 [8]  992.1464 1908.3581  947.1851 1171.9917 1049.0974 3140.2982 5484.2815
[15] 5555.2205 2942.4684 4969.7242 4827.8463


$beta
[1] 0.0008603493


$tau
 [1] 2309.9885 3560.9364 3736.7851 2782.6060 2571.7874 2727.6533 3950.6011
 [8]  992.1464 1908.3581  947.1851 1171.9917 1049.0974 3140.2982 5484.2816
[15] 5555.2206 2942.4684 4969.7243 4827.8464


$beta
[1] 0.0008603333


$tau
 [1] 2311.0003 3560.9364 3736.7851 2782.6060 2571.7874 2727.6533 3950.6012
 [8]  992.1464 1908.3582  947.1851 1171.9917 1049.0974 3140.2983 5484.2817
[15] 5555.2207 2942.4685 4969.7244 4827.8464


$beta
[1] 0.0008603174


$tau
 [1] 2312.0112 3560.9365 3736.7852 2782.6061 2571.7874 2727.6534 3950.6012
 [8]  992.1464 1908.3582  947.1851 1171.9917 1049.0974 3140.2983 5484.2818
[15] 5555.2207 2942.4685 4969.7245 4827.8465


$beta
[1] 0.0008603015


$tau
 [1] 2313.0213 3560.9365 3736.7852 2782.6061 2571.7875 2727.6534 3950.6013
 [8]  992.1465 1908.3582  947.1851 1171.9917 1049.0975 3140.2984 5484.2819
[15] 5555.2208 2942.4686 4969.7245 4827.8466


$beta
[1] 0.0008602856


$tau
 [1] 2314.0306 3560.9366 3736.7853 2782.6062 2571.7875 2727.6534 3950.6014
 [8]  992.1465 1908.3583  947.1852 1171.9918 1049.0975 3140.2984 5484.2820
[15] 5555.2209 2942.4686 4969.7246 4827.8467


$beta
[1] 0.0008602697


$tau
 [1] 2315.0391 3560.9366 3736.7854 2782.6062 2571.7876 2727.6535 3950.6014
 [8]  992.1465 1908.3583  947.1852 1171.9918 1049.0975 3140.2985 5484.2820
[15] 5555.2210 2942.4687 4969.7247 4827.8467


$beta
[1] 0.0008602538


$tau
 [1] 2316.0467 3560.9367 3736.7854 2782.6063 2571.7876 2727.6535 3950.6015
 [8]  992.1465 1908.3583  947.1852 1171.9918 1049.0975 3140.2985 5484.2821
[15] 5555.2211 2942.4687 4969.7248 4827.8468


$beta
[1] 0.0008602379


$tau
 [1] 2317.0535 3560.9367 3736.7855 2782.6063 2571.7877 2727.6536 3950.6015
 [8]  992.1465 1908.3584  947.1852 1171.9918 1049.0975 3140.2986 5484.2822
[15] 5555.2212 2942.4688 4969.7249 4827.8469


$beta
[1] 0.0008602221


$tau
 [1] 2318.0594 3560.9368 3736.7855 2782.6063 2571.7877 2727.6536 3950.6016
 [8]  992.1465 1908.3584  947.1852 1171.9918 1049.0975 3140.2986 5484.2823
[15] 5555.2213 2942.4688 4969.7249 4827.8470


$beta
[1] 0.0008602062


$tau
 [1] 2319.0645 3560.9369 3736.7856 2782.6064 2571.7877 2727.6537 3950.6017
 [8]  992.1465 1908.3584  947.1852 1171.9918 1049.0976 3140.2987 5484.2824
[15] 5555.2214 2942.4689 4969.7250 4827.8471


$beta
[1] 0.0008601904


$tau
 [1] 2320.0688 3560.9369 3736.7857 2782.6064 2571.7878 2727.6537 3950.6017
 [8]  992.1466 1908.3584  947.1852 1171.9919 1049.0976 3140.2987 5484.2825
[15] 5555.2215 2942.4689 4969.7251 4827.8471


$beta
[1] 0.0008601746


$tau
 [1] 2321.0723 3560.9370 3736.7857 2782.6065 2571.7878 2727.6537 3950.6018
 [8]  992.1466 1908.3585  947.1853 1171.9919 1049.0976 3140.2988 5484.2826
[15] 5555.2215 2942.4690 4969.7252 4827.8472


$beta
[1] 0.0008601588


$tau
 [1] 2322.0750 3560.9370 3736.7858 2782.6065 2571.7879 2727.6538 3950.6019
 [8]  992.1466 1908.3585  947.1853 1171.9919 1049.0976 3140.2988 5484.2826
[15] 5555.2216 2942.4690 4969.7252 4827.8473


$beta
[1] 0.000860143


$tau
 [1] 2323.0768 3560.9371 3736.7858 2782.6066 2571.7879 2727.6538 3950.6019
 [8]  992.1466 1908.3585  947.1853 1171.9919 1049.0976 3140.2989 5484.2827
[15] 5555.2217 2942.4691 4969.7253 4827.8474


$beta
[1] 0.0008601273


$tau
 [1] 2324.0778 3560.9371 3736.7859 2782.6066 2571.7879 2727.6539 3950.6020
 [8]  992.1466 1908.3586  947.1853 1171.9919 1049.0976 3140.2989 5484.2828
[15] 5555.2218 2942.4691 4969.7254 4827.8474


$beta
[1] 0.0008601115


$tau
 [1] 2325.0779 3560.9372 3736.7859 2782.6066 2571.7880 2727.6539 3950.6020
 [8]  992.1466 1908.3586  947.1853 1171.9920 1049.0977 3140.2990 5484.2829
[15] 5555.2219 2942.4691 4969.7255 4827.8475


$beta
[1] 0.0008600958


$tau
 [1] 2326.0773 3560.9373 3736.7860 2782.6067 2571.7880 2727.6540 3950.6021
 [8]  992.1467 1908.3586  947.1853 1171.9920 1049.0977 3140.2990 5484.2830
[15] 5555.2220 2942.4692 4969.7256 4827.8476
> round(param.temp$beta * param.temp$tau)
 [1] 2 3 3 2 2 2 3 1 2 1 1 1 3 5 5 3 4 4

MLE logic with x missing

> y <- c(3540, 3560, 3739, 2784, 2571, 2729, 3952, 993, 1908, 948,
+     1172, 1047, 3138, 5485, 5554, 2943, 4969, 4828)
> x <- c(3, 4, 1, 1, 3, 1, 2, 0, 2, 0, 1, 3, 5, 4, 6, 2, 5, 4)
> x.incomplete <- x[-1]
> beta.0 <- 3000
> tau.0 <- rep(0.01, 18)
> param <- list(beta = beta.0, tau = tau.0)
> em <- function(y, x, param) {
+     temp.beta <- param$beta
+     temp.tau <- param$tau
+     new.beta <- sum(y)/(temp.tau[1] + sum(x))
+     new.tau1 <- (temp.tau[1] + y[1])/(new.beta + 1)
+     new.tau2 <- (y[-1] + x)/(new.beta + 1)
+     new.tau <- c(new.tau1, new.tau2)
+     param.new <- list(beta = new.beta, tau = new.tau)
+     return(param.new)
+ }
> param.start <- param
> counter <- 1
> param.temp <- em(y, x.incomplete, param)
> while (abs(param.temp$beta - param.start$beta) > 1e-05) {
+     param.start <- param.temp
+     param.temp <- em(y, x.incomplete, param.temp)
+     counter <- counter + 1
+     print(param.temp)
+ }
$beta
[1] 1193.925


$tau
 [1] 2.9648604 2.9826131 3.1299027 2.3306896 2.1541095 2.2846616 3.3089934
 [8] 0.8310143 1.5984262 0.7933550 0.9816513 0.8787160 2.6302899 4.5935925
[15] 4.6530104 2.4645891 4.1626032 4.0437673


$beta
[1] 1189.4


$tau
 [1] 2.9762812 2.9939519 3.1418014 2.3395499 2.1622986 2.2933470 3.3215729
 [8] 0.8341735 1.6045028 0.7963710 0.9853832 0.8820565 2.6402892 4.6110555
[15] 4.6706993 2.4739586 4.1784278 4.0591401


$beta
[1] 1189.111


$tau
 [1] 2.9770140 2.9946793 3.1425647 2.3401184 2.1628240 2.2939042 3.3223799
 [8] 0.8343761 1.6048927 0.7965645 0.9856226 0.8822708 2.6409307 4.6121759
[15] 4.6718342 2.4745597 4.1794430 4.0601264


$beta
[1] 1189.092


$tau
 [1] 2.9770610 2.9947260 3.1426137 2.3401549 2.1628577 2.2939400 3.3224317
 [8] 0.8343891 1.6049177 0.7965770 0.9856379 0.8822846 2.6409719 4.6122478
[15] 4.6719070 2.4745982 4.1795082 4.0601897


$beta
[1] 1189.091


$tau
 [1] 2.9770640 2.9947290 3.1426168 2.3401572 2.1628598 2.2939422 3.3224350
 [8] 0.8343900 1.6049193 0.7965777 0.9856389 0.8822855 2.6409745 4.6122524
[15] 4.6719117 2.4746007 4.1795124 4.0601937


$beta
[1] 1189.091


$tau
 [1] 2.9770642 2.9947292 3.1426171 2.3401573 2.1628600 2.2939424 3.3224352
 [8] 0.8343900 1.6049194 0.7965778 0.9856390 0.8822855 2.6409747 4.6122527
[15] 4.6719120 2.4746009 4.1795126 4.0601940


$beta
[1] 1189.091


$tau
 [1] 2.9770642 2.9947292 3.1426171 2.3401574 2.1628600 2.2939424 3.3224353
 [8] 0.8343900 1.6049194 0.7965778 0.9856390 0.8822855 2.6409747 4.6122527
[15] 4.6719120 2.4746009 4.1795126 4.0601940
> k <- 1
> for (k in 1:18) {
+     while (abs(param.temp$tau[k] - param.start$tau[k]) > 0.01) {
+         param.start <- param.temp
+         param.temp <- em(y, x.incomplete, param.temp)
+         print(param.temp)
+     }
+ }
> round(param.temp$beta * param.temp$tau)
 [1] 3540 3561 3737 2783 2572 2728 3951  992 1908  947 1172 1049 3140 5484 5555
[16] 2943 4970 4828