diff --git a/cmdstanpy/model.py b/cmdstanpy/model.py index 8a3d067b..11519886 100644 --- a/cmdstanpy/model.py +++ b/cmdstanpy/model.py @@ -205,9 +205,9 @@ def __init__( self._compiler_options.add_include_path(path) # try to detect models w/out parameters, needed for sampler - if not cmdstan_version_before( - 2, 27 - ): # unknown end of version range + if (not cmdstan_version_before(2, 27)) and cmdstan_version_before( + 2, 36 + ): try: model_info = self.src_info() if 'parameters' in model_info: diff --git a/cmdstanpy/stanfit/mcmc.py b/cmdstanpy/stanfit/mcmc.py index 0bc1e599..8c93c8e3 100644 --- a/cmdstanpy/stanfit/mcmc.py +++ b/cmdstanpy/stanfit/mcmc.py @@ -412,7 +412,7 @@ def _assemble_draws(self) -> None: self._step_size[chain] = float(step_size.strip()) if self._metadata.cmdstan_config['metric'] != 'unit_e': line = fd.readline().strip() # metric type - line = fd.readline().lstrip(' #\t') + line = fd.readline().lstrip(' #\t').rstrip() num_unconstrained_params = len(line.split(',')) if chain == 0: # can't allocate w/o num params if self.metric_type == 'diag_e': @@ -429,18 +429,21 @@ def _assemble_draws(self) -> None: ), dtype=float, ) - if self.metric_type == 'diag_e': - xs = line.split(',') - self._metric[chain, :] = [float(x) for x in xs] - else: - xs = line.split(',') - self._metric[chain, 0, :] = [float(x) for x in xs] - for i in range(1, num_unconstrained_params): - line = fd.readline().lstrip(' #\t').strip() + if line: + if self.metric_type == 'diag_e': xs = line.split(',') - self._metric[chain, i, :] = [ + self._metric[chain, :] = [float(x) for x in xs] + else: + xs = line.strip().split(',') + self._metric[chain, 0, :] = [ float(x) for x in xs ] + for i in range(1, num_unconstrained_params): + line = fd.readline().lstrip(' #\t').rstrip() + xs = line.split(',') + self._metric[chain, i, :] = [ + float(x) for x in xs + ] else: # unit_e changed in 2.34 to have an extra line pos = fd.tell() line = fd.readline().strip() diff --git a/test/data/no_param_hmc_sample.csv b/test/data/no_param_hmc_sample.csv new file mode 100644 index 00000000..91999e46 --- /dev/null +++ b/test/data/no_param_hmc_sample.csv @@ -0,0 +1,157 @@ +# stan_version_major = 2 +# stan_version_minor = 35 +# stan_version_patch = 0 +# model = datagen_poisson_glm_model +# start_datetime = 2024-12-02 14:45:12 UTC +# method = sample (Default) +# sample +# num_samples = 100 +# num_warmup = 1000 (Default) +# save_warmup = false (Default) +# thin = 1 (Default) +# adapt +# engaged = true (Default) +# gamma = 0.05 (Default) +# delta = 0.8 (Default) +# kappa = 0.75 (Default) +# t0 = 10 (Default) +# init_buffer = 75 (Default) +# term_buffer = 50 (Default) +# window = 25 (Default) +# save_metric = false (Default) +# algorithm = hmc (Default) +# hmc +# engine = nuts (Default) +# nuts +# max_depth = 10 (Default) +# metric = diag_e (Default) +# metric_file = (Default) +# stepsize = 1 (Default) +# stepsize_jitter = 0 (Default) +# num_chains = 1 (Default) +# id = 1 (Default) +# data +# file = (Default) +# init = 2 (Default) +# random +# seed = 2399056448 (Default) +# output +# file = no_param_hmc_sample.csv +# diagnostic_file = (Default) +# refresh = 100 (Default) +# sig_figs = -1 (Default) +# profile_file = profile.csv (Default) +# save_cmdstan_config = false (Default) +# num_threads = 1 (Default) +# stanc_version = %%NAME%%3 %%VERSION%% +# stancflags = +lp__,accept_stat__,stepsize__,treedepth__,n_leapfrog__,divergent__,energy__,N,y_sim.1,y_sim.2,y_sim.3,y_sim.4,y_sim.5,y_sim.6,y_sim.7,y_sim.8,y_sim.9,y_sim.10,y_sim.11,y_sim.12,y_sim.13,y_sim.14,y_sim.15,y_sim.16,y_sim.17,y_sim.18,y_sim.19,y_sim.20,x_sim.1,x_sim.2,x_sim.3,x_sim.4,x_sim.5,x_sim.6,x_sim.7,x_sim.8,x_sim.9,x_sim.10,x_sim.11,x_sim.12,x_sim.13,x_sim.14,x_sim.15,x_sim.16,x_sim.17,x_sim.18,x_sim.19,x_sim.20,pop_sim.1,pop_sim.2,pop_sim.3,pop_sim.4,pop_sim.5,pop_sim.6,pop_sim.7,pop_sim.8,pop_sim.9,pop_sim.10,pop_sim.11,pop_sim.12,pop_sim.13,pop_sim.14,pop_sim.15,pop_sim.16,pop_sim.17,pop_sim.18,pop_sim.19,pop_sim.20,alpha_sim,beta_sim,eta.1,eta.2,eta.3,eta.4,eta.5,eta.6,eta.7,eta.8,eta.9,eta.10,eta.11,eta.12,eta.13,eta.14,eta.15,eta.16,eta.17,eta.18,eta.19,eta.20 +# Adaptation terminated +# Step size = nan +# Diagonal elements of inverse mass matrix: +# +0,1,nan,1,1,0,0,20,6750,1423,3205,2667,1966,2508,1669,6584,5527,4653,1713,1633,7137,5913,5501,6521,3586,4828,6227,1643,0.0358751,0.3613,0.0568551,0.328804,0.551602,0.907548,0.986684,0.0930657,0.171859,0.110039,0.73272,0.4307,0.0434005,0.228586,0.520925,0.0700279,0.758075,0.507096,0.237532,0.668033,2247.9,540.034,1081.91,976.307,806.367,1137.97,796.464,2221.16,1913.15,1565.87,774.843,619.084,2367.8,2142.01,2192.39,2205.31,1545.16,1948.16,2211.37,670.302,1.11747,-0.370597,8.82193,7.27521,8.08288,7.87939,7.60559,7.81814,7.43199,8.78877,8.61029,8.43289,7.49859,7.3861,8.8711,8.70226,8.61717,8.79014,8.17941,8.50418,8.73081,7.37763 +0,1,nan,1,1,0,0,20,185,112,151,255,170,353,296,380,381,271,97,125,354,78,270,100,384,339,285,285,0.767628,0.833944,0.890277,0.498764,0.675498,0.836085,0.242281,0.812561,0.794329,0.669046,0.215553,0.895611,0.557843,0.409189,0.124885,0.478168,0.805523,0.346138,0.709774,0.513694,1069.1,719.358,925.497,1534.46,1206.01,2233.36,2277.88,2458.2,2444.86,1863.5,747.707,683.77,2215.12,593.475,2040.06,697.22,2248.59,2212.21,1899.32,1974.05,-2.05261,0.252962,5.11614,4.73671,5.00293,5.40949,5.21334,5.87015,5.73968,5.96012,5.95007,5.64685,4.61893,4.70157,5.79157,4.43689,5.59972,4.61545,5.86922,5.7367,5.67619,5.66517 +0,1,nan,1,1,0,0,20,843,864,429,2615,676,711,1088,868,1123,1228,2468,1432,891,428,738,1797,1437,884,809,314,0.947415,0.743777,0.807928,0.136811,0.478371,0.535572,0.959296,0.612102,0.951738,0.964139,0.0390374,0.0933832,0.615627,0.720168,0.711152,0.583801,0.436258,0.323434,0.428534,0.876856,1605.16,1312.44,702.864,2145.46,776.288,912.863,2034.61,1119.84,2050.12,2338.09,1889.06,1125.96,1142.08,634.935,1112.62,2308.65,1683.7,860.648,930.225,562.119,0.319134,-0.988456,6.76363,6.76359,6.0757,7.85501,6.50081,6.60633,6.98897,6.73504,7.00403,7.12321,7.82438,7.25322,6.75122,6.0608,6.63067,7.48649,7.31666,6.75712,6.73097,5.78411 +0,1,nan,1,1,0,0,20,4626,7667,7150,6264,13074,3342,9266,5430,5295,4175,10641,3156,1896,7649,15203,6965,16358,12796,7438,10223,0.436206,0.422895,0.801021,0.310217,0.681982,0.212956,0.741143,0.629379,0.642764,0.405095,0.52686,0.248013,0.0521713,0.956044,0.823841,0.267075,0.852849,0.753056,0.229052,0.444614,928.396,1566.52,993.323,1407.32,2070.25,821.631,1398.95,906.779,879.288,861.596,1948.85,755.504,582.826,915.524,2074.02,1597.1,2195.01,1845.26,1811.89,2041.54,1.18941,0.967583,8.44493,8.9552,8.86552,8.73901,9.4847,8.10675,9.15,8.60828,8.59045,8.34016,9.27418,8.05676,7.60778,8.93395,9.62379,8.82377,9.70855,9.43843,8.91316,9.24107 +0,1,nan,1,1,0,0,20,4095,1360,2978,1819,8232,2063,2480,2967,5355,3163,5762,3549,2042,5731,2936,6234,8516,2470,5481,7486,0.501393,0.888994,0.777071,0.806781,0.0108512,0.230819,0.739926,0.126279,0.166629,0.361696,0.526674,0.0648142,0.621933,0.319603,0.493195,0.085164,0.0289425,0.520329,0.299834,0.253178,1598.59,728.386,1477.87,908.126,2110.19,649.029,1214.07,857.999,1585.71,1101.75,2267.72,985.286,869.934,1953.35,1168.12,1739.56,2247.73,951.516,1812.55,2416.79,1.3499,-0.824398,8.31343,7.20784,8.00764,7.49617,8.99548,7.63509,7.84164,8.0004,8.58131,8.05637,8.64224,8.1894,7.60559,8.66372,8.00646,8.74108,9.04371,7.779,8.6052,8.93137 +0,1,nan,1,1,0,0,20,6227,2649,4227,10507,2752,3279,6138,11051,6879,3860,3595,8573,3448,3422,7835,10183,5530,10221,10047,3074,0.297983,0.770367,0.361813,0.200616,0.745457,0.606017,0.154471,0.0559532,0.580136,0.0263179,0.395554,0.0530344,0.101761,0.375638,0.221684,0.500883,0.313314,0.574526,0.78758,0.915164,1345.19,543.862,889.649,2228.94,565.773,667.706,1337.69,2400.58,1472.33,816.038,738.181,1854.86,727.853,709.643,1687.62,2110.02,1145.61,2126.91,2041.24,635.18,1.53426,0.0448381,8.75191,7.8675,8.34131,9.25254,7.90588,8.06528,8.73989,9.32024,8.85487,8.2399,8.15619,9.06221,8.12892,8.11587,8.97528,9.21117,8.592,9.22245,9.19089,8.02921 +0,1,nan,1,1,0,0,20,4056,6060,5017,2441,4405,4766,4001,5555,2360,2897,8106,8261,5911,2784,7002,8574,8799,4092,4910,5228,0.791505,0.603934,0.248974,0.658333,0.107072,0.27714,0.708065,0.0564639,0.0153666,0.900236,0.718289,0.665758,0.713403,0.285957,0.185811,0.361335,0.293796,0.0505286,0.365573,0.0467961,1237.88,1659.8,1215.39,710.076,1001.08,1167.57,1161.7,1258,506.187,910.725,2421.66,2367.79,1747.43,685.453,1643.12,2185.1,2165.69,919.254,1252.32,1151.1,1.518,-0.401637,8.32126,8.68989,8.52082,7.81896,8.38383,8.46937,8.29125,8.63259,7.73873,7.97067,9.02171,9.02032,8.69737,7.93322,8.84772,9.06229,9.08049,8.32126,8.50393,8.54768 +0,1,nan,1,1,0,0,20,2071,1213,2032,1776,1252,1414,1463,492,2455,689,1257,1286,625,1041,419,2305,1432,1876,733,1217,0.405632,0.967693,0.559558,0.519434,0.0490468,0.67918,0.937888,0.601061,0.139926,0.277525,0.0996672,0.795317,0.7662,0.951606,0.576368,0.0787189,0.879105,0.0729796,0.36986,0.329068,2248.41,1916.69,2480.32,2182.48,1050.64,1896.5,2269.11,653.138,2328.57,712.835,1188.68,1886.91,931.899,1747.18,510.632,2114.26,2199.74,1738.29,782.672,1286.53,0.162614,-0.684445,7.60296,7.05863,7.59577,7.49531,7.0862,7.24552,7.24783,6.23301,7.81985,6.54191,7.175,7.16096,6.47542,6.97705,6.00377,7.7652,7.25701,7.57332,6.57218,7.09709 +0,1,nan,1,1,0,0,20,4972,2599,6229,9671,23949,6731,23704,8655,7080,25276,20314,7323,24184,20640,11404,2714,24434,26529,13240,13908,0.373493,0.288499,0.291193,0.52083,0.842037,0.923357,0.961181,0.556488,0.620886,0.929656,0.886698,0.374971,0.822285,0.87098,0.671873,0.096518,0.957934,0.984322,0.704954,0.601297,1379.08,864.215,2059.38,1917.75,2279.14,534.455,1771.25,1592.63,1129.64,2026.35,1782.29,2007.74,2471.76,1877.71,1621.71,1409.45,1862.24,1886.49,1743.01,2282.31,0.45743,2.22334,8.517,7.86068,8.73501,9.17432,10.0611,8.79161,10.0739,9.06783,8.86752,10.1384,9.91452,8.89588,10.0983,9.93172,9.34247,7.92298,10.1168,10.1884,9.48815,9.52726 +0,1,nan,1,1,0,0,20,1092,818,524,764,517,463,732,191,522,290,534,453,669,316,541,234,325,557,798,438,0.194648,0.334521,0.34465,0.258072,0.677127,0.954013,0.37216,0.628786,0.62251,0.882827,0.584254,0.921892,0.701107,0.346831,0.0817019,0.65102,0.96946,0.780851,0.312601,0.793243,1966.99,1815.2,1230.38,1459.5,1620.91,2354.98,1638.1,526.895,1645.48,1326.61,1537.25,1891.42,2333.9,706.415,805.43,818.488,1565.43,2052.85,1726.49,1723.14,-0.343843,-1.25755,6.99564,6.73943,6.33782,6.61747,6.19538,6.22073,6.58944,5.13243,6.27911,5.73635,6.25918,6.04192,6.52978,5.7802,6.24479,5.54493,5.79293,6.30118,6.71689,6.11052 +0,1,nan,1,1,0,0,20,292,742,365,338,350,612,521,279,479,204,232,258,589,686,572,228,752,620,645,273,0.541701,0.21662,0.745382,0.939855,0.766218,0.580384,0.640727,0.292165,0.925369,0.961324,0.045072,0.1397,0.95681,0.0267728,0.935204,0.106021,0.441185,0.402812,0.507985,0.511399,922.801,2291.18,1122.16,1216.76,1128.35,2058.7,1700.57,875.201,1598.69,668.713,715.675,750.986,2139.6,2031.93,1958.71,680.389,2480.46,1934.13,2119.55,857.835,-1.15281,-0.063023,5.64047,6.57036,5.82323,5.8919,5.82741,6.44044,6.24553,5.60323,6.16582,5.29196,5.41758,5.45978,6.45527,6.46225,6.36829,5.36318,6.63559,6.38922,6.47414,5.56937 +0,1,nan,1,1,0,0,20,156,346,176,122,48,278,107,140,136,118,42,136,378,412,155,90,265,76,144,67,0.714888,0.253292,0.67711,0.322324,0.943431,0.48986,0.535242,0.0333192,0.435112,0.561044,0.770841,0.472364,0.287827,0.0606392,0.259426,0.565634,0.316429,0.884418,0.817882,0.934713,2311.03,1849.82,2293.04,883.222,901.054,2443.87,980.612,531.624,1077.4,1102.79,574.836,1205.88,2105.79,1616.14,888.596,933.869,1738.21,1376.68,2223.04,1296.32,-1.30293,-1.73425,5.20273,5.78064,5.26043,4.92166,3.86449,5.64887,4.657,4.91522,4.92478,4.72968,3.71432,4.97284,5.85035,5.9797,5.0368,4.55546,5.60892,4.3907,4.98529,4.24333 +0,1,nan,1,1,0,0,20,2839,1975,4617,863,2897,2446,3210,3793,2159,1288,2038,2973,3522,5534,2599,4040,1673,1310,1246,1936,0.0386616,0.154354,0.36429,0.684252,0.730749,0.0350048,0.856691,0.525144,0.31466,0.45809,0.763649,0.472726,0.395567,0.069763,0.874876,0.446831,0.150085,0.474475,0.970144,0.810379,1243.07,926.709,2480.12,590.961,2016.59,1096.85,2493.87,2237.67,1137.15,714.038,1449.48,1806.92,1963.4,2474.51,1965.17,2313.66,791.255,757.806,1010.62,1404.44,0.839915,-0.635436,7.94069,7.57347,8.4245,6.78687,7.98473,7.81787,8.11713,8.21941,7.67625,7.11976,7.63362,8.03891,8.17099,8.60938,7.86732,8.30257,7.41817,7.16884,7.14177,7.57237 +0,1,nan,1,1,0,0,20,244,854,1165,1349,1509,834,626,619,962,433,1172,1178,218,1226,229,1900,538,184,1494,1117,0.845421,0.112016,0.69006,0.260768,0.0595716,0.164388,0.3695,0.341829,0.0229098,0.594879,0.352966,0.388701,0.950042,0.340168,0.757614,0.150975,0.743249,0.968631,0.242343,0.531475,622.797,1166.9,2386.22,2100.18,1846.37,1106.67,996.185,1027.67,1154.81,848.417,1954.96,1919.53,523.207,1882.22,576.586,2436.09,1334.96,521.402,2088.53,2108.84,-0.157622,-0.8724,5.53905,6.80676,7.01783,7.26466,7.31138,6.70807,6.42396,6.47921,6.87408,6.06678,7.11258,7.06311,5.27354,7.08582,5.53856,7.50882,6.39062,5.25387,7.27517,7.03261 +0,1,nan,1,1,0,0,20,541,106,352,315,609,507,1063,805,183,1089,238,229,738,905,926,1331,539,618,312,844,0.22406,0.846915,0.575589,0.727692,0.92987,0.652203,0.421833,0.254651,0.745783,0.0471186,0.687394,0.963717,0.609098,0.272508,0.337515,0.241766,0.152904,0.904961,0.668324,0.256438,1011.79,519.545,967.977,936.538,2325.92,1532.49,2351.02,1452.92,590.564,1542.63,678.284,923.361,2095.76,1769.39,1876.87,2414.52,874.854,2393.15,902.536,1616.96,-0.32661,-1.1461,6.33607,4.95569,5.88892,5.68157,6.35954,6.26055,6.95253,6.66287,5.19973,6.96063,5.40513,5.39689,6.62297,6.83946,6.82393,7.18556,6.2722,6.41658,5.71263,6.76779 +0,1,nan,1,1,0,0,20,354,897,933,860,1242,196,758,993,430,783,366,1210,287,804,235,204,419,579,287,289,0.813235,0.40656,0.648275,0.486813,0.315005,0.897693,0.476831,0.316242,0.722584,0.781615,0.292832,0.32573,0.939333,0.0583231,0.653324,0.940356,0.841331,0.51211,0.56978,0.898714,1055.27,1951.22,2379.73,1821.79,2268.25,658.245,1864.84,1862.51,1305.11,2174.89,657.133,2381.53,912.074,1179.56,558.203,739.777,1356.15,1359.06,665.041,888.043,-0.361622,-0.946519,5.83019,6.82977,6.79951,6.68518,7.06698,5.27827,6.71798,6.86873,6.12848,6.5833,5.84909,7.10557,5.565,6.65607,5.34472,5.35466,6.05445,6.36821,5.59892,5.57675 +0,1,nan,1,1,0,0,20,179,150,601,432,531,313,160,772,363,394,543,619,260,160,404,191,326,573,340,514,0.28998,0.457523,0.716736,0.806416,0.920665,0.173062,0.391975,0.628685,0.451822,0.448918,0.569642,0.401996,0.0892001,0.105935,0.215479,0.465604,0.244993,0.819135,0.392844,0.827421,635.804,596.576,1735.41,1195.45,1287.61,1466.83,604.314,2476.48,1315.2,1486.69,1781.38,2483.97,1263.58,740.947,1908.81,759.417,1390.62,1536.13,1348.19,1373.79,-1.70613,0.891911,5.0074,5.09315,6.39214,6.09941,6.27557,5.73909,5.04758,6.6692,5.8786,5.99858,6.28709,6.47003,5.51514,4.99629,6.0403,5.3417,5.74989,6.36149,5.85077,6.25719 +0,1,nan,1,1,0,0,20,1279,2328,1990,2442,660,1073,1032,2227,859,735,2212,3105,1832,1166,1730,880,2368,1009,2714,1936,0.616276,0.688323,0.651943,0.566577,0.671402,0.585238,0.616201,0.625799,0.648034,0.463534,0.22056,0.849019,0.0286537,0.594697,0.0925145,0.239034,0.277688,0.831966,0.659392,0.622903,1137.56,1972.49,1711.76,2128.8,538.755,954.28,898.973,1908.86,768.02,655.559,2421.86,2361.62,2186.42,955.489,1967.35,937.998,2449.46,843.129,2306.02,1705.62,-0.179653,0.523354,7.17951,7.76764,7.60682,7.78018,6.46099,6.98759,6.94409,7.70212,6.80331,6.54843,7.72807,8.03179,7.52536,6.99381,7.45321,6.78919,7.7693,6.99288,7.90872,7.58803 +0,1,nan,1,1,0,0,20,954,530,2028,649,1122,1199,1512,756,1390,1408,779,947,864,620,1091,1111,1242,2276,674,693,0.140401,0.0274482,0.19641,0.258911,0.0112119,0.243648,0.270094,0.640717,0.687124,0.647701,0.130999,0.0567381,0.562602,0.823509,0.402627,0.919274,0.862565,0.0912952,0.689497,0.610787,963.393,539.244,2172,666.217,1149.27,1314.13,1633.38,899.221,1632,1624.49,775.518,909.01,965.886,779.283,1234.18,1422.17,1545.24,2232.58,839.559,858.72,0.0339215,-0.307443,6.86122,6.31565,7.65694,6.45594,7.07736,7.13994,7.34929,6.63847,7.22023,7.22774,6.64718,6.82883,6.734,6.43911,7.0283,7.01124,7.11166,7.71677,6.55482,6.60158 +0,1,nan,1,1,0,0,20,453,527,579,777,884,194,299,478,586,399,953,679,868,560,558,322,679,382,606,464,0.0820672,0.730608,0.0272608,0.231752,0.365555,0.912126,0.505664,0.683578,0.615139,0.646075,0.468961,0.108413,0.214686,0.641089,0.742928,0.972599,0.699753,0.426643,0.476298,0.775739,926.394,1833.89,1126.67,1764.72,2100.68,715.735,789.299,1294.74,1752.66,1097.97,2411.08,1361.94,1879.53,1643.34,1670.92,1140.72,2095.83,882.175,1611.25,1452.84,-0.670528,-0.634026,6.10874,6.38044,6.33921,6.65828,6.74772,5.32447,5.68001,6.06213,6.40835,5.92107,6.81997,6.4774,6.73213,6.32749,6.27957,5.75223,6.53352,5.84136,6.41225,6.11891 +0,1,nan,1,1,0,0,20,1719,3527,2472,1129,3158,2188,1808,3364,3325,1184,1487,2563,2980,3744,1367,4530,1572,4767,4698,3615,0.762739,0.352392,0.5012,0.0826345,0.6361,0.558698,0.887798,0.227606,0.0452816,0.705276,0.604945,0.94324,0.707242,0.748455,0.825976,0.524052,0.646019,0.705226,0.644336,0.893124,765.164,1837.82,1251.37,629.851,1490.7,1007.06,759.353,1835.27,1939.49,536.022,705.244,1043.94,1327.51,1666.54,615.787,2163.69,744.76,2169.11,2198.67,1539.5,0.517627,0.389061,7.45447,8.17107,7.84462,6.99526,8.07211,7.64978,7.4955,8.12113,8.10543,7.0762,7.31153,7.83536,7.98385,8.22733,7.26188,8.40108,7.38203,8.47408,8.46392,8.20432 +0,1,nan,1,1,0,0,20,7431,5935,13356,4609,9503,9986,4982,5055,7918,8674,5831,7581,7777,11757,7783,5515,3609,5996,10486,10658,0.902545,0.944,0.760401,0.463714,0.961644,0.943944,0.759569,0.878302,0.465959,0.91497,0.153685,0.504429,0.983068,0.448239,0.0815461,0.565156,0.0375136,0.166647,0.919678,0.338781,890.844,691.613,1815,842.15,1079.65,1138.7,689.844,620.807,1415.69,1031.42,1406.64,1320.58,871.202,2193.3,1987.76,903.155,960.317,1416.14,1244.98,2140.24,1.28721,0.921374,8.91096,8.69601,9.49166,8.45042,9.15763,9.19457,8.52352,8.52747,8.9719,9.06893,8.67777,8.9378,8.96285,9.39336,8.9571,8.61382,8.18904,8.69644,9.26145,9.26803 +0,1,nan,1,1,0,0,20,1163,994,1755,2206,1667,1325,4120,504,1011,2053,1323,5634,1039,2258,1301,823,724,4054,778,308,0.30439,0.0807929,0.48143,0.670393,0.337415,0.87188,0.925827,0.0517862,0.362957,0.874023,0.615888,0.980122,0.838979,0.981717,0.321917,0.060808,0.645391,0.896799,0.180711,0.0804142,1698.12,2254.58,1916.28,1655.25,2259.29,687.945,1966.34,1110.53,1279.26,1044.93,1165.05,2406.84,571.935,947.205,1867.77,1791.63,576.178,2009.63,1379.34,659.376,-0.925382,1.81517,7.06441,6.94199,7.50664,7.7032,7.40989,7.19094,8.33908,6.18121,6.88748,7.61283,7.25308,8.63978,6.94653,7.71012,7.19145,6.67588,6.60253,8.30817,6.632,5.71188 +0,1,nan,1,1,0,0,20,2909,4586,3722,2594,4473,2766,3701,3297,5600,6943,4789,2952,4013,3460,2252,3731,2875,6539,4406,6638,0.913951,0.239153,0.744684,0.821258,0.537099,0.545808,0.608822,0.215829,0.200733,0.249445,0.719564,0.844075,0.396276,0.698819,0.849924,0.419868,0.503009,0.280438,0.58404,0.364778,1510.82,1449.43,1761.57,1278.55,1747.15,1098.15,1601.38,1077.58,1749.47,2241.09,2187.07,1508.8,1425.04,1534.09,1164.94,1350.33,1091.98,2162.42,1796.1,2349.61,1.32729,-0.763953,7.94948,8.42352,8.23235,7.85337,8.38272,7.9117,8.2408,8.14488,8.64101,8.85145,8.46789,8.00153,8.28651,8.12912,7.73841,8.21464,7.93877,8.79203,8.37448,8.81062 +0,1,nan,1,1,0,0,20,1864,844,1009,1023,1290,2035,969,712,679,1818,565,1297,943,863,1418,1756,420,2040,1026,1580,0.262762,0.567183,0.14184,0.510855,0.811743,0.325913,0.166878,0.371032,0.183555,0.681595,0.861223,0.952088,0.125828,0.0499009,0.598734,0.525057,0.743983,0.568187,0.515129,0.0634032,2086.53,971.069,1183.71,1171.43,1585.08,2304.62,1065.49,759.319,782.408,2165.72,675.938,1612.8,1021.99,1005.54,1727.03,2000.85,515.093,2370.31,1206.88,1725.35,-0.0967947,-0.115744,7.51605,6.71595,6.9632,6.91005,7.17764,7.60815,6.85508,6.49268,6.54434,7.50482,6.31963,7.17873,6.81815,6.81071,7.28807,7.44376,6.06144,7.60822,6.93938,7.34905 +0,1,nan,1,1,0,0,20,490,1421,1671,1872,1317,866,1022,620,1229,801,2144,1740,1464,600,1702,682,768,735,501,1366,0.936357,0.295183,0.446483,0.424811,0.872958,0.189493,0.871156,0.958609,0.203574,0.301544,0.848728,0.546557,0.129516,0.19631,0.236336,0.583419,0.537759,0.171939,0.346599,0.249338,534.295,1914.6,2159.45,2387.8,1473.5,1124.9,1104.55,627.27,1609.58,1054.42,2375.89,2075.97,2019.01,823.459,2312.78,791.671,902.271,953.84,631.5,1736.89,-0.375002,0.336916,6.22142,7.28171,7.45304,7.54625,7.2145,6.71429,6.9257,6.38935,7.07731,6.68734,7.68407,7.44733,7.279,6.40465,7.45083,6.49571,6.61109,6.54342,6.18987,7.16885 +0,1,nan,1,1,0,0,20,434,277,234,738,493,161,298,167,248,197,477,467,351,692,311,359,658,158,258,220,0.138259,0.86238,0.30999,0.17072,0.408031,0.877984,0.184979,0.697913,0.398311,0.705871,0.694715,0.806909,0.600338,0.267107,0.948973,0.461634,0.420584,0.772916,0.471577,0.569967,875.311,1773.37,682.942,1509.25,1528.14,766.309,675.937,852.473,744.036,953.138,2204.88,2489.62,1319.49,1789.75,2104.33,1230.69,2054.4,941.697,1007.74,805.115,-0.545602,-1.38947,6.03687,5.73679,5.55009,6.53656,6.21926,4.87605,5.71348,5.23281,5.51305,5.33337,6.18754,6.15311,5.80525,6.57309,5.78759,5.92831,6.49775,5.22814,5.71463,5.35343 +0,1,nan,1,1,0,0,20,1803,2413,801,2230,3876,2753,2613,880,4400,3183,4145,1619,3143,1960,1498,2057,2265,3282,2960,801,0.847495,0.629479,0.881672,0.700546,0.290247,0.451036,0.644779,0.64034,0.0834584,0.318525,0.0987166,0.698976,0.328824,0.827806,0.836058,0.423394,0.52304,0.401067,0.166001,0.841919,2127.8,2039.98,1028.18,2089.89,1884.57,1738.27,2202.51,751.156,1470.26,1573.05,1462.29,1509.54,1645.36,2160.25,1714.02,1212.45,1550.59,1840.09,1156.45,955.036,1.20729,-1.60001,7.51413,7.82081,6.73214,7.73127,8.28435,7.94627,7.87299,6.80435,8.36695,8.05842,8.3371,7.40848,8.08688,7.56077,7.31618,7.63025,7.71681,8.08315,7.99479,6.72196 +0,1,nan,1,1,0,0,20,7332,9842,5981,2682,4386,3846,4770,3866,5224,6296,4205,8465,3918,3813,4347,5066,4793,5698,8790,2876,0.0235321,0.52178,0.0784249,0.0194006,0.23858,0.95092,0.97727,0.957704,0.0547994,0.812137,0.326011,0.788304,0.929242,0.498798,0.457443,0.503751,0.0668417,0.231915,0.899186,0.0937347,2126.9,2442.32,1676.6,797.261,1180.39,840.712,1042.74,847.01,1481.42,1422.43,1130.58,1971.02,839.072,935.231,1110.26,1249.2,1408.44,1572.37,1950.08,800.573,1.23373,0.311099,8.90348,9.19676,8.68266,7.92095,8.38156,8.26381,8.48737,8.27339,8.55154,8.74651,8.36564,9.06528,8.25512,8.2297,8.3884,8.52071,8.50477,8.66622,9.0891,7.94822 +0,1,nan,1,1,0,0,20,1638,3356,943,2117,3330,1657,826,545,846,1284,2450,2642,375,1123,584,1107,543,1228,2596,2381,0.321009,0.0942679,0.784887,0.138525,0.0323874,0.481924,0.63726,0.60396,0.558934,0.592813,0.227964,0.215162,0.568763,0.722858,0.370023,0.578076,0.356787,0.214412,0.0776184,0.263955,1582.44,2429.82,1738.31,1676.69,2187.46,1951.08,1287.49,804.862,1186.11,1810.69,2147.79,2201.27,515.585,1872.99,600.512,1532.95,575.864,1074.7,1834.24,2232.86,0.454103,-1.35822,7.38483,8.12164,6.84873,7.69053,8.10061,7.37568,6.74902,6.32446,6.77339,7.1504,7.81667,7.85866,5.9269,7.0076,6.34931,7.0039,6.32538,7.14268,7.86307,7.80663 +0,1,nan,1,1,0,0,20,940,1202,2443,1402,1133,1327,1280,606,836,588,1591,1740,554,1594,2320,1902,1336,2148,945,1233,0.61096,0.500377,0.618116,0.122257,0.176022,0.582552,0.752646,0.359693,0.560225,0.366051,0.349189,0.338931,0.211304,0.106894,0.705834,0.557702,0.146824,0.532347,0.397039,0.357239,823.229,1128.7,2024.99,1633.91,1368.93,1176.87,1045.39,639.589,717.16,598.141,1632.9,1858.91,612.756,2033.91,1899.21,1652.23,1594.56,1960.16,916.879,1186.95,-0.284575,0.729644,6.87444,7.10935,7.77975,7.20336,7.06564,7.21109,7.21673,6.4387,6.69949,6.37634,7.36832,7.49047,6.28757,7.41113,7.77963,7.53223,7.19691,7.68463,6.8261,7.05523 +0,1,nan,1,1,0,0,20,204,379,370,381,382,146,496,754,337,180,386,563,656,116,620,352,865,455,383,347,0.472443,0.86189,0.438211,0.144334,0.568086,0.432983,0.214439,0.306439,0.89489,0.806109,0.737024,0.763331,0.287949,0.834367,0.662172,0.406248,0.310153,0.366073,0.639523,0.174057,704.345,1768.15,1101.16,935.491,1400.17,503.339,1317.24,2134.2,1593.71,765.774,1700.28,2211.32,1877.96,514.977,2459.21,1054.39,2398.39,1361.95,1448.33,828.524,-0.755033,-0.931869,5.36198,5.91949,5.84073,5.95154,5.95994,5.06275,6.22843,6.62525,5.78487,5.13467,5.99671,6.23499,6.51458,4.71157,6.4355,5.82712,6.7385,6.12051,5.92718,5.80242 +0,1,nan,1,1,0,0,20,1469,1349,387,1015,625,482,440,899,581,715,378,1234,658,812,578,1116,687,1210,521,1033,0.932571,0.944271,0.336979,0.608479,0.0902825,0.0350474,0.237855,0.402199,0.285721,0.597577,0.0962241,0.3888,0.204741,0.367236,0.0263303,0.234307,0.368627,0.782925,0.878237,0.392138,2368.63,2263.28,680.769,1799.75,1394.47,1072.79,874.526,1762.9,1150.45,1362.5,805.061,2372.53,1293.65,1588.19,1215.91,2247.64,1372.35,1972.71,936.319,2008.59,-0.758592,0.291365,7.28319,7.24111,5.86282,6.9141,6.50798,6.22964,6.08439,6.83331,6.37257,6.6326,5.96036,7.1264,6.46628,6.71876,6.35233,7.02731,6.57309,7.05669,6.33925,6.96085 +0,1,nan,1,1,0,0,20,1006,1110,984,737,976,864,1838,432,1218,1324,2444,2223,2462,831,2279,1517,858,1496,1432,485,0.490527,0.493183,0.561018,0.818101,0.584241,0.973441,0.494102,0.702409,0.71825,0.517364,0.267357,0.200204,0.0361623,0.683539,0.350271,0.494304,0.608827,0.664155,0.689347,0.51584,1234.49,1356.62,1335.53,1318.28,1324.87,1814.36,2222.05,686.737,1944.69,1648.35,2338.19,1891.11,1808.1,1228.17,2374.53,1854.36,1219.69,2158.83,2144.27,566.544,0.368885,-1.14431,6.92598,7.01728,6.92399,6.6168,6.8894,6.75845,7.50966,6.09706,7.11984,7.18439,7.82008,7.68471,7.82754,6.69998,7.74062,7.32854,6.77855,7.28621,7.25061,6.11816 +0,1,nan,1,1,0,0,20,1428,915,1356,1315,1331,893,384,576,1453,831,558,1119,327,1281,1575,1354,931,1481,1510,825,0.209359,0.915524,0.160794,0.819128,0.803807,0.321043,0.353384,0.586163,0.584678,0.0404712,0.881924,0.562587,0.814296,0.587597,0.0809581,0.553374,0.47173,0.113072,0.206671,0.901613,2210.57,1754.88,1979.85,2349.48,2419.37,1453.79,565.751,899,2434,1262.72,1032.85,1885.46,599.402,2193.04,2339.85,2378.59,1561.93,2244.82,2297.94,1575.33,-0.365682,-0.276357,7.27746,6.85146,7.18066,7.16989,7.20344,6.82753,5.87481,6.27361,7.27003,6.76416,6.33067,7.02077,5.80521,7.16498,7.36979,7.25565,6.85763,7.31945,7.31697,6.74737 +0,1,nan,1,1,0,0,20,930,949,1318,502,588,288,490,895,533,880,1236,1110,1003,1249,362,1301,312,656,435,922,0.953426,0.543728,0.700277,0.943923,0.586137,0.643898,0.534349,0.704539,0.110691,0.482031,0.841018,0.673418,0.754122,0.853147,0.962589,0.446439,0.097771,0.183526,0.197904,0.108985,1732.94,1717.16,2400.83,886.372,1122.03,560.123,994.559,1631.87,1106.47,1610.57,2180.15,2124.38,1941.61,2296.06,663.764,2469.17,566.253,1259.69,827.574,1852.81,-0.710343,0.12668,6.86801,6.80696,7.16194,6.19637,6.3868,5.69938,6.25965,6.77639,6.31261,6.73506,7.08335,7.0362,6.95646,7.13668,5.90952,7.15785,5.64108,6.45153,6.03323,6.82792 +0,1,nan,1,1,0,0,20,3335,2776,1528,2058,4293,2568,2232,1821,1955,2152,3390,1147,1406,3915,3144,3213,4338,2033,832,1431,0.203435,0.621399,0.560647,0.602122,0.111456,0.681307,0.348563,0.148257,0.0400774,0.310926,0.0738642,0.413991,0.0386209,0.0126074,0.389139,0.815109,0.312206,0.341728,0.492654,0.119535,1618.71,1721.28,944.799,1299.15,1959.9,1712.71,1187.61,867.824,841.716,1078.3,1502.1,592.683,602.976,1665.91,1663.15,2318.77,2232.31,1040.55,505.787,659.63,0.873779,-0.679425,8.12494,7.90241,7.34383,7.63415,8.3787,7.85671,7.71666,7.53904,7.58199,7.64567,8.13821,6.97716,7.24942,8.28334,8.02586,8.06876,8.37245,7.58911,6.76517,7.28424 +0,1,nan,1,1,0,0,20,938,1085,1951,2246,2398,1309,2486,631,1848,1328,495,2097,1231,2475,2487,810,1786,2416,634,1733,0.800639,0.40872,0.573967,0.455093,0.357115,0.80532,0.383635,0.343021,0.368833,0.247426,0.89945,0.46843,0.724922,0.399037,0.0198634,0.336686,0.0726758,0.263784,0.225744,0.58067,922.858,1000.82,1918.42,2120.09,2182.42,1331.31,2452.02,602.554,1781.61,1231.29,513.763,2046.08,1246.17,2340.54,2170.45,744.017,1568.15,2247.07,574.144,1729.35,0.14992,-0.242484,6.78325,6.95939,7.57,7.69878,7.75152,7.14856,7.86156,6.46792,7.54576,7.20574,6.17358,7.66002,7.10197,7.8113,7.82779,6.68034,7.48995,7.80334,6.44806,7.46462 +0,1,nan,1,1,0,0,20,393,127,281,469,238,443,189,249,335,378,471,198,331,144,499,284,176,235,463,498,0.13052,0.326804,0.732552,0.454428,0.156875,0.279897,0.561504,0.945473,0.810501,0.360894,0.417175,0.17159,0.739988,0.315672,0.312419,0.343763,0.222217,0.755943,0.679849,0.358285,1654.24,580.85,1634.56,2331.64,1088.16,2018.93,858.132,1405.17,2058.59,1852.83,2427.49,859.379,1921.21,645.912,2320.23,1342.39,786.26,1272.6,2398.52,2301.73,-1.44053,-0.339026,5.92632,4.81317,5.71025,6.15974,5.49853,6.07491,5.12387,5.48685,5.91447,5.96159,6.21265,5.25751,5.86931,4.92312,6.20298,5.64514,5.15142,5.45201,6.11159,6.17942 +0,1,nan,1,1,0,0,20,271,368,604,343,400,344,190,139,579,246,247,320,352,445,279,246,359,418,575,263,0.388268,0.871498,0.367487,0.979364,0.173296,0.907038,0.761908,0.236629,0.129154,0.434893,0.83487,0.951357,0.340338,0.468607,0.611264,0.966679,0.73414,0.395804,0.405275,0.98097,1099.08,2119.77,2474.86,2013.17,1466.5,1854.72,1087.07,522.751,1980.73,964.431,1500.3,1991.38,1362.04,1918.18,1376.23,1522.41,1861.14,1723.44,2283.14,1556.24,-1.21469,-0.616657,5.54811,5.90696,6.37263,5.78885,5.96908,5.75147,5.30671,4.8985,6.29689,5.38867,5.5839,5.79523,5.79218,6.05547,5.63547,5.51725,5.86154,5.99331,6.2687,5.53042 +0,1,nan,1,1,0,0,20,890,893,2578,856,1393,2174,842,1091,2259,1160,637,2155,1856,1756,1819,1583,604,2009,849,1586,0.538395,0.898348,0.0437401,0.589493,0.0229319,0.193809,0.967603,0.0657352,0.394098,0.523444,0.92974,0.311173,0.913195,0.98755,0.924894,0.150735,0.939381,0.367331,0.597148,0.513026,890.61,999.146,2207.28,862.149,1161.68,2009.9,927.229,897.565,2223.67,1181.01,734.819,2097.5,2068.53,2063.44,2162.6,1470.28,744.891,1912.23,903.967,1601.07,0.154158,-0.299485,6.78482,6.79202,7.84058,6.73704,7.20491,7.70196,6.69658,6.93416,7.74305,7.07152,6.47534,7.70947,7.51526,7.49053,7.55623,7.40222,6.48607,7.60017,6.78211,7.37894 +0,1,nan,1,1,0,0,20,2255,12543,2726,2326,3415,2688,1314,1070,5847,2464,1168,3498,2194,1943,6880,3982,4842,7747,1057,1854,0.255247,0.86216,0.153828,0.695444,0.0364358,0.852118,0.159694,0.214627,0.72145,0.336051,0.0867557,0.67457,0.425662,0.0415563,0.61684,0.384397,0.72708,0.799594,0.0552253,0.0150913,1128.85,2334.12,1609.87,565.653,2462.26,502.954,742.095,590.051,1356.29,1089.26,756.9,892.298,831.175,1336.72,1910.73,1627.44,1120.9,1597.81,732.536,1323.65,0.281034,1.61898,7.72323,9.43224,7.91399,7.74492,8.14886,7.88109,7.14905,7.00872,8.66156,7.81835,7.05072,8.16695,7.69301,7.54629,8.83493,8.29813,8.48004,8.95195,6.96696,7.49362 +0,1,nan,1,1,0,0,20,1769,1171,1262,592,1715,736,2160,714,1440,922,580,971,2021,326,2613,495,1245,3164,2316,689,0.581688,0.130793,0.18882,0.438954,0.36555,0.61216,0.23239,0.648676,0.169667,0.912111,0.733361,0.705218,0.435997,0.812778,0.0878077,0.716377,0.0531729,0.0485505,0.340922,0.20484,2439.3,848.047,994.247,610.549,1759.54,1019.85,1729.77,1043.76,1093.91,2087.37,1012.65,1528.31,2352.64,599.498,1780.23,850.972,820.493,2041.12,2225.83,573.567,0.49728,-1.39694,7.48417,7.05751,7.1355,6.29845,7.45944,6.56955,7.62839,6.5417,7.25778,6.86678,6.39315,6.84405,7.65151,5.75797,7.85912,6.24293,7.13291,8.05071,7.72892,6.56301 +0,1,nan,1,1,0,0,20,266,145,99,174,341,266,85,204,102,274,177,144,82,134,88,416,223,227,146,282,0.851922,0.439817,0.691636,0.378209,0.152311,0.809049,0.557437,0.601162,0.47726,0.455553,0.270824,0.352884,0.550973,0.538337,0.986857,0.0934999,0.9285,0.966322,0.34938,0.333169,1725.6,930.488,794.746,1094.22,1915.03,1673.73,614.964,1330.72,537.269,1478.71,911.339,940.864,548.11,917.781,606.427,2281.5,1594.3,1585.65,884.146,1638.32,-1.71192,-0.215313,5.55798,5.02909,4.81719,5.20445,5.81278,5.53669,4.58962,5.35212,4.47182,5.48892,5.04469,5.0589,4.47593,4.99413,4.48318,6.00054,5.46235,5.44877,4.99748,5.61777 +0,1,nan,1,1,0,0,20,1259,3225,3942,1951,4524,2397,2625,2627,3588,4139,3669,1424,2229,4237,4773,2481,4344,2907,2280,3080,0.834391,0.0668328,0.0480156,0.092594,0.0801366,0.809104,0.014034,0.83052,0.429682,0.750145,0.320479,0.364839,0.357439,0.369114,0.0753914,0.537096,0.725116,0.801124,0.61521,0.140444,741.058,1426.61,1765.52,860.495,2062.84,1323.48,1205.26,1548.09,1811.76,2288.4,1811.87,717.03,1123.52,2029.79,2174.61,1276.64,2483.31,1639.01,1223.3,1439.1,0.805065,-0.295862,7.16628,8.04835,8.26706,7.53518,8.4132,7.7537,7.89537,7.90412,8.17999,8.31873,8.21236,7.27224,7.72353,8.31155,8.46737,7.79814,8.40788,7.96989,7.73235,8.03528 +0,1,nan,1,1,0,0,20,3927,5913,1948,3229,5135,8074,2961,3195,4010,2110,1078,2532,617,1711,2298,2443,7122,2332,4606,5057,0.659355,0.62989,0.954099,0.224752,0.720005,0.858308,0.488143,0.833262,0.467075,0.611745,0.290778,0.371003,0.0911802,0.189974,0.131405,0.246728,0.781733,0.20823,0.568493,0.666773,1575.74,2422.8,523.175,2290.32,1892.3,2451.1,1450.26,1007.54,2125.71,897.356,704.257,1529.76,521.824,1295.35,1860.24,1753.86,2429.03,1669.76,2131.59,2046.59,0.0346177,1.33816,8.27942,8.67019,7.57127,8.07182,8.54365,8.98746,7.96733,8.06492,8.3215,7.65268,6.98087,7.86395,6.41396,7.45537,7.73892,7.83435,8.87595,7.7337,8.45998,8.5508 +0,1,nan,1,1,0,0,20,1037,547,697,685,315,475,988,416,726,990,792,1074,1301,565,1244,580,302,833,296,422,0.050952,0.939714,0.875461,0.841528,0.770046,0.976107,0.792871,0.735619,0.657713,0.371945,0.596124,0.165866,0.955511,0.177819,0.732014,0.967561,0.361598,0.681857,0.589335,0.15973,1956.01,973.591,1185.32,1219.5,579.756,903.926,1742.52,866.759,1330.1,1734.3,1395.4,1979.56,2313.03,1077.57,2255.11,1051.07,581.908,1573.78,546.382,770.876,-0.605037,0.00547661,6.9739,6.2811,6.47753,6.50576,5.76179,6.20706,6.86239,6.16375,6.59157,6.85536,6.63916,6.9865,7.14651,6.3784,7.11993,6.35783,5.76326,6.75993,5.70151,6.04337 +0,1,nan,1,1,0,0,20,4918,12236,12244,7225,1511,8366,6848,9145,2971,8252,3790,5964,9310,4220,12315,9586,7884,18344,3510,3345,0.716487,0.926442,0.759263,0.26822,0.138514,0.477458,0.876789,0.514861,0.179035,0.436451,0.366697,0.737101,0.609908,0.585134,0.968581,0.327877,0.342298,0.797343,0.507373,0.315676,575.278,960.185,1297.99,1936.42,529.969,1525.08,597.923,1531.8,923.431,1623.31,848.018,659.896,1339.62,625.999,890.526,2283.65,1816.5,1813.42,621.726,807.135,0.821789,1.85696,8.50713,9.40928,9.40028,8.88846,7.35182,9.03821,8.84341,9.11207,7.98235,9.02448,8.24563,8.68264,9.1545,8.34771,9.41222,9.16418,8.96209,9.80539,8.19646,8.10148 +0,1,nan,1,1,0,0,20,2773,666,5476,1823,1782,5958,1429,904,4843,1020,2908,1862,3214,2489,1486,3473,4626,901,2209,2632,0.455172,0.276031,0.665805,0.0960374,0.204903,0.901528,0.271204,0.366709,0.672626,0.124881,0.755282,0.127593,0.503242,0.635484,0.356902,0.517858,0.866161,0.313752,0.117754,0.418716,1775.41,568.744,2372.96,2121.52,1642.27,1765.6,1183.3,672.676,2060.27,1089.81,1088.49,1980.03,1840.94,1150.68,1107.86,1915.85,1460.39,709.163,2451.18,1744.15,-0.289506,1.66958,7.95223,6.51478,8.594,7.53072,7.45643,8.69191,7.23935,6.83401,8.46409,6.91275,7.96405,7.51439,8.06873,7.8196,7.31656,8.13302,8.44308,6.79841,7.71142,7.8736 +0,1,nan,1,1,0,0,20,267,647,948,282,192,326,411,228,258,272,351,364,412,557,311,487,352,263,239,59,0.482609,0.137508,0.102111,0.712786,0.54833,0.495283,0.378577,0.957883,0.819317,0.15206,0.59224,0.0156632,0.126134,0.167535,0.354497,0.428903,0.48478,0.660854,0.587248,0.902566,1306.1,1781.6,2310.89,1929.69,1013.46,1567.36,1540,2299.47,2200.97,727.777,2117.59,811.208,1091.6,1577.38,1377.91,2114.73,1702.33,1713.92,1487.26,507.436,-0.792539,-1.66397,5.57922,6.46392,6.78294,5.58652,5.21619,5.74048,5.91706,5.35401,5.5408,5.54443,5.88003,5.87992,5.99298,6.29221,5.84591,6.15046,5.84056,5.55436,5.53499,3.93499 +0,1,nan,1,1,0,0,20,880,890,769,1330,1358,899,500,1853,1085,1483,1290,2642,1636,943,1476,1337,1357,1281,1430,483,0.604493,0.304612,0.535654,0.770047,0.721813,0.706561,0.2934,0.856657,0.169375,0.411299,0.244081,0.952887,0.515311,0.604367,0.89777,0.680144,0.791823,0.38077,0.945254,0.462562,1123.78,1531.67,1059.48,1491.28,1613.01,1102.05,920.875,1902.82,2243.65,2347.74,2357.65,2476.5,2291.71,1217.98,1453.16,1666.34,1536.65,2023.15,1482.09,760.614,-0.86847,0.970019,6.74235,6.76112,6.61666,7.18588,7.21756,6.82184,6.24146,7.51359,7.01169,7.29171,7.13371,7.87045,7.36844,6.82273,7.28388,7.20967,7.23698,7.11329,7.34965,6.21435 +0,1,nan,1,1,0,0,20,1997,541,507,862,1441,607,560,1387,624,1955,1269,2112,1497,1223,473,480,536,1600,2130,2096,0.747692,0.725219,0.791756,0.943705,0.0307617,0.88913,0.348304,0.35561,0.216765,0.900877,0.483698,0.785728,0.0521466,0.219534,0.491915,0.977531,0.732091,0.55571,0.484304,0.292757,2118.27,581.365,516.028,908.101,1704.28,637.679,661.051,1568.48,715.146,2008.72,1389.97,2307.64,1720.64,1394.88,555.728,500.488,573.717,1751.07,2402.82,2354.01,-0.164993,0.114205,7.57875,6.28321,6.17159,6.75414,7.27942,6.39439,6.36861,7.23348,6.43225,7.54314,7.12728,7.66872,7.29141,7.10064,6.21147,6.16223,6.27075,7.36645,7.67471,7.63232 +0,1,nan,1,1,0,0,20,2757,2764,1734,555,219,1965,627,402,1489,540,2940,1126,2053,3312,761,1222,208,396,1626,1159,0.0189587,0.197048,0.1412,0.439175,0.981124,0.41347,0.774754,0.591979,0.435957,0.818953,0.249405,0.518241,0.475995,0.203617,0.544368,0.546047,0.899203,0.604869,0.2804,0.388554,1034.19,1694.81,897.828,609.776,981.43,2091.72,1596.73,711.126,1545.49,1555.94,2013.27,1588.77,2483.02,2034.28,1164.07,1840.63,744.547,695.316,1149,1090.48,1.02819,-2.58584,7.92054,7.95398,7.46304,6.30564,5.38017,7.60476,6.40051,6.06427,7.24397,6.26034,7.99078,7.05882,7.61457,8.11956,6.68021,7.13406,5.31577,6.00846,7.34976,7.01782 +0,1,nan,1,1,0,0,20,276,554,1170,642,504,414,1549,636,534,1017,1119,821,737,576,438,1358,920,912,738,1069,0.449243,0.196351,0.0428725,0.722847,0.857804,0.512721,0.0191991,0.389847,0.894317,0.0471885,0.172323,0.12685,0.629847,0.843444,0.28759,0.139849,0.289047,0.509158,0.545709,0.396757,660.776,910.769,1666.5,2260.34,1941.66,962.835,2080.51,1379.79,2298.48,1474.2,1847.31,1363.99,2161.88,2237.11,816.782,2199.71,1714.32,2233.36,2051.62,2369.09,-0.27402,-1.30462,5.63331,6.28411,7.08853,6.50621,6.17817,5.92696,7.3413,6.44706,6.29924,6.96029,7.02265,6.77866,6.58301,6.33855,6.05616,7.23961,6.79566,6.77299,6.64043,6.97863 +0,1,nan,1,1,0,0,20,436,693,113,682,325,143,150,539,509,394,463,165,347,560,434,304,454,499,100,650,0.624559,0.845676,0.165432,0.790934,0.615191,0.127033,0.545312,0.490399,0.875944,0.496787,0.810969,0.275675,0.210007,0.728553,0.670005,0.284696,0.700243,0.459326,0.285206,0.812957,1679.76,2410.27,742.961,2387.95,1385.11,832.544,584.422,2437.56,1680.84,1857.08,1671,809.047,2187.42,2124.72,1807.49,1899.19,1809.72,2199.86,507.717,2248.31,-1.96789,0.891268,6.01517,6.57333,4.7902,6.51524,5.81394,4.86982,4.88876,6.26794,6.23986,6.00164,6.17608,4.97367,5.90976,6.34285,6.12896,5.83504,6.15714,6.13764,4.51623,6.47461 +0,1,nan,1,1,0,0,20,5979,7118,4631,6993,4722,5437,4430,5932,4681,5814,2361,4666,7939,5795,3179,5666,5965,2875,3321,4299,0.861806,0.333975,0.852047,0.721153,0.684412,0.936233,0.20495,0.981528,0.434631,0.740947,0.404035,0.0869599,0.239591,0.46754,0.740458,0.586114,0.928437,0.2069,0.879387,0.660977,2285.76,2169.24,1702.97,2490.4,1655.36,2059.43,1294.85,2310.09,1492.7,2072.58,739.703,1274.22,2319.2,1876.06,1156.59,1880.82,2287.28,803.27,1243.09,1514.32,1.32516,-0.395791,8.71852,8.87511,8.42806,8.85993,8.46605,8.58479,8.41019,8.68172,8.46148,8.66845,7.7715,8.44084,8.97931,8.67704,8.08532,8.63265,8.69281,7.93196,8.10246,8.38628 +0,1,nan,1,1,0,0,20,4301,1569,859,2199,3318,869,1967,2256,1542,3797,694,749,2502,2109,3604,1072,4737,2764,1779,883,0.926991,0.596126,0.0899479,0.92268,0.977385,0.439769,0.550416,0.441229,0.51806,0.907132,0.114911,0.385171,0.699495,0.475903,0.876189,0.138562,0.723551,0.762604,0.300989,0.526599,1606.13,1032.37,1362.16,815.17,1115.77,722.391,1350.61,1890.34,1171.03,1449.45,1059.96,702.583,1365.42,1677.59,1455.98,1579.72,2457.87,1392.07,1940.84,651.57,-0.612766,1.74601,8.38735,7.36769,6.76111,7.70164,8.11105,6.73764,7.55658,7.70213,7.35741,8.25004,6.55386,6.61451,7.82777,7.64328,8.2005,6.99417,8.45761,7.9573,7.48364,6.78606 +0,1,nan,1,1,0,0,20,900,6173,1634,1050,3730,957,2410,3658,3378,4225,3058,8365,2289,2177,3456,3768,2295,1978,5377,4622,0.9081,0.345745,0.879187,0.77818,0.618751,0.649556,0.408293,0.659604,0.595257,0.264726,0.648071,0.144167,0.764808,0.78294,0.594,0.580476,0.32042,0.584046,0.300727,0.551905,859.21,2395.22,1480.96,830.833,2327.15,631.381,1044.38,2353.01,1904.78,1468.4,1999.75,2352.89,1712.04,1773.78,1953.64,2193.72,840.974,1099.71,1957.71,2493.94,1.49575,-1.62095,6.77977,8.71654,7.37108,6.95678,8.24518,6.89076,7.7851,8.19001,8.08299,8.35857,8.04604,9.02546,7.70147,7.70751,8.11035,8.24818,7.71092,7.55184,8.58781,8.42275 +0,1,nan,1,1,0,0,20,322,400,689,332,487,569,523,336,362,248,787,148,600,224,215,325,325,537,959,373,0.287443,0.365837,0.729732,0.556504,0.393446,0.425746,0.383036,0.740244,0.721807,0.558413,0.737721,0.125033,0.760587,0.352248,0.50268,0.320694,0.554022,0.972884,0.932419,0.0610205,1351.47,1492.22,1793.68,1034.69,1684.74,2174.7,2001.54,909.204,1000.14,815.667,2044.82,688.629,1614.99,839.168,767.509,1313.56,976.177,1147.62,2159.93,2030.24,-1.73865,1.01462,5.76194,5.94056,6.49378,5.76785,6.08992,6.37797,6.25166,5.82499,5.9016,5.53194,6.63293,4.92292,6.42015,5.35116,5.41453,5.76723,5.70712,6.29391,6.88524,5.93917 +0,1,nan,1,1,0,0,20,393,755,1216,607,620,495,359,345,540,452,572,468,905,301,445,782,899,206,509,1469,0.459275,0.839402,0.401578,0.26445,0.451178,0.402059,0.0829979,0.734222,0.54869,0.980772,0.743945,0.369032,0.857517,0.884212,0.236928,0.420386,0.016725,0.928641,0.0273498,0.0536907,830.436,2066.83,2299.46,1081.26,1179.2,906.144,518.701,871.677,1186.56,1321.68,1365.37,860.712,2399.73,871.751,828.344,1640.8,1274.05,628.691,755.416,2070.07,-0.341963,-0.767635,6.02743,6.64746,7.0902,6.44092,6.38428,6.1586,5.84565,5.86484,6.31565,6.09182,6.30614,6.13251,6.78289,5.74979,6.19559,6.73827,6.79516,5.38882,6.26431,7.25216 +0,1,nan,1,1,0,0,20,1239,1017,920,313,688,747,415,644,1346,441,1013,975,728,447,1114,325,857,1246,868,420,0.665968,0.648821,0.598085,0.939177,0.0746171,0.753718,0.259253,0.594773,0.687996,0.293386,0.427662,0.687119,0.584859,0.0439918,0.513896,0.442795,0.0524387,0.583998,0.501448,0.638023,2357.23,1951.91,1637,660.669,1269.03,1421.31,756.604,1263.91,2469.82,821.527,1878.71,1787.01,1405.76,892.611,2261.81,672.377,1595.98,2326.4,1602.39,797.222,-0.608008,-0.0641436,7.11452,6.92694,6.75425,5.825,6.53322,6.60298,6.0042,6.49581,7.15976,6.08434,6.9029,6.83622,6.60281,6.18332,7.08295,5.87441,6.76387,7.10661,6.73908,6.0322 +0,1,nan,1,1,0,0,20,1357,1300,2355,983,406,502,417,1125,755,1048,1340,596,1479,501,999,798,875,889,1287,1379,0.716763,0.524996,0.163899,0.5417,0.741099,0.616259,0.45033,0.328314,0.42033,0.088978,0.509539,0.0695294,0.676304,0.980129,0.731464,0.570085,0.667578,0.0354951,0.541628,0.385966,2102.75,1754.08,2449.79,1389.12,639.898,744.799,526.302,1321.41,942.712,1021.24,1836.16,563.89,2271.07,964.277,1606.64,1195.99,1251.68,847.822,1824.32,1631.44,0.0735644,-0.727158,7.20336,7.16151,7.75814,6.91609,5.99598,6.23856,6.01198,7.02128,6.61668,6.93763,7.21848,6.35786,7.30979,6.23223,6.92358,6.74575,6.72037,6.79042,7.18868,7.19012 +0,1,nan,1,1,0,0,20,155,701,295,498,731,370,491,177,240,367,174,242,597,139,248,441,745,428,206,224,0.294196,0.147028,0.316184,0.578899,0.222584,0.197085,0.0595372,0.363917,0.69153,0.790576,0.609452,0.346054,0.398417,0.47363,0.298412,0.45562,0.0463764,0.402425,0.915314,0.495019,583.781,1954.96,988.207,2397.96,2253.49,1077.69,1158.68,688.377,1470.55,1995.25,878.092,860.626,2143.59,586.559,749.074,1768.6,1858.73,1598.37,1478.53,972.537,-0.805489,-1.28741,5.18529,6.58335,5.68334,6.2316,6.62819,5.92335,6.1729,5.26034,5.59762,5.77524,5.18765,5.50666,6.35182,4.95903,5.42917,6.08588,6.66245,6.05316,5.31493,5.43713 +0,1,nan,1,1,0,0,20,1147,2173,1028,1335,4228,1680,2713,3246,3876,3975,2673,2371,1087,1657,3491,901,3518,1136,3121,2611,0.576117,0.232933,0.443324,0.726366,0.964935,0.989595,0.973579,0.609228,0.563186,0.495564,0.269662,0.743246,0.556633,0.628241,0.483128,0.438453,0.46888,0.1315,0.662748,0.0499491,613.351,1664.27,670.537,652.688,1627.93,627.848,1029.73,1676.69,2126.27,2300.86,1929.41,1088.99,605.436,890.539,2037.64,555.688,2137.97,965.991,1632.38,2361.41,0.0753678,0.920481,7.02461,7.70692,6.99152,7.22507,8.35864,7.42857,7.90858,8.06072,8.25589,8.27257,7.88855,7.75251,6.99369,7.44548,8.13963,6.79916,8.17458,7.06957,8.08321,7.88836 +0,1,nan,1,1,0,0,20,1805,1690,1874,1828,880,841,1329,2265,1976,826,609,2119,1509,2235,2173,639,433,1046,1044,973,0.981781,0.0154689,0.520448,0.0255511,0.307084,0.09069,0.428129,0.177018,0.340316,0.275551,0.195267,0.178243,0.211375,0.241748,0.153967,0.77614,0.812879,0.247368,0.703974,0.396234,2202.43,1406.54,1868.95,1427.2,772.055,724.911,1294.44,1925.29,1807.06,762.983,565.503,1819.24,1340.51,1989.91,1908.8,701.431,503.938,902.297,1156.8,926.854,0.212245,-0.447421,7.47029,7.45421,7.51252,7.46428,6.72391,6.75772,7.18652,7.69588,7.55944,6.72619,6.46259,7.63867,7.31847,7.69992,7.69759,6.41811,6.071,6.90651,6.95068,6.86676 +0,1,nan,1,1,0,0,20,7070,2211,2258,5778,1454,4449,6251,5145,1486,7772,7129,4788,4499,3123,6175,3490,5211,4541,7566,6014,0.702547,0.0802025,0.907944,0.813971,0.799993,0.552059,0.178871,0.0583901,0.661903,0.409324,0.315107,0.169594,0.586037,0.407915,0.395885,0.174647,0.195603,0.316398,0.341669,0.734988,2483.51,642.667,838.829,2097.73,518.754,1504.78,1838.47,1440.59,523.354,2494.51,2179.26,1408.37,1506.96,1016.38,1936.59,1020.16,1516.85,1416.67,2363.92,2189.63,1.29019,-0.358963,8.85543,7.72703,7.69628,8.64662,7.25445,8.40843,8.74267,8.54204,7.31285,8.96511,8.86382,8.4795,8.39768,8.06777,8.71677,8.15521,8.54437,8.43268,8.93562,8.71785 +0,1,nan,1,1,0,0,20,441,637,1299,823,1046,1163,988,1064,616,1286,1040,1090,1450,896,346,1378,677,358,1274,1681,0.679757,0.716009,0.330443,0.789177,0.924543,0.533125,0.804684,0.443649,0.832082,0.582593,0.870992,0.860963,0.436355,0.420154,0.756392,0.24236,0.877775,0.517232,0.303937,0.05068,891.65,1241.62,1946.86,1579.16,2411.69,1950.7,1998.19,1671.63,1357.71,2249.11,2242.7,2335.96,2284.21,1472.95,660.457,2067.33,1526.35,608.107,1874.73,2148.3,-0.162978,-0.692241,6.15954,6.46554,7.18225,6.65537,6.9851,7.04391,6.87998,6.95146,6.47457,7.15202,6.94952,6.99721,7.26873,6.8412,5.80635,7.30326,6.56002,5.88932,7.16284,7.47437 +0,1,nan,1,1,0,0,20,823,348,871,676,720,351,421,592,848,311,812,573,1053,776,1061,637,725,1208,524,441,0.148811,0.0153398,0.0666096,0.190806,0.196127,0.545081,0.595958,0.29808,0.716506,0.493668,0.887178,0.575338,0.321742,0.433945,0.193928,0.920517,0.186178,0.182789,0.690256,0.819047,1740.36,749.286,1810.42,1388.62,1597.54,812.395,901.891,1315.03,2020.3,611.681,1943.91,1341.88,2380.43,1677.25,2302.5,1494.74,1547.35,2443.76,1249.87,1115.07,-0.722564,-0.166943,6.71444,5.894,6.76763,6.48165,6.62091,5.88643,5.98244,6.40929,6.76882,5.61123,6.70178,6.38321,6.99876,6.6299,6.98681,6.43347,6.59065,7.04822,6.293,6.15738 +0,1,nan,1,1,0,0,20,3056,5463,5879,10890,6194,6790,6481,3040,6197,9117,6353,5569,13034,7894,6354,2457,9975,7336,7187,2272,0.720784,0.36404,0.442426,0.638817,0.564133,0.488917,0.504577,0.743072,0.309944,0.738598,0.663616,0.27674,0.895615,0.864772,0.720202,0.544584,0.532854,0.978004,0.4205,0.334628,633.706,1546.93,1529.03,2407.69,1469.54,1719.16,1592.16,629.407,1786.72,1818.43,1358.41,1672.63,2292.88,1445.39,1322.76,583.838,2387.04,1211.31,1891.68,665.588,0.971908,0.844136,8.03193,8.62324,8.67776,9.29758,8.74082,8.83421,8.77069,8.04394,8.72168,9.10111,8.74616,8.62766,9.46549,8.97803,8.76733,7.80123,9.19952,8.89694,8.87209,7.75505 +0,1,nan,1,1,0,0,20,3910,4574,3317,7366,5649,2761,2758,3486,4271,1878,4549,3727,2606,7290,2027,3146,4163,3024,1982,4122,0.25919,0.306778,0.300721,0.700001,0.733291,0.105456,0.327271,0.0705159,0.424375,0.123746,0.798671,0.178104,0.425166,0.606696,0.477799,0.0598027,0.419543,0.901112,0.666778,0.335234,1802.77,1981.93,1458.39,1767.41,1281.16,1586.91,1196.35,2217.65,1550.28,1097.38,944.226,1956.52,958.625,2092.94,682.31,2029.5,1568.58,546.067,511.41,1749.53,0.366811,1.50246,8.25331,8.41956,8.10372,8.89581,8.62407,7.8948,7.94555,8.17696,8.35061,7.55341,8.41715,8.21333,7.8711,8.92467,7.61017,8.07221,8.35508,8.02343,7.60579,8.33759 +0,1,nan,1,1,0,0,20,1326,414,1015,299,459,385,598,611,114,169,272,150,450,842,1156,492,536,127,66,1093,0.273072,0.643342,0.382165,0.875952,0.179883,0.673926,0.481009,0.395601,0.829372,0.880625,0.819436,0.842249,0.289016,0.113232,0.289215,0.547142,0.320325,0.928827,0.798903,0.296432,1936.67,1676.34,2055.69,2375.34,556.578,1705.48,1594.16,1325.03,848.033,1611.64,1785.14,993.303,747.739,873.763,1990.47,1770.95,984.88,1164.61,506.075,1734.14,0.307002,-2.7124,7.13505,5.98637,6.89878,5.70397,6.1409,5.92065,6.37642,6.42317,4.80034,5.3034,5.57162,4.92352,6.14013,6.77268,7.11866,6.3022,6.33067,4.84779,4.36674,6.96123 +0,1,nan,1,1,0,0,20,676,944,955,1892,1170,1131,394,410,722,2033,2229,932,1234,1864,1648,722,962,746,1812,1344,0.729291,0.696522,0.225057,0.28201,0.501611,0.133899,0.558823,0.53544,0.938566,0.103158,0.0186089,0.624619,0.506183,0.310939,0.490883,0.072591,0.826483,0.716329,0.183153,0.585504,1115.47,1581.49,1176.89,2478.64,1768.68,1262.53,620.869,614.464,1544.05,2210.39,2334.3,1577.91,1829.81,2391.02,2364.66,791.176,1710.21,1259.12,2072.03,2100.45,-0.010293,-0.757523,6.45428,6.8282,6.88986,7.59154,7.08771,7.02915,5.99751,6.00485,6.62088,7.61248,7.73108,6.8804,7.11823,7.53364,7.38624,6.60824,6.808,6.58524,7.48725,7.19608 +0,1,nan,1,1,0,0,20,915,458,711,523,1210,485,706,465,463,638,386,350,262,879,744,295,301,754,417,1048,0.714084,0.977042,0.115814,0.185564,0.225641,0.696142,0.916782,0.968571,0.777952,0.52012,0.875689,0.667238,0.29458,0.282867,0.201273,0.533997,0.938047,0.320209,0.0790342,0.625087,2436.99,1328.73,1423.23,1061.43,2428.18,1171.73,1896.02,1290.97,1252.62,1488.01,1046.5,954.31,513.984,1939.44,1464.56,730.127,830.177,1646.88,827.997,2465.02,-0.640757,-0.398822,6.87297,6.16156,6.57374,6.25261,7.06415,6.14784,6.54112,6.1361,6.18197,6.45701,5.96321,5.95412,5.48395,6.81658,6.56828,5.73949,5.70677,6.63818,6.04673,6.9199 +0,1,nan,1,1,0,0,20,689,1955,665,1276,2082,1740,1054,1945,713,2852,1544,857,3022,1699,2559,1381,1466,1096,1279,1037,0.627684,0.90328,0.253809,0.977633,0.0170037,0.629557,0.546779,0.244321,0.0140161,0.808896,0.365161,0.337534,0.704362,0.195574,0.482791,0.607629,0.210785,0.898315,0.644844,0.328176,573.284,1394.92,689.215,867.388,2390.31,1458.2,916.771,1982.48,866.692,2135.66,1516.9,822.09,2313.52,1717.09,2366.58,1136.74,1570.07,819.873,1104.62,914.005,-0.158642,0.5569,6.5423,7.58499,6.51826,7.15129,7.63001,7.47691,6.96672,7.56953,6.61385,7.95836,7.36914,6.74118,7.98014,7.39866,7.87942,7.21567,7.31762,7.05078,7.20773,6.84196 +0,1,nan,1,1,0,0,20,4537,6842,11872,3744,4635,2603,1845,7095,2195,13020,5223,1459,6689,2373,3369,15492,2533,2482,2129,3183,0.279044,0.810962,0.979848,0.525383,0.932881,0.68995,0.172076,0.591649,0.66393,0.883299,0.344467,0.299704,0.549758,0.161355,0.306888,0.950523,0.656736,0.288452,0.314014,0.036234,2291.31,1353.86,1731.73,1241.61,731.131,659.06,1129.36,2030.63,554.647,2290.63,2390.91,697.432,2103.87,1529.71,1627.9,2383.6,651.354,1327.47,1004.5,2416.43,0.176053,1.77557,8.40839,8.82669,9.37272,8.23307,8.42704,7.89192,7.51099,8.84267,7.67324,9.48099,8.56711,7.2556,8.80372,7.79539,8.116,9.64014,7.82119,7.87925,7.64585,8.03043 +0,1,nan,1,1,0,0,20,608,458,905,152,515,194,222,313,287,116,308,525,71,365,128,427,244,530,90,369,0.0298502,0.499326,0.109213,0.658442,0.296609,0.547729,0.92089,0.0116529,0.482388,0.816753,0.440909,0.467446,0.866671,0.634558,0.691072,0.427298,0.935731,0.430601,0.81849,0.169301,1620.8,2336.74,2452.98,931.138,1891.97,1024.51,1768.83,857.058,1446.18,902.426,1395.79,2396.71,537.803,2279.5,757.552,1929.87,2134.99,2188.53,645.062,1121.65,-0.895124,-1.39743,6.45384,6.16362,6.75732,5.02116,6.23576,5.27143,5.29607,5.8421,5.70745,4.76861,5.72995,6.2335,4.18125,5.94984,4.76924,6.07297,5.46348,6.19412,4.43044,5.89085 +0,1,nan,1,1,0,0,20,4492,5135,6323,2663,4374,4929,4794,1739,2943,4992,3294,3445,5537,4759,1673,3343,1737,4109,5474,3127,0.172988,0.952797,0.742638,0.46981,0.44627,0.529484,0.984434,0.759271,0.839205,0.743241,0.836229,0.295936,0.11531,0.164388,0.756778,0.696234,0.90351,0.341663,0.646758,0.0527071,1932.39,1895.44,2410.99,1082.1,1762.29,1970.48,1699.24,694.644,1057.98,1859.18,1202.88,1457.46,2358.5,1988.57,620.923,1261.24,635.716,1730.35,2085.31,1326.71,0.831038,0.184106,8.4294,8.55366,8.75555,7.90419,8.38757,8.51455,8.45021,7.51422,7.94966,8.49576,8.07746,8.16997,8.61805,8.45647,7.40157,8.09907,7.45213,8.35002,8.59278,8.0312 +0,1,nan,1,1,0,0,20,5146,1004,675,2213,717,1731,1856,1198,974,1636,2582,4053,2155,2393,1282,2196,1403,2336,2413,585,0.824977,0.0458041,0.288838,0.319523,0.108416,0.8528,0.699042,0.36712,0.0682025,0.508292,0.344697,0.736322,0.23704,0.637906,0.252615,0.282704,0.372836,0.496234,0.980536,0.413387,2478.61,1546.33,695.998,2235.11,914.829,827.355,1075.12,1115.22,1380.06,1309.28,2498.11,2245.17,2457.75,1523.81,1425.61,2372.45,1236.92,1852.75,961.693,521.523,-0.461903,1.41966,8.52473,6.94676,6.49349,7.70376,6.51075,7.46702,7.51068,7.07609,6.8648,7.43693,7.85074,8.29996,7.68161,7.77268,7.15908,7.71112,7.18778,7.76701,7.79882,6.38172 +0,1,nan,1,1,0,0,20,1479,2147,4737,2900,1161,2234,1261,1150,760,2668,1702,3901,1820,5165,5446,1783,3342,1720,2377,4842,0.95382,0.474546,0.921795,0.191653,0.686125,0.834299,0.495949,0.89178,0.114565,0.424292,0.443455,0.525687,0.220823,0.815289,0.939018,0.870685,0.328747,0.511552,0.672397,0.984354,603.545,1325.39,1979.28,2137.86,576.055,1016.9,768.24,518.077,608.116,1625.78,1040.74,2304.49,1313.13,2421.65,2274.35,796.516,2160.4,953.476,1205.65,1916.78,0.166302,0.74748,7.28208,7.71048,8.44581,7.97712,7.03537,7.71443,7.18112,7.08301,6.6623,7.87719,7.44546,8.30186,7.51153,8.56792,8.59765,7.49737,8.09008,7.40879,7.76368,8.46049 +0,1,nan,1,1,0,0,20,1765,924,1640,1720,948,764,567,1329,647,659,2033,1776,2458,1414,1267,1722,1635,582,2250,1489,0.0467935,0.215945,0.750858,0.390883,0.894789,0.508232,0.91646,0.0562382,0.524791,0.657736,0.0593365,0.198216,0.134067,0.459427,0.322037,0.443921,0.267538,0.966506,0.406767,0.647066,1404.04,929.855,2347.04,1792.46,1459.54,928.862,927.775,1131.76,744.137,814.739,1650.16,1686.16,2144.98,1660.56,1291.29,1912.45,1589.98,948.674,2285.51,1971.19,0.235887,-0.756195,7.44761,6.90762,7.429,7.43165,6.84513,6.68553,6.37565,7.22489,6.45127,6.44138,7.59964,7.5162,7.80539,7.30338,7.15576,7.45634,7.40505,6.36009,7.66264,7.33297 +0,1,nan,1,1,0,0,20,149,64,244,224,266,299,263,346,382,192,185,298,285,561,102,245,340,102,185,108,0.21682,0.584615,0.511398,0.530665,0.185481,0.0786197,0.135563,0.582272,0.183515,0.590027,0.312072,0.0207253,0.619418,0.0391593,0.202558,0.0606131,0.611385,0.914254,0.952532,0.0970036,715.39,617.098,1538.99,1358.28,1343.2,1384.48,1239.5,2185.12,1875.44,1456.06,1037.8,1315.51,2071.14,2433.3,634.371,1182.44,2265.91,871.389,1607.68,533.926,-1.49082,-0.692784,4.9318,4.5292,5.49377,5.35552,5.58349,5.68779,5.53773,5.79522,5.91864,5.38391,5.23784,5.6768,5.71591,6.27906,4.82148,5.54252,5.81135,4.64589,5.23183,4.72223 +0,1,nan,1,1,0,0,20,2146,1487,1923,1315,1035,1988,1652,1906,1230,1584,685,624,726,560,1950,1511,2330,917,679,458,0.848775,0.543537,0.17801,0.808686,0.844417,0.862205,0.180735,0.150892,0.57558,0.0304566,0.233496,0.348817,0.243188,0.136192,0.307159,0.0138192,0.979439,0.931254,0.403637,0.372795,2126.65,1575.36,2272.65,1247.43,973.239,1922.51,2134.17,2478.04,1335.68,2164.01,810.41,733.496,886.887,688.988,2401.31,2040.82,2099.42,851.56,772.801,561.859,-0.31948,0.439978,7.71627,7.2819,7.48754,7.16516,6.93267,7.62126,7.42587,7.56213,7.13096,7.37364,6.48079,6.43181,6.57523,6.27566,7.59943,7.30771,7.76087,6.83732,6.50813,6.17579 +0,1,nan,1,1,0,0,20,2967,1981,4368,5007,6951,1844,3791,2840,7734,6571,4552,2382,3550,2454,5630,4336,3769,1214,3215,5371,0.679074,0.3279,0.167523,0.423438,0.930644,0.658298,0.150177,0.716386,0.94829,0.785376,0.265909,0.577489,0.779453,0.136752,0.78119,0.175054,0.203145,0.24452,0.882693,0.561541,1135.66,912.801,2344.87,2237.39,2158.75,717.04,2139.67,1033.43,2401.93,2309.17,2289,948.833,1217.24,1397.67,1972.65,2308.72,1974.1,616.555,1058,2139.04,0.50844,0.692397,8.0136,7.552,8.38442,8.51469,8.8301,7.53938,8.28083,7.9451,8.94906,8.79688,8.42843,7.76352,8.15247,7.84569,8.63647,8.3741,8.23697,7.10189,8.08375,8.56536 +0,1,nan,1,1,0,0,20,752,965,400,1145,1219,873,422,482,1098,1189,644,451,916,413,367,899,330,547,1192,518,0.190536,0.573375,0.215679,0.172067,0.567356,0.11119,0.335458,0.680709,0.964136,0.185081,0.481325,0.911961,0.89776,0.756624,0.731985,0.811164,0.91009,0.377246,0.678176,0.177795,1515.22,1931.91,808.196,2222.53,2459.92,1739.45,876.931,1034.4,2248.4,2457.47,1318.89,946.467,1859.05,902.316,715.347,1820.29,631.135,1112.8,2494.73,1036.38,-0.721988,0.0163225,6.60443,6.85364,5.97634,6.98722,7.09516,6.74115,6.05992,6.2307,7.01172,7.08792,6.47041,6.14563,6.82049,6.09533,5.86273,6.79801,5.74039,6.29881,7.11102,6.22441 +0,1,nan,1,1,0,0,20,76,158,195,145,135,74,211,328,370,102,37,67,421,106,152,113,58,258,92,137,0.527471,0.0355812,0.664252,0.0604114,0.130286,0.963852,0.20497,0.192953,0.177696,0.814573,0.895395,0.956962,0.234995,0.549003,0.490054,0.529999,0.557164,0.293524,0.857873,0.628627,663.586,756.428,2291.94,639.055,835.314,1406.89,1217.44,1767.28,2340.97,1290.82,572.633,1392.92,2431.59,984.448,1271.05,1149.48,612.924,1850.16,1268,1527.35,-1.469,-1.49686,4.23911,5.10635,5.27387,4.90057,5.06379,4.33739,5.3287,5.71938,6.02334,4.47474,3.54097,4.33772,5.97555,4.60131,4.94506,4.78474,4.11525,5.61467,4.39208,4.92133 +0,1,nan,1,1,0,0,20,29570,11579,14350,3126,9504,34673,11395,3255,5953,4915,14203,5728,21234,7246,16495,11126,25178,14051,3594,29465,0.920584,0.677503,0.228205,0.0984071,0.355647,0.79791,0.160366,0.114083,0.0422038,0.207321,0.248263,0.0739405,0.967086,0.0753921,0.717811,0.521126,0.486138,0.437107,0.0859789,0.617516,1344.59,823.827,2269.14,642.573,1202.04,1973.19,2057.88,636.458,1306.99,806.06,2200.21,1215.15,887.506,1506.6,1093.2,1037.37,2495.18,1555.81,740.965,2299.79,1.42418,1.80747,10.292,9.36271,9.56381,8.06753,9.15878,10.4538,9.34347,8.0863,8.67595,8.49106,9.56921,8.66045,9.96057,8.87806,9.71847,9.31055,10.125,9.56399,8.18754,10.2809 +0,1,nan,1,1,0,0,20,309,931,626,463,470,1131,648,411,1259,471,786,446,1169,522,831,1216,736,792,804,700,0.327688,0.592699,0.342164,0.881933,0.226538,0.676143,0.654325,0.034139,0.53232,0.257018,0.0208108,0.368868,0.495116,0.210359,0.64942,0.174005,0.938466,0.187066,0.728207,0.958149,583.365,1739.01,1170.21,889.75,934.183,2142.36,1213.43,797.835,2474.02,867.127,1579.2,899.007,2297.28,1059.61,1683.29,2358.21,1470.62,1575,1583.35,1180.81,-0.671246,0.0369214,5.70967,6.81171,6.40632,6.15226,6.17679,7.02338,6.45411,6.01192,7.16201,6.10343,6.6942,6.14366,7.08652,6.30217,6.78123,7.10084,6.65684,6.69767,6.72294,6.43809 +0,1,nan,1,1,0,0,20,25419,13997,15526,9672,31929,19558,9373,14884,14667,11058,13747,29186,12971,23167,4315,15917,12195,8051,14630,9871,0.983799,0.407075,0.376921,0.722454,0.868013,0.63028,0.147605,0.873715,0.727887,0.699624,0.472524,0.770674,0.256959,0.818913,0.177803,0.587193,0.765471,0.133271,0.441941,0.49319,1630.69,1937.61,2196.18,863.909,2343.68,1981.05,1763.16,1077.59,1287.28,1007.97,1699.63,2464.9,2139.16,1820.93,803.456,1715.56,1031.76,1586.71,1880.97,1189.46,1.45753,1.32687,10.1597,9.56688,9.65213,9.1776,10.3687,9.88521,9.12825,9.59932,9.58363,9.30153,9.52268,10.29,9.46665,10.0512,8.38238,9.68416,9.41223,9.00378,9.58347,9.19319 +0,1,nan,1,1,0,0,20,819,1688,1074,1284,1824,2432,4334,2685,1728,4238,2626,5326,2394,2428,3583,2278,2929,1223,1370,1394,0.286494,0.719722,0.637121,0.473142,0.388404,0.635843,0.75863,0.345745,0.203765,0.815624,0.693377,0.829008,0.157765,0.637516,0.648527,0.290258,0.788491,0.634903,0.804101,0.636749,658.689,819.222,584.716,829.108,1322.29,1313.91,2111.45,1983.55,1463.86,1961.37,1368.83,2385.23,2128.61,1336.33,1935.03,1819.66,1375.61,640.854,615.688,741.638,-0.0292189,1.00876,6.75003,7.40516,6.98461,7.16842,7.5497,7.79295,8.39119,7.9122,7.46516,8.37494,7.89194,8.5841,7.79315,7.81156,8.19287,7.76999,7.99283,7.07404,7.20466,7.22197 +0,1,nan,1,1,0,0,20,2306,977,2547,2715,1078,1247,819,1953,1390,1902,1293,665,588,2636,581,603,732,1495,1347,420,0.153743,0.773233,0.0529484,0.105242,0.397993,0.607416,0.320483,0.0217379,0.399692,0.512486,0.555134,0.419535,0.493898,0.0124404,0.668962,0.813537,0.201808,0.329773,0.278239,0.790225,1796.38,1761.82,1924.81,2106.85,1207.15,1779.23,863.87,1315.05,1630.11,2333.64,1848.56,736.181,766.004,1834.2,864.052,1077.44,679.507,1589.45,1215.48,837.915,0.373779,-1.20049,7.68274,6.91962,7.8728,7.90039,6.992,7.12852,6.75046,7.52931,7.29035,7.51373,7.2295,6.47161,6.42204,7.87321,6.33233,6.37948,6.65288,7.34903,7.14265,6.15603 +0,1,nan,1,1,0,0,20,365,201,281,396,117,302,336,82,246,255,390,306,377,440,220,359,130,187,176,408,0.843113,0.861921,0.848235,0.580952,0.184685,0.540008,0.691623,0.709864,0.529617,0.4014,0.878034,0.960413,0.302791,0.192189,0.313439,0.0381922,0.108852,0.638911,0.798223,0.471598,2150.51,1301.53,1608.82,2295.56,710.561,1941.83,1848.71,554.934,1561.98,1813.56,2390.44,2146.35,2241.64,2410.31,1321.09,1817.97,776.361,1096.71,1104.06,2212.16,-1.6965,-0.184846,5.82112,5.31547,5.52996,5.93485,4.83542,5.77507,5.6979,4.49113,5.55931,5.73235,5.92043,5.7975,5.96249,6.05548,5.43178,5.80192,4.938,5.18547,5.1627,5.91805 +0,1,nan,1,1,0,0,20,581,1103,1129,283,1230,4192,1839,844,1238,2477,1282,2452,1746,2629,812,2847,902,916,941,739,0.39395,0.0314425,0.756811,0.0111658,0.461562,0.978681,0.341407,0.0971613,0.208556,0.825811,0.173667,0.759563,0.271476,0.648194,0.244068,0.627351,0.618649,0.148487,0.233893,0.108801,717.728,2185.49,780.291,561.816,1377.49,2157.12,2317.09,1419.69,1972.64,1580.02,2113.71,1715.22,2450.42,2125.18,1236.29,2469.85,794.081,1557.28,1438.47,1300.99,-0.756643,1.46926,6.39826,6.97915,7.01498,5.59094,7.14953,8.35783,7.49304,6.64431,7.13691,7.82188,7.15472,7.80665,7.44624,7.85734,6.72183,7.97701,6.8295,6.81222,6.85834,6.5741 +0,1,nan,1,1,0,0,20,1373,1299,1899,2210,1382,1223,2204,1655,1675,1925,843,839,1476,2194,2426,1460,2501,1446,2581,847,0.643077,0.381267,0.946283,0.703033,0.095629,0.307798,0.945593,0.099031,0.14563,0.932409,0.472038,0.861157,0.0255042,0.456745,0.708671,0.95804,0.676515,0.591462,0.987416,0.813072,1154.36,1494.78,1217.4,1797.98,2064.03,1514.27,1428.15,2407.78,2352.82,1324.34,838.053,590.124,2328.28,2319.21,1986.88,966.32,2081.28,1297.31,1599.48,617.481,-0.489505,0.982602,7.19368,7.19486,7.54479,7.69571,7.23688,7.13563,7.70377,7.39426,7.41696,7.61535,6.7054,6.737,7.28844,7.70828,7.80116,7.32536,7.81598,7.25972,7.85817,6.73507 +0,1,nan,1,1,0,0,20,310,87,123,296,28,121,161,227,159,516,256,399,193,196,69,391,118,79,66,110,0.390018,0.70064,0.32598,0.406197,0.791016,0.774308,0.767359,0.543944,0.686644,0.164031,0.359296,0.188146,0.693511,0.72835,0.539191,0.0909529,0.89735,0.803991,0.95464,0.710277,2239.05,1210.79,872.915,2086.95,517.958,1913.32,2309.62,2166.36,2153.96,2384.95,1611.39,2038.3,2387.84,2396.68,726.365,1501.82,2084.77,1242.07,1815.03,1550.1,-1.25184,-1.92127,5.71264,4.50107,4.8937,5.61121,3.4783,4.8171,5.0187,5.3839,5.104,6.20995,5.44271,6.00655,5.19388,5.13064,4.30028,5.88785,4.66652,4.32801,4.4179,4.7296 +0,1,nan,1,1,0,0,20,4024,7503,6741,2845,4875,3461,15294,1400,10925,4005,3092,5590,5670,1376,4621,11754,13340,4829,6414,8023,0.642358,0.387393,0.559959,0.29301,0.445057,0.407473,0.887511,0.162736,0.67213,0.186717,0.483146,0.729516,0.642176,0.149805,0.0550886,0.655179,0.917189,0.0867304,0.458533,0.445872,792.519,2122.29,1489.28,928.016,1282.23,948.487,2156.62,552.692,2076.28,1445.07,764.768,979.45,1117.53,545.479,2042.41,2305.25,1812.93,2030.84,1653.67,2091.94,0.743484,1.36561,8.29591,8.93276,8.81422,7.97667,8.50761,8.1548,9.63177,7.28052,9.29968,8.27438,8.04284,8.62671,8.63932,7.24972,8.4406,9.38115,9.49871,8.47813,8.78041,8.99822 +0,1,nan,1,1,0,0,20,289,595,260,848,529,344,206,761,322,348,638,569,368,221,804,202,199,544,298,627,0.0243731,0.853609,0.242478,0.614938,0.221023,0.251662,0.538739,0.515453,0.342744,0.400408,0.539081,0.175395,0.497547,0.0628076,0.886734,0.621373,0.527338,0.0938144,0.396867,0.159362,874.374,1823.12,854.117,2455.26,1727.48,997.449,607.616,2334.29,964.389,1077.37,1901.59,1856.4,1168.64,630.898,2240.53,573.214,616.172,1821.26,922.987,1941.08,-1.17363,0.121402,5.60284,6.4383,5.60588,6.70701,6.30762,5.76212,5.30132,6.64441,5.73947,5.85726,6.44226,6.37406,5.95036,5.28114,6.64849,5.25306,5.31391,6.34504,5.70216,6.41671 +0,1,nan,1,1,0,0,20,275,320,890,865,1267,484,676,728,349,371,1509,202,1198,484,1080,644,199,303,490,1519,0.153452,0.504129,0.416003,0.416251,0.633733,0.945822,0.29985,0.326243,0.0282175,0.139667,0.789661,0.133303,0.863943,0.50078,0.693518,0.87538,0.0279258,0.324871,0.810333,0.919699,888.302,664.919,2343.93,2264.05,2335.06,619.574,1905.42,2049.68,1365.36,1313.12,2368.45,773.413,1665.49,977.414,1924.97,892.521,836.037,810.874,687.176,1965.04,-1.4587,1.25338,5.52295,5.67283,6.8223,6.78793,7.0914,6.15581,6.46959,6.57565,5.79584,5.89652,7.30104,5.3592,7.04202,6.05388,6.97321,6.43253,5.30498,5.6466,6.08955,7.2773 +0,1,nan,1,1,0,0,20,247,777,818,1110,295,358,913,536,660,1116,512,1199,528,528,395,484,930,448,565,494,0.321914,0.479373,0.145062,0.739594,0.0331412,0.6056,0.492783,0.539423,0.0227995,0.401186,0.10937,0.725921,0.184788,0.715367,0.544223,0.640477,0.263407,0.882402,0.586073,0.423947,548.724,1608.33,1828.02,2172.98,615.184,760.302,1849.36,1092.44,1444.75,2343.3,1064.36,2380.22,1086.08,1083.24,778.04,1022.33,2019.35,812.007,1190.58,1031.14,-0.806579,0.169506,5.55558,6.65763,6.729,7.00264,5.62096,5.92979,6.79955,6.28102,6.47298,7.02074,6.18209,7.09142,6.21508,6.30239,5.94245,6.23183,6.8486,6.0425,6.37496,6.2037 +0,1,nan,1,1,0,0,20,391,547,910,1086,794,1030,1570,963,711,661,744,2127,1248,594,1363,1020,594,871,469,647,0.955185,0.375997,0.294382,0.623574,0.628965,0.0151622,0.156802,0.634918,0.756372,0.605548,0.150361,0.1057,0.219394,0.762905,0.371702,0.210037,0.70016,0.790868,0.586307,0.735489,1129.96,826.025,1219.34,2152.35,1695.8,1028.88,1879.42,1962.07,1714.5,1271.05,883.67,2439.81,1638.45,1501.95,2135.28,1255.73,1314.92,2138.92,848.72,1520.56,-0.00149578,-1.14881,5.93112,6.28318,6.76638,6.95645,6.71186,6.91731,7.35709,6.85086,6.57646,6.45045,6.60985,7.67675,7.14797,6.4366,7.23784,6.89268,6.37569,6.75801,6.06868,6.4804 +0,1,nan,1,1,0,0,20,440,809,634,404,465,576,1447,1383,1556,761,430,1047,1207,644,1436,847,1426,1466,697,995,0.727459,0.966522,0.662259,0.822898,0.708831,0.410689,0.451146,0.925878,0.60364,0.407197,0.942708,0.403372,0.631485,0.854492,0.93927,0.123669,0.923239,0.750658,0.737335,0.771862,702.291,1166.76,964.548,619.941,688.07,965.222,2328.8,2041.58,2401.14,1227.56,623.957,1657.76,1848.88,935.851,2014.09,1442.56,2153.12,2293.33,1053.01,1523.97,-0.547,0.167672,6.12932,6.67705,6.4357,6.0206,6.10574,6.39422,7.28175,7.22972,7.33791,6.63406,6.04715,6.93386,7.08122,6.43773,7.21841,6.74791,7.28247,7.31662,6.53603,6.91149 +# +# Elapsed Time: 0 seconds (Warm-up) +# 0.002 seconds (Sampling) +# 0.002 seconds (Total) +# diff --git a/test/test_model.py b/test/test_model.py index 55c2bfe7..29995ff0 100644 --- a/test/test_model.py +++ b/test/test_model.py @@ -95,10 +95,13 @@ def test_exe_only() -> None: assert not model2._fixed_param -def test_fixed_param() -> None: +def test_legacy_fixed_param() -> None: stan = os.path.join(DATAFILES_PATH, 'datagen_poisson_glm.stan') model = CmdStanModel(stan_file=stan) - assert model._fixed_param + if cmdstan_version_before(2, 36): + assert model._fixed_param + else: + assert not model._fixed_param def test_model_pedantic(caplog: pytest.LogCaptureFixture) -> None: diff --git a/test/test_sample.py b/test/test_sample.py index 369132dc..944ce7a6 100644 --- a/test/test_sample.py +++ b/test/test_sample.py @@ -642,13 +642,18 @@ def test_fixed_param_good() -> None: assert datagen_fit.step_size is None -def test_fixed_param_unspecified() -> None: +def test_sample_no_params() -> None: stan = os.path.join(DATAFILES_PATH, 'datagen_poisson_glm.stan') datagen_model = CmdStanModel(stan_file=stan) datagen_fit = datagen_model.sample(iter_sampling=100, show_progress=False) - assert datagen_fit.step_size is None summary = datagen_fit.summary() - assert 'lp__' not in list(summary.index) + + if cmdstan_version_before(2, 36): + assert 'lp__' not in list(summary.index) + assert datagen_fit.step_size is None + else: + assert 'lp__' in list(summary.index) + assert np.isnan(datagen_fit.step_size).all() exe_only = os.path.join(DATAFILES_PATH, 'exe_only') shutil.copyfile(datagen_model.exe_file, exe_only) @@ -656,9 +661,14 @@ def test_fixed_param_unspecified() -> None: datagen2_model = CmdStanModel(exe_file=exe_only) datagen2_fit = datagen2_model.sample(iter_sampling=200, show_console=True) assert datagen2_fit.chains == 4 - assert datagen2_fit.step_size is None summary = datagen2_fit.summary() - assert 'lp__' not in list(summary.index) + + if cmdstan_version_before(2, 36): + assert datagen2_fit.step_size is None + assert 'lp__' not in list(summary.index) + else: + assert np.isnan(datagen2_fit.step_size).all() + assert 'lp__' in list(summary.index) def test_index_bounds_error() -> None: @@ -823,7 +833,7 @@ def test_validate_good_run() -> None: assert 'Treedepth satisfactory for all transitions.' in diagnostics assert 'No divergent transitions found.' in diagnostics assert 'E-BFMI satisfactory' in diagnostics - assert 'Effective sample size satisfactory.' in diagnostics + assert 'effective sample size satisfactory' in diagnostics.lower() def test_validate_big_run() -> None: @@ -986,6 +996,12 @@ def test_from_csv_fixed_param() -> None: assert fixed_param_sample.draws_pd().shape == (100, 88) +def test_from_csv_no_param_hmc() -> None: + csv_path = os.path.join(DATAFILES_PATH, 'no_param_hmc_sample.csv') + no_parameters_sample = from_csv(path=csv_path) + assert no_parameters_sample.draws_pd().shape == (100, 93) + + def test_custom_metric() -> None: stan = os.path.join(DATAFILES_PATH, 'bernoulli.stan') jdata = os.path.join(DATAFILES_PATH, 'bernoulli.data.json') @@ -1621,33 +1637,18 @@ def test_validate_sample_sig_figs(stanfile='bernoulli.stan'): def test_validate_summary_sig_figs() -> None: - # construct CmdStanMCMC from logistic model output, config - exe = os.path.join(DATAFILES_PATH, 'logistic' + EXTENSION) - rdata = os.path.join(DATAFILES_PATH, 'logistic.data.R') - sampler_args = SamplerArgs(iter_sampling=100) - cmdstan_args = CmdStanArgs( - model_name='logistic', - model_exe=exe, - chain_ids=[1, 2, 3, 4], - seed=12345, - data=rdata, - output_dir=DATAFILES_PATH, - sig_figs=17, - method_args=sampler_args, + # construct CmdStanMCMC from logistic model output + fit = from_csv( + [ + os.path.join(DATAFILES_PATH, 'logistic_output_1.csv'), + os.path.join(DATAFILES_PATH, 'logistic_output_2.csv'), + os.path.join(DATAFILES_PATH, 'logistic_output_3.csv'), + os.path.join(DATAFILES_PATH, 'logistic_output_4.csv'), + ] ) - runset = RunSet(args=cmdstan_args, chains=4) - runset._csv_files = [ - os.path.join(DATAFILES_PATH, 'logistic_output_1.csv'), - os.path.join(DATAFILES_PATH, 'logistic_output_2.csv'), - os.path.join(DATAFILES_PATH, 'logistic_output_3.csv'), - os.path.join(DATAFILES_PATH, 'logistic_output_4.csv'), - ] - retcodes = runset._retcodes - for i in range(len(retcodes)): - runset._set_retcode(i, 0) - fit = CmdStanMCMC(runset) sum_default = fit.summary() + beta1_default = format(sum_default.iloc[1, 0], '.18g') assert beta1_default.startswith('1.3')