Ilustrasi: data menyebar Binomial dengan prior sebaran Beta(a,b)

library(Bolstad)
simulasi.post.pi <- binogcp(7, 20, density = "beta", params = c(2, 5))

## hasil simulasi berdasarkan MCMC 
simulasi.post.pi
## $name
## [1] "pi"
## 
## $param.x
##    [1] 0.0005 0.0015 0.0025 0.0035 0.0045 0.0055 0.0065 0.0075 0.0085 0.0095
##   [11] 0.0105 0.0115 0.0125 0.0135 0.0145 0.0155 0.0165 0.0175 0.0185 0.0195
##   [21] 0.0205 0.0215 0.0225 0.0235 0.0245 0.0255 0.0265 0.0275 0.0285 0.0295
##   [31] 0.0305 0.0315 0.0325 0.0335 0.0345 0.0355 0.0365 0.0375 0.0385 0.0395
##   [41] 0.0405 0.0415 0.0425 0.0435 0.0445 0.0455 0.0465 0.0475 0.0485 0.0495
##   [51] 0.0505 0.0515 0.0525 0.0535 0.0545 0.0555 0.0565 0.0575 0.0585 0.0595
##   [61] 0.0605 0.0615 0.0625 0.0635 0.0645 0.0655 0.0665 0.0675 0.0685 0.0695
##   [71] 0.0705 0.0715 0.0725 0.0735 0.0745 0.0755 0.0765 0.0775 0.0785 0.0795
##   [81] 0.0805 0.0815 0.0825 0.0835 0.0845 0.0855 0.0865 0.0875 0.0885 0.0895
##   [91] 0.0905 0.0915 0.0925 0.0935 0.0945 0.0955 0.0965 0.0975 0.0985 0.0995
##  [101] 0.1005 0.1015 0.1025 0.1035 0.1045 0.1055 0.1065 0.1075 0.1085 0.1095
##  [111] 0.1105 0.1115 0.1125 0.1135 0.1145 0.1155 0.1165 0.1175 0.1185 0.1195
##  [121] 0.1205 0.1215 0.1225 0.1235 0.1245 0.1255 0.1265 0.1275 0.1285 0.1295
##  [131] 0.1305 0.1315 0.1325 0.1335 0.1345 0.1355 0.1365 0.1375 0.1385 0.1395
##  [141] 0.1405 0.1415 0.1425 0.1435 0.1445 0.1455 0.1465 0.1475 0.1485 0.1495
##  [151] 0.1505 0.1515 0.1525 0.1535 0.1545 0.1555 0.1565 0.1575 0.1585 0.1595
##  [161] 0.1605 0.1615 0.1625 0.1635 0.1645 0.1655 0.1665 0.1675 0.1685 0.1695
##  [171] 0.1705 0.1715 0.1725 0.1735 0.1745 0.1755 0.1765 0.1775 0.1785 0.1795
##  [181] 0.1805 0.1815 0.1825 0.1835 0.1845 0.1855 0.1865 0.1875 0.1885 0.1895
##  [191] 0.1905 0.1915 0.1925 0.1935 0.1945 0.1955 0.1965 0.1975 0.1985 0.1995
##  [201] 0.2005 0.2015 0.2025 0.2035 0.2045 0.2055 0.2065 0.2075 0.2085 0.2095
##  [211] 0.2105 0.2115 0.2125 0.2135 0.2145 0.2155 0.2165 0.2175 0.2185 0.2195
##  [221] 0.2205 0.2215 0.2225 0.2235 0.2245 0.2255 0.2265 0.2275 0.2285 0.2295
##  [231] 0.2305 0.2315 0.2325 0.2335 0.2345 0.2355 0.2365 0.2375 0.2385 0.2395
##  [241] 0.2405 0.2415 0.2425 0.2435 0.2445 0.2455 0.2465 0.2475 0.2485 0.2495
##  [251] 0.2505 0.2515 0.2525 0.2535 0.2545 0.2555 0.2565 0.2575 0.2585 0.2595
##  [261] 0.2605 0.2615 0.2625 0.2635 0.2645 0.2655 0.2665 0.2675 0.2685 0.2695
##  [271] 0.2705 0.2715 0.2725 0.2735 0.2745 0.2755 0.2765 0.2775 0.2785 0.2795
##  [281] 0.2805 0.2815 0.2825 0.2835 0.2845 0.2855 0.2865 0.2875 0.2885 0.2895
##  [291] 0.2905 0.2915 0.2925 0.2935 0.2945 0.2955 0.2965 0.2975 0.2985 0.2995
##  [301] 0.3005 0.3015 0.3025 0.3035 0.3045 0.3055 0.3065 0.3075 0.3085 0.3095
##  [311] 0.3105 0.3115 0.3125 0.3135 0.3145 0.3155 0.3165 0.3175 0.3185 0.3195
##  [321] 0.3205 0.3215 0.3225 0.3235 0.3245 0.3255 0.3265 0.3275 0.3285 0.3295
##  [331] 0.3305 0.3315 0.3325 0.3335 0.3345 0.3355 0.3365 0.3375 0.3385 0.3395
##  [341] 0.3405 0.3415 0.3425 0.3435 0.3445 0.3455 0.3465 0.3475 0.3485 0.3495
##  [351] 0.3505 0.3515 0.3525 0.3535 0.3545 0.3555 0.3565 0.3575 0.3585 0.3595
##  [361] 0.3605 0.3615 0.3625 0.3635 0.3645 0.3655 0.3665 0.3675 0.3685 0.3695
##  [371] 0.3705 0.3715 0.3725 0.3735 0.3745 0.3755 0.3765 0.3775 0.3785 0.3795
##  [381] 0.3805 0.3815 0.3825 0.3835 0.3845 0.3855 0.3865 0.3875 0.3885 0.3895
##  [391] 0.3905 0.3915 0.3925 0.3935 0.3945 0.3955 0.3965 0.3975 0.3985 0.3995
##  [401] 0.4005 0.4015 0.4025 0.4035 0.4045 0.4055 0.4065 0.4075 0.4085 0.4095
##  [411] 0.4105 0.4115 0.4125 0.4135 0.4145 0.4155 0.4165 0.4175 0.4185 0.4195
##  [421] 0.4205 0.4215 0.4225 0.4235 0.4245 0.4255 0.4265 0.4275 0.4285 0.4295
##  [431] 0.4305 0.4315 0.4325 0.4335 0.4345 0.4355 0.4365 0.4375 0.4385 0.4395
##  [441] 0.4405 0.4415 0.4425 0.4435 0.4445 0.4455 0.4465 0.4475 0.4485 0.4495
##  [451] 0.4505 0.4515 0.4525 0.4535 0.4545 0.4555 0.4565 0.4575 0.4585 0.4595
##  [461] 0.4605 0.4615 0.4625 0.4635 0.4645 0.4655 0.4665 0.4675 0.4685 0.4695
##  [471] 0.4705 0.4715 0.4725 0.4735 0.4745 0.4755 0.4765 0.4775 0.4785 0.4795
##  [481] 0.4805 0.4815 0.4825 0.4835 0.4845 0.4855 0.4865 0.4875 0.4885 0.4895
##  [491] 0.4905 0.4915 0.4925 0.4935 0.4945 0.4955 0.4965 0.4975 0.4985 0.4995
##  [501] 0.5005 0.5015 0.5025 0.5035 0.5045 0.5055 0.5065 0.5075 0.5085 0.5095
##  [511] 0.5105 0.5115 0.5125 0.5135 0.5145 0.5155 0.5165 0.5175 0.5185 0.5195
##  [521] 0.5205 0.5215 0.5225 0.5235 0.5245 0.5255 0.5265 0.5275 0.5285 0.5295
##  [531] 0.5305 0.5315 0.5325 0.5335 0.5345 0.5355 0.5365 0.5375 0.5385 0.5395
##  [541] 0.5405 0.5415 0.5425 0.5435 0.5445 0.5455 0.5465 0.5475 0.5485 0.5495
##  [551] 0.5505 0.5515 0.5525 0.5535 0.5545 0.5555 0.5565 0.5575 0.5585 0.5595
##  [561] 0.5605 0.5615 0.5625 0.5635 0.5645 0.5655 0.5665 0.5675 0.5685 0.5695
##  [571] 0.5705 0.5715 0.5725 0.5735 0.5745 0.5755 0.5765 0.5775 0.5785 0.5795
##  [581] 0.5805 0.5815 0.5825 0.5835 0.5845 0.5855 0.5865 0.5875 0.5885 0.5895
##  [591] 0.5905 0.5915 0.5925 0.5935 0.5945 0.5955 0.5965 0.5975 0.5985 0.5995
##  [601] 0.6005 0.6015 0.6025 0.6035 0.6045 0.6055 0.6065 0.6075 0.6085 0.6095
##  [611] 0.6105 0.6115 0.6125 0.6135 0.6145 0.6155 0.6165 0.6175 0.6185 0.6195
##  [621] 0.6205 0.6215 0.6225 0.6235 0.6245 0.6255 0.6265 0.6275 0.6285 0.6295
##  [631] 0.6305 0.6315 0.6325 0.6335 0.6345 0.6355 0.6365 0.6375 0.6385 0.6395
##  [641] 0.6405 0.6415 0.6425 0.6435 0.6445 0.6455 0.6465 0.6475 0.6485 0.6495
##  [651] 0.6505 0.6515 0.6525 0.6535 0.6545 0.6555 0.6565 0.6575 0.6585 0.6595
##  [661] 0.6605 0.6615 0.6625 0.6635 0.6645 0.6655 0.6665 0.6675 0.6685 0.6695
##  [671] 0.6705 0.6715 0.6725 0.6735 0.6745 0.6755 0.6765 0.6775 0.6785 0.6795
##  [681] 0.6805 0.6815 0.6825 0.6835 0.6845 0.6855 0.6865 0.6875 0.6885 0.6895
##  [691] 0.6905 0.6915 0.6925 0.6935 0.6945 0.6955 0.6965 0.6975 0.6985 0.6995
##  [701] 0.7005 0.7015 0.7025 0.7035 0.7045 0.7055 0.7065 0.7075 0.7085 0.7095
##  [711] 0.7105 0.7115 0.7125 0.7135 0.7145 0.7155 0.7165 0.7175 0.7185 0.7195
##  [721] 0.7205 0.7215 0.7225 0.7235 0.7245 0.7255 0.7265 0.7275 0.7285 0.7295
##  [731] 0.7305 0.7315 0.7325 0.7335 0.7345 0.7355 0.7365 0.7375 0.7385 0.7395
##  [741] 0.7405 0.7415 0.7425 0.7435 0.7445 0.7455 0.7465 0.7475 0.7485 0.7495
##  [751] 0.7505 0.7515 0.7525 0.7535 0.7545 0.7555 0.7565 0.7575 0.7585 0.7595
##  [761] 0.7605 0.7615 0.7625 0.7635 0.7645 0.7655 0.7665 0.7675 0.7685 0.7695
##  [771] 0.7705 0.7715 0.7725 0.7735 0.7745 0.7755 0.7765 0.7775 0.7785 0.7795
##  [781] 0.7805 0.7815 0.7825 0.7835 0.7845 0.7855 0.7865 0.7875 0.7885 0.7895
##  [791] 0.7905 0.7915 0.7925 0.7935 0.7945 0.7955 0.7965 0.7975 0.7985 0.7995
##  [801] 0.8005 0.8015 0.8025 0.8035 0.8045 0.8055 0.8065 0.8075 0.8085 0.8095
##  [811] 0.8105 0.8115 0.8125 0.8135 0.8145 0.8155 0.8165 0.8175 0.8185 0.8195
##  [821] 0.8205 0.8215 0.8225 0.8235 0.8245 0.8255 0.8265 0.8275 0.8285 0.8295
##  [831] 0.8305 0.8315 0.8325 0.8335 0.8345 0.8355 0.8365 0.8375 0.8385 0.8395
##  [841] 0.8405 0.8415 0.8425 0.8435 0.8445 0.8455 0.8465 0.8475 0.8485 0.8495
##  [851] 0.8505 0.8515 0.8525 0.8535 0.8545 0.8555 0.8565 0.8575 0.8585 0.8595
##  [861] 0.8605 0.8615 0.8625 0.8635 0.8645 0.8655 0.8665 0.8675 0.8685 0.8695
##  [871] 0.8705 0.8715 0.8725 0.8735 0.8745 0.8755 0.8765 0.8775 0.8785 0.8795
##  [881] 0.8805 0.8815 0.8825 0.8835 0.8845 0.8855 0.8865 0.8875 0.8885 0.8895
##  [891] 0.8905 0.8915 0.8925 0.8935 0.8945 0.8955 0.8965 0.8975 0.8985 0.8995
##  [901] 0.9005 0.9015 0.9025 0.9035 0.9045 0.9055 0.9065 0.9075 0.9085 0.9095
##  [911] 0.9105 0.9115 0.9125 0.9135 0.9145 0.9155 0.9165 0.9175 0.9185 0.9195
##  [921] 0.9205 0.9215 0.9225 0.9235 0.9245 0.9255 0.9265 0.9275 0.9285 0.9295
##  [931] 0.9305 0.9315 0.9325 0.9335 0.9345 0.9355 0.9365 0.9375 0.9385 0.9395
##  [941] 0.9405 0.9415 0.9425 0.9435 0.9445 0.9455 0.9465 0.9475 0.9485 0.9495
##  [951] 0.9505 0.9515 0.9525 0.9535 0.9545 0.9555 0.9565 0.9575 0.9585 0.9595
##  [961] 0.9605 0.9615 0.9625 0.9635 0.9645 0.9655 0.9665 0.9675 0.9685 0.9695
##  [971] 0.9705 0.9715 0.9725 0.9735 0.9745 0.9755 0.9765 0.9775 0.9785 0.9795
##  [981] 0.9805 0.9815 0.9825 0.9835 0.9845 0.9855 0.9865 0.9875 0.9885 0.9895
##  [991] 0.9905 0.9915 0.9925 0.9935 0.9945 0.9955 0.9965 0.9975 0.9985 0.9995
## 
## $prior
##    [1] 1.497002e-02 4.473061e-02 7.425281e-02 1.035377e-01 1.325864e-01
##    [6] 1.613998e-01 1.899792e-01 2.183256e-01 2.464399e-01 2.743234e-01
##   [11] 3.019769e-01 3.294017e-01 3.565986e-01 3.835689e-01 4.103135e-01
##   [16] 4.368334e-01 4.631297e-01 4.892035e-01 5.150557e-01 5.406874e-01
##   [21] 5.660996e-01 5.912934e-01 6.162697e-01 6.410296e-01 6.655741e-01
##   [26] 6.899042e-01 7.140209e-01 7.379253e-01 7.616182e-01 7.851008e-01
##   [31] 8.083740e-01 8.314388e-01 8.542963e-01 8.769473e-01 8.993929e-01
##   [36] 9.216341e-01 9.436718e-01 9.655071e-01 9.871409e-01 1.008574e+00
##   [41] 1.029808e+00 1.050843e+00 1.071680e+00 1.092321e+00 1.112766e+00
##   [46] 1.133017e+00 1.153074e+00 1.172937e+00 1.192609e+00 1.212090e+00
##   [51] 1.231381e+00 1.250483e+00 1.269397e+00 1.288124e+00 1.306664e+00
##   [56] 1.325019e+00 1.343190e+00 1.361177e+00 1.378981e+00 1.396604e+00
##   [61] 1.414047e+00 1.431309e+00 1.448393e+00 1.465299e+00 1.482027e+00
##   [66] 1.498579e+00 1.514957e+00 1.531159e+00 1.547189e+00 1.563045e+00
##   [71] 1.578730e+00 1.594245e+00 1.609589e+00 1.624764e+00 1.639771e+00
##   [76] 1.654611e+00 1.669284e+00 1.683792e+00 1.698135e+00 1.712315e+00
##   [81] 1.726331e+00 1.740185e+00 1.753878e+00 1.767411e+00 1.780784e+00
##   [86] 1.793999e+00 1.807056e+00 1.819956e+00 1.832699e+00 1.845288e+00
##   [91] 1.857722e+00 1.870002e+00 1.882130e+00 1.894105e+00 1.905930e+00
##   [96] 1.917604e+00 1.929129e+00 1.940505e+00 1.951733e+00 1.962814e+00
##  [101] 1.973749e+00 1.984539e+00 1.995184e+00 2.005685e+00 2.016044e+00
##  [106] 2.026260e+00 2.036334e+00 2.046268e+00 2.056063e+00 2.065718e+00
##  [111] 2.075235e+00 2.084615e+00 2.093858e+00 2.102965e+00 2.111937e+00
##  [116] 2.120775e+00 2.129479e+00 2.138051e+00 2.146490e+00 2.154798e+00
##  [121] 2.162976e+00 2.171024e+00 2.178943e+00 2.186734e+00 2.194397e+00
##  [126] 2.201934e+00 2.209344e+00 2.216630e+00 2.223791e+00 2.230828e+00
##  [131] 2.237743e+00 2.244535e+00 2.251205e+00 2.257755e+00 2.264185e+00
##  [136] 2.270495e+00 2.276687e+00 2.282761e+00 2.288717e+00 2.294558e+00
##  [141] 2.300282e+00 2.305892e+00 2.311387e+00 2.316768e+00 2.322037e+00
##  [146] 2.327194e+00 2.332239e+00 2.337173e+00 2.341997e+00 2.346712e+00
##  [151] 2.351318e+00 2.355815e+00 2.360206e+00 2.364490e+00 2.368668e+00
##  [156] 2.372741e+00 2.376709e+00 2.380573e+00 2.384334e+00 2.387992e+00
##  [161] 2.391548e+00 2.395003e+00 2.398357e+00 2.401612e+00 2.404767e+00
##  [166] 2.407823e+00 2.410782e+00 2.413643e+00 2.416407e+00 2.419076e+00
##  [171] 2.421649e+00 2.424127e+00 2.426512e+00 2.428802e+00 2.431000e+00
##  [176] 2.433106e+00 2.435120e+00 2.437043e+00 2.438876e+00 2.440619e+00
##  [181] 2.442273e+00 2.443839e+00 2.445317e+00 2.446708e+00 2.448012e+00
##  [186] 2.449230e+00 2.450362e+00 2.451410e+00 2.452374e+00 2.453254e+00
##  [191] 2.454051e+00 2.454766e+00 2.455399e+00 2.455951e+00 2.456422e+00
##  [196] 2.456814e+00 2.457125e+00 2.457358e+00 2.457513e+00 2.457590e+00
##  [201] 2.457590e+00 2.457514e+00 2.457361e+00 2.457134e+00 2.456831e+00
##  [206] 2.456454e+00 2.456004e+00 2.455480e+00 2.454884e+00 2.454216e+00
##  [211] 2.453476e+00 2.452666e+00 2.451785e+00 2.450835e+00 2.449815e+00
##  [216] 2.448727e+00 2.447570e+00 2.446346e+00 2.445055e+00 2.443697e+00
##  [221] 2.442274e+00 2.440785e+00 2.439231e+00 2.437612e+00 2.435930e+00
##  [226] 2.434184e+00 2.432376e+00 2.430505e+00 2.428573e+00 2.426579e+00
##  [231] 2.424525e+00 2.422410e+00 2.420236e+00 2.418003e+00 2.415710e+00
##  [236] 2.413360e+00 2.410952e+00 2.408487e+00 2.405965e+00 2.403387e+00
##  [241] 2.400753e+00 2.398064e+00 2.395320e+00 2.392522e+00 2.389670e+00
##  [246] 2.386765e+00 2.383808e+00 2.380798e+00 2.377736e+00 2.374623e+00
##  [251] 2.371459e+00 2.368244e+00 2.364980e+00 2.361666e+00 2.358303e+00
##  [256] 2.354892e+00 2.351432e+00 2.347926e+00 2.344371e+00 2.340771e+00
##  [261] 2.337124e+00 2.333431e+00 2.329693e+00 2.325910e+00 2.322082e+00
##  [266] 2.318211e+00 2.314296e+00 2.310338e+00 2.306338e+00 2.302295e+00
##  [271] 2.298210e+00 2.294084e+00 2.289917e+00 2.285710e+00 2.281462e+00
##  [276] 2.277175e+00 2.272849e+00 2.268484e+00 2.264080e+00 2.259639e+00
##  [281] 2.255160e+00 2.250644e+00 2.246091e+00 2.241502e+00 2.236877e+00
##  [286] 2.232216e+00 2.227521e+00 2.222791e+00 2.218026e+00 2.213228e+00
##  [291] 2.208396e+00 2.203531e+00 2.198634e+00 2.193704e+00 2.188742e+00
##  [296] 2.183749e+00 2.178725e+00 2.173670e+00 2.168584e+00 2.163469e+00
##  [301] 2.158324e+00 2.153150e+00 2.147947e+00 2.142715e+00 2.137456e+00
##  [306] 2.132168e+00 2.126854e+00 2.121512e+00 2.116144e+00 2.110749e+00
##  [311] 2.105329e+00 2.099883e+00 2.094412e+00 2.088916e+00 2.083396e+00
##  [316] 2.077851e+00 2.072283e+00 2.066691e+00 2.061077e+00 2.055439e+00
##  [321] 2.049779e+00 2.044098e+00 2.038394e+00 2.032669e+00 2.026924e+00
##  [326] 2.021157e+00 2.015370e+00 2.009563e+00 2.003737e+00 1.997891e+00
##  [331] 1.992026e+00 1.986143e+00 1.980241e+00 1.974321e+00 1.968383e+00
##  [336] 1.962428e+00 1.956456e+00 1.950467e+00 1.944461e+00 1.938440e+00
##  [341] 1.932402e+00 1.926349e+00 1.920281e+00 1.914198e+00 1.908101e+00
##  [346] 1.901989e+00 1.895863e+00 1.889723e+00 1.883570e+00 1.877404e+00
##  [351] 1.871225e+00 1.865033e+00 1.858829e+00 1.852613e+00 1.846386e+00
##  [356] 1.840147e+00 1.833897e+00 1.827637e+00 1.821365e+00 1.815084e+00
##  [361] 1.808792e+00 1.802491e+00 1.796181e+00 1.789861e+00 1.783532e+00
##  [366] 1.777195e+00 1.770850e+00 1.764496e+00 1.758135e+00 1.751766e+00
##  [371] 1.745390e+00 1.739006e+00 1.732617e+00 1.726220e+00 1.719817e+00
##  [376] 1.713409e+00 1.706994e+00 1.700574e+00 1.694149e+00 1.687719e+00
##  [381] 1.681284e+00 1.674845e+00 1.668401e+00 1.661954e+00 1.655502e+00
##  [386] 1.649048e+00 1.642589e+00 1.636128e+00 1.629664e+00 1.623198e+00
##  [391] 1.616729e+00 1.610258e+00 1.603785e+00 1.597310e+00 1.590834e+00
##  [396] 1.584357e+00 1.577878e+00 1.571399e+00 1.564920e+00 1.558440e+00
##  [401] 1.551960e+00 1.545480e+00 1.539001e+00 1.532522e+00 1.526043e+00
##  [406] 1.519566e+00 1.513090e+00 1.506615e+00 1.500142e+00 1.493670e+00
##  [411] 1.487201e+00 1.480734e+00 1.474269e+00 1.467807e+00 1.461347e+00
##  [416] 1.454891e+00 1.448437e+00 1.441987e+00 1.435541e+00 1.429098e+00
##  [421] 1.422660e+00 1.416225e+00 1.409795e+00 1.403369e+00 1.396948e+00
##  [426] 1.390532e+00 1.384121e+00 1.377715e+00 1.371314e+00 1.364919e+00
##  [431] 1.358530e+00 1.352147e+00 1.345770e+00 1.339399e+00 1.333035e+00
##  [436] 1.326677e+00 1.320326e+00 1.313982e+00 1.307645e+00 1.301315e+00
##  [441] 1.294993e+00 1.288679e+00 1.282372e+00 1.276073e+00 1.269782e+00
##  [446] 1.263500e+00 1.257226e+00 1.250960e+00 1.244703e+00 1.238455e+00
##  [451] 1.232216e+00 1.225986e+00 1.219765e+00 1.213554e+00 1.207353e+00
##  [456] 1.201161e+00 1.194979e+00 1.188807e+00 1.182645e+00 1.176494e+00
##  [461] 1.170352e+00 1.164222e+00 1.158102e+00 1.151993e+00 1.145895e+00
##  [466] 1.139808e+00 1.133733e+00 1.127668e+00 1.121615e+00 1.115574e+00
##  [471] 1.109545e+00 1.103527e+00 1.097521e+00 1.091528e+00 1.085546e+00
##  [476] 1.079577e+00 1.073621e+00 1.067677e+00 1.061746e+00 1.055827e+00
##  [481] 1.049922e+00 1.044029e+00 1.038150e+00 1.032284e+00 1.026431e+00
##  [486] 1.020592e+00 1.014766e+00 1.008954e+00 1.003156e+00 9.973715e-01
##  [491] 9.916011e-01 9.858449e-01 9.801029e-01 9.743752e-01 9.686618e-01
##  [496] 9.629630e-01 9.572787e-01 9.516091e-01 9.459543e-01 9.403144e-01
##  [501] 9.346894e-01 9.290794e-01 9.234846e-01 9.179050e-01 9.123407e-01
##  [506] 9.067918e-01 9.012584e-01 8.957406e-01 8.902384e-01 8.847519e-01
##  [511] 8.792812e-01 8.738264e-01 8.683876e-01 8.629648e-01 8.575581e-01
##  [516] 8.521677e-01 8.467935e-01 8.414356e-01 8.360941e-01 8.307692e-01
##  [521] 8.254608e-01 8.201690e-01 8.148939e-01 8.096355e-01 8.043940e-01
##  [526] 7.991694e-01 7.939617e-01 7.887711e-01 7.835975e-01 7.784411e-01
##  [531] 7.733018e-01 7.681798e-01 7.630752e-01 7.579879e-01 7.529180e-01
##  [536] 7.478657e-01 7.428309e-01 7.378136e-01 7.328140e-01 7.278322e-01
##  [541] 7.228680e-01 7.179217e-01 7.129932e-01 7.080826e-01 7.031900e-01
##  [546] 6.983153e-01 6.934587e-01 6.886201e-01 6.837997e-01 6.789974e-01
##  [551] 6.742134e-01 6.694476e-01 6.647000e-01 6.599708e-01 6.552600e-01
##  [556] 6.505675e-01 6.458935e-01 6.412380e-01 6.366009e-01 6.319824e-01
##  [561] 6.273825e-01 6.228012e-01 6.182384e-01 6.136944e-01 6.091691e-01
##  [566] 6.046624e-01 6.001746e-01 5.957054e-01 5.912551e-01 5.868236e-01
##  [571] 5.824110e-01 5.780172e-01 5.736424e-01 5.692864e-01 5.649494e-01
##  [576] 5.606313e-01 5.563322e-01 5.520521e-01 5.477911e-01 5.435490e-01
##  [581] 5.393260e-01 5.351220e-01 5.309371e-01 5.267713e-01 5.226246e-01
##  [586] 5.184971e-01 5.143886e-01 5.102993e-01 5.062291e-01 5.021780e-01
##  [591] 4.981462e-01 4.941334e-01 4.901399e-01 4.861655e-01 4.822104e-01
##  [596] 4.782744e-01 4.743576e-01 4.704600e-01 4.665816e-01 4.627224e-01
##  [601] 4.588824e-01 4.550616e-01 4.512600e-01 4.474776e-01 4.437144e-01
##  [606] 4.399704e-01 4.362455e-01 4.325399e-01 4.288534e-01 4.251862e-01
##  [611] 4.215380e-01 4.179091e-01 4.142993e-01 4.107086e-01 4.071371e-01
##  [616] 4.035847e-01 4.000514e-01 3.965372e-01 3.930422e-01 3.895661e-01
##  [621] 3.861092e-01 3.826713e-01 3.792525e-01 3.758527e-01 3.724719e-01
##  [626] 3.691100e-01 3.657672e-01 3.624433e-01 3.591384e-01 3.558524e-01
##  [631] 3.525852e-01 3.493370e-01 3.461076e-01 3.428971e-01 3.397054e-01
##  [636] 3.365324e-01 3.333783e-01 3.302429e-01 3.271262e-01 3.240283e-01
##  [641] 3.209490e-01 3.178883e-01 3.148463e-01 3.118229e-01 3.088180e-01
##  [646] 3.058317e-01 3.028639e-01 2.999145e-01 2.969837e-01 2.940712e-01
##  [651] 2.911771e-01 2.883014e-01 2.854441e-01 2.826050e-01 2.797842e-01
##  [656] 2.769816e-01 2.741972e-01 2.714309e-01 2.686828e-01 2.659528e-01
##  [661] 2.632408e-01 2.605468e-01 2.578708e-01 2.552128e-01 2.525727e-01
##  [666] 2.499504e-01 2.473459e-01 2.447593e-01 2.421903e-01 2.396391e-01
##  [671] 2.371055e-01 2.345896e-01 2.320912e-01 2.296104e-01 2.271471e-01
##  [676] 2.247012e-01 2.222727e-01 2.198616e-01 2.174678e-01 2.150913e-01
##  [681] 2.127320e-01 2.103899e-01 2.080649e-01 2.057570e-01 2.034661e-01
##  [686] 2.011923e-01 1.989354e-01 1.966953e-01 1.944721e-01 1.922658e-01
##  [691] 1.900761e-01 1.879032e-01 1.857469e-01 1.836072e-01 1.814840e-01
##  [696] 1.793773e-01 1.772871e-01 1.752133e-01 1.731558e-01 1.711145e-01
##  [701] 1.690895e-01 1.670807e-01 1.650880e-01 1.631113e-01 1.611507e-01
##  [706] 1.592060e-01 1.572772e-01 1.553643e-01 1.534671e-01 1.515857e-01
##  [711] 1.497199e-01 1.478698e-01 1.460352e-01 1.442161e-01 1.424124e-01
##  [716] 1.406242e-01 1.388512e-01 1.370935e-01 1.353510e-01 1.336237e-01
##  [721] 1.319114e-01 1.302142e-01 1.285319e-01 1.268646e-01 1.252120e-01
##  [726] 1.235743e-01 1.219512e-01 1.203428e-01 1.187491e-01 1.171698e-01
##  [731] 1.156050e-01 1.140546e-01 1.125185e-01 1.109967e-01 1.094892e-01
##  [736] 1.079957e-01 1.065164e-01 1.050511e-01 1.035997e-01 1.021622e-01
##  [741] 1.007386e-01 9.932867e-02 9.793246e-02 9.654987e-02 9.518085e-02
##  [746] 9.382532e-02 9.248322e-02 9.115447e-02 8.983902e-02 8.853680e-02
##  [751] 8.724773e-02 8.597176e-02 8.470881e-02 8.345881e-02 8.222171e-02
##  [756] 8.099742e-02 7.978588e-02 7.858703e-02 7.740079e-02 7.622709e-02
##  [761] 7.506587e-02 7.391706e-02 7.278059e-02 7.165638e-02 7.054438e-02
##  [766] 6.944450e-02 6.835668e-02 6.728085e-02 6.621694e-02 6.516488e-02
##  [771] 6.412460e-02 6.309603e-02 6.207909e-02 6.107372e-02 6.007985e-02
##  [776] 5.909741e-02 5.812632e-02 5.716651e-02 5.621792e-02 5.528047e-02
##  [781] 5.435409e-02 5.343871e-02 5.253426e-02 5.164067e-02 5.075786e-02
##  [786] 4.988576e-02 4.902431e-02 4.817343e-02 4.733304e-02 4.650309e-02
##  [791] 4.568349e-02 4.487417e-02 4.407506e-02 4.328610e-02 4.250720e-02
##  [796] 4.173829e-02 4.097931e-02 4.023019e-02 3.949084e-02 3.876120e-02
##  [801] 3.804120e-02 3.733076e-02 3.662981e-02 3.593829e-02 3.525611e-02
##  [806] 3.458320e-02 3.391950e-02 3.326494e-02 3.261943e-02 3.198291e-02
##  [811] 3.135531e-02 3.073656e-02 3.012657e-02 2.952529e-02 2.893264e-02
##  [816] 2.834854e-02 2.777293e-02 2.720574e-02 2.664689e-02 2.609631e-02
##  [821] 2.555393e-02 2.501968e-02 2.449349e-02 2.397528e-02 2.346499e-02
##  [826] 2.296255e-02 2.246788e-02 2.198091e-02 2.150157e-02 2.102980e-02
##  [831] 2.056552e-02 2.010866e-02 1.965915e-02 1.921692e-02 1.878190e-02
##  [836] 1.835402e-02 1.793322e-02 1.751941e-02 1.711254e-02 1.671253e-02
##  [841] 1.631931e-02 1.593281e-02 1.555297e-02 1.517972e-02 1.481298e-02
##  [846] 1.445269e-02 1.409879e-02 1.375119e-02 1.340984e-02 1.307467e-02
##  [851] 1.274560e-02 1.242258e-02 1.210553e-02 1.179438e-02 1.148908e-02
##  [856] 1.118955e-02 1.089573e-02 1.060754e-02 1.032493e-02 1.004783e-02
##  [861] 9.776175e-03 9.509893e-03 9.248923e-03 8.993199e-03 8.742657e-03
##  [866] 8.497232e-03 8.256859e-03 8.021475e-03 7.791016e-03 7.565417e-03
##  [871] 7.344616e-03 7.128549e-03 6.917154e-03 6.710367e-03 6.508126e-03
##  [876] 6.310370e-03 6.117036e-03 5.928062e-03 5.743388e-03 5.562952e-03
##  [881] 5.386693e-03 5.214552e-03 5.046467e-03 4.882379e-03 4.722229e-03
##  [886] 4.565956e-03 4.413503e-03 4.264810e-03 4.119819e-03 3.978472e-03
##  [891] 3.840711e-03 3.706479e-03 3.575719e-03 3.448373e-03 3.324386e-03
##  [896] 3.203701e-03 3.086263e-03 2.972016e-03 2.860905e-03 2.752876e-03
##  [901] 2.647874e-03 2.545845e-03 2.446735e-03 2.350491e-03 2.257061e-03
##  [906] 2.166392e-03 2.078430e-03 1.993126e-03 1.910428e-03 1.830283e-03
##  [911] 1.752643e-03 1.677456e-03 1.604672e-03 1.534243e-03 1.466120e-03
##  [916] 1.400253e-03 1.336594e-03 1.275097e-03 1.215712e-03 1.158395e-03
##  [921] 1.103097e-03 1.049773e-03 9.983767e-04 9.488638e-04 9.011892e-04
##  [926] 8.553085e-04 8.111778e-04 7.687537e-04 7.279932e-04 6.888539e-04
##  [931] 6.512937e-04 6.152709e-04 5.807446e-04 5.476740e-04 5.160191e-04
##  [936] 4.857400e-04 4.567977e-04 4.291534e-04 4.027690e-04 3.776066e-04
##  [941] 3.536290e-04 3.307997e-04 3.090822e-04 2.884410e-04 2.688408e-04
##  [946] 2.502470e-04 2.326253e-04 2.159423e-04 2.001647e-04 1.852600e-04
##  [951] 1.711962e-04 1.579418e-04 1.454657e-04 1.337377e-04 1.227278e-04
##  [956] 1.124066e-04 1.027456e-04 9.371643e-05 8.529150e-05 7.744374e-05
##  [961] 7.014666e-05 6.337434e-05 5.710144e-05 5.130320e-05 4.595544e-05
##  [966] 4.103457e-05 3.651761e-05 3.238215e-05 2.860639e-05 2.516914e-05
##  [971] 2.204981e-05 1.922842e-05 1.668559e-05 1.440259e-05 1.236129e-05
##  [976] 1.054418e-05 8.934391e-06 7.515677e-06 6.272431e-06 5.189687e-06
##  [981] 4.253117e-06 3.449042e-06 2.764433e-06 2.186912e-06 1.704762e-06
##  [986] 1.306923e-06 9.829998e-07 7.232666e-07 5.186678e-07 3.608230e-07
##  [991] 2.420305e-07 1.552708e-07 9.420996e-08 5.320379e-08 2.730089e-08
##  [996] 1.224652e-08 4.486118e-09 1.168945e-09 1.516472e-10 1.874062e-12
## 
## $likelihood
##    [1] 7.761871e-24 1.675574e-20 5.908100e-19 6.147255e-18 3.523904e-17
##    [6] 1.417102e-16 4.503766e-16 1.210391e-15 2.869048e-15 6.168413e-15
##   [11] 1.226678e-14 2.288662e-14 4.049044e-14 6.848549e-14 1.114582e-13
##   [16] 1.754390e-13 2.681973e-13 3.995670e-13 5.818031e-13 8.299674e-13
##   [21] 1.162344e-12 1.600885e-12 2.171679e-12 2.905456e-12 3.838108e-12
##   [26] 5.011240e-12 6.472746e-12 8.277397e-12 1.048746e-11 1.317335e-11
##   [31] 1.641423e-11 2.029874e-11 2.492563e-11 3.040450e-11 3.685643e-11
##   [36] 4.441478e-11 5.322582e-11 6.344954e-11 7.526028e-11 8.884755e-11
##   [41] 1.044167e-10 1.221896e-10 1.424055e-10 1.653215e-10 1.912137e-10
##   [46] 2.203773e-10 2.531277e-10 2.898011e-10 3.307550e-10 3.763691e-10
##   [51] 4.270456e-10 4.832101e-10 5.453119e-10 6.138249e-10 6.892479e-10
##   [56] 7.721050e-10 8.629465e-10 9.623490e-10 1.070916e-09 1.189277e-09
##   [61] 1.318092e-09 1.458047e-09 1.609855e-09 1.774260e-09 1.952034e-09
##   [66] 2.143977e-09 2.350919e-09 2.573720e-09 2.813268e-09 3.070482e-09
##   [71] 3.346308e-09 3.641725e-09 3.957737e-09 4.295383e-09 4.655726e-09
##   [76] 5.039861e-09 5.448911e-09 5.884029e-09 6.346394e-09 6.837215e-09
##   [81] 7.357727e-09 7.909196e-09 8.492911e-09 9.110190e-09 9.762375e-09
##   [86] 1.045084e-08 1.117697e-08 1.194218e-08 1.274793e-08 1.359567e-08
##   [91] 1.448689e-08 1.542310e-08 1.640584e-08 1.743664e-08 1.851708e-08
##   [96] 1.964875e-08 2.083326e-08 2.207222e-08 2.336728e-08 2.472008e-08
##  [101] 2.613230e-08 2.760561e-08 2.914171e-08 3.074229e-08 3.240908e-08
##  [106] 3.414379e-08 3.594815e-08 3.782389e-08 3.977276e-08 4.179649e-08
##  [111] 4.389685e-08 4.607557e-08 4.833440e-08 5.067510e-08 5.309942e-08
##  [116] 5.560909e-08 5.820586e-08 6.089145e-08 6.366761e-08 6.653605e-08
##  [121] 6.949848e-08 7.255660e-08 7.571209e-08 7.896664e-08 8.232191e-08
##  [126] 8.577953e-08 8.934114e-08 9.300835e-08 9.678276e-08 1.006659e-07
##  [131] 1.046594e-07 1.087647e-07 1.129833e-07 1.173168e-07 1.217665e-07
##  [136] 1.263339e-07 1.310204e-07 1.358274e-07 1.407561e-07 1.458078e-07
##  [141] 1.509839e-07 1.562856e-07 1.617140e-07 1.672704e-07 1.729558e-07
##  [146] 1.787713e-07 1.847180e-07 1.907968e-07 1.970088e-07 2.033548e-07
##  [151] 2.098357e-07 2.164524e-07 2.232056e-07 2.300961e-07 2.371245e-07
##  [156] 2.442914e-07 2.515976e-07 2.590435e-07 2.666296e-07 2.743564e-07
##  [161] 2.822243e-07 2.902335e-07 2.983844e-07 3.066773e-07 3.151122e-07
##  [166] 3.236894e-07 3.324088e-07 3.412706e-07 3.502746e-07 3.594208e-07
##  [171] 3.687090e-07 3.781391e-07 3.877106e-07 3.974235e-07 4.072771e-07
##  [176] 4.172712e-07 4.274052e-07 4.376786e-07 4.480907e-07 4.586409e-07
##  [181] 4.693285e-07 4.801526e-07 4.911125e-07 5.022072e-07 5.134358e-07
##  [186] 5.247972e-07 5.362905e-07 5.479144e-07 5.596678e-07 5.715494e-07
##  [191] 5.835579e-07 5.956920e-07 6.079502e-07 6.203312e-07 6.328333e-07
##  [196] 6.454550e-07 6.581946e-07 6.710505e-07 6.840210e-07 6.971042e-07
##  [201] 7.102982e-07 7.236013e-07 7.370114e-07 7.505265e-07 7.641447e-07
##  [206] 7.778637e-07 7.916815e-07 8.055959e-07 8.196046e-07 8.337053e-07
##  [211] 8.478957e-07 8.621735e-07 8.765362e-07 8.909813e-07 9.055064e-07
##  [216] 9.201089e-07 9.347862e-07 9.495358e-07 9.643548e-07 9.792407e-07
##  [221] 9.941907e-07 1.009202e-06 1.024272e-06 1.039398e-06 1.054576e-06
##  [226] 1.069804e-06 1.085080e-06 1.100399e-06 1.115759e-06 1.131158e-06
##  [231] 1.146591e-06 1.162057e-06 1.177551e-06 1.193071e-06 1.208614e-06
##  [236] 1.224176e-06 1.239754e-06 1.255346e-06 1.270947e-06 1.286555e-06
##  [241] 1.302167e-06 1.317779e-06 1.333388e-06 1.348990e-06 1.364584e-06
##  [246] 1.380164e-06 1.395728e-06 1.411273e-06 1.426795e-06 1.442291e-06
##  [251] 1.457758e-06 1.473192e-06 1.488590e-06 1.503950e-06 1.519266e-06
##  [256] 1.534537e-06 1.549759e-06 1.564929e-06 1.580043e-06 1.595098e-06
##  [261] 1.610091e-06 1.625019e-06 1.639878e-06 1.654665e-06 1.669377e-06
##  [266] 1.684012e-06 1.698564e-06 1.713033e-06 1.727414e-06 1.741704e-06
##  [271] 1.755900e-06 1.769999e-06 1.783999e-06 1.797895e-06 1.811686e-06
##  [276] 1.825367e-06 1.838937e-06 1.852392e-06 1.865730e-06 1.878947e-06
##  [281] 1.892041e-06 1.905008e-06 1.917847e-06 1.930554e-06 1.943127e-06
##  [286] 1.955563e-06 1.967860e-06 1.980014e-06 1.992024e-06 2.003886e-06
##  [291] 2.015599e-06 2.027160e-06 2.038566e-06 2.049816e-06 2.060907e-06
##  [296] 2.071836e-06 2.082602e-06 2.093202e-06 2.103634e-06 2.113896e-06
##  [301] 2.123986e-06 2.133902e-06 2.143643e-06 2.153205e-06 2.162588e-06
##  [306] 2.171790e-06 2.180808e-06 2.189641e-06 2.198288e-06 2.206747e-06
##  [311] 2.215015e-06 2.223093e-06 2.230978e-06 2.238669e-06 2.246164e-06
##  [316] 2.253463e-06 2.260564e-06 2.267465e-06 2.274166e-06 2.280666e-06
##  [321] 2.286963e-06 2.293057e-06 2.298947e-06 2.304631e-06 2.310109e-06
##  [326] 2.315380e-06 2.320444e-06 2.325299e-06 2.329946e-06 2.334383e-06
##  [331] 2.338610e-06 2.342627e-06 2.346432e-06 2.350027e-06 2.353410e-06
##  [336] 2.356582e-06 2.359541e-06 2.362288e-06 2.364824e-06 2.367146e-06
##  [341] 2.369257e-06 2.371156e-06 2.372842e-06 2.374317e-06 2.375580e-06
##  [346] 2.376631e-06 2.377472e-06 2.378101e-06 2.378521e-06 2.378730e-06
##  [351] 2.378730e-06 2.378521e-06 2.378104e-06 2.377480e-06 2.376648e-06
##  [356] 2.375610e-06 2.374367e-06 2.372920e-06 2.371268e-06 2.369414e-06
##  [361] 2.367358e-06 2.365102e-06 2.362645e-06 2.359990e-06 2.357138e-06
##  [366] 2.354089e-06 2.350844e-06 2.347406e-06 2.343776e-06 2.339954e-06
##  [371] 2.335942e-06 2.331741e-06 2.327354e-06 2.322781e-06 2.318024e-06
##  [376] 2.313084e-06 2.307963e-06 2.302663e-06 2.297185e-06 2.291531e-06
##  [381] 2.285702e-06 2.279701e-06 2.273529e-06 2.267189e-06 2.260680e-06
##  [386] 2.254007e-06 2.247170e-06 2.240171e-06 2.233013e-06 2.225696e-06
##  [391] 2.218224e-06 2.210599e-06 2.202821e-06 2.194894e-06 2.186820e-06
##  [396] 2.178599e-06 2.170236e-06 2.161731e-06 2.153087e-06 2.144305e-06
##  [401] 2.135389e-06 2.126341e-06 2.117162e-06 2.107854e-06 2.098421e-06
##  [406] 2.088864e-06 2.079186e-06 2.069388e-06 2.059473e-06 2.049444e-06
##  [411] 2.039303e-06 2.029051e-06 2.018692e-06 2.008227e-06 1.997659e-06
##  [416] 1.986991e-06 1.976224e-06 1.965362e-06 1.954405e-06 1.943357e-06
##  [421] 1.932221e-06 1.920997e-06 1.909690e-06 1.898300e-06 1.886832e-06
##  [426] 1.875285e-06 1.863665e-06 1.851971e-06 1.840208e-06 1.828376e-06
##  [431] 1.816480e-06 1.804520e-06 1.792499e-06 1.780421e-06 1.768286e-06
##  [436] 1.756097e-06 1.743856e-06 1.731567e-06 1.719230e-06 1.706849e-06
##  [441] 1.694426e-06 1.681962e-06 1.669461e-06 1.656924e-06 1.644353e-06
##  [446] 1.631752e-06 1.619121e-06 1.606464e-06 1.593782e-06 1.581077e-06
##  [451] 1.568352e-06 1.555609e-06 1.542850e-06 1.530077e-06 1.517292e-06
##  [456] 1.504496e-06 1.491694e-06 1.478885e-06 1.466072e-06 1.453258e-06
##  [461] 1.440444e-06 1.427632e-06 1.414823e-06 1.402021e-06 1.389227e-06
##  [466] 1.376442e-06 1.363669e-06 1.350910e-06 1.338165e-06 1.325437e-06
##  [471] 1.312729e-06 1.300040e-06 1.287374e-06 1.274732e-06 1.262115e-06
##  [476] 1.249525e-06 1.236965e-06 1.224434e-06 1.211936e-06 1.199471e-06
##  [481] 1.187042e-06 1.174649e-06 1.162294e-06 1.149978e-06 1.137704e-06
##  [486] 1.125472e-06 1.113283e-06 1.101140e-06 1.089044e-06 1.076995e-06
##  [491] 1.064995e-06 1.053045e-06 1.041147e-06 1.029302e-06 1.017511e-06
##  [496] 1.005775e-06 9.940951e-07 9.824726e-07 9.709085e-07 9.594040e-07
##  [501] 9.479599e-07 9.365774e-07 9.252574e-07 9.140009e-07 9.028087e-07
##  [506] 8.916819e-07 8.806213e-07 8.696277e-07 8.587020e-07 8.478450e-07
##  [511] 8.370575e-07 8.263401e-07 8.156938e-07 8.051191e-07 7.946167e-07
##  [516] 7.841874e-07 7.738317e-07 7.635502e-07 7.533436e-07 7.432124e-07
##  [521] 7.331572e-07 7.231784e-07 7.132766e-07 7.034522e-07 6.937056e-07
##  [526] 6.840374e-07 6.744479e-07 6.649375e-07 6.555065e-07 6.461553e-07
##  [531] 6.368842e-07 6.276935e-07 6.185834e-07 6.095543e-07 6.006063e-07
##  [536] 5.917396e-07 5.829545e-07 5.742511e-07 5.656296e-07 5.570900e-07
##  [541] 5.486326e-07 5.402573e-07 5.319643e-07 5.237536e-07 5.156252e-07
##  [546] 5.075792e-07 4.996156e-07 4.917342e-07 4.839352e-07 4.762184e-07
##  [551] 4.685838e-07 4.610313e-07 4.535608e-07 4.461721e-07 4.388651e-07
##  [556] 4.316397e-07 4.244957e-07 4.174329e-07 4.104511e-07 4.035501e-07
##  [561] 3.967297e-07 3.899897e-07 3.833297e-07 3.767495e-07 3.702488e-07
##  [566] 3.638274e-07 3.574849e-07 3.512210e-07 3.450354e-07 3.389277e-07
##  [571] 3.328976e-07 3.269448e-07 3.210687e-07 3.152691e-07 3.095456e-07
##  [576] 3.038976e-07 2.983249e-07 2.928271e-07 2.874035e-07 2.820539e-07
##  [581] 2.767777e-07 2.715745e-07 2.664439e-07 2.613853e-07 2.563982e-07
##  [586] 2.514822e-07 2.466368e-07 2.418614e-07 2.371555e-07 2.325186e-07
##  [591] 2.279503e-07 2.234498e-07 2.190168e-07 2.146507e-07 2.103508e-07
##  [596] 2.061167e-07 2.019478e-07 1.978435e-07 1.938033e-07 1.898266e-07
##  [601] 1.859128e-07 1.820613e-07 1.782716e-07 1.745430e-07 1.708750e-07
##  [606] 1.672669e-07 1.637183e-07 1.602284e-07 1.567967e-07 1.534226e-07
##  [611] 1.501055e-07 1.468448e-07 1.436398e-07 1.404900e-07 1.373948e-07
##  [616] 1.343535e-07 1.313655e-07 1.284303e-07 1.255472e-07 1.227156e-07
##  [621] 1.199349e-07 1.172045e-07 1.145238e-07 1.118922e-07 1.093091e-07
##  [626] 1.067738e-07 1.042859e-07 1.018445e-07 9.944927e-08 9.709948e-08
##  [631] 9.479455e-08 9.253390e-08 9.031693e-08 8.814305e-08 8.601166e-08
##  [636] 8.392219e-08 8.187405e-08 7.986665e-08 7.789942e-08 7.597179e-08
##  [641] 7.408317e-08 7.223301e-08 7.042073e-08 6.864577e-08 6.690758e-08
##  [646] 6.520559e-08 6.353926e-08 6.190803e-08 6.031136e-08 5.874872e-08
##  [651] 5.721955e-08 5.572334e-08 5.425954e-08 5.282764e-08 5.142712e-08
##  [656] 5.005746e-08 4.871815e-08 4.740867e-08 4.612853e-08 4.487724e-08
##  [661] 4.365428e-08 4.245919e-08 4.129146e-08 4.015063e-08 3.903622e-08
##  [666] 3.794775e-08 3.688476e-08 3.584680e-08 3.483340e-08 3.384411e-08
##  [671] 3.287849e-08 3.193611e-08 3.101651e-08 3.011928e-08 2.924399e-08
##  [676] 2.839021e-08 2.755753e-08 2.674554e-08 2.595384e-08 2.518203e-08
##  [681] 2.442970e-08 2.369647e-08 2.298196e-08 2.228578e-08 2.160757e-08
##  [686] 2.094694e-08 2.030355e-08 1.967701e-08 1.906700e-08 1.847314e-08
##  [691] 1.789511e-08 1.733255e-08 1.678514e-08 1.625255e-08 1.573445e-08
##  [696] 1.523052e-08 1.474045e-08 1.426394e-08 1.380066e-08 1.335034e-08
##  [701] 1.291266e-08 1.248734e-08 1.207410e-08 1.167266e-08 1.128273e-08
##  [706] 1.090405e-08 1.053636e-08 1.017938e-08 9.832869e-09 9.496566e-09
##  [711] 9.170224e-09 8.853601e-09 8.546457e-09 8.248557e-09 7.959671e-09
##  [716] 7.679571e-09 7.408038e-09 7.144852e-09 6.889800e-09 6.642672e-09
##  [721] 6.403265e-09 6.171375e-09 5.946806e-09 5.729365e-09 5.518862e-09
##  [726] 5.315112e-09 5.117933e-09 4.927148e-09 4.742581e-09 4.564062e-09
##  [731] 4.391426e-09 4.224507e-09 4.063148e-09 3.907191e-09 3.756483e-09
##  [736] 3.610876e-09 3.470222e-09 3.334381e-09 3.203211e-09 3.076576e-09
##  [741] 2.954345e-09 2.836385e-09 2.722571e-09 2.612779e-09 2.506887e-09
##  [746] 2.404777e-09 2.306334e-09 2.211446e-09 2.120003e-09 2.031898e-09
##  [751] 1.947027e-09 1.865288e-09 1.786583e-09 1.710814e-09 1.637888e-09
##  [756] 1.567713e-09 1.500200e-09 1.435262e-09 1.372814e-09 1.312776e-09
##  [761] 1.255065e-09 1.199605e-09 1.146320e-09 1.095136e-09 1.045982e-09
##  [766] 9.987881e-10 9.534867e-10 9.100120e-10 8.683005e-10 8.282901e-10
##  [771] 7.899207e-10 7.531339e-10 7.178730e-10 6.840832e-10 6.517110e-10
##  [776] 6.207048e-10 5.910144e-10 5.625913e-10 5.353882e-10 5.093597e-10
##  [781] 4.844614e-10 4.606506e-10 4.378858e-10 4.161268e-10 3.953349e-10
##  [786] 3.754725e-10 3.565033e-10 3.383920e-10 3.211047e-10 3.046087e-10
##  [791] 2.888720e-10 2.738642e-10 2.595555e-10 2.459174e-10 2.329223e-10
##  [796] 2.205435e-10 2.087555e-10 1.975332e-10 1.868530e-10 1.766917e-10
##  [801] 1.670271e-10 1.578379e-10 1.491034e-10 1.408038e-10 1.329200e-10
##  [806] 1.254335e-10 1.183268e-10 1.115828e-10 1.051851e-10 9.911804e-11
##  [811] 9.336647e-11 8.791587e-11 8.275231e-11 7.786240e-11 7.323328e-11
##  [816] 6.885264e-11 6.470863e-11 6.078993e-11 5.708565e-11 5.358539e-11
##  [821] 5.027916e-11 4.715740e-11 4.421096e-11 4.143110e-11 3.880942e-11
##  [826] 3.633793e-11 3.400895e-11 3.181517e-11 2.974961e-11 2.780556e-11
##  [831] 2.597668e-11 2.425686e-11 2.264030e-11 2.112148e-11 1.969511e-11
##  [836] 1.835616e-11 1.709986e-11 1.592163e-11 1.481714e-11 1.378226e-11
##  [841] 1.281306e-11 1.190582e-11 1.105698e-11 1.026318e-11 9.521222e-12
##  [846] 8.828066e-12 8.180836e-12 7.576802e-12 7.013376e-12 6.488110e-12
##  [851] 5.998680e-12 5.542893e-12 5.118668e-12 4.724043e-12 4.357159e-12
##  [856] 4.016263e-12 3.699698e-12 3.405902e-12 3.133400e-12 2.880802e-12
##  [861] 2.646798e-12 2.430155e-12 2.229712e-12 2.044376e-12 1.873121e-12
##  [866] 1.714981e-12 1.569049e-12 1.434474e-12 1.310459e-12 1.196255e-12
##  [871] 1.091160e-12 9.945174e-13 9.057127e-13 8.241707e-13 7.493542e-13
##  [876] 6.807609e-13 6.179223e-13 5.604012e-13 5.077898e-13 4.597085e-13
##  [881] 4.158036e-13 3.757464e-13 3.392308e-13 3.059729e-13 2.757088e-13
##  [886] 2.481940e-13 2.232016e-13 2.005217e-13 1.799597e-13 1.613360e-13
##  [891] 1.444844e-13 1.292516e-13 1.154962e-13 1.030877e-13 9.190611e-14
##  [896] 8.184100e-14 7.279082e-14 6.466235e-14 5.737004e-14 5.083549e-14
##  [901] 4.498690e-14 3.975860e-14 3.509057e-14 3.092802e-14 2.722098e-14
##  [906] 2.392394e-14 2.099547e-14 1.839793e-14 1.609714e-14 1.406211e-14
##  [911] 1.226476e-14 1.067970e-14 9.283993e-15 8.056924e-15 6.979840e-15
##  [916] 6.035950e-15 5.210165e-15 4.488940e-15 3.860140e-15 3.312902e-15
##  [921] 2.837522e-15 2.425342e-15 2.068649e-15 1.760585e-15 1.495061e-15
##  [926] 1.266679e-15 1.070664e-15 9.027999e-16 7.593673e-16 6.370951e-16
##  [931] 5.331103e-16 4.448949e-16 3.702465e-16 3.072428e-16 2.542097e-16
##  [936] 2.096925e-16 1.724300e-16 1.413316e-16 1.154560e-16 9.399340e-17
##  [941] 7.624865e-17 6.162665e-17 4.961940e-17 3.979451e-17 3.178502e-17
##  [946] 2.528048e-17 2.001905e-17 1.578067e-17 1.238099e-17 9.666131e-18
##  [951] 7.508141e-18 5.801017e-18 4.457299e-18 3.405127e-18 2.585709e-18
##  [956] 1.951169e-18 1.462700e-18 1.088999e-18 8.049522e-19 5.905139e-19
##  [961] 4.297777e-19 3.101948e-19 2.219278e-19 1.573147e-19 1.104293e-19
##  [966] 7.672099e-20 5.272244e-20 3.581298e-20 2.402894e-20 1.591229e-20
##  [971] 1.039094e-20 6.684714e-21 4.232065e-21 2.633589e-21 1.608775e-21
##  [976] 9.632744e-22 5.644012e-22 3.229908e-22 1.801465e-22 9.768576e-23
##  [981] 5.135505e-23 2.608943e-23 1.275917e-23 5.980147e-24 2.671912e-24
##  [986] 1.130786e-24 4.497991e-25 1.665674e-25 5.674252e-26 1.751344e-26
##  [991] 4.801625e-27 1.138925e-27 2.253772e-28 3.532228e-29 4.054604e-30
##  [996] 3.006433e-31 1.154053e-32 1.464234e-34 1.925852e-37 1.216437e-43
## 
## $posterior
##    [1] 1.089175e-19 7.025500e-16 4.112155e-14 5.966072e-13 4.379571e-12
##    [6] 2.143941e-11 8.020301e-11 2.477073e-10 6.627620e-10 1.586153e-09
##   [11] 3.472266e-09 7.066692e-09 1.353446e-08 2.462355e-08 4.286834e-08
##   [16] 7.183744e-08 1.164302e-07 1.832263e-07 2.808917e-07 4.206453e-07
##   [21] 6.167886e-07 8.873026e-07 1.254513e-06 1.745827e-06 2.394541e-06
##   [26] 3.240729e-06 4.332197e-06 5.725520e-06 7.487150e-06 9.694609e-06
##   [31] 1.243774e-05 1.582006e-05 1.996013e-05 2.499309e-05 3.107216e-05
##   [36] 3.837026e-05 4.708170e-05 5.742390e-05 6.963919e-05 8.399666e-05
##   [41] 1.007940e-04 1.203596e-04 1.430542e-04 1.692734e-04 1.994490e-04
##   [46] 2.340518e-04 2.735933e-04 3.186278e-04 3.697544e-04 4.276196e-04
##   [51] 4.929189e-04 5.663991e-04 6.488602e-04 7.411579e-04 8.442053e-04
##   [56] 9.589748e-04 1.086501e-03 1.227880e-03 1.384275e-03 1.556916e-03
##   [61] 1.747102e-03 1.956202e-03 2.185656e-03 2.436981e-03 2.711766e-03
##   [66] 3.011678e-03 3.338463e-03 3.693945e-03 4.080028e-03 4.498698e-03
##   [71] 4.952023e-03 5.442153e-03 5.971324e-03 6.541855e-03 7.156149e-03
##   [76] 7.816695e-03 8.526067e-03 9.286926e-03 1.010202e-02 1.097417e-02
##   [81] 1.190629e-02 1.290139e-02 1.396255e-02 1.509294e-02 1.629580e-02
##   [86] 1.757446e-02 1.893234e-02 2.037292e-02 2.189977e-02 2.351654e-02
##   [91] 2.522694e-02 2.703476e-02 2.894387e-02 3.095819e-02 3.308173e-02
##   [96] 3.531853e-02 3.767274e-02 4.014852e-02 4.275011e-02 4.548182e-02
##  [101] 4.834798e-02 5.135299e-02 5.450128e-02 5.779733e-02 6.124566e-02
##  [106] 6.485083e-02 6.861741e-02 7.255002e-02 7.665329e-02 8.093188e-02
##  [111] 8.539047e-02 9.003374e-02 9.486638e-02 9.989309e-02 1.051186e-01
##  [116] 1.105475e-01 1.161847e-01 1.220346e-01 1.281021e-01 1.343917e-01
##  [121] 1.409080e-01 1.476557e-01 1.546393e-01 1.618633e-01 1.693321e-01
##  [126] 1.770503e-01 1.850221e-01 1.932519e-01 2.017440e-01 2.105025e-01
##  [131] 2.195316e-01 2.288353e-01 2.384176e-01 2.482823e-01 2.584333e-01
##  [136] 2.688743e-01 2.796089e-01 2.906406e-01 3.019729e-01 3.136090e-01
##  [141] 3.255521e-01 3.378054e-01 3.503718e-01 3.632540e-01 3.764549e-01
##  [146] 3.899771e-01 4.038229e-01 4.179947e-01 4.324946e-01 4.473248e-01
##  [151] 4.624870e-01 4.779830e-01 4.938144e-01 5.099826e-01 5.264890e-01
##  [156] 5.433345e-01 5.605201e-01 5.780467e-01 5.959148e-01 6.141249e-01
##  [161] 6.326772e-01 6.515719e-01 6.708088e-01 6.903878e-01 7.103084e-01
##  [166] 7.305699e-01 7.511716e-01 7.721125e-01 7.933915e-01 8.150072e-01
##  [171] 8.369580e-01 8.592424e-01 8.818583e-01 9.048038e-01 9.280765e-01
##  [176] 9.516739e-01 9.755936e-01 9.998325e-01 1.024388e+00 1.049256e+00
##  [181] 1.074434e+00 1.099919e+00 1.125706e+00 1.151791e+00 1.178171e+00
##  [186] 1.204841e+00 1.231797e+00 1.259034e+00 1.286547e+00 1.314332e+00
##  [191] 1.342383e+00 1.370694e+00 1.399262e+00 1.428078e+00 1.457139e+00
##  [196] 1.486438e+00 1.515969e+00 1.545726e+00 1.575702e+00 1.605890e+00
##  [201] 1.636285e+00 1.666879e+00 1.697665e+00 1.728636e+00 1.759785e+00
##  [206] 1.791104e+00 1.822587e+00 1.854225e+00 1.886010e+00 1.917936e+00
##  [211] 1.949993e+00 1.982174e+00 2.014471e+00 2.046875e+00 2.079379e+00
##  [216] 2.111973e+00 2.144649e+00 2.177399e+00 2.210213e+00 2.243084e+00
##  [221] 2.276003e+00 2.308960e+00 2.341946e+00 2.374953e+00 2.407972e+00
##  [226] 2.440993e+00 2.474008e+00 2.507007e+00 2.539981e+00 2.572921e+00
##  [231] 2.605818e+00 2.638662e+00 2.671445e+00 2.704156e+00 2.736788e+00
##  [236] 2.769329e+00 2.801772e+00 2.834107e+00 2.866325e+00 2.898417e+00
##  [241] 2.930373e+00 2.962184e+00 2.993841e+00 3.025336e+00 3.056658e+00
##  [246] 3.087800e+00 3.118751e+00 3.149504e+00 3.180049e+00 3.210378e+00
##  [251] 3.240482e+00 3.270353e+00 3.299981e+00 3.329358e+00 3.358476e+00
##  [256] 3.387327e+00 3.415903e+00 3.444194e+00 3.472194e+00 3.499895e+00
##  [261] 3.527287e+00 3.554365e+00 3.581120e+00 3.607545e+00 3.633632e+00
##  [266] 3.659374e+00 3.684764e+00 3.709796e+00 3.734461e+00 3.758754e+00
##  [271] 3.782668e+00 3.806196e+00 3.829332e+00 3.852070e+00 3.874404e+00
##  [276] 3.896327e+00 3.917835e+00 3.938922e+00 3.959581e+00 3.979809e+00
##  [281] 3.999599e+00 4.018947e+00 4.037848e+00 4.056298e+00 4.074291e+00
##  [286] 4.091823e+00 4.108891e+00 4.125490e+00 4.141617e+00 4.157267e+00
##  [291] 4.172438e+00 4.187126e+00 4.201327e+00 4.215040e+00 4.228260e+00
##  [296] 4.240986e+00 4.253215e+00 4.264944e+00 4.276172e+00 4.286896e+00
##  [301] 4.297115e+00 4.306828e+00 4.316032e+00 4.324726e+00 4.332910e+00
##  [306] 4.340582e+00 4.347742e+00 4.354388e+00 4.360522e+00 4.366142e+00
##  [311] 4.371248e+00 4.375840e+00 4.379919e+00 4.383485e+00 4.386539e+00
##  [316] 4.389081e+00 4.391112e+00 4.392634e+00 4.393646e+00 4.394152e+00
##  [321] 4.394152e+00 4.393648e+00 4.392642e+00 4.391136e+00 4.389132e+00
##  [326] 4.386631e+00 4.383638e+00 4.380153e+00 4.376181e+00 4.371723e+00
##  [331] 4.366783e+00 4.361363e+00 4.355467e+00 4.349099e+00 4.342262e+00
##  [336] 4.334958e+00 4.327193e+00 4.318970e+00 4.310293e+00 4.301166e+00
##  [341] 4.291593e+00 4.281578e+00 4.271127e+00 4.260243e+00 4.248931e+00
##  [346] 4.237195e+00 4.225041e+00 4.212474e+00 4.199498e+00 4.186119e+00
##  [351] 4.172341e+00 4.158170e+00 4.143612e+00 4.128672e+00 4.113354e+00
##  [356] 4.097665e+00 4.081611e+00 4.065197e+00 4.048429e+00 4.031312e+00
##  [361] 4.013853e+00 3.996057e+00 3.977931e+00 3.959481e+00 3.940712e+00
##  [366] 3.921631e+00 3.902243e+00 3.882556e+00 3.862575e+00 3.842307e+00
##  [371] 3.821758e+00 3.800934e+00 3.779842e+00 3.758488e+00 3.736878e+00
##  [376] 3.715019e+00 3.692918e+00 3.670580e+00 3.648013e+00 3.625222e+00
##  [381] 3.602214e+00 3.578997e+00 3.555575e+00 3.531956e+00 3.508146e+00
##  [386] 3.484152e+00 3.459980e+00 3.435637e+00 3.411128e+00 3.386461e+00
##  [391] 3.361641e+00 3.336676e+00 3.311571e+00 3.286333e+00 3.260968e+00
##  [396] 3.235483e+00 3.209883e+00 3.184175e+00 3.158365e+00 3.132459e+00
##  [401] 3.106464e+00 3.080385e+00 3.054229e+00 3.028000e+00 3.001706e+00
##  [406] 2.975353e+00 2.948945e+00 2.922489e+00 2.895991e+00 2.869456e+00
##  [411] 2.842890e+00 2.816298e+00 2.789687e+00 2.763060e+00 2.736425e+00
##  [416] 2.709786e+00 2.683148e+00 2.656517e+00 2.629898e+00 2.603296e+00
##  [421] 2.576716e+00 2.550162e+00 2.523641e+00 2.497156e+00 2.470712e+00
##  [426] 2.444314e+00 2.417967e+00 2.391676e+00 2.365444e+00 2.339275e+00
##  [431] 2.313176e+00 2.287148e+00 2.261198e+00 2.235328e+00 2.209544e+00
##  [436] 2.183848e+00 2.158244e+00 2.132737e+00 2.107331e+00 2.082028e+00
##  [441] 2.056832e+00 2.031747e+00 2.006777e+00 1.981923e+00 1.957191e+00
##  [446] 1.932582e+00 1.908101e+00 1.883749e+00 1.859531e+00 1.835448e+00
##  [451] 1.811504e+00 1.787701e+00 1.764041e+00 1.740529e+00 1.717165e+00
##  [456] 1.693952e+00 1.670893e+00 1.647989e+00 1.625244e+00 1.602658e+00
##  [461] 1.580235e+00 1.557976e+00 1.535882e+00 1.513956e+00 1.492199e+00
##  [466] 1.470613e+00 1.449200e+00 1.427961e+00 1.406897e+00 1.386010e+00
##  [471] 1.365301e+00 1.344771e+00 1.324422e+00 1.304254e+00 1.284269e+00
##  [476] 1.264467e+00 1.244850e+00 1.225417e+00 1.206171e+00 1.187111e+00
##  [481] 1.168238e+00 1.149554e+00 1.131057e+00 1.112749e+00 1.094631e+00
##  [486] 1.076701e+00 1.058962e+00 1.041412e+00 1.024053e+00 1.006884e+00
##  [491] 9.899045e-01 9.731156e-01 9.565169e-01 9.401083e-01 9.238896e-01
##  [496] 9.078606e-01 8.920210e-01 8.763706e-01 8.609090e-01 8.456357e-01
##  [501] 8.305504e-01 8.156527e-01 8.009418e-01 7.864173e-01 7.720787e-01
##  [506] 7.579251e-01 7.439560e-01 7.301706e-01 7.165682e-01 7.031479e-01
##  [511] 6.899090e-01 6.768505e-01 6.639716e-01 6.512713e-01 6.387487e-01
##  [516] 6.264027e-01 6.142324e-01 6.022367e-01 5.904145e-01 5.787647e-01
##  [521] 5.672862e-01 5.559779e-01 5.448384e-01 5.338667e-01 5.230615e-01
##  [526] 5.124216e-01 5.019457e-01 4.916324e-01 4.814806e-01 4.714888e-01
##  [531] 4.616557e-01 4.519800e-01 4.424603e-01 4.330952e-01 4.238833e-01
##  [536] 4.148232e-01 4.059134e-01 3.971525e-01 3.885390e-01 3.800716e-01
##  [541] 3.717486e-01 3.635687e-01 3.555303e-01 3.476319e-01 3.398721e-01
##  [546] 3.322493e-01 3.247620e-01 3.174088e-01 3.101879e-01 3.030980e-01
##  [551] 2.961375e-01 2.893049e-01 2.825985e-01 2.760170e-01 2.695588e-01
##  [556] 2.632222e-01 2.570058e-01 2.509081e-01 2.449274e-01 2.390624e-01
##  [561] 2.333114e-01 2.276728e-01 2.221453e-01 2.167273e-01 2.114172e-01
##  [566] 2.062135e-01 2.011148e-01 1.961195e-01 1.912262e-01 1.864333e-01
##  [571] 1.817394e-01 1.771430e-01 1.726426e-01 1.682368e-01 1.639241e-01
##  [576] 1.597031e-01 1.555724e-01 1.515305e-01 1.475760e-01 1.437075e-01
##  [581] 1.399237e-01 1.362231e-01 1.326043e-01 1.290660e-01 1.256069e-01
##  [586] 1.222256e-01 1.189208e-01 1.156912e-01 1.125354e-01 1.094521e-01
##  [591] 1.064402e-01 1.034983e-01 1.006251e-01 9.781943e-02 9.508006e-02
##  [596] 9.240576e-02 8.979533e-02 8.724756e-02 8.476129e-02 8.233535e-02
##  [601] 7.996858e-02 7.765986e-02 7.540804e-02 7.321203e-02 7.107072e-02
##  [606] 6.898303e-02 6.694790e-02 6.496426e-02 6.303107e-02 6.114730e-02
##  [611] 5.931195e-02 5.752400e-02 5.578247e-02 5.408639e-02 5.243480e-02
##  [616] 5.082674e-02 4.926130e-02 4.773755e-02 4.625459e-02 4.481152e-02
##  [621] 4.340748e-02 4.204158e-02 4.071300e-02 3.942088e-02 3.816441e-02
##  [626] 3.694277e-02 3.575517e-02 3.460083e-02 3.347897e-02 3.238885e-02
##  [631] 3.132970e-02 3.030081e-02 2.930145e-02 2.833091e-02 2.738851e-02
##  [636] 2.647357e-02 2.558541e-02 2.472337e-02 2.388682e-02 2.307512e-02
##  [641] 2.228765e-02 2.152380e-02 2.078298e-02 2.006460e-02 1.936808e-02
##  [646] 1.869287e-02 1.803841e-02 1.740417e-02 1.678960e-02 1.619420e-02
##  [651] 1.561746e-02 1.505888e-02 1.451797e-02 1.399426e-02 1.348727e-02
##  [656] 1.299656e-02 1.252168e-02 1.206218e-02 1.161765e-02 1.118766e-02
##  [661] 1.077181e-02 1.036970e-02 9.980936e-03 9.605137e-03 9.241932e-03
##  [666] 8.890958e-03 8.551858e-03 8.224286e-03 7.907903e-03 7.602379e-03
##  [671] 7.307391e-03 7.022625e-03 6.747773e-03 6.482536e-03 6.226622e-03
##  [676] 5.979746e-03 5.741631e-03 5.512005e-03 5.290606e-03 5.077176e-03
##  [681] 4.871466e-03 4.673231e-03 4.482235e-03 4.298246e-03 4.121040e-03
##  [686] 3.950397e-03 3.786104e-03 3.627956e-03 3.475749e-03 3.329288e-03
##  [691] 3.188384e-03 3.052849e-03 2.922505e-03 2.797176e-03 2.676693e-03
##  [696] 2.560891e-03 2.449609e-03 2.342692e-03 2.239988e-03 2.141351e-03
##  [701] 2.046639e-03 1.955713e-03 1.868440e-03 1.784689e-03 1.704336e-03
##  [706] 1.627257e-03 1.553335e-03 1.482454e-03 1.414504e-03 1.349378e-03
##  [711] 1.286970e-03 1.227180e-03 1.169910e-03 1.115066e-03 1.062556e-03
##  [716] 1.012292e-03 9.641878e-04 9.181612e-04 8.741319e-04 8.320225e-04
##  [721] 7.917584e-04 7.532672e-04 7.164792e-04 6.813270e-04 6.477454e-04
##  [726] 6.156718e-04 5.850454e-04 5.558077e-04 5.279024e-04 5.012749e-04
##  [731] 4.758727e-04 4.516454e-04 4.285440e-04 4.065215e-04 3.855328e-04
##  [736] 3.655341e-04 3.464835e-04 3.283405e-04 3.110662e-04 2.946231e-04
##  [741] 2.789752e-04 2.640879e-04 2.499278e-04 2.364629e-04 2.236624e-04
##  [746] 2.114967e-04 1.999373e-04 1.889570e-04 1.785296e-04 1.686299e-04
##  [751] 1.592337e-04 1.503178e-04 1.418602e-04 1.338393e-04 1.262349e-04
##  [756] 1.190273e-04 1.121977e-04 1.057282e-04 9.960152e-05 9.380124e-05
##  [761] 8.831156e-05 8.311737e-05 7.820422e-05 7.355832e-05 6.916644e-05
##  [766] 6.501596e-05 6.109481e-05 5.739147e-05 5.389493e-05 5.059467e-05
##  [771] 4.748067e-05 4.454335e-05 4.177358e-05 3.916265e-05 3.670225e-05
##  [776] 3.438447e-05 3.220176e-05 3.014695e-05 2.821320e-05 2.639399e-05
##  [781] 2.468313e-05 2.307472e-05 2.156315e-05 2.014310e-05 1.880950e-05
##  [786] 1.755754e-05 1.638263e-05 1.528046e-05 1.424688e-05 1.327800e-05
##  [791] 1.237011e-05 1.151968e-05 1.072338e-05 9.978066e-06 9.280732e-06
##  [796] 8.628549e-06 8.018835e-06 7.449052e-06 6.916800e-06 6.419809e-06
##  [801] 5.955935e-06 5.523151e-06 5.119541e-06 4.743299e-06 4.392718e-06
##  [806] 4.066190e-06 3.762196e-06 3.479307e-06 3.216173e-06 2.971526e-06
##  [811] 2.744169e-06 2.532977e-06 2.336892e-06 2.154918e-06 1.986119e-06
##  [816] 1.829616e-06 1.684584e-06 1.550247e-06 1.425877e-06 1.310793e-06
##  [821] 1.204355e-06 1.105962e-06 1.015054e-06 9.311055e-07 8.536234e-07
##  [826] 7.821480e-07 7.162489e-07 6.555242e-07 5.995980e-07 5.481200e-07
##  [831] 5.007628e-07 4.572212e-07 4.172110e-07 3.804669e-07 3.467422e-07
##  [836] 3.158072e-07 2.874481e-07 2.614663e-07 2.376772e-07 2.159092e-07
##  [841] 1.960033e-07 1.778118e-07 1.611977e-07 1.460342e-07 1.322038e-07
##  [846] 1.195978e-07 1.081155e-07 9.766411e-08 8.815755e-08 7.951655e-08
##  [851] 7.166792e-08 6.454415e-08 5.808305e-08 5.222733e-08 4.692426e-08
##  [856] 4.212535e-08 3.778603e-08 3.386536e-08 3.032577e-08 2.713280e-08
##  [861] 2.425485e-08 2.166299e-08 1.933075e-08 1.723391e-08 1.535034e-08
##  [866] 1.365983e-08 1.214395e-08 1.078589e-08 9.570318e-09 8.483311e-09
##  [871] 7.512185e-09 6.645418e-09 5.872549e-09 5.184086e-09 4.571428e-09
##  [876] 4.026782e-09 3.543102e-09 3.114014e-09 2.733763e-09 2.397157e-09
##  [881] 2.099516e-09 1.836625e-09 1.604691e-09 1.400307e-09 1.220412e-09
##  [886] 1.062262e-09 9.233993e-10 8.016223e-10 6.949638e-10 6.016672e-10
##  [891] 5.201653e-10 4.490620e-10 3.871147e-10 3.332190e-10 2.863944e-10
##  [896] 2.457716e-10 2.105806e-10 1.801405e-10 1.538499e-10 1.311784e-10
##  [901] 1.116586e-10 9.487936e-11 8.047966e-11 6.814274e-11 5.759117e-11
##  [906] 4.858233e-11 4.090438e-11 3.437261e-11 2.882624e-11 2.412556e-11
##  [911] 2.014936e-11 1.679264e-11 1.396464e-11 1.158703e-11 9.592319e-12
##  [916] 7.922474e-12 6.527695e-12 5.365321e-12 4.398885e-12 3.597277e-12
##  [921] 2.934010e-12 2.386585e-12 1.935931e-12 1.565920e-12 1.262942e-12
##  [926] 1.015542e-12 8.141007e-13 6.505601e-13 5.181887e-13 4.113771e-13
##  [931] 3.254638e-13 2.565858e-13 2.015510e-13 1.577294e-13 1.229608e-13
##  [936] 9.547628e-14 7.383216e-14 5.685393e-14 4.358943e-14 3.326945e-14
##  [941] 2.527486e-14 1.910919e-14 1.437587e-14 1.075942e-14 8.009887e-15
##  [946] 5.930111e-15 4.365251e-15 3.194272e-15 2.323013e-15 1.678584e-15
##  [951] 1.204857e-15 8.588353e-16 6.077728e-16 4.268701e-16 2.974618e-16
##  [956] 2.055870e-16 1.408728e-16 9.566478e-17 6.435534e-17 4.286721e-17
##  [961] 2.825918e-17 1.842709e-17 1.187866e-17 7.565235e-18 4.756964e-18
##  [966] 2.951025e-18 1.804706e-18 1.087064e-18 6.443271e-19 3.754134e-19
##  [971] 2.147675e-19 1.204856e-19 6.619157e-20 3.555474e-20 1.864094e-20
##  [976] 9.520760e-21 4.726739e-21 2.275448e-21 1.059182e-21 4.752052e-22
##  [981] 2.047383e-22 8.434742e-23 3.306262e-23 1.225891e-23 4.269673e-24
##  [986] 1.385285e-24 4.144582e-25 1.129268e-25 2.758714e-26 5.923446e-27
##  [991] 1.089349e-27 1.657652e-28 1.990286e-29 1.761570e-30 1.037610e-31
##  [996] 3.451222e-33 4.852943e-35 1.604403e-37 2.737574e-41 2.136891e-49
## 
## $pi
##    [1] 0.0005 0.0015 0.0025 0.0035 0.0045 0.0055 0.0065 0.0075 0.0085 0.0095
##   [11] 0.0105 0.0115 0.0125 0.0135 0.0145 0.0155 0.0165 0.0175 0.0185 0.0195
##   [21] 0.0205 0.0215 0.0225 0.0235 0.0245 0.0255 0.0265 0.0275 0.0285 0.0295
##   [31] 0.0305 0.0315 0.0325 0.0335 0.0345 0.0355 0.0365 0.0375 0.0385 0.0395
##   [41] 0.0405 0.0415 0.0425 0.0435 0.0445 0.0455 0.0465 0.0475 0.0485 0.0495
##   [51] 0.0505 0.0515 0.0525 0.0535 0.0545 0.0555 0.0565 0.0575 0.0585 0.0595
##   [61] 0.0605 0.0615 0.0625 0.0635 0.0645 0.0655 0.0665 0.0675 0.0685 0.0695
##   [71] 0.0705 0.0715 0.0725 0.0735 0.0745 0.0755 0.0765 0.0775 0.0785 0.0795
##   [81] 0.0805 0.0815 0.0825 0.0835 0.0845 0.0855 0.0865 0.0875 0.0885 0.0895
##   [91] 0.0905 0.0915 0.0925 0.0935 0.0945 0.0955 0.0965 0.0975 0.0985 0.0995
##  [101] 0.1005 0.1015 0.1025 0.1035 0.1045 0.1055 0.1065 0.1075 0.1085 0.1095
##  [111] 0.1105 0.1115 0.1125 0.1135 0.1145 0.1155 0.1165 0.1175 0.1185 0.1195
##  [121] 0.1205 0.1215 0.1225 0.1235 0.1245 0.1255 0.1265 0.1275 0.1285 0.1295
##  [131] 0.1305 0.1315 0.1325 0.1335 0.1345 0.1355 0.1365 0.1375 0.1385 0.1395
##  [141] 0.1405 0.1415 0.1425 0.1435 0.1445 0.1455 0.1465 0.1475 0.1485 0.1495
##  [151] 0.1505 0.1515 0.1525 0.1535 0.1545 0.1555 0.1565 0.1575 0.1585 0.1595
##  [161] 0.1605 0.1615 0.1625 0.1635 0.1645 0.1655 0.1665 0.1675 0.1685 0.1695
##  [171] 0.1705 0.1715 0.1725 0.1735 0.1745 0.1755 0.1765 0.1775 0.1785 0.1795
##  [181] 0.1805 0.1815 0.1825 0.1835 0.1845 0.1855 0.1865 0.1875 0.1885 0.1895
##  [191] 0.1905 0.1915 0.1925 0.1935 0.1945 0.1955 0.1965 0.1975 0.1985 0.1995
##  [201] 0.2005 0.2015 0.2025 0.2035 0.2045 0.2055 0.2065 0.2075 0.2085 0.2095
##  [211] 0.2105 0.2115 0.2125 0.2135 0.2145 0.2155 0.2165 0.2175 0.2185 0.2195
##  [221] 0.2205 0.2215 0.2225 0.2235 0.2245 0.2255 0.2265 0.2275 0.2285 0.2295
##  [231] 0.2305 0.2315 0.2325 0.2335 0.2345 0.2355 0.2365 0.2375 0.2385 0.2395
##  [241] 0.2405 0.2415 0.2425 0.2435 0.2445 0.2455 0.2465 0.2475 0.2485 0.2495
##  [251] 0.2505 0.2515 0.2525 0.2535 0.2545 0.2555 0.2565 0.2575 0.2585 0.2595
##  [261] 0.2605 0.2615 0.2625 0.2635 0.2645 0.2655 0.2665 0.2675 0.2685 0.2695
##  [271] 0.2705 0.2715 0.2725 0.2735 0.2745 0.2755 0.2765 0.2775 0.2785 0.2795
##  [281] 0.2805 0.2815 0.2825 0.2835 0.2845 0.2855 0.2865 0.2875 0.2885 0.2895
##  [291] 0.2905 0.2915 0.2925 0.2935 0.2945 0.2955 0.2965 0.2975 0.2985 0.2995
##  [301] 0.3005 0.3015 0.3025 0.3035 0.3045 0.3055 0.3065 0.3075 0.3085 0.3095
##  [311] 0.3105 0.3115 0.3125 0.3135 0.3145 0.3155 0.3165 0.3175 0.3185 0.3195
##  [321] 0.3205 0.3215 0.3225 0.3235 0.3245 0.3255 0.3265 0.3275 0.3285 0.3295
##  [331] 0.3305 0.3315 0.3325 0.3335 0.3345 0.3355 0.3365 0.3375 0.3385 0.3395
##  [341] 0.3405 0.3415 0.3425 0.3435 0.3445 0.3455 0.3465 0.3475 0.3485 0.3495
##  [351] 0.3505 0.3515 0.3525 0.3535 0.3545 0.3555 0.3565 0.3575 0.3585 0.3595
##  [361] 0.3605 0.3615 0.3625 0.3635 0.3645 0.3655 0.3665 0.3675 0.3685 0.3695
##  [371] 0.3705 0.3715 0.3725 0.3735 0.3745 0.3755 0.3765 0.3775 0.3785 0.3795
##  [381] 0.3805 0.3815 0.3825 0.3835 0.3845 0.3855 0.3865 0.3875 0.3885 0.3895
##  [391] 0.3905 0.3915 0.3925 0.3935 0.3945 0.3955 0.3965 0.3975 0.3985 0.3995
##  [401] 0.4005 0.4015 0.4025 0.4035 0.4045 0.4055 0.4065 0.4075 0.4085 0.4095
##  [411] 0.4105 0.4115 0.4125 0.4135 0.4145 0.4155 0.4165 0.4175 0.4185 0.4195
##  [421] 0.4205 0.4215 0.4225 0.4235 0.4245 0.4255 0.4265 0.4275 0.4285 0.4295
##  [431] 0.4305 0.4315 0.4325 0.4335 0.4345 0.4355 0.4365 0.4375 0.4385 0.4395
##  [441] 0.4405 0.4415 0.4425 0.4435 0.4445 0.4455 0.4465 0.4475 0.4485 0.4495
##  [451] 0.4505 0.4515 0.4525 0.4535 0.4545 0.4555 0.4565 0.4575 0.4585 0.4595
##  [461] 0.4605 0.4615 0.4625 0.4635 0.4645 0.4655 0.4665 0.4675 0.4685 0.4695
##  [471] 0.4705 0.4715 0.4725 0.4735 0.4745 0.4755 0.4765 0.4775 0.4785 0.4795
##  [481] 0.4805 0.4815 0.4825 0.4835 0.4845 0.4855 0.4865 0.4875 0.4885 0.4895
##  [491] 0.4905 0.4915 0.4925 0.4935 0.4945 0.4955 0.4965 0.4975 0.4985 0.4995
##  [501] 0.5005 0.5015 0.5025 0.5035 0.5045 0.5055 0.5065 0.5075 0.5085 0.5095
##  [511] 0.5105 0.5115 0.5125 0.5135 0.5145 0.5155 0.5165 0.5175 0.5185 0.5195
##  [521] 0.5205 0.5215 0.5225 0.5235 0.5245 0.5255 0.5265 0.5275 0.5285 0.5295
##  [531] 0.5305 0.5315 0.5325 0.5335 0.5345 0.5355 0.5365 0.5375 0.5385 0.5395
##  [541] 0.5405 0.5415 0.5425 0.5435 0.5445 0.5455 0.5465 0.5475 0.5485 0.5495
##  [551] 0.5505 0.5515 0.5525 0.5535 0.5545 0.5555 0.5565 0.5575 0.5585 0.5595
##  [561] 0.5605 0.5615 0.5625 0.5635 0.5645 0.5655 0.5665 0.5675 0.5685 0.5695
##  [571] 0.5705 0.5715 0.5725 0.5735 0.5745 0.5755 0.5765 0.5775 0.5785 0.5795
##  [581] 0.5805 0.5815 0.5825 0.5835 0.5845 0.5855 0.5865 0.5875 0.5885 0.5895
##  [591] 0.5905 0.5915 0.5925 0.5935 0.5945 0.5955 0.5965 0.5975 0.5985 0.5995
##  [601] 0.6005 0.6015 0.6025 0.6035 0.6045 0.6055 0.6065 0.6075 0.6085 0.6095
##  [611] 0.6105 0.6115 0.6125 0.6135 0.6145 0.6155 0.6165 0.6175 0.6185 0.6195
##  [621] 0.6205 0.6215 0.6225 0.6235 0.6245 0.6255 0.6265 0.6275 0.6285 0.6295
##  [631] 0.6305 0.6315 0.6325 0.6335 0.6345 0.6355 0.6365 0.6375 0.6385 0.6395
##  [641] 0.6405 0.6415 0.6425 0.6435 0.6445 0.6455 0.6465 0.6475 0.6485 0.6495
##  [651] 0.6505 0.6515 0.6525 0.6535 0.6545 0.6555 0.6565 0.6575 0.6585 0.6595
##  [661] 0.6605 0.6615 0.6625 0.6635 0.6645 0.6655 0.6665 0.6675 0.6685 0.6695
##  [671] 0.6705 0.6715 0.6725 0.6735 0.6745 0.6755 0.6765 0.6775 0.6785 0.6795
##  [681] 0.6805 0.6815 0.6825 0.6835 0.6845 0.6855 0.6865 0.6875 0.6885 0.6895
##  [691] 0.6905 0.6915 0.6925 0.6935 0.6945 0.6955 0.6965 0.6975 0.6985 0.6995
##  [701] 0.7005 0.7015 0.7025 0.7035 0.7045 0.7055 0.7065 0.7075 0.7085 0.7095
##  [711] 0.7105 0.7115 0.7125 0.7135 0.7145 0.7155 0.7165 0.7175 0.7185 0.7195
##  [721] 0.7205 0.7215 0.7225 0.7235 0.7245 0.7255 0.7265 0.7275 0.7285 0.7295
##  [731] 0.7305 0.7315 0.7325 0.7335 0.7345 0.7355 0.7365 0.7375 0.7385 0.7395
##  [741] 0.7405 0.7415 0.7425 0.7435 0.7445 0.7455 0.7465 0.7475 0.7485 0.7495
##  [751] 0.7505 0.7515 0.7525 0.7535 0.7545 0.7555 0.7565 0.7575 0.7585 0.7595
##  [761] 0.7605 0.7615 0.7625 0.7635 0.7645 0.7655 0.7665 0.7675 0.7685 0.7695
##  [771] 0.7705 0.7715 0.7725 0.7735 0.7745 0.7755 0.7765 0.7775 0.7785 0.7795
##  [781] 0.7805 0.7815 0.7825 0.7835 0.7845 0.7855 0.7865 0.7875 0.7885 0.7895
##  [791] 0.7905 0.7915 0.7925 0.7935 0.7945 0.7955 0.7965 0.7975 0.7985 0.7995
##  [801] 0.8005 0.8015 0.8025 0.8035 0.8045 0.8055 0.8065 0.8075 0.8085 0.8095
##  [811] 0.8105 0.8115 0.8125 0.8135 0.8145 0.8155 0.8165 0.8175 0.8185 0.8195
##  [821] 0.8205 0.8215 0.8225 0.8235 0.8245 0.8255 0.8265 0.8275 0.8285 0.8295
##  [831] 0.8305 0.8315 0.8325 0.8335 0.8345 0.8355 0.8365 0.8375 0.8385 0.8395
##  [841] 0.8405 0.8415 0.8425 0.8435 0.8445 0.8455 0.8465 0.8475 0.8485 0.8495
##  [851] 0.8505 0.8515 0.8525 0.8535 0.8545 0.8555 0.8565 0.8575 0.8585 0.8595
##  [861] 0.8605 0.8615 0.8625 0.8635 0.8645 0.8655 0.8665 0.8675 0.8685 0.8695
##  [871] 0.8705 0.8715 0.8725 0.8735 0.8745 0.8755 0.8765 0.8775 0.8785 0.8795
##  [881] 0.8805 0.8815 0.8825 0.8835 0.8845 0.8855 0.8865 0.8875 0.8885 0.8895
##  [891] 0.8905 0.8915 0.8925 0.8935 0.8945 0.8955 0.8965 0.8975 0.8985 0.8995
##  [901] 0.9005 0.9015 0.9025 0.9035 0.9045 0.9055 0.9065 0.9075 0.9085 0.9095
##  [911] 0.9105 0.9115 0.9125 0.9135 0.9145 0.9155 0.9165 0.9175 0.9185 0.9195
##  [921] 0.9205 0.9215 0.9225 0.9235 0.9245 0.9255 0.9265 0.9275 0.9285 0.9295
##  [931] 0.9305 0.9315 0.9325 0.9335 0.9345 0.9355 0.9365 0.9375 0.9385 0.9395
##  [941] 0.9405 0.9415 0.9425 0.9435 0.9445 0.9455 0.9465 0.9475 0.9485 0.9495
##  [951] 0.9505 0.9515 0.9525 0.9535 0.9545 0.9555 0.9565 0.9575 0.9585 0.9595
##  [961] 0.9605 0.9615 0.9625 0.9635 0.9645 0.9655 0.9665 0.9675 0.9685 0.9695
##  [971] 0.9705 0.9715 0.9725 0.9735 0.9745 0.9755 0.9765 0.9775 0.9785 0.9795
##  [981] 0.9805 0.9815 0.9825 0.9835 0.9845 0.9855 0.9865 0.9875 0.9885 0.9895
##  [991] 0.9905 0.9915 0.9925 0.9935 0.9945 0.9955 0.9965 0.9975 0.9985 0.9995
## 
## $pi.prior
##    [1] 1.497002e-02 4.473061e-02 7.425281e-02 1.035377e-01 1.325864e-01
##    [6] 1.613998e-01 1.899792e-01 2.183256e-01 2.464399e-01 2.743234e-01
##   [11] 3.019769e-01 3.294017e-01 3.565986e-01 3.835689e-01 4.103135e-01
##   [16] 4.368334e-01 4.631297e-01 4.892035e-01 5.150557e-01 5.406874e-01
##   [21] 5.660996e-01 5.912934e-01 6.162697e-01 6.410296e-01 6.655741e-01
##   [26] 6.899042e-01 7.140209e-01 7.379253e-01 7.616182e-01 7.851008e-01
##   [31] 8.083740e-01 8.314388e-01 8.542963e-01 8.769473e-01 8.993929e-01
##   [36] 9.216341e-01 9.436718e-01 9.655071e-01 9.871409e-01 1.008574e+00
##   [41] 1.029808e+00 1.050843e+00 1.071680e+00 1.092321e+00 1.112766e+00
##   [46] 1.133017e+00 1.153074e+00 1.172937e+00 1.192609e+00 1.212090e+00
##   [51] 1.231381e+00 1.250483e+00 1.269397e+00 1.288124e+00 1.306664e+00
##   [56] 1.325019e+00 1.343190e+00 1.361177e+00 1.378981e+00 1.396604e+00
##   [61] 1.414047e+00 1.431309e+00 1.448393e+00 1.465299e+00 1.482027e+00
##   [66] 1.498579e+00 1.514957e+00 1.531159e+00 1.547189e+00 1.563045e+00
##   [71] 1.578730e+00 1.594245e+00 1.609589e+00 1.624764e+00 1.639771e+00
##   [76] 1.654611e+00 1.669284e+00 1.683792e+00 1.698135e+00 1.712315e+00
##   [81] 1.726331e+00 1.740185e+00 1.753878e+00 1.767411e+00 1.780784e+00
##   [86] 1.793999e+00 1.807056e+00 1.819956e+00 1.832699e+00 1.845288e+00
##   [91] 1.857722e+00 1.870002e+00 1.882130e+00 1.894105e+00 1.905930e+00
##   [96] 1.917604e+00 1.929129e+00 1.940505e+00 1.951733e+00 1.962814e+00
##  [101] 1.973749e+00 1.984539e+00 1.995184e+00 2.005685e+00 2.016044e+00
##  [106] 2.026260e+00 2.036334e+00 2.046268e+00 2.056063e+00 2.065718e+00
##  [111] 2.075235e+00 2.084615e+00 2.093858e+00 2.102965e+00 2.111937e+00
##  [116] 2.120775e+00 2.129479e+00 2.138051e+00 2.146490e+00 2.154798e+00
##  [121] 2.162976e+00 2.171024e+00 2.178943e+00 2.186734e+00 2.194397e+00
##  [126] 2.201934e+00 2.209344e+00 2.216630e+00 2.223791e+00 2.230828e+00
##  [131] 2.237743e+00 2.244535e+00 2.251205e+00 2.257755e+00 2.264185e+00
##  [136] 2.270495e+00 2.276687e+00 2.282761e+00 2.288717e+00 2.294558e+00
##  [141] 2.300282e+00 2.305892e+00 2.311387e+00 2.316768e+00 2.322037e+00
##  [146] 2.327194e+00 2.332239e+00 2.337173e+00 2.341997e+00 2.346712e+00
##  [151] 2.351318e+00 2.355815e+00 2.360206e+00 2.364490e+00 2.368668e+00
##  [156] 2.372741e+00 2.376709e+00 2.380573e+00 2.384334e+00 2.387992e+00
##  [161] 2.391548e+00 2.395003e+00 2.398357e+00 2.401612e+00 2.404767e+00
##  [166] 2.407823e+00 2.410782e+00 2.413643e+00 2.416407e+00 2.419076e+00
##  [171] 2.421649e+00 2.424127e+00 2.426512e+00 2.428802e+00 2.431000e+00
##  [176] 2.433106e+00 2.435120e+00 2.437043e+00 2.438876e+00 2.440619e+00
##  [181] 2.442273e+00 2.443839e+00 2.445317e+00 2.446708e+00 2.448012e+00
##  [186] 2.449230e+00 2.450362e+00 2.451410e+00 2.452374e+00 2.453254e+00
##  [191] 2.454051e+00 2.454766e+00 2.455399e+00 2.455951e+00 2.456422e+00
##  [196] 2.456814e+00 2.457125e+00 2.457358e+00 2.457513e+00 2.457590e+00
##  [201] 2.457590e+00 2.457514e+00 2.457361e+00 2.457134e+00 2.456831e+00
##  [206] 2.456454e+00 2.456004e+00 2.455480e+00 2.454884e+00 2.454216e+00
##  [211] 2.453476e+00 2.452666e+00 2.451785e+00 2.450835e+00 2.449815e+00
##  [216] 2.448727e+00 2.447570e+00 2.446346e+00 2.445055e+00 2.443697e+00
##  [221] 2.442274e+00 2.440785e+00 2.439231e+00 2.437612e+00 2.435930e+00
##  [226] 2.434184e+00 2.432376e+00 2.430505e+00 2.428573e+00 2.426579e+00
##  [231] 2.424525e+00 2.422410e+00 2.420236e+00 2.418003e+00 2.415710e+00
##  [236] 2.413360e+00 2.410952e+00 2.408487e+00 2.405965e+00 2.403387e+00
##  [241] 2.400753e+00 2.398064e+00 2.395320e+00 2.392522e+00 2.389670e+00
##  [246] 2.386765e+00 2.383808e+00 2.380798e+00 2.377736e+00 2.374623e+00
##  [251] 2.371459e+00 2.368244e+00 2.364980e+00 2.361666e+00 2.358303e+00
##  [256] 2.354892e+00 2.351432e+00 2.347926e+00 2.344371e+00 2.340771e+00
##  [261] 2.337124e+00 2.333431e+00 2.329693e+00 2.325910e+00 2.322082e+00
##  [266] 2.318211e+00 2.314296e+00 2.310338e+00 2.306338e+00 2.302295e+00
##  [271] 2.298210e+00 2.294084e+00 2.289917e+00 2.285710e+00 2.281462e+00
##  [276] 2.277175e+00 2.272849e+00 2.268484e+00 2.264080e+00 2.259639e+00
##  [281] 2.255160e+00 2.250644e+00 2.246091e+00 2.241502e+00 2.236877e+00
##  [286] 2.232216e+00 2.227521e+00 2.222791e+00 2.218026e+00 2.213228e+00
##  [291] 2.208396e+00 2.203531e+00 2.198634e+00 2.193704e+00 2.188742e+00
##  [296] 2.183749e+00 2.178725e+00 2.173670e+00 2.168584e+00 2.163469e+00
##  [301] 2.158324e+00 2.153150e+00 2.147947e+00 2.142715e+00 2.137456e+00
##  [306] 2.132168e+00 2.126854e+00 2.121512e+00 2.116144e+00 2.110749e+00
##  [311] 2.105329e+00 2.099883e+00 2.094412e+00 2.088916e+00 2.083396e+00
##  [316] 2.077851e+00 2.072283e+00 2.066691e+00 2.061077e+00 2.055439e+00
##  [321] 2.049779e+00 2.044098e+00 2.038394e+00 2.032669e+00 2.026924e+00
##  [326] 2.021157e+00 2.015370e+00 2.009563e+00 2.003737e+00 1.997891e+00
##  [331] 1.992026e+00 1.986143e+00 1.980241e+00 1.974321e+00 1.968383e+00
##  [336] 1.962428e+00 1.956456e+00 1.950467e+00 1.944461e+00 1.938440e+00
##  [341] 1.932402e+00 1.926349e+00 1.920281e+00 1.914198e+00 1.908101e+00
##  [346] 1.901989e+00 1.895863e+00 1.889723e+00 1.883570e+00 1.877404e+00
##  [351] 1.871225e+00 1.865033e+00 1.858829e+00 1.852613e+00 1.846386e+00
##  [356] 1.840147e+00 1.833897e+00 1.827637e+00 1.821365e+00 1.815084e+00
##  [361] 1.808792e+00 1.802491e+00 1.796181e+00 1.789861e+00 1.783532e+00
##  [366] 1.777195e+00 1.770850e+00 1.764496e+00 1.758135e+00 1.751766e+00
##  [371] 1.745390e+00 1.739006e+00 1.732617e+00 1.726220e+00 1.719817e+00
##  [376] 1.713409e+00 1.706994e+00 1.700574e+00 1.694149e+00 1.687719e+00
##  [381] 1.681284e+00 1.674845e+00 1.668401e+00 1.661954e+00 1.655502e+00
##  [386] 1.649048e+00 1.642589e+00 1.636128e+00 1.629664e+00 1.623198e+00
##  [391] 1.616729e+00 1.610258e+00 1.603785e+00 1.597310e+00 1.590834e+00
##  [396] 1.584357e+00 1.577878e+00 1.571399e+00 1.564920e+00 1.558440e+00
##  [401] 1.551960e+00 1.545480e+00 1.539001e+00 1.532522e+00 1.526043e+00
##  [406] 1.519566e+00 1.513090e+00 1.506615e+00 1.500142e+00 1.493670e+00
##  [411] 1.487201e+00 1.480734e+00 1.474269e+00 1.467807e+00 1.461347e+00
##  [416] 1.454891e+00 1.448437e+00 1.441987e+00 1.435541e+00 1.429098e+00
##  [421] 1.422660e+00 1.416225e+00 1.409795e+00 1.403369e+00 1.396948e+00
##  [426] 1.390532e+00 1.384121e+00 1.377715e+00 1.371314e+00 1.364919e+00
##  [431] 1.358530e+00 1.352147e+00 1.345770e+00 1.339399e+00 1.333035e+00
##  [436] 1.326677e+00 1.320326e+00 1.313982e+00 1.307645e+00 1.301315e+00
##  [441] 1.294993e+00 1.288679e+00 1.282372e+00 1.276073e+00 1.269782e+00
##  [446] 1.263500e+00 1.257226e+00 1.250960e+00 1.244703e+00 1.238455e+00
##  [451] 1.232216e+00 1.225986e+00 1.219765e+00 1.213554e+00 1.207353e+00
##  [456] 1.201161e+00 1.194979e+00 1.188807e+00 1.182645e+00 1.176494e+00
##  [461] 1.170352e+00 1.164222e+00 1.158102e+00 1.151993e+00 1.145895e+00
##  [466] 1.139808e+00 1.133733e+00 1.127668e+00 1.121615e+00 1.115574e+00
##  [471] 1.109545e+00 1.103527e+00 1.097521e+00 1.091528e+00 1.085546e+00
##  [476] 1.079577e+00 1.073621e+00 1.067677e+00 1.061746e+00 1.055827e+00
##  [481] 1.049922e+00 1.044029e+00 1.038150e+00 1.032284e+00 1.026431e+00
##  [486] 1.020592e+00 1.014766e+00 1.008954e+00 1.003156e+00 9.973715e-01
##  [491] 9.916011e-01 9.858449e-01 9.801029e-01 9.743752e-01 9.686618e-01
##  [496] 9.629630e-01 9.572787e-01 9.516091e-01 9.459543e-01 9.403144e-01
##  [501] 9.346894e-01 9.290794e-01 9.234846e-01 9.179050e-01 9.123407e-01
##  [506] 9.067918e-01 9.012584e-01 8.957406e-01 8.902384e-01 8.847519e-01
##  [511] 8.792812e-01 8.738264e-01 8.683876e-01 8.629648e-01 8.575581e-01
##  [516] 8.521677e-01 8.467935e-01 8.414356e-01 8.360941e-01 8.307692e-01
##  [521] 8.254608e-01 8.201690e-01 8.148939e-01 8.096355e-01 8.043940e-01
##  [526] 7.991694e-01 7.939617e-01 7.887711e-01 7.835975e-01 7.784411e-01
##  [531] 7.733018e-01 7.681798e-01 7.630752e-01 7.579879e-01 7.529180e-01
##  [536] 7.478657e-01 7.428309e-01 7.378136e-01 7.328140e-01 7.278322e-01
##  [541] 7.228680e-01 7.179217e-01 7.129932e-01 7.080826e-01 7.031900e-01
##  [546] 6.983153e-01 6.934587e-01 6.886201e-01 6.837997e-01 6.789974e-01
##  [551] 6.742134e-01 6.694476e-01 6.647000e-01 6.599708e-01 6.552600e-01
##  [556] 6.505675e-01 6.458935e-01 6.412380e-01 6.366009e-01 6.319824e-01
##  [561] 6.273825e-01 6.228012e-01 6.182384e-01 6.136944e-01 6.091691e-01
##  [566] 6.046624e-01 6.001746e-01 5.957054e-01 5.912551e-01 5.868236e-01
##  [571] 5.824110e-01 5.780172e-01 5.736424e-01 5.692864e-01 5.649494e-01
##  [576] 5.606313e-01 5.563322e-01 5.520521e-01 5.477911e-01 5.435490e-01
##  [581] 5.393260e-01 5.351220e-01 5.309371e-01 5.267713e-01 5.226246e-01
##  [586] 5.184971e-01 5.143886e-01 5.102993e-01 5.062291e-01 5.021780e-01
##  [591] 4.981462e-01 4.941334e-01 4.901399e-01 4.861655e-01 4.822104e-01
##  [596] 4.782744e-01 4.743576e-01 4.704600e-01 4.665816e-01 4.627224e-01
##  [601] 4.588824e-01 4.550616e-01 4.512600e-01 4.474776e-01 4.437144e-01
##  [606] 4.399704e-01 4.362455e-01 4.325399e-01 4.288534e-01 4.251862e-01
##  [611] 4.215380e-01 4.179091e-01 4.142993e-01 4.107086e-01 4.071371e-01
##  [616] 4.035847e-01 4.000514e-01 3.965372e-01 3.930422e-01 3.895661e-01
##  [621] 3.861092e-01 3.826713e-01 3.792525e-01 3.758527e-01 3.724719e-01
##  [626] 3.691100e-01 3.657672e-01 3.624433e-01 3.591384e-01 3.558524e-01
##  [631] 3.525852e-01 3.493370e-01 3.461076e-01 3.428971e-01 3.397054e-01
##  [636] 3.365324e-01 3.333783e-01 3.302429e-01 3.271262e-01 3.240283e-01
##  [641] 3.209490e-01 3.178883e-01 3.148463e-01 3.118229e-01 3.088180e-01
##  [646] 3.058317e-01 3.028639e-01 2.999145e-01 2.969837e-01 2.940712e-01
##  [651] 2.911771e-01 2.883014e-01 2.854441e-01 2.826050e-01 2.797842e-01
##  [656] 2.769816e-01 2.741972e-01 2.714309e-01 2.686828e-01 2.659528e-01
##  [661] 2.632408e-01 2.605468e-01 2.578708e-01 2.552128e-01 2.525727e-01
##  [666] 2.499504e-01 2.473459e-01 2.447593e-01 2.421903e-01 2.396391e-01
##  [671] 2.371055e-01 2.345896e-01 2.320912e-01 2.296104e-01 2.271471e-01
##  [676] 2.247012e-01 2.222727e-01 2.198616e-01 2.174678e-01 2.150913e-01
##  [681] 2.127320e-01 2.103899e-01 2.080649e-01 2.057570e-01 2.034661e-01
##  [686] 2.011923e-01 1.989354e-01 1.966953e-01 1.944721e-01 1.922658e-01
##  [691] 1.900761e-01 1.879032e-01 1.857469e-01 1.836072e-01 1.814840e-01
##  [696] 1.793773e-01 1.772871e-01 1.752133e-01 1.731558e-01 1.711145e-01
##  [701] 1.690895e-01 1.670807e-01 1.650880e-01 1.631113e-01 1.611507e-01
##  [706] 1.592060e-01 1.572772e-01 1.553643e-01 1.534671e-01 1.515857e-01
##  [711] 1.497199e-01 1.478698e-01 1.460352e-01 1.442161e-01 1.424124e-01
##  [716] 1.406242e-01 1.388512e-01 1.370935e-01 1.353510e-01 1.336237e-01
##  [721] 1.319114e-01 1.302142e-01 1.285319e-01 1.268646e-01 1.252120e-01
##  [726] 1.235743e-01 1.219512e-01 1.203428e-01 1.187491e-01 1.171698e-01
##  [731] 1.156050e-01 1.140546e-01 1.125185e-01 1.109967e-01 1.094892e-01
##  [736] 1.079957e-01 1.065164e-01 1.050511e-01 1.035997e-01 1.021622e-01
##  [741] 1.007386e-01 9.932867e-02 9.793246e-02 9.654987e-02 9.518085e-02
##  [746] 9.382532e-02 9.248322e-02 9.115447e-02 8.983902e-02 8.853680e-02
##  [751] 8.724773e-02 8.597176e-02 8.470881e-02 8.345881e-02 8.222171e-02
##  [756] 8.099742e-02 7.978588e-02 7.858703e-02 7.740079e-02 7.622709e-02
##  [761] 7.506587e-02 7.391706e-02 7.278059e-02 7.165638e-02 7.054438e-02
##  [766] 6.944450e-02 6.835668e-02 6.728085e-02 6.621694e-02 6.516488e-02
##  [771] 6.412460e-02 6.309603e-02 6.207909e-02 6.107372e-02 6.007985e-02
##  [776] 5.909741e-02 5.812632e-02 5.716651e-02 5.621792e-02 5.528047e-02
##  [781] 5.435409e-02 5.343871e-02 5.253426e-02 5.164067e-02 5.075786e-02
##  [786] 4.988576e-02 4.902431e-02 4.817343e-02 4.733304e-02 4.650309e-02
##  [791] 4.568349e-02 4.487417e-02 4.407506e-02 4.328610e-02 4.250720e-02
##  [796] 4.173829e-02 4.097931e-02 4.023019e-02 3.949084e-02 3.876120e-02
##  [801] 3.804120e-02 3.733076e-02 3.662981e-02 3.593829e-02 3.525611e-02
##  [806] 3.458320e-02 3.391950e-02 3.326494e-02 3.261943e-02 3.198291e-02
##  [811] 3.135531e-02 3.073656e-02 3.012657e-02 2.952529e-02 2.893264e-02
##  [816] 2.834854e-02 2.777293e-02 2.720574e-02 2.664689e-02 2.609631e-02
##  [821] 2.555393e-02 2.501968e-02 2.449349e-02 2.397528e-02 2.346499e-02
##  [826] 2.296255e-02 2.246788e-02 2.198091e-02 2.150157e-02 2.102980e-02
##  [831] 2.056552e-02 2.010866e-02 1.965915e-02 1.921692e-02 1.878190e-02
##  [836] 1.835402e-02 1.793322e-02 1.751941e-02 1.711254e-02 1.671253e-02
##  [841] 1.631931e-02 1.593281e-02 1.555297e-02 1.517972e-02 1.481298e-02
##  [846] 1.445269e-02 1.409879e-02 1.375119e-02 1.340984e-02 1.307467e-02
##  [851] 1.274560e-02 1.242258e-02 1.210553e-02 1.179438e-02 1.148908e-02
##  [856] 1.118955e-02 1.089573e-02 1.060754e-02 1.032493e-02 1.004783e-02
##  [861] 9.776175e-03 9.509893e-03 9.248923e-03 8.993199e-03 8.742657e-03
##  [866] 8.497232e-03 8.256859e-03 8.021475e-03 7.791016e-03 7.565417e-03
##  [871] 7.344616e-03 7.128549e-03 6.917154e-03 6.710367e-03 6.508126e-03
##  [876] 6.310370e-03 6.117036e-03 5.928062e-03 5.743388e-03 5.562952e-03
##  [881] 5.386693e-03 5.214552e-03 5.046467e-03 4.882379e-03 4.722229e-03
##  [886] 4.565956e-03 4.413503e-03 4.264810e-03 4.119819e-03 3.978472e-03
##  [891] 3.840711e-03 3.706479e-03 3.575719e-03 3.448373e-03 3.324386e-03
##  [896] 3.203701e-03 3.086263e-03 2.972016e-03 2.860905e-03 2.752876e-03
##  [901] 2.647874e-03 2.545845e-03 2.446735e-03 2.350491e-03 2.257061e-03
##  [906] 2.166392e-03 2.078430e-03 1.993126e-03 1.910428e-03 1.830283e-03
##  [911] 1.752643e-03 1.677456e-03 1.604672e-03 1.534243e-03 1.466120e-03
##  [916] 1.400253e-03 1.336594e-03 1.275097e-03 1.215712e-03 1.158395e-03
##  [921] 1.103097e-03 1.049773e-03 9.983767e-04 9.488638e-04 9.011892e-04
##  [926] 8.553085e-04 8.111778e-04 7.687537e-04 7.279932e-04 6.888539e-04
##  [931] 6.512937e-04 6.152709e-04 5.807446e-04 5.476740e-04 5.160191e-04
##  [936] 4.857400e-04 4.567977e-04 4.291534e-04 4.027690e-04 3.776066e-04
##  [941] 3.536290e-04 3.307997e-04 3.090822e-04 2.884410e-04 2.688408e-04
##  [946] 2.502470e-04 2.326253e-04 2.159423e-04 2.001647e-04 1.852600e-04
##  [951] 1.711962e-04 1.579418e-04 1.454657e-04 1.337377e-04 1.227278e-04
##  [956] 1.124066e-04 1.027456e-04 9.371643e-05 8.529150e-05 7.744374e-05
##  [961] 7.014666e-05 6.337434e-05 5.710144e-05 5.130320e-05 4.595544e-05
##  [966] 4.103457e-05 3.651761e-05 3.238215e-05 2.860639e-05 2.516914e-05
##  [971] 2.204981e-05 1.922842e-05 1.668559e-05 1.440259e-05 1.236129e-05
##  [976] 1.054418e-05 8.934391e-06 7.515677e-06 6.272431e-06 5.189687e-06
##  [981] 4.253117e-06 3.449042e-06 2.764433e-06 2.186912e-06 1.704762e-06
##  [986] 1.306923e-06 9.829998e-07 7.232666e-07 5.186678e-07 3.608230e-07
##  [991] 2.420305e-07 1.552708e-07 9.420996e-08 5.320379e-08 2.730089e-08
##  [996] 1.224652e-08 4.486118e-09 1.168945e-09 1.516472e-10 1.874062e-12
## 
## attr(,"class")
## [1] "Bolstad"

Mencari sebaran kumulatif posterior (CDF) menggunakan Simpson’s rule

cdf.post = cdf(simulasi.post.pi)
plot(cdf.post, type = "l", xlab = expression(pi[0]),
     ylab = expression(Pr(pi <= pi[0])))

Inferensi Bayesian dari posterior: nilai harapan dan ragam

round(mean(simulasi.post.pi),5)
## [1] 0.33333
round(var(simulasi.post.pi),5)
## [1] 0.00794

Menggunakan fungsi quantile untuk menentukan 95% credible interval

nilai.kuantil <- quantile(simulasi.post.pi, probs = c(0.025, 0.975))
cat(paste("Approximate 95% credible interval : ["
          , round(nilai.kuantil[1], 5), " ", round(nilai.kuantil[2], 5), "]\n", sep = ""))
## Approximate 95% credible interval : [0.17164 0.5174]

.