simple_poisson_lognormal_model.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "import pyro\n",
    "import pyro.distributions as dist\n",
    "import pyro.poutine as poutine\n",
    "import torch\n",
    "import pandas as pd\n",
    "import plotly.express as px\n",
    "\n",
    "from pyro.infer import MCMC, NUTS, HMC\n",
    "from pyro.infer.mcmc.util import initialize_model, summary"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [],
   "source": [
    "def generate(num_predictors, num_samples):\n",
    "    \"\"\"Generates samples from a simple poisson lognormal model (without noise)\"\"\"\n",
    "    betas = pyro.sample('betas', dist.Normal(torch.zeros(num_predictors+1), 0.5 * torch.ones(num_predictors+1)))\n",
    "    with pyro.plate('gen_data_plate', num_samples):\n",
    "        noise = pyro.sample('noise', dist.Normal(0, 0.025))\n",
    "        pred = pyro.sample('gen_predictors', dist.Uniform(torch.zeros((num_predictors, num_samples)), torch.ones((num_predictors, num_samples)) * 5))\n",
    "        X = torch.cat((torch.ones(1, num_samples), pred), 0)\n",
    "        thetas = betas @ X + noise\n",
    "        accidents = pyro.sample('gen_accidents', dist.Poisson(torch.exp(thetas)))\n",
    "    return betas, X, accidents\n",
    "\n",
    "def prelim_model(num_predictors, num_observations, predictors, data):\n",
    "    \"\"\"Specifies the prior for a simple poisson lognormal model (without noise)\"\"\"\n",
    "    sigma_squared = 1.0/pyro.sample('tau_squared', dist.Chi2(10))\n",
    "    betas = pyro.sample('betas', dist.Normal(torch.zeros(num_predictors+1), 10 * torch.ones(num_predictors+1)))\n",
    "    with pyro.plate('observation_plate', num_observations):\n",
    "        noise = pyro.sample('noise', dist.Normal(0, 1))\n",
    "        thetas = betas @ predictors + noise\n",
    "        accidents = pyro.sample('accidents', dist.Poisson(torch.exp(thetas)), obs=data)\n",
    "    return accidents\n",
    "    "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [],
   "source": [
    "num_samples = 1000\n",
    "num_predictors = 10"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Warmup:   0%|          | 0/2000 [00:00, ?it/s]/home/asifmallik/.local/lib/python3.8/site-packages/pyro/poutine/subsample_messenger.py:58: TracerWarning: torch.tensor results are registered as constants in the trace. You can safely ignore this warning if you use this function to create tensors out of constant variables that would be the same every time you call this function. In any other case, this might cause the trace to be incorrect.\n",
      "  result = torch.tensor(0., device=self.device)\n",
      "Sample: 100%|██████████| 2000/2000 [06:14,  5.34it/s, step size=6.49e-02, acc. prob=0.922]\n"
     ]
    }
   ],
   "source": [
    "betas, X, obs = generate(num_predictors, num_samples)\n",
    "nuts_kernel = NUTS(prelim_model, jit_compile=True)\n",
    "mcmc = MCMC(nuts_kernel, num_samples=1000, warmup_steps=1000, num_chains=1, mp_context = \"spawn\")\n",
    "mcmc.run(num_predictors, num_samples, X, obs)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "tensor([ 0.6999, -0.1162,  0.3481, -0.4646,  0.2658, -0.0985, -0.0074, -0.7256,\n",
      "         0.4661,  0.0450, -0.1655])\n",
      "\n",
      "                   mean       std    median      2.5%     97.5%     n_eff     r_hat\n",
      "  tau_squared      9.99      4.32      9.51      2.30     18.19   1599.79      1.00\n",
      "     betas[0]      0.01      0.27      0.00     -0.56      0.48    306.27      1.00\n",
      "     betas[1]     -0.09      0.04     -0.09     -0.15     -0.02    262.11      1.00\n",
      "     betas[2]      0.40      0.04      0.40      0.32      0.46    328.25      1.00\n",
      "     betas[3]     -0.50      0.04     -0.50     -0.57     -0.43    311.73      1.00\n",
      "     betas[4]      0.33      0.03      0.33      0.26      0.39    411.64      1.00\n",
      "     betas[5]     -0.11      0.04     -0.11     -0.18     -0.04    405.07      1.00\n",
      "     betas[6]      0.06      0.03      0.06     -0.01      0.12    346.84      1.02\n",
      "     betas[7]     -0.80      0.04     -0.80     -0.87     -0.72    305.95      1.00\n",
      "     betas[8]      0.53      0.03      0.53      0.46      0.60    347.71      1.00\n",
      "     betas[9]      0.05      0.03      0.05     -0.02      0.11    307.79      1.00\n",
      "    betas[10]     -0.19      0.03     -0.18     -0.25     -0.11    386.80      1.00\n",
      "     noise[0]     -0.01      1.02      0.02     -1.96      1.99   1574.64      1.00\n",
      "     noise[1]     -0.64      0.83     -0.60     -2.13      0.97   1847.06      1.00\n",
      "     noise[2]     -0.05      0.97     -0.01     -2.03      1.94   2154.68      1.00\n",
      "     noise[3]      0.78      0.91      0.77     -0.98      2.49   2171.07      1.00\n",
      "     noise[4]     -0.13      0.97     -0.11     -2.12      1.69   2259.40      1.00\n",
      "     noise[5]      0.24      0.81      0.27     -1.31      1.92   1667.15      1.00\n",
      "     noise[6]     -0.11      0.92     -0.11     -1.89      1.70   2261.33      1.00\n",
      "     noise[7]      0.17      0.79      0.20     -1.42      1.51   1775.25      1.00\n",
      "     noise[8]      1.10      0.74      1.16     -0.33      2.54   2145.99      1.00\n",
      "     noise[9]     -0.27      0.90     -0.27     -1.95      1.53   1814.07      1.00\n",
      "    noise[10]      0.18      0.74      0.20     -1.33      1.52   1759.47      1.00\n",
      "    noise[11]     -0.23      0.88     -0.25     -1.95      1.54   1403.38      1.00\n",
      "    noise[12]     -0.05      0.98     -0.02     -1.86      1.80   2198.52      1.00\n",
      "    noise[13]     -0.47      0.88     -0.41     -2.26      1.13   1553.62      1.00\n",
      "    noise[14]      0.01      0.60      0.01     -1.06      1.16   1362.55      1.00\n",
      "    noise[15]     -0.40      0.56     -0.36     -1.52      0.60   1497.64      1.00\n",
      "    noise[16]     -0.49      0.55     -0.46     -1.49      0.62   1483.74      1.00\n",
      "    noise[17]      0.61      0.88      0.63     -1.09      2.28   1921.50      1.00\n",
      "    noise[18]      0.28      0.33      0.29     -0.35      0.92    864.85      1.00\n",
      "    noise[19]      0.21      0.75      0.24     -1.25      1.61   1703.80      1.00\n",
      "    noise[20]      0.22      0.54      0.27     -0.86      1.32   1300.90      1.00\n",
      "    noise[21]     -0.27      0.86     -0.23     -2.00      1.47   1457.09      1.00\n",
      "    noise[22]     -0.08      0.36     -0.06     -0.74      0.63    910.32      1.00\n",
      "    noise[23]     -0.34      0.91     -0.34     -2.16      1.42   1503.53      1.00\n",
      "    noise[24]     -0.16      0.68     -0.13     -1.41      1.09   1493.68      1.00\n",
      "    noise[25]      0.28      0.81      0.34     -1.25      1.80   1794.56      1.00\n",
      "    noise[26]     -0.03      0.98     -0.02     -1.75      2.06   1744.39      1.00\n",
      "    noise[27]     -0.21      0.84     -0.19     -1.88      1.43   1678.10      1.00\n",
      "    noise[28]      0.40      0.82      0.45     -1.36      1.85   1579.84      1.00\n",
      "    noise[29]     -0.07      0.97     -0.05     -1.84      1.78   2511.53      1.00\n",
      "    noise[30]     -0.11      0.88     -0.10     -1.82      1.65   2422.85      1.00\n",
      "    noise[31]      0.44      0.80      0.45     -1.18      1.86   3127.97      1.00\n",
      "    noise[32]      0.88      0.51      0.90     -0.17      1.83   1341.01      1.00\n",
      "    noise[33]     -0.09      0.95     -0.07     -1.77      1.82   1922.12      1.00\n",
      "    noise[34]      0.38      0.85      0.45     -1.30      2.03   2053.75      1.00\n",
      "    noise[35]     -0.30      0.90     -0.31     -1.99      1.41   1736.19      1.00\n",
      "    noise[36]      0.91      0.38      0.93      0.13      1.59   1041.22      1.00\n",
      "    noise[37]     -0.27      0.39     -0.28     -1.04      0.49    952.66      1.00\n",
      "    noise[38]      0.12      0.41      0.14     -0.64      0.91    747.45      1.00\n",
      "    noise[39]     -0.06      0.96     -0.02     -1.92      1.80   2699.34      1.00\n",
      "    noise[40]      0.09      0.75      0.11     -1.39      1.63   2681.35      1.00\n",
      "    noise[41]     -0.33      0.87     -0.29     -2.08      1.24   1733.37      1.00\n",
      "    noise[42]      0.30      0.65      0.31     -0.94      1.60   1385.82      1.00\n",
      "    noise[43]     -0.18      0.94     -0.14     -1.91      1.75   1835.11      1.00\n",
      "    noise[44]      0.03      0.49      0.06     -0.95      0.96   1353.51      1.00\n",
      "    noise[45]     -0.03      0.96     -0.01     -1.73      1.88   1851.26      1.00\n",
      "    noise[46]     -0.48      0.82     -0.46     -2.10      0.99   1472.51      1.00\n",
      "    noise[47]     -0.50      0.84     -0.48     -2.23      0.99   1741.65      1.00\n",
      "    noise[48]      0.23      0.78      0.25     -1.23      1.74   1728.34      1.00\n",
      "    noise[49]     -0.25      0.98     -0.23     -2.17      1.54   1552.69      1.00\n",
      "    noise[50]     -0.17      0.92     -0.16     -1.97      1.71   1716.41      1.00\n",
      "    noise[51]     -0.10      1.01     -0.09     -1.92      1.81   2005.63      1.00\n",
      "    noise[52]      0.24      0.60      0.27     -1.06      1.22   1345.87      1.00\n",
      "    noise[53]     -0.20      0.95     -0.16     -2.04      1.65   2190.79      1.00\n",
      "    noise[54]      0.45      0.80      0.51     -1.03      2.02   2310.05      1.00\n",
      "    noise[55]      0.81      0.97      0.77     -0.99      2.64   2038.93      1.00\n",
      "    noise[56]     -0.23      0.89     -0.22     -1.94      1.47   1788.89      1.00\n",
      "    noise[57]     -0.21      0.91     -0.20     -2.08      1.53   2261.35      1.00\n",
      "    noise[58]     -0.19      0.96     -0.18     -1.97      1.77   1866.42      1.00\n",
      "    noise[59]     -0.04      1.01     -0.04     -1.84      2.03   1778.33      1.00\n",
      "    noise[60]     -0.09      0.97     -0.09     -1.97      1.74   1718.14      1.00\n",
      "    noise[61]      1.03      0.48      1.05      0.12      1.99   1094.20      1.00\n",
      "    noise[62]     -0.57      0.84     -0.52     -2.18      0.95   2005.97      1.00\n",
      "    noise[63]     -0.08      0.99     -0.07     -1.96      1.75   2061.63      1.00\n",
      "    noise[64]     -0.91      0.79     -0.85     -2.49      0.50   1762.48      1.00\n",
      "    noise[65]     -0.12      0.93     -0.12     -2.02      1.69   1871.06      1.00\n",
      "    noise[66]     -0.13      0.32     -0.11     -0.77      0.48    745.49      1.00\n",
      "    noise[67]     -0.38      0.92     -0.35     -2.07      1.45   1511.38      1.00\n",
      "    noise[68]     -0.51      0.83     -0.46     -2.07      1.01   2780.12      1.00\n",
      "    noise[69]     -0.36      0.88     -0.35     -2.09      1.21   1838.40      1.00\n",
      "    noise[70]     -0.32      0.72     -0.30     -1.81      0.99   1299.68      1.00\n",
      "    noise[71]     -0.12      0.91     -0.12     -1.76      1.86   1859.63      1.00\n",
      "    noise[72]     -0.22      0.86     -0.22     -1.86      1.47   1669.70      1.00\n",
      "    noise[73]     -0.18      0.91     -0.15     -2.01      1.44   1892.56      1.00\n",
      "    noise[74]      0.29      0.62      0.34     -0.87      1.49   1314.96      1.00\n",
      "    noise[75]      0.29      0.78      0.30     -1.22      1.76   1619.32      1.00\n",
      "    noise[76]     -0.06      1.04     -0.04     -2.07      1.95   1547.50      1.00\n",
      "    noise[77]      0.36      0.64      0.37     -0.88      1.55   1439.09      1.00\n",
      "    noise[78]     -0.08      0.71     -0.06     -1.60      1.13   1719.08      1.00\n",
      "    noise[79]     -0.03      0.98     -0.02     -1.88      1.85   1401.17      1.00\n",
      "    noise[80]     -0.46      0.42     -0.44     -1.29      0.29   1556.17      1.00\n",
      "    noise[81]      0.24      0.77      0.28     -1.21      1.80   1347.18      1.00\n",
      "    noise[82]     -0.02      0.70      0.01     -1.27      1.42   1614.25      1.00\n",
      "    noise[83]      0.28      0.76      0.31     -1.21      1.68   1453.13      1.00\n",
      "    noise[84]     -0.22      0.91     -0.21     -1.99      1.54   2180.45      1.00\n",
      "    noise[85]     -0.00      0.73      0.04     -1.54      1.27   1555.64      1.00\n",
      "    noise[86]     -0.14      1.01     -0.13     -2.02      1.88   2144.11      1.00\n",
      "    noise[87]     -0.29      0.68     -0.25     -1.62      1.05   1677.53      1.00\n",
      "    noise[88]     -0.22      0.99     -0.19     -2.08      1.81   2196.32      1.00\n",
      "    noise[89]      0.11      0.50      0.13     -0.93      0.99   1420.30      1.00\n",
      "    noise[90]      0.61      0.34      0.63     -0.00      1.30    974.27      1.00\n",
      "    noise[91]      0.48      0.83      0.50     -1.06      2.17   1398.10      1.00\n",
      "    noise[92]     -0.08      0.96     -0.07     -1.97      1.67   2170.11      1.00\n",
      "    noise[93]     -0.74      0.77     -0.75     -2.14      0.83   1917.66      1.00\n",
      "    noise[94]      0.40      0.40      0.41     -0.44      1.13   1160.88      1.00\n",
      "    noise[95]     -0.07      0.96     -0.09     -1.80      1.81   1502.26      1.00\n",
      "    noise[96]     -0.03      0.30     -0.03     -0.64      0.52    707.44      1.00\n",
      "    noise[97]      0.12      0.76      0.16     -1.41      1.49   1633.43      1.00\n",
      "    noise[98]     -0.06      0.49     -0.02     -1.07      0.83   1502.19      1.00\n",
      "    noise[99]     -0.21      0.96     -0.20     -2.00      1.77   2261.84      1.00\n",
      "   noise[100]     -0.07      1.03     -0.07     -2.27      1.74   1945.92      1.00\n",
      "   noise[101]     -0.19      0.28     -0.18     -0.69      0.40    699.70      1.00\n",
      "   noise[102]      0.07      0.78      0.10     -1.52      1.45   1975.18      1.00\n",
      "   noise[103]      0.36      0.44      0.39     -0.51      1.22   1072.02      1.00\n",
      "   noise[104]     -0.30      0.87     -0.32     -2.03      1.33   1977.12      1.00\n",
      "   noise[105]      0.12      0.54      0.13     -0.92      1.04   1270.89      1.00\n",
      "   noise[106]     -0.10      0.99     -0.09     -2.08      1.63   1571.33      1.00\n",
      "   noise[107]     -0.17      0.88     -0.16     -1.72      1.59   2042.05      1.00\n",
      "   noise[108]      0.37      0.77      0.43     -1.01      1.99   1435.52      1.00\n",
      "   noise[109]     -0.56      0.83     -0.54     -2.33      0.89   1963.59      1.00\n",
      "   noise[110]      0.43      0.43      0.42     -0.33      1.30   1380.91      1.00\n",
      "   noise[111]     -0.48      0.90     -0.46     -2.21      1.22   1657.77      1.00\n",
      "   noise[112]     -0.49      0.84     -0.43     -2.16      1.04   2103.43      1.00\n",
      "   noise[113]      0.00      0.94      0.01     -1.81      1.78   1875.82      1.00\n",
      "   noise[114]      0.57      0.58      0.61     -0.51      1.74   1433.31      1.00\n",
      "   noise[115]      0.56      0.90      0.59     -1.25      2.16   1623.70      1.00\n",
      "   noise[116]     -0.09      0.55     -0.06     -1.12      0.96   1198.76      1.00\n",
      "   noise[117]     -0.25      0.94     -0.23     -2.03      1.47   2142.35      1.00\n",
      "   noise[118]      0.52      0.82      0.53     -1.09      2.05   1564.96      1.00\n",
      "   noise[119]     -0.06      1.00     -0.02     -1.80      2.18   2094.60      1.00\n",
      "   noise[120]     -0.45      0.84     -0.43     -2.13      1.05   2144.73      1.00\n",
      "   noise[121]     -0.19      0.96     -0.17     -2.10      1.58   1880.32      1.00\n",
      "   noise[122]      0.06      0.51      0.11     -0.97      1.01   1651.45      1.00\n",
      "   noise[123]     -0.03      1.06     -0.04     -1.93      2.12   2199.59      1.00\n",
      "   noise[124]     -0.09      0.28     -0.08     -0.69      0.47    555.94      1.00\n",
      "   noise[125]     -0.12      0.94     -0.13     -1.84      1.74   1508.66      1.00\n",
      "   noise[126]     -0.19      0.92     -0.16     -1.95      1.73   1780.69      1.00\n",
      "   noise[127]     -0.04      1.01     -0.05     -2.03      1.89   1720.24      1.00\n",
      "   noise[128]     -0.10      0.94     -0.08     -2.01      1.62   2256.84      1.00\n",
      "   noise[129]     -0.05      0.32     -0.04     -0.66      0.61    670.33      1.00\n",
      "   noise[130]     -0.31      0.51     -0.27     -1.35      0.58   1509.16      1.00\n",
      "   noise[131]     -0.06      1.01     -0.08     -2.14      1.80   1583.92      1.00\n",
      "   noise[132]     -0.17      0.94     -0.17     -1.95      1.68   1776.04      1.00\n",
      "   noise[133]     -0.25      0.90     -0.23     -2.00      1.59   1743.76      1.00\n",
      "   noise[134]      0.73      0.87      0.78     -1.03      2.33   2090.92      1.00\n",
      "   noise[135]     -0.09      0.86     -0.08     -1.76      1.49   1980.61      1.00\n",
      "   noise[136]     -0.12      0.97     -0.10     -1.82      1.95   2243.15      1.00\n",
      "   noise[137]     -0.38      0.84     -0.35     -2.01      1.19   1667.19      1.00\n",
      "   noise[138]      0.88      0.91      0.86     -0.79      2.70   2449.07      1.00\n",
      "   noise[139]     -0.21      0.90     -0.21     -1.87      1.54   1765.06      1.00\n",
      "   noise[140]     -0.42      0.84     -0.37     -2.15      1.05   1880.79      1.00\n",
      "   noise[141]     -0.40      0.87     -0.35     -2.09      1.30   1555.28      1.00\n",
      "   noise[142]      0.35      0.32      0.37     -0.31      0.92    708.94      1.00\n",
      "   noise[143]      0.46      0.52      0.49     -0.56      1.38   2154.70      1.00\n",
      "   noise[144]     -0.42      0.84     -0.40     -2.06      1.25   1381.21      1.00\n",
      "   noise[145]     -0.35      0.90     -0.32     -2.04      1.37   2279.32      1.00\n",
      "   noise[146]     -0.03      0.99     -0.03     -1.99      1.94   2143.75      1.00\n",
      "   noise[147]     -0.25      0.92     -0.22     -2.17      1.42   3207.92      1.00\n",
      "   noise[148]      0.02      0.47      0.04     -0.93      0.91   1285.76      1.00\n",
      "   noise[149]      0.01      0.56      0.04     -1.12      1.01   1916.73      1.00\n",
      "   noise[150]     -0.44      0.91     -0.40     -2.21      1.29   1699.40      1.00\n",
      "   noise[151]     -0.15      0.88     -0.13     -1.82      1.55   1795.71      1.00\n",
      "   noise[152]     -0.12      0.95     -0.12     -2.03      1.59   1632.18      1.00\n",
      "   noise[153]      0.45      0.82      0.47     -1.06      2.18   1712.40      1.00\n",
      "   noise[154]     -0.27      0.92     -0.27     -2.12      1.49   1974.10      1.00\n",
      "   noise[155]     -0.36      0.92     -0.37     -2.15      1.38   2161.93      1.00\n",
      "   noise[156]      0.15      0.47      0.18     -0.80      1.03   1565.08      1.00\n",
      "   noise[157]      0.08      0.43      0.10     -0.83      0.85   1434.82      1.00\n",
      "   noise[158]     -0.23      0.95     -0.21     -2.05      1.59   1381.63      1.00\n",
      "   noise[159]     -0.30      0.89     -0.28     -1.99      1.48   1526.61      1.00\n",
      "   noise[160]     -0.07      0.99     -0.07     -1.97      1.92   1889.24      1.00\n",
      "   noise[161]      0.70      0.67      0.73     -0.60      2.00   1647.02      1.00\n",
      "   noise[162]      0.05      0.40      0.06     -0.72      0.78   1248.65      1.00\n",
      "   noise[163]     -0.80      0.62     -0.78     -2.00      0.39   1367.10      1.00\n",
      "   noise[164]     -0.08      0.30     -0.08     -0.63      0.57    759.76      1.00\n",
      "   noise[165]     -0.12      0.92     -0.16     -1.94      1.55   1650.09      1.00\n",
      "   noise[166]      0.53      0.57      0.54     -0.74      1.52   1429.01      1.00\n",
      "   noise[167]      0.62      0.85      0.63     -0.92      2.29   2472.11      1.00\n",
      "   noise[168]     -0.07      0.44     -0.05     -0.93      0.77   1147.89      1.00\n",
      "   noise[169]     -0.50      0.46     -0.46     -1.43      0.34   1149.56      1.00\n",
      "   noise[170]     -0.10      0.37     -0.09     -0.81      0.60    879.63      1.00\n",
      "   noise[171]      0.28      0.83      0.30     -1.31      1.88   1753.71      1.00\n",
      "   noise[172]      0.41      0.83      0.44     -1.17      2.02   1669.00      1.00\n",
      "   noise[173]     -0.05      1.02      0.02     -1.95      1.90   1672.05      1.00\n",
      "   noise[174]      0.07      0.75      0.10     -1.34      1.54   1337.82      1.00\n",
      "   noise[175]     -0.16      0.88     -0.14     -1.82      1.57   2001.76      1.00\n",
      "   noise[176]      0.37      0.58      0.40     -0.63      1.57   1497.11      1.00\n",
      "   noise[177]     -0.02      0.94      0.02     -1.88      1.81   1646.44      1.00\n",
      "   noise[178]     -0.15      0.95     -0.09     -1.89      1.73   1797.53      1.00\n",
      "   noise[179]      0.33      0.66      0.36     -1.06      1.51   1377.27      1.00\n",
      "   noise[180]     -0.04      0.94     -0.06     -1.82      1.67   2096.71      1.00\n",
      "   noise[181]      0.20      0.80      0.21     -1.35      1.63   2504.07      1.00\n",
      "   noise[182]      0.24      0.32      0.24     -0.42      0.83   1022.39      1.00\n",
      "   noise[183]      0.01      0.47      0.04     -0.94      0.94   1694.11      1.00\n",
      "   noise[184]     -0.21      0.97     -0.21     -1.98      1.59   1951.80      1.00\n",
      "   noise[185]     -0.07      0.94     -0.10     -1.93      1.68   2351.57      1.00\n",
      "   noise[186]     -0.16      0.94     -0.17     -2.02      1.58   1984.92      1.00\n",
      "   noise[187]     -0.16      0.74     -0.14     -1.58      1.30   1514.47      1.00\n",
      "   noise[188]     -0.03      0.38     -0.01     -0.82      0.68   1115.60      1.00\n",
      "   noise[189]     -0.24      1.00     -0.25     -2.13      1.59   1646.42      1.00\n",
      "   noise[190]      0.55      0.84      0.55     -1.03      2.22   1387.63      1.00\n",
      "   noise[191]     -0.09      0.47     -0.06     -0.98      0.78   1452.07      1.00\n",
      "   noise[192]      0.41      0.84      0.44     -1.10      2.00   1938.32      1.00\n",
      "   noise[193]     -0.93      0.73     -0.89     -2.39      0.52   1648.85      1.00\n",
      "   noise[194]     -0.68      0.82     -0.64     -2.22      0.84   1777.84      1.00\n",
      "   noise[195]     -0.42      0.88     -0.37     -2.12      1.27   1428.67      1.00\n",
      "   noise[196]      0.12      0.63      0.14     -1.22      1.27   1891.34      1.00\n",
      "   noise[197]     -0.95      0.74     -0.94     -2.26      0.62   1786.34      1.00\n",
      "   noise[198]     -0.09      0.95     -0.06     -2.14      1.54   1625.49      1.00\n",
      "   noise[199]     -0.27      0.88     -0.22     -1.91      1.53   2385.16      1.00\n",
      "   noise[200]      0.56      0.90      0.57     -1.15      2.28   2750.00      1.00\n",
      "   noise[201]     -0.32      0.70     -0.25     -1.74      0.97   1219.00      1.00\n",
      "   noise[202]      0.39      0.46      0.40     -0.46      1.26   1263.62      1.00\n",
      "   noise[203]     -0.08      0.96     -0.07     -1.87      1.65   2084.15      1.00\n",
      "   noise[204]      0.45      0.79      0.48     -1.16      1.87   2133.58      1.00\n",
      "   noise[205]      0.01      0.29      0.01     -0.49      0.60    612.50      1.00\n",
      "   noise[206]      0.43      0.44      0.46     -0.52      1.28    650.80      1.00\n",
      "   noise[207]     -0.31      0.88     -0.31     -2.04      1.34   1655.20      1.00\n",
      "   noise[208]     -0.75      0.61     -0.71     -1.87      0.48   1566.63      1.00\n",
      "   noise[209]     -0.10      0.97     -0.10     -1.99      1.82   1520.54      1.00\n",
      "   noise[210]     -0.35      0.86     -0.31     -2.01      1.30   1720.95      1.00\n",
      "   noise[211]     -0.09      0.93     -0.08     -1.95      1.61   1729.74      1.00\n",
      "   noise[212]      0.34      0.37      0.35     -0.39      1.03    894.79      1.00\n",
      "   noise[213]      0.37      0.51      0.39     -0.67      1.28   1521.95      1.00\n",
      "   noise[214]     -0.32      0.48     -0.32     -1.16      0.68   1545.34      1.00\n",
      "   noise[215]      0.13      0.81      0.18     -1.51      1.67   2522.82      1.00\n",
      "   noise[216]     -0.06      0.43     -0.04     -0.87      0.75   1007.47      1.00\n",
      "   noise[217]     -0.06      0.98     -0.07     -1.87      1.92   1991.15      1.00\n",
      "   noise[218]      0.06      0.72      0.11     -1.44      1.32   1204.77      1.00\n",
      "   noise[219]     -0.29      0.86     -0.25     -1.81      1.46   1695.75      1.00\n",
      "   noise[220]     -0.61      0.84     -0.59     -2.23      0.93   1328.55      1.00\n",
      "   noise[221]     -0.02      0.30     -0.03     -0.59      0.59    627.41      1.00\n",
      "   noise[222]      0.16      0.74      0.20     -1.40      1.48   2729.38      1.00\n",
      "   noise[223]      0.58      0.84      0.59     -1.03      2.28   1954.67      1.00\n",
      "   noise[224]      0.37      0.42      0.39     -0.40      1.28   1599.71      1.00\n",
      "   noise[225]      0.37      0.74      0.39     -0.99      1.76   1708.96      1.00\n",
      "   noise[226]     -0.23      0.88     -0.22     -1.94      1.39   1507.99      1.00\n",
      "   noise[227]      0.31      0.84      0.36     -1.24      2.05   1602.87      1.00\n",
      "   noise[228]      0.89      0.57      0.92     -0.26      2.01   1698.66      1.00\n",
      "   noise[229]     -0.06      0.98     -0.09     -2.00      1.85   1812.58      1.00\n",
      "   noise[230]     -0.27      0.91     -0.20     -2.07      1.48   1488.78      1.00\n",
      "   noise[231]     -0.57      0.79     -0.54     -1.96      1.02   2350.05      1.00\n",
      "   noise[232]     -0.15      0.96     -0.14     -1.87      1.87   2108.47      1.00\n",
      "   noise[233]      0.23      0.43      0.23     -0.63      1.01   1163.37      1.00\n",
      "   noise[234]      0.44      0.53      0.46     -0.66      1.39   1674.70      1.00\n",
      "   noise[235]     -0.13      0.97     -0.13     -1.82      1.83   1479.98      1.00\n",
      "   noise[236]     -0.10      0.93     -0.13     -1.98      1.61   1538.63      1.00\n",
      "   noise[237]      0.61      0.82      0.64     -0.91      2.26   2012.34      1.00\n",
      "   noise[238]     -0.16      0.93     -0.13     -2.02      1.63   1938.81      1.00\n",
      "   noise[239]     -0.21      0.98     -0.21     -2.05      1.66   1654.30      1.00\n",
      "   noise[240]      0.88      0.71      0.90     -0.60      2.12   1442.35      1.00\n",
      "   noise[241]     -0.10      0.96     -0.11     -1.86      1.74   1535.17      1.00\n",
      "   noise[242]      0.38      0.46      0.39     -0.53      1.23   1452.15      1.00\n",
      "   noise[243]     -0.32      0.83     -0.31     -1.87      1.38   1521.49      1.00\n",
      "   noise[244]     -0.06      1.00     -0.05     -1.90      1.85   1899.63      1.00\n",
      "   noise[245]      0.84      0.94      0.86     -1.13      2.53   1558.99      1.00\n",
      "   noise[246]      0.01      1.00      0.04     -2.06      1.79   1906.72      1.00\n",
      "   noise[247]     -0.59      0.84     -0.56     -2.06      1.13   2314.68      1.00\n",
      "   noise[248]     -0.40      0.49     -0.38     -1.38      0.53   1533.61      1.00\n",
      "   noise[249]      0.22      0.61      0.26     -1.03      1.30   1267.18      1.00\n",
      "   noise[250]     -0.26      0.91     -0.25     -2.07      1.45   2080.37      1.00\n",
      "   noise[251]      0.16      0.33      0.17     -0.54      0.76    749.43      1.00\n",
      "   noise[252]     -0.04      0.70      0.01     -1.41      1.38   1750.39      1.00\n",
      "   noise[253]     -0.13      0.92     -0.11     -1.88      1.67   1884.50      1.00\n",
      "   noise[254]     -0.25      0.81     -0.25     -1.83      1.27   1674.20      1.00\n",
      "   noise[255]     -0.26      0.87     -0.27     -2.06      1.43   1290.66      1.00\n",
      "   noise[256]     -0.43      0.83     -0.40     -2.06      1.02   1854.30      1.00\n",
      "   noise[257]      0.49      0.39      0.52     -0.33      1.16   1180.03      1.00\n",
      "   noise[258]     -0.31      0.84     -0.27     -1.91      1.33   2053.18      1.00\n",
      "   noise[259]      0.58      0.54      0.63     -0.58      1.53   1149.95      1.00\n",
      "   noise[260]      0.17      0.45      0.18     -0.85      0.95   1925.44      1.00\n",
      "   noise[261]     -0.21      0.87     -0.19     -1.90      1.60   1990.21      1.00\n",
      "   noise[262]      0.85      0.72      0.87     -0.47      2.23   1734.28      1.00\n",
      "   noise[263]      0.47      0.53      0.46     -0.58      1.49   1358.83      1.00\n",
      "   noise[264]     -0.14      0.96     -0.11     -2.09      1.67   1738.89      1.00\n",
      "   noise[265]      0.15      0.38      0.20     -0.56      0.89    893.59      1.00\n",
      "   noise[266]     -0.05      0.93     -0.02     -1.73      1.70   2146.66      1.00\n",
      "   noise[267]     -0.12      0.92     -0.08     -1.83      1.69   1611.08      1.00\n",
      "   noise[268]     -0.09      0.76     -0.04     -1.56      1.34   1297.62      1.00\n",
      "   noise[269]     -0.00      1.01      0.00     -1.98      1.78   2611.84      1.00\n",
      "   noise[270]     -0.52      0.85     -0.47     -2.13      1.12   1884.30      1.00\n",
      "   noise[271]     -0.01      0.45      0.02     -0.87      0.83   1653.58      1.00\n",
      "   noise[272]     -0.52      0.50     -0.48     -1.56      0.39   1578.37      1.00\n",
      "   noise[273]     -0.24      0.93     -0.20     -2.04      1.47   1677.01      1.00\n",
      "   noise[274]     -0.59      0.79     -0.58     -2.25      0.77   2114.55      1.00\n",
      "   noise[275]     -0.47      0.66     -0.41     -1.79      0.75   1386.51      1.00\n",
      "   noise[276]      0.21      0.63      0.24     -1.07      1.43   1648.65      1.00\n",
      "   noise[277]     -0.01      0.95     -0.02     -1.99      1.80   1710.71      1.00\n",
      "   noise[278]      0.40      0.76      0.43     -1.16      1.75   2353.39      1.00\n",
      "   noise[279]     -0.09      0.69     -0.08     -1.34      1.30   2153.15      1.00\n",
      "   noise[280]     -0.21      0.94     -0.18     -2.08      1.51   1513.49      1.00\n",
      "   noise[281]     -0.49      0.67     -0.47     -1.90      0.70   1460.03      1.00\n",
      "   noise[282]      0.68      0.93      0.72     -0.94      2.63   1895.17      1.00\n",
      "   noise[283]      1.16      0.67      1.20     -0.12      2.41   1314.44      1.00\n",
      "   noise[284]      0.15      0.76      0.19     -1.28      1.58   1566.88      1.00\n",
      "   noise[285]     -0.16      0.88     -0.14     -2.13      1.28   1614.37      1.00\n",
      "   noise[286]     -0.69      0.76     -0.69     -2.10      0.83   1342.91      1.00\n",
      "   noise[287]     -0.25      0.35     -0.24     -0.96      0.38   1049.10      1.00\n",
      "   noise[288]     -0.28      0.89     -0.27     -1.99      1.48   1533.21      1.00\n",
      "   noise[289]      0.09      0.23      0.09     -0.36      0.52    427.82      1.00\n",
      "   noise[290]     -0.15      0.94     -0.13     -2.02      1.66   1796.67      1.00\n",
      "   noise[291]      0.31      0.77      0.34     -1.19      1.80   1881.24      1.00\n",
      "   noise[292]     -0.23      0.91     -0.22     -2.08      1.56   2575.91      1.00\n",
      "   noise[293]     -0.19      0.46     -0.15     -1.08      0.69   1393.92      1.00\n",
      "   noise[294]     -0.32      0.91     -0.31     -2.06      1.50   1385.42      1.00\n",
      "   noise[295]      0.23      0.54      0.24     -0.86      1.16   1698.55      1.00\n",
      "   noise[296]      0.19      0.55      0.21     -0.85      1.31   1840.08      1.00\n",
      "   noise[297]     -0.43      0.81     -0.43     -1.90      1.14   1626.51      1.00\n",
      "   noise[298]     -0.01      0.51      0.04     -1.03      0.92   1659.04      1.00\n",
      "   noise[299]      0.48      0.80      0.52     -1.11      2.01   1796.02      1.00\n",
      "   noise[300]     -0.02      0.95     -0.02     -1.83      1.77   2002.01      1.00\n",
      "   noise[301]     -0.32      0.50     -0.27     -1.32      0.58   1833.12      1.00\n",
      "   noise[302]      0.17      0.47      0.20     -0.77      1.01   1152.71      1.00\n",
      "   noise[303]     -0.47      0.81     -0.45     -2.00      1.10   1160.10      1.00\n",
      "   noise[304]      0.66      0.86      0.67     -1.01      2.22   1958.78      1.00\n",
      "   noise[305]     -0.04      1.02     -0.03     -1.78      2.09   2010.87      1.00\n",
      "   noise[306]      0.29      0.82      0.35     -1.30      1.90   1902.07      1.00\n",
      "   noise[307]     -0.17      0.91     -0.13     -2.05      1.47   2237.46      1.00\n",
      "   noise[308]     -0.11      0.92     -0.10     -1.89      1.61   1782.62      1.00\n",
      "   noise[309]     -0.28      0.70     -0.23     -1.77      0.98   1388.77      1.00\n",
      "   noise[310]      0.35      0.83      0.38     -1.23      1.98   1769.48      1.00\n",
      "   noise[311]     -0.45      0.37     -0.43     -1.16      0.30    856.41      1.00\n",
      "   noise[312]      0.73      0.34      0.75      0.08      1.36    733.23      1.00\n",
      "   noise[313]     -0.09      0.97     -0.10     -1.89      1.87   2473.49      1.00\n",
      "   noise[314]     -0.20      0.27     -0.20     -0.75      0.28    597.34      1.00\n",
      "   noise[315]     -0.22      0.94     -0.19     -1.85      1.66   1840.20      1.00\n",
      "   noise[316]      0.16      0.31      0.16     -0.45      0.73    825.41      1.00\n",
      "   noise[317]      0.14      0.64      0.16     -1.07      1.34   1176.54      1.00\n",
      "   noise[318]      0.46      0.54      0.48     -0.63      1.46   1308.05      1.00\n",
      "   noise[319]     -0.12      1.00     -0.11     -2.04      1.68   1864.90      1.00\n",
      "   noise[320]     -0.39      0.90     -0.38     -2.08      1.33   1507.90      1.00\n",
      "   noise[321]      0.88      0.70      0.93     -0.53      2.26   1232.23      1.00\n",
      "   noise[322]     -0.26      0.92     -0.24     -1.98      1.60   1556.68      1.00\n",
      "   noise[323]     -0.27      0.87     -0.27     -2.06      1.34   1576.33      1.00\n",
      "   noise[324]      0.01      0.62      0.06     -1.14      1.15   1439.22      1.00\n",
      "   noise[325]     -0.22      0.98     -0.18     -2.19      1.62   3251.92      1.00\n",
      "   noise[326]      0.05      0.35      0.07     -0.64      0.72   1140.15      1.00\n",
      "   noise[327]     -0.14      0.70     -0.10     -1.57      1.15   1364.67      1.00\n",
      "   noise[328]     -0.05      0.93     -0.06     -1.73      1.86   1508.51      1.00\n",
      "   noise[329]     -0.73      0.63     -0.69     -1.98      0.51   1839.47      1.00\n",
      "   noise[330]     -0.34      0.88     -0.35     -2.10      1.38   1510.30      1.00\n",
      "   noise[331]      0.56      0.83      0.58     -1.13      1.99   2112.88      1.00\n",
      "   noise[332]     -0.54      0.83     -0.48     -2.20      1.00   1964.44      1.00\n",
      "   noise[333]     -0.07      0.47     -0.06     -0.95      0.80   1317.43      1.00\n",
      "   noise[334]     -0.50      0.83     -0.43     -2.23      1.01   1884.26      1.00\n",
      "   noise[335]     -0.36      0.28     -0.36     -0.90      0.18    570.15      1.00\n",
      "   noise[336]     -0.41      0.86     -0.38     -2.11      1.10   1323.40      1.00\n",
      "   noise[337]     -0.03      0.70      0.00     -1.32      1.37   1842.28      1.00\n",
      "   noise[338]      0.34      0.79      0.36     -1.21      1.80   1771.07      1.00\n",
      "   noise[339]      0.18      0.76      0.21     -1.30      1.58   1409.75      1.00\n",
      "   noise[340]      0.61      0.43      0.62     -0.18      1.42   1175.57      1.00\n",
      "   noise[341]     -0.54      0.81     -0.51     -2.15      0.92   1447.33      1.00\n",
      "   noise[342]     -0.35      0.90     -0.33     -2.18      1.33   1960.57      1.00\n",
      "   noise[343]     -0.03      1.03     -0.03     -2.10      2.04   1882.22      1.00\n",
      "   noise[344]      0.01      1.00      0.00     -1.85      1.94   1629.31      1.00\n",
      "   noise[345]     -0.71      0.61     -0.68     -1.92      0.38   1441.45      1.00\n",
      "   noise[346]      0.30      0.64      0.36     -0.99      1.56   1806.04      1.00\n",
      "   noise[347]     -0.33      0.88     -0.30     -2.08      1.38   1451.79      1.00\n",
      "   noise[348]     -0.33      0.57     -0.29     -1.44      0.76   1221.94      1.00\n",
      "   noise[349]      0.41      0.50      0.43     -0.53      1.40   1392.52      1.00\n",
      "   noise[350]     -0.17      0.96     -0.19     -1.92      1.79   1833.78      1.00\n",
      "   noise[351]     -0.43      0.84     -0.41     -2.21      1.05   1483.44      1.00\n",
      "   noise[352]     -0.08      0.89     -0.06     -1.89      1.61   1252.36      1.00\n",
      "   noise[353]      0.50      0.32      0.51     -0.08      1.16    963.84      1.00\n",
      "   noise[354]      0.53      0.53      0.58     -0.47      1.61   1320.90      1.00\n",
      "   noise[355]     -0.08      0.92     -0.09     -1.81      1.67   2031.09      1.00\n",
      "   noise[356]     -0.12      0.99     -0.12     -1.94      1.82   1952.10      1.00\n",
      "   noise[357]     -0.04      1.03     -0.05     -2.09      1.85   1519.98      1.00\n",
      "   noise[358]     -0.22      0.92     -0.20     -1.94      1.62   1743.55      1.00\n",
      "   noise[359]     -0.30      0.85     -0.29     -2.09      1.26   1748.99      1.00\n",
      "   noise[360]      0.09      0.54      0.12     -0.99      1.13   1435.59      1.00\n",
      "   noise[361]     -0.41      0.90     -0.36     -2.17      1.42   2563.21      1.00\n",
      "   noise[362]     -0.16      0.94     -0.12     -1.99      1.60   1447.99      1.00\n",
      "   noise[363]      0.76      0.72      0.81     -0.57      2.22   1932.11      1.00\n",
      "   noise[364]     -0.27      0.97     -0.27     -2.14      1.51   2259.77      1.00\n",
      "   noise[365]     -0.02      1.03     -0.02     -1.90      2.12   2222.01      1.00\n",
      "   noise[366]     -0.23      0.94     -0.22     -1.91      1.64   2256.60      1.00\n",
      "   noise[367]     -1.04      0.69     -0.99     -2.53      0.18   1632.24      1.00\n",
      "   noise[368]     -0.23      0.86     -0.22     -1.76      1.60   1494.76      1.00\n",
      "   noise[369]     -0.14      0.67     -0.12     -1.49      1.14   2235.98      1.00\n",
      "   noise[370]     -0.02      0.98     -0.01     -1.81      2.07   2637.98      1.00\n",
      "   noise[371]     -0.15      0.99     -0.17     -2.05      1.78   2354.69      1.00\n",
      "   noise[372]     -0.61      0.26     -0.60     -1.13     -0.10    487.46      1.00\n",
      "   noise[373]     -0.95      0.73     -0.89     -2.36      0.55   1891.67      1.00\n",
      "   noise[374]     -0.34      0.89     -0.26     -2.20      1.40   1626.72      1.00\n",
      "   noise[375]      0.64      0.57      0.65     -0.44      1.74   1193.27      1.00\n",
      "   noise[376]     -0.05      0.50     -0.02     -0.97      0.94   1238.80      1.00\n",
      "   noise[377]      0.23      0.64      0.29     -1.17      1.33   1474.96      1.00\n",
      "   noise[378]     -0.11      0.33     -0.11     -0.74      0.50    508.50      1.00\n",
      "   noise[379]     -0.07      1.05     -0.09     -1.94      2.01   1577.22      1.00\n",
      "   noise[380]      0.06      0.39      0.07     -0.67      0.83    832.89      1.00\n",
      "   noise[381]      0.34      0.82      0.39     -1.57      1.70   2594.02      1.00\n",
      "   noise[382]     -0.27      0.93     -0.25     -1.96      1.54   1468.40      1.00\n",
      "   noise[383]     -0.66      0.87     -0.65     -2.31      0.98   1922.82      1.00\n",
      "   noise[384]     -0.22      0.89     -0.20     -1.93      1.51   1520.89      1.00\n",
      "   noise[385]      0.63      0.37      0.64     -0.03      1.40    711.71      1.00\n",
      "   noise[386]     -0.23      0.81     -0.19     -2.01      1.08   1856.45      1.00\n",
      "   noise[387]      0.09      0.74      0.11     -1.35      1.44   1489.97      1.00\n",
      "   noise[388]      0.48      0.79      0.51     -1.04      1.96   2151.51      1.00\n",
      "   noise[389]     -0.24      0.86     -0.23     -2.00      1.31   1650.73      1.00\n",
      "   noise[390]     -0.06      0.98     -0.03     -1.87      1.88   2022.64      1.00\n",
      "   noise[391]      0.14      0.73      0.15     -1.44      1.36   2075.79      1.00\n",
      "   noise[392]     -0.28      0.96     -0.26     -2.17      1.56   2095.42      1.00\n",
      "   noise[393]      0.43      0.62      0.47     -0.93      1.49   1560.35      1.00\n",
      "   noise[394]     -0.12      0.96     -0.11     -1.92      1.88   2127.79      1.00\n",
      "   noise[395]     -0.12      0.25     -0.11     -0.66      0.31    457.64      1.00\n",
      "   noise[396]      1.30      0.66      1.36     -0.01      2.49   1730.44      1.00\n",
      "   noise[397]     -0.09      1.03     -0.11     -2.03      1.97   4366.03      1.00\n",
      "   noise[398]     -0.03      0.94     -0.03     -1.95      1.66   2459.00      1.00\n",
      "   noise[399]     -0.14      0.89     -0.16     -1.82      1.58   1923.97      1.00\n",
      "   noise[400]      0.36      0.77      0.38     -1.10      1.90   1755.39      1.00\n",
      "   noise[401]     -0.07      0.95     -0.05     -1.91      1.80   1677.18      1.00\n",
      "   noise[402]     -0.02      1.05     -0.00     -1.97      2.12   1772.31      1.00\n",
      "   noise[403]     -0.18      0.91     -0.17     -2.03      1.47   1459.21      1.00\n",
      "   noise[404]     -0.09      0.68     -0.06     -1.40      1.20   1748.17      1.00\n",
      "   noise[405]     -0.31      0.68     -0.25     -1.68      1.01   1830.88      1.00\n",
      "   noise[406]     -0.37      0.82     -0.33     -1.95      1.08   1579.38      1.00\n",
      "   noise[407]      0.34      0.40      0.37     -0.47      1.09   1405.93      1.00\n",
      "   noise[408]     -0.05      0.96     -0.05     -1.69      1.98   1579.25      1.00\n",
      "   noise[409]      0.59      0.38      0.59     -0.25      1.22    836.48      1.00\n",
      "   noise[410]      0.85      0.97      0.86     -1.01      2.61   2237.40      1.00\n",
      "   noise[411]     -0.26      0.98     -0.23     -2.24      1.56   2054.76      1.00\n",
      "   noise[412]     -0.23      0.88     -0.22     -1.81      1.57   2003.81      1.00\n",
      "   noise[413]      0.37      0.25      0.36     -0.06      0.89    423.23      1.00\n",
      "   noise[414]      0.69      0.41      0.72     -0.13      1.52   1089.97      1.00\n",
      "   noise[415]     -0.40      0.84     -0.35     -2.05      1.16   1758.45      1.00\n",
      "   noise[416]      0.18      0.46      0.22     -0.76      1.04   1156.86      1.00\n",
      "   noise[417]      0.01      1.04      0.04     -1.98      1.88   2186.28      1.00\n",
      "   noise[418]      0.48      0.80      0.52     -1.08      2.06   2083.17      1.00\n",
      "   noise[419]     -0.03      1.01     -0.04     -1.95      1.92   1834.31      1.00\n",
      "   noise[420]      0.13      0.56      0.16     -1.00      1.12   1314.16      1.00\n",
      "   noise[421]     -0.27      0.84     -0.25     -1.87      1.45   1886.49      1.00\n",
      "   noise[422]     -0.02      1.00     -0.03     -1.96      1.78   1619.29      1.00\n",
      "   noise[423]     -0.19      0.85     -0.20     -1.89      1.41   2013.79      1.00\n",
      "   noise[424]     -0.40      0.90     -0.37     -2.25      1.28   2468.14      1.00\n",
      "   noise[425]     -0.07      0.93     -0.07     -1.91      1.63   1742.17      1.00\n",
      "   noise[426]     -0.68      0.79     -0.62     -2.25      0.77   1755.37      1.00\n",
      "   noise[427]      0.85      1.00      0.86     -1.12      2.81   1869.71      1.00\n",
      "   noise[428]      0.45      0.33      0.46     -0.21      1.07    863.83      1.00\n",
      "   noise[429]     -0.17      0.87     -0.14     -2.04      1.37   1911.48      1.00\n",
      "   noise[430]      0.54      0.69      0.60     -0.83      1.80   2029.32      1.00\n",
      "   noise[431]      0.48      0.49      0.50     -0.48      1.39   1963.41      1.00\n",
      "   noise[432]      0.14      0.78      0.17     -1.43      1.55   1311.78      1.00\n",
      "   noise[433]     -0.32      0.58     -0.29     -1.46      0.78   1191.98      1.00\n",
      "   noise[434]      0.49      0.54      0.52     -0.52      1.55   1361.16      1.00\n",
      "   noise[435]     -0.16      0.93     -0.14     -1.83      1.61   2783.71      1.00\n",
      "   noise[436]      0.19      0.81      0.27     -1.37      1.67   1961.79      1.00\n",
      "   noise[437]     -1.01      0.75     -0.97     -2.41      0.46   1499.29      1.00\n",
      "   noise[438]     -0.06      0.96     -0.07     -2.00      1.69   2795.51      1.00\n",
      "   noise[439]     -0.26      0.93     -0.23     -2.01      1.58   1501.03      1.00\n",
      "   noise[440]      0.28      0.61      0.33     -0.96      1.44   1172.16      1.00\n",
      "   noise[441]      0.71      0.90      0.72     -0.98      2.41   1775.08      1.00\n",
      "   noise[442]     -0.09      0.59     -0.05     -1.26      0.98   1341.81      1.00\n",
      "   noise[443]      0.39      0.65      0.42     -0.91      1.60   1674.39      1.00\n",
      "   noise[444]      0.26      0.36      0.27     -0.46      0.92    775.89      1.00\n",
      "   noise[445]     -0.74      0.84     -0.71     -2.53      0.77   1931.90      1.00\n",
      "   noise[446]     -0.11      0.93     -0.12     -1.92      1.72   2053.30      1.00\n",
      "   noise[447]      0.09      0.37      0.11     -0.64      0.77   1186.85      1.01\n",
      "   noise[448]     -0.77      0.78     -0.77     -2.29      0.71   1777.79      1.00\n",
      "   noise[449]     -0.02      0.61      0.03     -1.10      1.28   1415.26      1.00\n",
      "   noise[450]     -0.61      0.82     -0.59     -2.10      1.05   1812.16      1.00\n",
      "   noise[451]     -0.13      0.96     -0.11     -2.08      1.50   2152.40      1.00\n",
      "   noise[452]      0.73      0.53      0.77     -0.29      1.69   1683.30      1.00\n",
      "   noise[453]     -0.05      1.05     -0.03     -2.07      2.08   1619.76      1.00\n",
      "   noise[454]      0.09      0.62      0.10     -1.12      1.23   1853.05      1.00\n",
      "   noise[455]     -0.52      0.77     -0.47     -1.99      0.95   1382.90      1.00\n",
      "   noise[456]      0.21      0.79      0.24     -1.42      1.63   2360.17      1.00\n",
      "   noise[457]      0.09      0.76      0.13     -1.43      1.53   2068.79      1.00\n",
      "   noise[458]     -0.14      0.93     -0.12     -2.02      1.68   2321.84      1.00\n",
      "   noise[459]      0.57      0.90      0.59     -1.24      2.28   1724.24      1.00\n",
      "   noise[460]      0.35      0.65      0.38     -1.10      1.44   2554.64      1.00\n",
      "   noise[461]      0.32      0.59      0.34     -0.80      1.45   1144.19      1.00\n",
      "   noise[462]     -0.08      0.95     -0.06     -1.77      1.96   3201.33      1.00\n",
      "   noise[463]      0.28      0.54      0.31     -0.83      1.27   1529.83      1.00\n",
      "   noise[464]      0.69      0.70      0.73     -0.81      1.84   2060.69      1.00\n",
      "   noise[465]      0.59      0.86      0.61     -1.01      2.27   1608.80      1.00\n",
      "   noise[466]      0.60      0.28      0.60      0.10      1.13    685.43      1.00\n",
      "   noise[467]     -0.19      0.72     -0.15     -1.47      1.40   1517.30      1.00\n",
      "   noise[468]     -0.04      0.97     -0.06     -1.78      1.96   2496.06      1.00\n",
      "   noise[469]     -0.54      0.85     -0.51     -2.23      1.02   2026.46      1.00\n",
      "   noise[470]      1.21      0.66      1.25      0.05      2.58   1168.43      1.00\n",
      "   noise[471]      0.02      0.54      0.06     -1.09      0.99   1362.44      1.00\n",
      "   noise[472]      0.37      0.72      0.36     -1.00      1.77   1516.14      1.00\n",
      "   noise[473]     -0.18      0.89     -0.16     -1.91      1.60   1365.87      1.00\n",
      "   noise[474]     -0.28      0.89     -0.23     -1.98      1.41   2204.54      1.00\n",
      "   noise[475]      0.11      0.75      0.14     -1.44      1.53   1487.38      1.00\n",
      "   noise[476]      0.03      0.54      0.05     -1.04      1.02   1330.24      1.00\n",
      "   noise[477]     -0.11      0.72     -0.05     -1.52      1.30   1043.59      1.00\n",
      "   noise[478]     -0.57      0.84     -0.54     -2.27      0.95   1537.43      1.00\n",
      "   noise[479]     -0.07      0.44     -0.05     -0.89      0.76   1065.63      1.00\n",
      "   noise[480]      0.21      0.22      0.22     -0.19      0.70    386.20      1.00\n",
      "   noise[481]      0.15      0.33      0.15     -0.46      0.83   1019.20      1.00\n",
      "   noise[482]     -0.69      0.48     -0.67     -1.61      0.24   1222.86      1.00\n",
      "   noise[483]     -0.14      0.95     -0.15     -2.09      1.51   1903.05      1.00\n",
      "   noise[484]     -0.18      0.93     -0.15     -1.97      1.47   1824.54      1.00\n",
      "   noise[485]      0.37      0.75      0.40     -1.09      1.76   2073.93      1.00\n",
      "   noise[486]     -0.08      0.99     -0.06     -2.06      1.75   1806.12      1.00\n",
      "   noise[487]      0.43      0.25      0.43     -0.07      0.88    545.87      1.00\n",
      "   noise[488]      0.03      0.52      0.06     -0.96      1.06   1441.09      1.00\n",
      "   noise[489]     -0.15      0.95     -0.15     -1.98      1.77   1883.93      1.00\n",
      "   noise[490]      0.52      0.88      0.53     -1.25      2.18   1809.03      1.00\n",
      "   noise[491]      0.54      0.30      0.54     -0.05      1.09    523.75      1.00\n",
      "   noise[492]     -0.22      0.83     -0.23     -1.81      1.38   1779.58      1.00\n",
      "   noise[493]      1.06      0.63      1.10     -0.27      2.20   1639.40      1.00\n",
      "   noise[494]     -0.12      0.99     -0.10     -2.01      1.65   2081.77      1.00\n",
      "   noise[495]      0.64      0.86      0.66     -0.84      2.33   1379.64      1.00\n",
      "   noise[496]     -0.09      0.23     -0.08     -0.52      0.35    264.18      1.00\n",
      "   noise[497]      0.58      0.67      0.62     -0.59      1.94   1717.01      1.00\n",
      "   noise[498]      0.32      0.80      0.36     -1.37      1.79   2683.68      1.00\n",
      "   noise[499]     -0.07      1.00     -0.05     -2.00      1.80   1958.33      1.00\n",
      "   noise[500]     -0.11      0.61     -0.06     -1.31      0.97   1230.17      1.00\n",
      "   noise[501]      0.07      0.41      0.08     -0.78      0.79   1079.69      1.00\n",
      "   noise[502]     -0.08      1.00     -0.05     -2.08      1.63   1169.33      1.00\n",
      "   noise[503]      0.11      0.77      0.16     -1.42      1.55   1750.28      1.00\n",
      "   noise[504]      0.14      0.53      0.16     -0.89      1.15   1534.32      1.00\n",
      "   noise[505]     -0.20      0.90     -0.18     -1.92      1.54   1987.51      1.00\n",
      "   noise[506]     -0.70      0.79     -0.67     -2.29      0.74   2721.99      1.00\n",
      "   noise[507]     -0.78      0.82     -0.76     -2.33      0.85   1364.28      1.00\n",
      "   noise[508]      0.41      0.47      0.44     -0.52      1.30   1508.40      1.00\n",
      "   noise[509]      0.78      0.70      0.79     -0.59      2.04   1656.87      1.00\n",
      "   noise[510]     -0.43      0.84     -0.41     -2.20      0.98   2469.90      1.00\n",
      "   noise[511]      0.89      0.68      0.90     -0.45      2.17   1707.21      1.00\n",
      "   noise[512]     -0.25      0.98     -0.23     -2.12      1.74   1637.88      1.00\n",
      "   noise[513]      0.05      0.63      0.10     -1.26      1.12   1455.10      1.00\n",
      "   noise[514]     -0.24      0.88     -0.21     -2.02      1.34   1621.69      1.00\n",
      "   noise[515]     -0.41      0.52     -0.36     -1.50      0.51   1283.07      1.00\n",
      "   noise[516]     -0.21      0.93     -0.23     -1.97      1.63   1748.97      1.00\n",
      "   noise[517]     -0.28      0.92     -0.26     -1.90      1.68   1598.22      1.00\n",
      "   noise[518]     -0.46      0.85     -0.43     -2.20      1.07   1740.98      1.00\n",
      "   noise[519]      0.26      0.81      0.32     -1.37      1.75   1309.37      1.00\n",
      "   noise[520]     -0.19      0.88     -0.20     -1.83      1.66   2135.65      1.00\n",
      "   noise[521]     -0.48      0.85     -0.41     -2.03      1.18   1684.99      1.00\n",
      "   noise[522]     -0.26      0.44     -0.24     -1.18      0.56   1108.18      1.00\n",
      "   noise[523]      0.09      0.49      0.08     -0.83      0.99   1130.80      1.00\n",
      "   noise[524]      0.74      0.94      0.75     -1.25      2.47   2636.81      1.00\n",
      "   noise[525]      0.72      0.59      0.74     -0.53      1.79   1982.61      1.00\n",
      "   noise[526]     -0.55      0.83     -0.52     -2.07      1.01   1539.72      1.00\n",
      "   noise[527]     -0.24      0.90     -0.20     -2.04      1.47   1178.77      1.00\n",
      "   noise[528]     -0.41      0.81     -0.37     -2.09      0.97   1962.24      1.00\n",
      "   noise[529]      0.16      0.71      0.20     -1.22      1.53   3064.82      1.00\n",
      "   noise[530]      0.69      0.48      0.71     -0.26      1.61    833.40      1.00\n",
      "   noise[531]      0.17      0.49      0.21     -0.81      1.09   1658.09      1.00\n",
      "   noise[532]     -0.03      0.98     -0.05     -1.74      2.21   2227.80      1.00\n",
      "   noise[533]     -0.10      0.99     -0.04     -2.07      1.80   2255.81      1.00\n",
      "   noise[534]     -0.51      0.55     -0.49     -1.55      0.55   1328.75      1.00\n",
      "   noise[535]     -0.09      0.94     -0.08     -1.93      1.80   1170.06      1.00\n",
      "   noise[536]     -0.02      1.01     -0.02     -2.06      1.78   1667.30      1.00\n",
      "   noise[537]      0.89      0.73      0.92     -0.49      2.33   1964.47      1.00\n",
      "   noise[538]     -0.08      1.02     -0.07     -2.15      1.73   2123.43      1.00\n",
      "   noise[539]      0.57      0.88      0.59     -1.18      2.23   1739.09      1.00\n",
      "   noise[540]      0.53      0.91      0.55     -1.18      2.28   1956.79      1.00\n",
      "   noise[541]     -0.05      1.00     -0.07     -1.94      1.90   2831.22      1.00\n",
      "   noise[542]      0.48      0.67      0.52     -0.82      1.75   1243.21      1.00\n",
      "   noise[543]      0.11      0.54      0.13     -0.94      1.09   1479.15      1.00\n",
      "   noise[544]     -0.18      0.88     -0.17     -1.87      1.42   1602.11      1.00\n",
      "   noise[545]     -0.03      0.98      0.00     -2.10      1.76   1649.15      1.00\n",
      "   noise[546]     -0.17      0.89     -0.14     -1.94      1.45   1720.19      1.00\n",
      "   noise[547]      0.70      0.47      0.71     -0.21      1.58   1094.62      1.00\n",
      "   noise[548]     -0.98      0.74     -0.95     -2.31      0.45   1009.33      1.00\n",
      "   noise[549]      0.05      0.56      0.09     -1.09      1.03   1369.48      1.00\n",
      "   noise[550]     -0.07      0.95     -0.05     -1.87      1.87   1599.29      1.00\n",
      "   noise[551]      0.68      0.67      0.69     -0.77      1.82   1418.78      1.00\n",
      "   noise[552]      0.39      0.70      0.45     -0.91      1.79   1447.30      1.00\n",
      "   noise[553]     -0.11      0.95     -0.08     -1.84      1.80   2070.70      1.00\n",
      "   noise[554]     -0.30      0.93     -0.28     -1.93      1.59   1616.87      1.00\n",
      "   noise[555]      0.52      0.86      0.54     -1.16      2.14   1748.75      1.00\n",
      "   noise[556]     -0.93      0.74     -0.90     -2.44      0.45   1383.45      1.00\n",
      "   noise[557]      0.38      0.77      0.43     -1.16      1.82   1842.33      1.00\n",
      "   noise[558]     -0.51      0.66     -0.47     -1.76      0.78   1934.70      1.00\n",
      "   noise[559]     -0.45      0.85     -0.43     -2.17      1.12   1918.94      1.00\n",
      "   noise[560]      0.72      0.90      0.76     -1.19      2.38   1615.44      1.00\n",
      "   noise[561]     -0.23      0.90     -0.23     -1.91      1.48   1745.74      1.00\n",
      "   noise[562]      0.45      0.82      0.47     -1.21      1.89   1241.39      1.00\n",
      "   noise[563]     -0.18      0.91     -0.17     -2.02      1.51   1682.02      1.00\n",
      "   noise[564]     -0.17      0.88     -0.16     -1.81      1.55   1487.02      1.00\n",
      "   noise[565]      0.79      0.89      0.82     -0.95      2.44   1638.45      1.00\n",
      "   noise[566]     -0.52      0.83     -0.48     -2.15      1.07   1809.62      1.00\n",
      "   noise[567]     -0.62      0.83     -0.61     -2.29      0.89   1987.99      1.00\n",
      "   noise[568]     -0.08      0.98     -0.08     -1.91      1.93   2217.32      1.00\n",
      "   noise[569]      0.09      0.79      0.13     -1.55      1.47   1454.12      1.00\n",
      "   noise[570]      1.18      0.78      1.23     -0.43      2.64   1461.90      1.00\n",
      "   noise[571]     -0.13      0.94     -0.13     -1.84      1.69   1483.82      1.00\n",
      "   noise[572]     -0.31      0.93     -0.29     -2.25      1.42   1760.17      1.00\n",
      "   noise[573]     -0.71      0.77     -0.64     -2.21      0.76   1515.86      1.00\n",
      "   noise[574]      0.64      0.50      0.67     -0.39      1.61   1528.10      1.00\n",
      "   noise[575]      0.38      0.83      0.41     -1.22      1.99   1956.92      1.00\n",
      "   noise[576]     -0.14      0.96     -0.11     -1.97      1.73   1836.18      1.00\n",
      "   noise[577]      0.22      0.36      0.24     -0.42      0.94    880.99      1.00\n",
      "   noise[578]     -0.10      1.02     -0.08     -2.14      1.91   2772.60      1.00\n",
      "   noise[579]      0.62      0.54      0.63     -0.36      1.76    816.92      1.00\n",
      "   noise[580]     -0.00      0.75      0.05     -1.47      1.36   1679.45      1.00\n",
      "   noise[581]     -0.21      0.91     -0.20     -2.03      1.46   1593.76      1.00\n",
      "   noise[582]     -0.21      0.70     -0.16     -1.62      0.99   1643.21      1.00\n",
      "   noise[583]      0.75      0.51      0.77     -0.18      1.80   1566.98      1.00\n",
      "   noise[584]     -0.05      0.99     -0.01     -2.15      1.74   1930.86      1.00\n",
      "   noise[585]      0.74      0.92      0.78     -0.99      2.58   1840.79      1.00\n",
      "   noise[586]     -0.40      0.87     -0.39     -1.89      1.43   2019.78      1.00\n",
      "   noise[587]     -1.14      0.72     -1.11     -2.64      0.18   1641.25      1.00\n",
      "   noise[588]     -0.61      0.77     -0.58     -2.12      0.87   1126.20      1.00\n",
      "   noise[589]     -0.64      0.81     -0.60     -2.21      0.80   1769.77      1.00\n",
      "   noise[590]      0.27      0.62      0.29     -0.89      1.45   2217.53      1.00\n",
      "   noise[591]      0.32      0.57      0.37     -0.88      1.43   1177.19      1.00\n",
      "   noise[592]     -0.04      0.97     -0.06     -2.03      1.67   2085.55      1.00\n",
      "   noise[593]      0.01      0.99     -0.00     -1.90      1.82   1032.82      1.00\n",
      "   noise[594]      0.33      0.40      0.37     -0.49      0.99   1149.12      1.00\n",
      "   noise[595]      0.01      0.98      0.04     -1.76      1.85   1637.90      1.00\n",
      "   noise[596]     -0.18      0.60     -0.15     -1.46      0.85   1239.10      1.00\n",
      "   noise[597]     -0.20      0.89     -0.17     -1.73      1.70   2107.67      1.00\n",
      "   noise[598]     -0.11      1.03     -0.09     -2.31      1.64   2146.18      1.00\n",
      "   noise[599]     -0.12      0.97     -0.13     -2.08      1.61   1733.90      1.00\n",
      "   noise[600]      0.23      0.38      0.25     -0.54      0.90    827.19      1.00\n",
      "   noise[601]     -0.20      0.92     -0.19     -1.94      1.51   1880.78      1.00\n",
      "   noise[602]     -0.07      0.96     -0.06     -2.08      1.65   2116.77      1.00\n",
      "   noise[603]     -0.14      0.70     -0.10     -1.57      1.13   1517.30      1.00\n",
      "   noise[604]     -0.49      0.87     -0.43     -2.32      1.15   1791.76      1.00\n",
      "   noise[605]     -0.02      0.99      0.01     -1.76      2.04   1689.01      1.00\n",
      "   noise[606]      0.24      0.31      0.26     -0.39      0.82   1136.05      1.00\n",
      "   noise[607]     -0.12      0.94     -0.11     -1.82      1.76   1544.41      1.00\n",
      "   noise[608]     -0.11      0.69     -0.05     -1.40      1.19   1894.67      1.00\n",
      "   noise[609]     -0.20      0.44     -0.18     -1.05      0.61   1280.53      1.00\n",
      "   noise[610]     -0.43      0.68     -0.40     -1.69      0.78   1015.01      1.00\n",
      "   noise[611]     -0.03      0.52     -0.00     -0.98      0.94   1671.90      1.00\n",
      "   noise[612]      1.34      0.80      1.35     -0.16      2.98   1438.41      1.00\n",
      "   noise[613]     -0.20      0.89     -0.18     -1.96      1.55   1568.05      1.00\n",
      "   noise[614]      0.13      0.50      0.14     -0.82      1.11   1278.18      1.00\n",
      "   noise[615]      0.20      0.67      0.25     -1.15      1.45   1581.75      1.00\n",
      "   noise[616]      0.66      0.41      0.68     -0.13      1.40   1379.31      1.00\n",
      "   noise[617]      0.96      0.50      1.00     -0.08      1.85   1435.55      1.00\n",
      "   noise[618]      0.43      0.42      0.44     -0.40      1.22   1005.90      1.00\n",
      "   noise[619]     -0.01      0.60      0.00     -1.25      1.07   1421.01      1.00\n",
      "   noise[620]     -0.11      0.95     -0.10     -1.91      1.66   1940.12      1.00\n",
      "   noise[621]     -0.25      0.90     -0.24     -1.91      1.52   2134.07      1.00\n",
      "   noise[622]     -0.01      1.02     -0.02     -2.16      1.95   2175.83      1.00\n",
      "   noise[623]      0.17      0.74      0.22     -1.29      1.58   1627.70      1.00\n",
      "   noise[624]      0.04      0.22      0.04     -0.36      0.44    249.42      1.00\n",
      "   noise[625]     -0.46      0.80     -0.43     -2.02      0.98   1455.37      1.00\n",
      "   noise[626]     -0.34      0.87     -0.31     -2.04      1.22   2675.43      1.00\n",
      "   noise[627]     -0.07      0.74     -0.05     -1.53      1.28   1502.62      1.00\n",
      "   noise[628]      0.22      0.61      0.26     -1.00      1.32   1937.07      1.00\n",
      "   noise[629]     -0.26      0.57     -0.24     -1.32      0.85   1677.46      1.00\n",
      "   noise[630]     -0.25      0.88     -0.27     -1.77      1.44   2223.01      1.00\n",
      "   noise[631]      0.50      0.82      0.52     -1.27      2.01   2423.38      1.00\n",
      "   noise[632]     -0.71      0.78     -0.70     -2.41      0.68   1557.45      1.00\n",
      "   noise[633]     -0.71      0.82     -0.66     -2.19      0.87   1630.26      1.00\n",
      "   noise[634]     -0.16      0.32     -0.14     -0.77      0.47    720.79      1.00\n",
      "   noise[635]     -0.41      0.68     -0.39     -1.71      0.94   1986.88      1.00\n",
      "   noise[636]      0.05      0.51      0.08     -0.93      1.00   1544.09      1.00\n",
      "   noise[637]     -0.08      0.97     -0.07     -1.96      1.73   1801.32      1.00\n",
      "   noise[638]     -0.07      0.89     -0.05     -1.86      1.57   2178.16      1.00\n",
      "   noise[639]     -0.34      0.57     -0.29     -1.50      0.71   1277.71      1.00\n",
      "   noise[640]      0.05      0.75      0.06     -1.51      1.45   1467.24      1.00\n",
      "   noise[641]      0.62      0.88      0.65     -1.00      2.35   1439.44      1.00\n",
      "   noise[642]      0.57      0.53      0.61     -0.40      1.58   1091.21      1.00\n",
      "   noise[643]     -0.05      0.98     -0.02     -1.90      1.69   2059.44      1.00\n",
      "   noise[644]     -0.15      0.97     -0.14     -2.06      1.73   1979.14      1.00\n",
      "   noise[645]     -0.00      0.55      0.03     -1.07      1.03   1186.24      1.00\n",
      "   noise[646]     -0.30      0.90     -0.28     -2.13      1.40   1478.01      1.00\n",
      "   noise[647]     -0.30      0.81     -0.32     -1.87      1.28   1684.38      1.00\n",
      "   noise[648]     -0.13      0.94     -0.14     -1.99      1.57   2220.33      1.00\n",
      "   noise[649]     -0.21      0.93     -0.20     -2.09      1.51   2247.75      1.00\n",
      "   noise[650]      1.04      0.69      1.09     -0.34      2.29   1351.11      1.00\n",
      "   noise[651]     -0.03      0.99     -0.04     -1.84      2.00   1538.43      1.00\n",
      "   noise[652]      0.13      0.78      0.11     -1.47      1.55   1289.03      1.00\n",
      "   noise[653]      0.27      0.55      0.28     -0.85      1.23   1053.68      1.00\n",
      "   noise[654]     -0.59      0.86     -0.54     -2.38      1.01   1352.74      1.00\n",
      "   noise[655]     -0.38      0.86     -0.37     -2.22      1.16   1645.95      1.00\n",
      "   noise[656]     -0.09      0.93     -0.10     -1.75      1.79   1555.02      1.00\n",
      "   noise[657]      0.28      0.41      0.29     -0.59      0.97    932.49      1.00\n",
      "   noise[658]     -0.85      0.76     -0.80     -2.38      0.48   1375.48      1.00\n",
      "   noise[659]      1.18      0.64      1.21     -0.14      2.40   1296.60      1.00\n",
      "   noise[660]     -0.19      0.92     -0.21     -1.84      1.66   1810.22      1.00\n",
      "   noise[661]      1.37      0.38      1.37      0.66      2.16   1027.38      1.00\n",
      "   noise[662]     -0.31      0.91     -0.27     -2.10      1.40   2500.49      1.00\n",
      "   noise[663]      1.01      0.76      1.06     -0.67      2.35   1759.60      1.00\n",
      "   noise[664]     -0.57      0.85     -0.53     -2.12      1.17   1780.37      1.00\n",
      "   noise[665]      0.44      0.66      0.47     -0.91      1.72   1071.98      1.00\n",
      "   noise[666]     -0.44      0.39     -0.43     -1.17      0.30    821.48      1.00\n",
      "   noise[667]     -0.17      0.90     -0.12     -1.87      1.64   2085.24      1.00\n",
      "   noise[668]     -0.11      0.94     -0.07     -2.01      1.55   1502.34      1.00\n",
      "   noise[669]     -0.45      0.85     -0.39     -2.10      1.21   1699.15      1.00\n",
      "   noise[670]     -0.16      0.89     -0.13     -1.77      1.55   1901.56      1.00\n",
      "   noise[671]      0.14      0.48      0.16     -0.81      1.11   1006.28      1.00\n",
      "   noise[672]      0.46      0.54      0.50     -0.67      1.42   1262.62      1.00\n",
      "   noise[673]      0.40      0.82      0.41     -1.37      1.88   1457.64      1.00\n",
      "   noise[674]     -0.29      0.90     -0.27     -2.08      1.39   1591.54      1.00\n",
      "   noise[675]      0.14      0.63      0.19     -1.02      1.39   1366.38      1.00\n",
      "   noise[676]      0.32      0.40      0.33     -0.48      1.07    851.08      1.00\n",
      "   noise[677]     -0.80      0.79     -0.75     -2.40      0.75   1631.61      1.00\n",
      "   noise[678]     -0.43      0.83     -0.42     -2.10      1.13   1504.14      1.00\n",
      "   noise[679]     -0.14      0.95     -0.13     -1.90      1.66   1857.54      1.00\n",
      "   noise[680]     -0.41      0.92     -0.37     -2.26      1.22   1886.65      1.00\n",
      "   noise[681]     -0.83      0.75     -0.79     -2.28      0.58   2180.12      1.00\n",
      "   noise[682]      0.66      0.45      0.68     -0.31      1.45   1124.44      1.00\n",
      "   noise[683]      0.92      0.99      0.91     -0.96      2.79   1862.71      1.00\n",
      "   noise[684]      0.28      0.75      0.30     -1.25      1.56   2130.98      1.00\n",
      "   noise[685]      0.26      0.55      0.27     -0.92      1.29   1401.31      1.00\n",
      "   noise[686]      0.17      0.62      0.18     -1.19      1.23   1714.79      1.00\n",
      "   noise[687]     -0.34      0.81     -0.32     -1.83      1.42   1872.37      1.00\n",
      "   noise[688]     -0.31      0.86     -0.26     -2.00      1.29   1684.53      1.00\n",
      "   noise[689]     -0.24      0.93     -0.20     -2.25      1.41   2046.75      1.00\n",
      "   noise[690]     -0.33      0.70     -0.30     -1.77      0.93   1660.29      1.00\n",
      "   noise[691]     -0.34      0.89     -0.32     -2.20      1.31   2119.88      1.00\n",
      "   noise[692]     -0.51      0.83     -0.46     -2.01      1.11   1773.60      1.00\n",
      "   noise[693]     -0.29      0.70     -0.27     -1.62      1.01   1600.54      1.00\n",
      "   noise[694]     -0.12      0.90     -0.10     -1.83      1.60   1823.86      1.00\n",
      "   noise[695]     -0.02      0.74     -0.02     -1.62      1.23   2307.77      1.00\n",
      "   noise[696]      0.18      0.34      0.21     -0.51      0.80    982.20      1.00\n",
      "   noise[697]      0.25      0.27      0.26     -0.27      0.80    603.31      1.00\n",
      "   noise[698]     -0.34      0.44     -0.31     -1.18      0.50   1004.25      1.00\n",
      "   noise[699]      0.31      0.82      0.35     -1.47      1.80   2761.66      1.00\n",
      "   noise[700]      0.39      0.80      0.43     -1.09      1.92   1872.25      1.00\n",
      "   noise[701]      0.07      0.66      0.11     -1.22      1.26   1647.73      1.00\n",
      "   noise[702]     -0.68      0.81     -0.61     -2.25      0.80   1680.65      1.00\n",
      "   noise[703]     -0.04      1.01     -0.04     -1.99      1.82   1717.55      1.00\n",
      "   noise[704]      0.32      0.79      0.34     -1.19      1.84   1308.87      1.00\n",
      "   noise[705]     -0.07      0.94     -0.09     -1.88      1.76   1842.47      1.00\n",
      "   noise[706]      0.62      0.50      0.64     -0.37      1.57   1212.85      1.00\n",
      "   noise[707]      0.70      0.71      0.73     -0.59      2.11   1697.43      1.00\n",
      "   noise[708]     -0.16      0.90     -0.14     -1.96      1.53   2248.19      1.00\n",
      "   noise[709]     -0.55      0.80     -0.52     -2.07      0.99   2168.08      1.00\n",
      "   noise[710]     -0.06      1.01     -0.06     -2.02      1.97   1617.57      1.00\n",
      "   noise[711]     -0.41      0.88     -0.41     -2.23      1.20   1781.62      1.00\n",
      "   noise[712]     -0.60      0.64     -0.56     -1.85      0.63   1508.26      1.00\n",
      "   noise[713]      0.35      0.50      0.38     -0.63      1.32    957.37      1.00\n",
      "   noise[714]     -0.18      0.89     -0.16     -1.86      1.64   2111.72      1.00\n",
      "   noise[715]     -0.04      0.96     -0.07     -2.08      1.69   2225.38      1.00\n",
      "   noise[716]      0.35      0.32      0.36     -0.25      0.98    803.92      1.00\n",
      "   noise[717]     -0.13      0.98     -0.10     -2.09      1.72   1907.79      1.00\n",
      "   noise[718]     -0.06      0.95     -0.07     -1.89      1.77   1744.04      1.00\n",
      "   noise[719]      0.66      0.85      0.68     -1.09      2.21   2063.00      1.00\n",
      "   noise[720]      0.90      0.75      0.94     -0.55      2.29   1548.57      1.00\n",
      "   noise[721]      0.27      0.53      0.30     -0.73      1.37   1120.09      1.00\n",
      "   noise[722]     -0.22      0.86     -0.18     -1.90      1.41   1966.74      1.00\n",
      "   noise[723]      0.74      0.47      0.77     -0.22      1.57   1453.85      1.00\n",
      "   noise[724]      0.19      0.79      0.24     -1.47      1.60   1536.24      1.00\n",
      "   noise[725]      0.28      0.55      0.31     -0.84      1.33   1230.86      1.00\n",
      "   noise[726]     -0.08      0.92     -0.06     -2.02      1.65   1799.36      1.00\n",
      "   noise[727]     -0.45      0.81     -0.42     -1.91      1.23   1521.80      1.00\n",
      "   noise[728]     -0.03      1.03     -0.03     -1.93      1.97   2272.33      1.00\n",
      "   noise[729]     -0.69      0.49     -0.69     -1.59      0.31   1532.45      1.00\n",
      "   noise[730]     -0.59      0.79     -0.54     -2.27      0.82   1733.45      1.00\n",
      "   noise[731]     -0.02      1.09      0.00     -2.24      1.92   2029.85      1.00\n",
      "   noise[732]      0.21      0.73      0.27     -1.18      1.63   1572.26      1.00\n",
      "   noise[733]     -0.09      0.95     -0.08     -1.81      1.74   1851.43      1.00\n",
      "   noise[734]     -0.02      0.47     -0.00     -0.91      0.84   1440.54      1.00\n",
      "   noise[735]     -0.29      0.68     -0.26     -1.64      1.00   1605.12      1.00\n",
      "   noise[736]      0.18      0.36      0.20     -0.53      0.90    853.76      1.00\n",
      "   noise[737]      0.66      0.67      0.71     -0.65      1.90   1321.34      1.00\n",
      "   noise[738]      0.65      0.53      0.67     -0.34      1.65   1597.49      1.00\n",
      "   noise[739]     -0.09      0.91     -0.11     -2.01      1.55   1818.20      1.00\n",
      "   noise[740]      0.32      0.78      0.37     -1.28      1.78   1464.64      1.00\n",
      "   noise[741]     -0.12      0.95     -0.13     -2.29      1.60   1376.61      1.00\n",
      "   noise[742]     -0.08      0.94     -0.10     -1.91      1.61   1596.08      1.00\n",
      "   noise[743]      0.27      0.78      0.32     -1.35      1.68   1602.70      1.00\n",
      "   noise[744]     -0.03      0.91     -0.02     -1.78      1.80   1736.16      1.00\n",
      "   noise[745]     -0.56      0.44     -0.55     -1.42      0.31   1273.63      1.00\n",
      "   noise[746]      0.56      0.30      0.57     -0.05      1.15    568.88      1.00\n",
      "   noise[747]      0.91      0.68      0.95     -0.47      2.14   1437.60      1.00\n",
      "   noise[748]     -0.33      0.69     -0.29     -1.64      1.07   1266.17      1.00\n",
      "   noise[749]      0.83      0.66      0.84     -0.38      2.01   1460.12      1.00\n",
      "   noise[750]     -0.04      0.40     -0.02     -0.76      0.81   1116.68      1.00\n",
      "   noise[751]     -0.20      0.97     -0.17     -2.07      1.76   1926.40      1.00\n",
      "   noise[752]     -0.37      0.87     -0.36     -1.92      1.41   1680.82      1.00\n",
      "   noise[753]     -0.45      0.80     -0.43     -1.98      1.10   1476.80      1.00\n",
      "   noise[754]      0.17      0.38      0.17     -0.62      0.81    758.92      1.00\n",
      "   noise[755]     -0.14      0.92     -0.13     -1.83      1.68   2433.04      1.00\n",
      "   noise[756]      0.79      0.70      0.83     -0.74      2.01   1971.73      1.00\n",
      "   noise[757]      0.19      0.69      0.23     -1.11      1.46   1844.47      1.00\n",
      "   noise[758]     -0.28      0.92     -0.24     -2.10      1.45   2082.94      1.00\n",
      "   noise[759]      0.07      0.56      0.10     -1.02      1.14   1196.59      1.00\n",
      "   noise[760]      0.26      0.75      0.25     -1.12      1.81   1512.78      1.00\n",
      "   noise[761]      0.22      0.62      0.28     -1.04      1.38   2254.71      1.00\n",
      "   noise[762]      0.52      0.85      0.58     -1.20      2.06   2701.28      1.00\n",
      "   noise[763]     -0.00      1.00      0.01     -2.04      1.87   2063.14      1.00\n",
      "   noise[764]      0.25      0.45      0.27     -0.71      1.02   1075.35      1.01\n",
      "   noise[765]     -0.26      0.60     -0.23     -1.36      0.87   1117.66      1.00\n",
      "   noise[766]      0.62      0.65      0.67     -0.64      1.78   1336.62      1.00\n",
      "   noise[767]     -0.28      0.85     -0.25     -1.85      1.37   2201.77      1.00\n",
      "   noise[768]      0.37      0.77      0.40     -0.96      1.99   1564.12      1.00\n",
      "   noise[769]     -0.09      0.74     -0.08     -1.54      1.34   1493.98      1.00\n",
      "   noise[770]      0.74      0.74      0.78     -0.88      2.10   1708.31      1.00\n",
      "   noise[771]      0.01      1.01      0.01     -1.95      1.83   1704.44      1.00\n",
      "   noise[772]      0.17      0.39      0.17     -0.60      0.90    878.26      1.00\n",
      "   noise[773]      1.05      0.56      1.09     -0.09      2.09   1867.66      1.00\n",
      "   noise[774]      1.02      0.78      1.03     -0.43      2.51   1591.73      1.00\n",
      "   noise[775]     -0.83      0.77     -0.81     -2.27      0.67   1122.58      1.00\n",
      "   noise[776]     -0.40      0.85     -0.37     -1.99      1.40   1895.85      1.00\n",
      "   noise[777]     -0.19      0.67     -0.16     -1.53      1.12   1421.20      1.00\n",
      "   noise[778]     -0.19      0.89     -0.19     -1.83      1.56   2415.71      1.00\n",
      "   noise[779]     -0.08      0.97     -0.07     -2.05      1.59   1569.14      1.00\n",
      "   noise[780]     -0.03      0.88     -0.04     -1.72      1.68   2330.79      1.00\n",
      "   noise[781]      0.69      0.82      0.72     -0.99      2.15   1510.18      1.00\n",
      "   noise[782]      0.25      0.41      0.28     -0.51      1.09    915.84      1.00\n",
      "   noise[783]      0.47      0.59      0.51     -0.72      1.62   1921.83      1.00\n",
      "   noise[784]     -0.18      0.50     -0.17     -1.28      0.67   1233.41      1.00\n",
      "   noise[785]     -0.16      0.93     -0.14     -1.87      1.77   2027.65      1.00\n",
      "   noise[786]     -0.19      0.91     -0.17     -1.88      1.58   1112.60      1.00\n",
      "   noise[787]     -0.32      0.84     -0.31     -1.92      1.20   1723.10      1.00\n",
      "   noise[788]     -0.12      0.94     -0.13     -1.92      1.86   2291.70      1.00\n",
      "   noise[789]      0.08      0.76      0.08     -1.39      1.54   1159.33      1.00\n",
      "   noise[790]     -0.40      0.66     -0.38     -1.63      0.95   1697.94      1.00\n",
      "   noise[791]      0.74      0.85      0.75     -1.00      2.30   1687.00      1.00\n",
      "   noise[792]     -0.19      0.50     -0.18     -1.13      0.79   1113.90      1.00\n",
      "   noise[793]     -0.22      0.92     -0.18     -1.93      1.69   2262.51      1.00\n",
      "   noise[794]     -0.09      0.95     -0.08     -1.96      1.68   1847.09      1.00\n",
      "   noise[795]      1.07      0.56      1.13     -0.03      2.15   1253.01      1.00\n",
      "   noise[796]      0.01      0.70      0.06     -1.41      1.22   1802.82      1.00\n",
      "   noise[797]     -0.74      0.61     -0.71     -1.97      0.38   1618.47      1.00\n",
      "   noise[798]     -0.35      0.93     -0.33     -2.14      1.48   2001.70      1.00\n",
      "   noise[799]      0.63      0.89      0.67     -0.93      2.31   1290.03      1.00\n",
      "   noise[800]     -0.13      0.93     -0.10     -2.00      1.65   2474.13      1.00\n",
      "   noise[801]     -0.02      1.00     -0.05     -1.98      1.79   2307.54      1.00\n",
      "   noise[802]     -0.04      1.00     -0.06     -2.03      1.91   2076.61      1.00\n",
      "   noise[803]     -0.49      0.81     -0.45     -2.08      1.08   1699.66      1.00\n",
      "   noise[804]     -0.34      0.72     -0.29     -1.88      0.82   1594.99      1.00\n",
      "   noise[805]      0.35      0.61      0.37     -0.98      1.40   1222.05      1.00\n",
      "   noise[806]     -0.28      0.88     -0.25     -1.95      1.33   1966.85      1.00\n",
      "   noise[807]     -0.03      0.91     -0.02     -1.69      1.78   1686.26      1.00\n",
      "   noise[808]     -0.18      0.90     -0.15     -1.93      1.54   1662.15      1.00\n",
      "   noise[809]      0.71      0.65      0.73     -0.63      1.94   1410.95      1.00\n",
      "   noise[810]     -0.38      0.92     -0.38     -2.07      1.39   1818.12      1.00\n",
      "   noise[811]     -0.11      0.95     -0.09     -2.14      1.49   1996.25      1.00\n",
      "   noise[812]      0.18      0.37      0.21     -0.54      0.87    953.20      1.00\n",
      "   noise[813]     -0.32      0.90     -0.28     -2.01      1.44   1469.40      1.00\n",
      "   noise[814]     -0.27      0.48     -0.24     -1.24      0.65   1306.92      1.00\n",
      "   noise[815]     -0.21      0.91     -0.21     -1.91      1.52   1909.65      1.00\n",
      "   noise[816]     -0.57      0.56     -0.55     -1.69      0.52   1326.50      1.00\n",
      "   noise[817]      0.77      0.49      0.79     -0.18      1.71   1154.65      1.00\n",
      "   noise[818]     -0.32      0.87     -0.31     -1.98      1.34   1788.11      1.00\n",
      "   noise[819]      0.75      0.87      0.78     -0.84      2.53   2519.81      1.00\n",
      "   noise[820]     -0.11      0.73     -0.11     -1.47      1.26   1554.36      1.00\n",
      "   noise[821]      1.29      0.65      1.31      0.04      2.53   1454.90      1.00\n",
      "   noise[822]      0.29      0.82      0.29     -1.39      1.88   1785.91      1.00\n",
      "   noise[823]     -0.02      0.96     -0.03     -1.89      1.87   2420.56      1.00\n",
      "   noise[824]     -0.36      0.46     -0.35     -1.26      0.50   1201.56      1.00\n",
      "   noise[825]     -0.20      0.52     -0.19     -1.22      0.79   1835.89      1.00\n",
      "   noise[826]      0.42      0.37      0.44     -0.30      1.15   1116.92      1.00\n",
      "   noise[827]     -0.06      0.87     -0.03     -1.70      1.68   1798.51      1.00\n",
      "   noise[828]     -0.05      1.00     -0.03     -1.99      1.86   1344.52      1.00\n",
      "   noise[829]     -0.24      0.72     -0.20     -1.73      1.07   1247.86      1.00\n",
      "   noise[830]      0.37      0.81      0.43     -1.14      1.93   1535.35      1.00\n",
      "   noise[831]      0.28      0.78      0.32     -1.33      1.62   1652.18      1.00\n",
      "   noise[832]     -0.04      0.94     -0.06     -1.93      1.72   1986.48      1.00\n",
      "   noise[833]     -0.04      0.98     -0.01     -1.98      1.77   2301.23      1.00\n",
      "   noise[834]     -0.48      0.36     -0.47     -1.16      0.24   1020.77      1.00\n",
      "   noise[835]      0.75      0.54      0.79     -0.39      1.74   1105.97      1.00\n",
      "   noise[836]     -0.13      0.96     -0.12     -2.14      1.60   1987.66      1.00\n",
      "   noise[837]     -0.22      0.94     -0.22     -2.11      1.49   1377.18      1.00\n",
      "   noise[838]      0.16      0.28      0.16     -0.32      0.76    592.37      1.00\n",
      "   noise[839]      0.53      0.88      0.54     -1.19      2.21   2103.27      1.00\n",
      "   noise[840]     -0.22      0.95     -0.22     -1.86      1.82   1438.03      1.00\n",
      "   noise[841]     -0.01      0.61      0.02     -1.27      1.05   1713.20      1.00\n",
      "   noise[842]     -0.10      0.92     -0.08     -1.84      1.72   2381.77      1.00\n",
      "   noise[843]     -0.64      0.83     -0.62     -2.36      0.81   1417.63      1.00\n",
      "   noise[844]     -0.69      0.80     -0.67     -2.28      0.78   1715.43      1.00\n",
      "   noise[845]     -0.14      0.93     -0.13     -2.04      1.72   1337.82      1.00\n",
      "   noise[846]      0.32      0.80      0.33     -1.24      1.86   2191.01      1.00\n",
      "   noise[847]     -0.20      0.93     -0.20     -2.10      1.53   1729.96      1.00\n",
      "   noise[848]     -0.16      0.92     -0.14     -1.80      1.73   2052.58      1.00\n",
      "   noise[849]     -0.17      0.71     -0.12     -1.61      1.13   1790.00      1.00\n",
      "   noise[850]      0.44      0.83      0.46     -1.22      1.96   1613.59      1.00\n",
      "   noise[851]     -1.14      0.72     -1.09     -2.59      0.26   1464.40      1.00\n",
      "   noise[852]      1.30      0.40      1.33      0.47      2.04   1176.74      1.00\n",
      "   noise[853]     -0.49      0.82     -0.45     -2.16      1.03   2484.94      1.00\n",
      "   noise[854]      0.05      0.76      0.10     -1.44      1.52   1388.29      1.00\n",
      "   noise[855]      0.14      0.43      0.14     -0.72      0.91   1249.86      1.00\n",
      "   noise[856]     -0.07      0.98     -0.08     -2.05      1.73   2515.07      1.00\n",
      "   noise[857]     -0.15      0.96     -0.15     -1.97      1.71   2650.77      1.00\n",
      "   noise[858]     -0.12      0.93     -0.11     -2.01      1.62   2328.72      1.00\n",
      "   noise[859]     -0.12      0.89     -0.13     -1.98      1.48   2574.35      1.00\n",
      "   noise[860]     -0.26      0.90     -0.24     -1.96      1.47   1738.01      1.00\n",
      "   noise[861]     -0.36      0.83     -0.36     -1.95      1.18   2100.02      1.00\n",
      "   noise[862]     -0.66      0.78     -0.66     -2.14      0.75   2046.97      1.00\n",
      "   noise[863]     -0.11      0.87     -0.12     -1.76      1.62   1907.53      1.00\n",
      "   noise[864]     -0.10      0.91     -0.08     -2.00      1.61   1710.52      1.00\n",
      "   noise[865]     -0.59      0.85     -0.54     -2.24      1.04   2196.11      1.00\n",
      "   noise[866]      0.15      0.39      0.17     -0.61      0.90   1097.26      1.00\n",
      "   noise[867]     -0.45      0.85     -0.43     -2.06      1.15   2239.51      1.00\n",
      "   noise[868]     -0.03      1.02     -0.00     -1.97      2.01   2122.15      1.00\n",
      "   noise[869]      0.09      0.79      0.12     -1.58      1.51   1652.59      1.00\n",
      "   noise[870]     -0.37      0.87     -0.35     -2.13      1.34   1658.91      1.00\n",
      "   noise[871]      0.90      0.47      0.92     -0.04      1.87    847.46      1.00\n",
      "   noise[872]     -0.27      0.90     -0.27     -2.10      1.38   2227.80      1.00\n",
      "   noise[873]      0.59      0.84      0.60     -1.01      2.13   1889.75      1.00\n",
      "   noise[874]      0.12      0.62      0.16     -1.05      1.25   1867.46      1.00\n",
      "   noise[875]      0.91      1.02      0.90     -0.91      2.99   2117.81      1.00\n",
      "   noise[876]     -0.39      0.81     -0.37     -1.83      1.24   1823.67      1.00\n",
      "   noise[877]     -0.62      0.79     -0.62     -2.16      0.83   2083.79      1.00\n",
      "   noise[878]     -0.03      1.00     -0.03     -1.90      1.86   3093.65      1.00\n",
      "   noise[879]     -0.05      0.94     -0.04     -1.92      1.71   1946.19      1.00\n",
      "   noise[880]      0.00      0.91      0.01     -1.81      1.67   2129.33      1.00\n",
      "   noise[881]     -0.07      0.99     -0.11     -1.98      1.94   2264.95      1.00\n",
      "   noise[882]     -0.45      0.47     -0.41     -1.45      0.37   1198.07      1.00\n",
      "   noise[883]      0.28      0.63      0.32     -0.96      1.47   1849.99      1.00\n",
      "   noise[884]      0.47      0.55      0.51     -0.62      1.53   1340.00      1.00\n",
      "   noise[885]     -0.01      0.94     -0.01     -1.93      1.75   1685.63      1.00\n",
      "   noise[886]     -0.06      0.99     -0.06     -2.00      1.86   2798.13      1.00\n",
      "   noise[887]     -0.17      0.90     -0.16     -1.85      1.53   2126.43      1.00\n",
      "   noise[888]     -0.23      0.34     -0.24     -0.84      0.44    716.73      1.00\n",
      "   noise[889]     -0.06      1.01     -0.04     -2.19      1.75   1165.45      1.00\n",
      "   noise[890]      0.35      0.80      0.36     -1.21      1.93   1744.87      1.00\n",
      "   noise[891]      0.34      0.86      0.36     -1.31      1.93   2660.46      1.00\n",
      "   noise[892]      0.49      0.80      0.53     -1.27      1.85   2006.07      1.00\n",
      "   noise[893]     -0.45      0.84     -0.43     -2.09      1.13   1659.46      1.00\n",
      "   noise[894]     -0.06      0.46     -0.04     -0.95      0.78   1344.00      1.00\n",
      "   noise[895]      0.58      0.49      0.59     -0.39      1.53   1389.38      1.00\n",
      "   noise[896]     -0.19      0.89     -0.16     -1.83      1.56   1470.21      1.00\n",
      "   noise[897]      0.17      0.39      0.18     -0.54      0.95   1076.40      1.00\n",
      "   noise[898]      0.18      0.51      0.18     -0.84      1.07    878.40      1.00\n",
      "   noise[899]     -0.08      0.60     -0.04     -1.24      1.07   1290.05      1.00\n",
      "   noise[900]     -0.06      1.00     -0.05     -1.95      1.85   2097.67      1.00\n",
      "   noise[901]      0.16      0.35      0.17     -0.51      0.82    849.11      1.00\n",
      "   noise[902]     -0.54      0.82     -0.49     -2.06      1.11   2269.42      1.00\n",
      "   noise[903]     -0.04      0.94     -0.06     -1.84      1.98   1615.32      1.00\n",
      "   noise[904]     -0.15      0.94     -0.12     -1.84      1.68   2204.99      1.00\n",
      "   noise[905]     -0.69      0.51     -0.67     -1.78      0.24   1697.22      1.00\n",
      "   noise[906]     -1.13      0.53     -1.10     -2.25     -0.17   1168.83      1.00\n",
      "   noise[907]     -0.15      0.88     -0.15     -1.87      1.53   1794.33      1.00\n",
      "   noise[908]      0.69      0.70      0.73     -0.61      2.08   1391.42      1.00\n",
      "   noise[909]     -0.48      0.85     -0.43     -1.96      1.31   1888.77      1.00\n",
      "   noise[910]     -0.72      0.65     -0.68     -2.01      0.49   1515.28      1.00\n",
      "   noise[911]     -0.11      0.61     -0.06     -1.23      1.09   1542.83      1.00\n",
      "   noise[912]     -0.07      0.93     -0.09     -1.93      1.68   1873.11      1.00\n",
      "   noise[913]     -0.27      0.88     -0.24     -1.98      1.49   1928.54      1.00\n",
      "   noise[914]     -0.04      0.89     -0.04     -1.76      1.67   1913.09      1.00\n",
      "   noise[915]      0.23      0.65      0.27     -0.99      1.46   1347.29      1.00\n",
      "   noise[916]      0.16      0.80      0.19     -1.47      1.62   1547.02      1.00\n",
      "   noise[917]     -0.23      0.91     -0.20     -2.08      1.59   2270.13      1.00\n",
      "   noise[918]     -0.26      0.85     -0.29     -1.99      1.26   1729.81      1.00\n",
      "   noise[919]     -0.13      0.68     -0.10     -1.39      1.15   2149.42      1.00\n",
      "   noise[920]     -0.28      0.84     -0.27     -1.97      1.32   1761.78      1.00\n",
      "   noise[921]     -0.63      0.81     -0.60     -2.28      0.92   1471.00      1.00\n",
      "   noise[922]     -0.10      0.96     -0.11     -1.92      1.79   1561.97      1.00\n",
      "   noise[923]      0.22      0.28      0.23     -0.32      0.75    834.51      1.00\n",
      "   noise[924]     -0.71      0.78     -0.64     -2.13      0.77   1058.11      1.00\n",
      "   noise[925]     -0.05      1.02     -0.04     -1.94      1.99   2053.98      1.00\n",
      "   noise[926]     -0.15      0.93     -0.14     -2.05      1.54   1776.74      1.00\n",
      "   noise[927]      0.81      0.54      0.85     -0.21      1.82   1079.15      1.00\n",
      "   noise[928]      0.51      0.40      0.53     -0.29      1.27   1428.69      1.00\n",
      "   noise[929]     -0.06      0.98     -0.04     -1.87      1.97   1273.00      1.00\n",
      "   noise[930]     -0.46      0.82     -0.43     -2.14      1.12   1879.36      1.00\n",
      "   noise[931]     -0.13      0.98     -0.13     -1.93      1.83   1910.32      1.00\n",
      "   noise[932]      0.46      0.60      0.48     -0.84      1.47   1071.15      1.00\n",
      "   noise[933]      0.20      0.35      0.22     -0.53      0.84    788.86      1.00\n",
      "   noise[934]      1.07      0.59      1.10     -0.04      2.28   1644.46      1.00\n",
      "   noise[935]      0.84      0.59      0.88     -0.37      1.95   2422.93      1.00\n",
      "   noise[936]     -0.44      0.89     -0.39     -2.22      1.18   1754.14      1.00\n",
      "   noise[937]     -0.14      0.98     -0.13     -2.05      1.66   2660.81      1.00\n",
      "   noise[938]     -0.06      0.77      0.00     -1.60      1.34   1895.48      1.00\n",
      "   noise[939]      1.21      0.56      1.23      0.10      2.25   1225.45      1.00\n",
      "   noise[940]     -0.26      0.94     -0.22     -2.23      1.48   1287.52      1.00\n",
      "   noise[941]     -0.29      0.23     -0.28     -0.75      0.14    280.36      1.00\n",
      "   noise[942]      0.26      0.84      0.29     -1.28      2.01   1884.16      1.00\n",
      "   noise[943]      0.32      0.44      0.33     -0.66      1.07   1540.61      1.00\n",
      "   noise[944]     -0.24      0.92     -0.23     -1.95      1.66   1752.57      1.00\n",
      "   noise[945]      0.47      0.25      0.49      0.02      0.96    656.19      1.00\n",
      "   noise[946]      0.22      0.79      0.28     -1.38      1.77   1808.58      1.00\n",
      "   noise[947]     -0.14      0.98     -0.11     -2.12      1.69   1825.11      1.00\n",
      "   noise[948]     -0.29      0.94     -0.25     -2.09      1.38   1784.43      1.00\n",
      "   noise[949]     -0.38      0.84     -0.34     -2.09      1.14   2147.40      1.00\n",
      "   noise[950]     -0.15      0.95     -0.15     -1.91      1.65   1531.06      1.00\n",
      "   noise[951]      0.24      0.56      0.26     -0.86      1.27   1033.47      1.00\n",
      "   noise[952]     -0.32      0.88     -0.25     -2.05      1.27   1663.34      1.00\n",
      "   noise[953]     -0.02      1.01      0.02     -2.14      1.84   1563.79      1.00\n",
      "   noise[954]      0.01      0.98     -0.00     -1.86      1.90   2288.54      1.00\n",
      "   noise[955]      0.35      0.53      0.38     -0.65      1.33   1225.05      1.00\n",
      "   noise[956]     -0.06      0.90     -0.01     -1.68      1.73   1665.44      1.00\n",
      "   noise[957]     -0.08      0.94     -0.07     -2.18      1.50   1299.44      1.00\n",
      "   noise[958]      0.13      0.75      0.18     -1.39      1.52   2549.69      1.00\n",
      "   noise[959]     -0.13      0.91     -0.11     -1.83      1.75   1573.17      1.00\n",
      "   noise[960]     -0.59      0.80     -0.57     -2.05      1.04   1569.11      1.00\n",
      "   noise[961]      0.70      0.91      0.73     -1.21      2.31   1815.14      1.00\n",
      "   noise[962]      0.25      0.32      0.28     -0.41      0.82    756.90      1.00\n",
      "   noise[963]     -0.00      0.63      0.05     -1.21      1.20   1200.66      1.00\n",
      "   noise[964]     -0.04      0.68      0.00     -1.51      1.11   1666.54      1.00\n",
      "   noise[965]     -0.13      0.90     -0.10     -1.75      1.63   2743.72      1.00\n",
      "   noise[966]     -0.43      0.44     -0.40     -1.29      0.42    696.01      1.00\n",
      "   noise[967]     -0.15      0.73     -0.12     -1.58      1.25   2009.46      1.00\n",
      "   noise[968]     -0.13      0.92     -0.11     -2.15      1.44   2180.97      1.00\n",
      "   noise[969]     -0.39      0.89     -0.30     -2.16      1.31   1920.63      1.00\n",
      "   noise[970]      0.27      0.50      0.31     -0.85      1.11   1357.70      1.00\n",
      "   noise[971]      0.64      0.85      0.67     -1.06      2.28   1932.41      1.00\n",
      "   noise[972]      0.24      0.77      0.28     -1.26      1.80   1292.13      1.00\n",
      "   noise[973]      0.41      0.41      0.42     -0.40      1.19   1027.10      1.00\n",
      "   noise[974]     -0.28      0.92     -0.29     -1.92      1.83   1679.00      1.00\n",
      "   noise[975]     -0.07      0.94     -0.09     -1.92      1.68   1991.14      1.00\n",
      "   noise[976]      0.32      0.83      0.35     -1.32      2.01   1460.60      1.00\n",
      "   noise[977]      0.26      0.67      0.30     -1.11      1.44   1109.30      1.00\n",
      "   noise[978]      0.37      0.66      0.39     -0.88      1.62   1322.08      1.00\n",
      "   noise[979]     -0.08      0.94     -0.10     -1.71      1.85   2309.07      1.00\n",
      "   noise[980]     -0.01      0.61      0.02     -1.19      1.18   1977.40      1.00\n",
      "   noise[981]     -0.21      0.91     -0.18     -1.92      1.53   1479.39      1.00\n",
      "   noise[982]     -0.16      0.91     -0.17     -1.76      1.61   1602.28      1.00\n",
      "   noise[983]     -0.06      1.00     -0.06     -1.89      1.91   2005.29      1.00\n",
      "   noise[984]      0.47      0.78      0.52     -1.09      1.92   1523.01      1.00\n",
      "   noise[985]      0.38      0.44      0.41     -0.38      1.29   1470.35      1.00\n",
      "   noise[986]     -0.56      0.63     -0.51     -1.74      0.65   1186.51      1.00\n",
      "   noise[987]      1.09      0.48      1.13      0.12      1.99    781.68      1.00\n",
      "   noise[988]     -0.51      0.85     -0.52     -2.22      1.01   1813.35      1.00\n",
      "   noise[989]     -0.09      0.97     -0.06     -2.04      1.68   1762.16      1.00\n",
      "   noise[990]      0.12      0.50      0.16     -0.85      1.02   1325.62      1.00\n",
      "   noise[991]      0.71      0.90      0.72     -0.89      2.55   2161.58      1.00\n",
      "   noise[992]     -0.48      0.49     -0.44     -1.51      0.42   1375.95      1.00\n",
      "   noise[993]     -0.04      0.97     -0.07     -2.00      1.78   1907.20      1.00\n",
      "   noise[994]     -0.25      0.91     -0.24     -1.98      1.47   2172.85      1.00\n",
      "   noise[995]     -0.15      1.04     -0.15     -2.09      1.87   2083.47      1.00\n",
      "   noise[996]      0.49      0.54      0.51     -0.53      1.55    905.31      1.00\n",
      "   noise[997]     -0.10      0.69     -0.04     -1.50      1.16   1469.40      1.00\n",
      "   noise[998]     -0.15      0.95     -0.16     -2.02      1.66   2189.48      1.00\n",
      "   noise[999]     -0.23      0.94     -0.24     -2.07      1.60   2846.75      1.00\n",
      "\n",
      "Number of divergences: 0\n",
      "tensor([ 0.,  0.,  0.,  1.,  0.,  1.,  0.,  1.,  2.,  0.,  1.,  0.,  0.,  0.,\n",
      "         2.,  2.,  2.,  1., 11.,  1.,  3.,  0.,  7.,  0.,  1.,  1.,  0.,  0.,\n",
      "         1.,  0.,  0.,  1.,  4.,  0.,  1.,  0.,  8.,  8.,  6.,  0.,  1.,  0.,\n",
      "         2.,  0.,  4.,  0.,  0.,  0.,  1.,  0.,  0.,  0.,  3.,  0.,  1.,  1.,\n",
      "         0.,  0.,  0.,  0.,  0.,  5.,  0.,  0.,  0.,  0., 12.,  0.,  0.,  0.,\n",
      "         1.,  0.,  0.,  0.,  2.,  1.,  0.,  2.,  1.,  0.,  5.,  1.,  1.,  1.,\n",
      "         0.,  1.,  0.,  1.,  0.,  4.,  9.,  1.,  0.,  0.,  6.,  0., 15.,  1.,\n",
      "         4.,  0.,  0., 14.,  1.,  5.,  0.,  3.,  0.,  0.,  1.,  0.,  5.,  0.,\n",
      "         0.,  0.,  3.,  1.,  3.,  0.,  1.,  0.,  0.,  0.,  3.,  0., 17.,  0.,\n",
      "         0.,  0.,  0., 13.,  3.,  0.,  0.,  0.,  1.,  0.,  0.,  0.,  1.,  0.,\n",
      "         0.,  0., 10.,  4.,  0.,  0.,  0.,  0.,  4.,  3.,  0.,  0.,  0.,  1.,\n",
      "         0.,  0.,  4.,  5.,  0.,  0.,  0.,  2.,  7.,  1., 13.,  0.,  3.,  1.,\n",
      "         4.,  4.,  7.,  1.,  1.,  0.,  1.,  0.,  3.,  0.,  0.,  2.,  0.,  1.,\n",
      "        10.,  4.,  0.,  0.,  0.,  1.,  8.,  0.,  1.,  4.,  1.,  0.,  0.,  0.,\n",
      "         2.,  0.,  0.,  0.,  1.,  1.,  4.,  0.,  1., 16.,  6.,  0.,  1.,  0.,\n",
      "         0.,  0.,  8.,  4.,  4.,  1.,  5.,  0.,  1.,  0.,  0., 13.,  1.,  1.,\n",
      "         6.,  1.,  0.,  1.,  3.,  0.,  0.,  0.,  0.,  6.,  3.,  0.,  0.,  1.,\n",
      "         0.,  0.,  2.,  0.,  5.,  0.,  0.,  1.,  0.,  0.,  3.,  2.,  0., 12.,\n",
      "         1.,  0.,  0.,  0.,  0.,  8.,  0.,  3.,  4.,  0.,  2.,  3.,  0.,  7.,\n",
      "         0.,  0.,  1.,  0.,  0.,  4.,  3.,  0.,  0.,  1.,  2.,  0.,  1.,  1.,\n",
      "         0.,  1.,  1.,  3.,  1.,  0.,  0.,  8.,  0., 27.,  0.,  1.,  0.,  4.,\n",
      "         0.,  3.,  3.,  0.,  3.,  1.,  0.,  3.,  5.,  0.,  1.,  0.,  1.,  0.,\n",
      "         0.,  1.,  1.,  7., 11.,  0., 20.,  0., 12.,  2.,  3.,  0.,  0.,  2.,\n",
      "         0.,  0.,  2.,  0.,  8.,  1.,  0.,  1.,  0.,  1.,  0.,  5.,  0., 18.,\n",
      "         0.,  1.,  1.,  1.,  6.,  0.,  0.,  0.,  0.,  1.,  2.,  0.,  2.,  4.,\n",
      "         0.,  0.,  0., 12.,  4.,  0.,  0.,  0.,  0.,  0.,  3.,  0.,  0.,  2.,\n",
      "         0.,  0.,  0.,  0.,  0.,  1.,  0.,  0., 19.,  0.,  0.,  3.,  4.,  2.,\n",
      "        12.,  0.,  7.,  1.,  0.,  0.,  0.,  8.,  0.,  1.,  1.,  0.,  0.,  1.,\n",
      "         0.,  2.,  0., 29.,  3.,  0.,  0.,  0.,  1.,  0.,  0.,  0.,  1.,  1.,\n",
      "         0.,  6.,  0.,  7.,  1.,  0.,  0., 29.,  6.,  0.,  4.,  0.,  1.,  0.,\n",
      "         3.,  0.,  0.,  0.,  0.,  0.,  0.,  1., 11.,  0.,  2.,  4.,  1.,  2.,\n",
      "         3.,  0.,  1.,  0.,  0.,  0.,  2.,  1.,  2.,  2., 12.,  0.,  0.,  8.,\n",
      "         0.,  2.,  0.,  0.,  3.,  0.,  2.,  0.,  1.,  1.,  0.,  1.,  2.,  2.,\n",
      "         0.,  3.,  2.,  1., 17.,  1.,  0.,  0.,  3.,  3.,  1.,  0.,  0.,  1.,\n",
      "         3.,  1.,  0.,  5., 36., 10.,  3.,  0.,  0.,  1.,  0., 23.,  3.,  0.,\n",
      "         1., 18.,  0.,  3.,  0.,  1., 69.,  2.,  1.,  0.,  2.,  6.,  0.,  1.,\n",
      "         3.,  0.,  0.,  0.,  4.,  2.,  0.,  2.,  0.,  2.,  0.,  3.,  0.,  0.,\n",
      "         0.,  1.,  0.,  0.,  4.,  5.,  1.,  3.,  0.,  0.,  0.,  1.,  5.,  4.,\n",
      "         0.,  0.,  2.,  0.,  0.,  2.,  0.,  1.,  1.,  0.,  2.,  3.,  0.,  0.,\n",
      "         0.,  5.,  0.,  3.,  0.,  2.,  2.,  0.,  0.,  1.,  0.,  1.,  1.,  0.,\n",
      "         1.,  0.,  1.,  0.,  0.,  1.,  0.,  0.,  0.,  1.,  2.,  0.,  0.,  0.,\n",
      "         4.,  1.,  0.,  8.,  0.,  4.,  1.,  0.,  1.,  4.,  0.,  1.,  0.,  0.,\n",
      "         0.,  0.,  2.,  3.,  0.,  0.,  7.,  0.,  2.,  0.,  0.,  0.,  9.,  0.,\n",
      "         0.,  1.,  0.,  0., 12.,  0.,  1.,  4.,  1.,  3.,  2.,  0.,  4.,  2.,\n",
      "         6.,  5.,  6.,  2.,  0.,  0.,  0.,  1., 91.,  0.,  0.,  1.,  2.,  2.,\n",
      "         0.,  1.,  0.,  0.,  9.,  1.,  3.,  0.,  0.,  2.,  1.,  1.,  4.,  0.,\n",
      "         0.,  3.,  0.,  0.,  0.,  0.,  2.,  0.,  1.,  3.,  0.,  0.,  0.,  7.,\n",
      "         0.,  3.,  0.,  9.,  0.,  2.,  0.,  2.,  8.,  0.,  0.,  0.,  0.,  4.,\n",
      "         3.,  1.,  0.,  2.,  7.,  0.,  0.,  0.,  0.,  0.,  5.,  1.,  1.,  3.,\n",
      "         2.,  0.,  0.,  0.,  1.,  0.,  0.,  1.,  0.,  1., 11., 20.,  5.,  1.,\n",
      "         1.,  2.,  0.,  0.,  1.,  0.,  4.,  2.,  0.,  0.,  0.,  0.,  1.,  4.,\n",
      "         0.,  0., 13.,  0.,  0.,  1.,  2.,  4.,  0.,  5.,  1.,  3.,  0.,  0.,\n",
      "         0.,  3.,  0.,  0.,  1.,  0.,  4.,  1.,  9.,  2.,  3.,  0.,  1.,  0.,\n",
      "         0.,  1.,  0.,  4., 14.,  2.,  1.,  2.,  7.,  0.,  0.,  0.,  9.,  0.,\n",
      "         2.,  1.,  0.,  3.,  1.,  2.,  1.,  0.,  5.,  2.,  2.,  0.,  1.,  1.,\n",
      "         2.,  0.,  6.,  4.,  2.,  0.,  0.,  1.,  0.,  0.,  0.,  1.,  7.,  3.,\n",
      "         3.,  0.,  0.,  0.,  0.,  1.,  1.,  1.,  3.,  0.,  0.,  4.,  1.,  1.,\n",
      "         0.,  1.,  0.,  0.,  0.,  0.,  1.,  2.,  0.,  0.,  0.,  2.,  0.,  0.,\n",
      "         9.,  0.,  3.,  0.,  2.,  4.,  0.,  1.,  1.,  3.,  1.,  0.,  3.,  3.,\n",
      "         7.,  0.,  0.,  1.,  1.,  1.,  0.,  0.,  8.,  4.,  0.,  0., 20.,  1.,\n",
      "         0.,  2.,  0.,  0.,  0.,  0.,  1.,  0.,  0.,  1.,  1.,  0.,  7.,  0.,\n",
      "         1.,  6.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  7.,  0.,\n",
      "         0.,  1.,  0.,  6.,  0.,  1.,  2.,  1.,  0.,  0.,  0.,  0.,  0.,  0.,\n",
      "         4.,  2.,  4.,  0.,  0.,  0., 11.,  0.,  1.,  1.,  1.,  0.,  4.,  4.,\n",
      "         0.,  7.,  3.,  2.,  0.,  8.,  0.,  0.,  0.,  2.,  2.,  0.,  2.,  0.,\n",
      "         1.,  2.,  0.,  0.,  0.,  2.,  1.,  0.,  0.,  1.,  0.,  0.,  0., 15.,\n",
      "         0.,  0.,  0.,  4.,  6.,  0.,  0.,  0.,  3., 11.,  3.,  3.,  0.,  0.,\n",
      "         1.,  4.,  0., 50.,  1.,  5.,  0., 25.,  1.,  0.,  0.,  0.,  0.,  3.,\n",
      "         0.,  0.,  0.,  3.,  0.,  0.,  1.,  0.,  0.,  1., 12.,  2.,  1.,  0.,\n",
      "         5.,  1.,  0.,  0.,  4.,  1.,  1.,  5.,  0.,  0.,  1.,  2.,  2.,  0.,\n",
      "         2.,  0.,  0.,  0.,  1.,  5.,  1.,  5.,  0.,  0.,  3.,  1.,  3.,  0.,\n",
      "         0.,  0.,  3.,  1.,  0.,  0.])\n"
     ]
    }
   ],
   "source": [
    "print(betas)\n",
    "mcmc.summary(prob=0.95)\n",
    "print(obs)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.8.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}