Skip to content

Commit

Permalink
updated
Browse files Browse the repository at this point in the history
  • Loading branch information
Jakob Runge committed Sep 26, 2017
1 parent a4b11b9 commit 201a4e6
Show file tree
Hide file tree
Showing 2 changed files with 40 additions and 4 deletions.
4 changes: 2 additions & 2 deletions run_pcmci_parallel.py
Expand Up @@ -217,8 +217,8 @@ def run_mci_parallel(j, pcmci_of_j, all_parents):
print("\n\n## Resulting condition sets:")
for j in [var for var in all_parents.keys()]:
pcmci_objects[j]._print_parents_single(j, all_parents[j],
pcmci_objects[j].test_statistic_values[j],
pcmci_objects[j].p_max[j])
pcmci_objects[j].val_matrix[j],
pcmci_objects[j].p_matrix[j])

if verbosity > -1:
print("\n##\n## Running Parallelized Tigramite MCI algorithm\n##"
Expand Down
40 changes: 38 additions & 2 deletions tigramite_tutorial.ipynb
Expand Up @@ -3945,7 +3945,7 @@
},
{
"cell_type": "code",
"execution_count": null,
"execution_count": 73,
"metadata": {},
"outputs": [
{
Expand Down Expand Up @@ -3982,8 +3982,44 @@
"## Variable 2\n",
"Null distribution for GPDC not available for deg. of freed. = 380.\n",
"Generating null distribution for df = 380. \n",
"For faster computations, run function generate_and_save_nulldists(...) to precompute null distribution and load *.npz file with argument null_dist_filename\n"
"For faster computations, run function generate_and_save_nulldists(...) to precompute null distribution and load *.npz file with argument null_dist_filename\n",
"\n",
"## Resulting condition sets:\n",
"\n",
" Variable 0 has 0 parent(s):\n",
"\n",
" Variable 1 has 0 parent(s):\n",
"\n",
" Variable 2 has 4 parent(s):\n",
" (1 -2): max_pval = 0.00000, min_val = 0.266\n",
" (0 -2): max_pval = 0.00000, min_val = 0.112\n",
" (2 -2): max_pval = 0.00000, min_val = 0.050\n",
" (0 -4): max_pval = 0.13600, min_val = 0.010\n",
"\n",
"##\n",
"## Predicting target 2\n",
"##\n"
]
},
{
"data": {
"text/plain": [
"<matplotlib.text.Text at 0x7f1c6af91a90>"
]
},
"execution_count": 73,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYQAAAEWCAYAAABmE+CbAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl8VPW9//HXhxAgIhJU3EAQl+JVUbFYWrl6tdqi1brV\nWq2/WkVxtyIWBXdbUCwu9br1oqj1itYdFVFcsHrFFUQ2kVZtLQYXrAQFI4bk8/tjZsLMZJYzycyc\nmcn7+XjkQc4y53wG8XzOdzd3R0REpFPYAYiISGlQQhAREUAJQUREopQQREQEUEIQEZEoJQQREQGU\nEEREJEoJQUqamf3TzD4zs+5x+042s7/GbbuZrTGz1WZWZ2bXmVlV3PG/Rs/ZLenaj0b37xvdrjWz\nO8zsEzP7ysz+ZmZj09wn9nN+nr9v12gMX0bjGJ3h3P3MbKGZ1ZvZv6Pfp09briUCSghSHqqAc7Kc\ns5u7bwj8F/ALYETS8b8Bx8c2zGwT4AfAirhzrgc2BP4D6AkcCryX6j5xP3/I9ctkcTmwA9Af2A84\n38wOTHPuO8BPgF7AVsDfgVvbeC0RJQQpC5OA35pZbbYT3f09YDawe9KhqcAv4koOxwKPAt/GnbMn\ncK+7r3T3Znd/190fan/4Ofk18PtoDEuAycAJqU5090/dfZmvn26gCdi+LdcSASUEKQ9zgL8Cv812\nopntCOxN6zf75UTeqH8c3T4euDvpnNeACWZ2opnt0J6AzWxstCon5U+az/QCtgTmx+2eD+yc4T79\notdrIPL384e2XktECUHKxaXA2WbWO83xt8xsDbCESPK4JcU5dwPHR5NGrbu/mnT8bCIlibOAd8zs\nPTM7KMV94h/uw1MF4+4T3b023U+a77Bh9M9Vcfu+BHqkOR93/1f0epsCFwPvtvVaIkoIUhbcfREw\nHRib5pQ9iDwEfwEMBbqnOOcR4IdEHvj/m+IeDe5+pbt/F9gEeAB40Mw2jr9P0sN9Zpu/VGuro39u\nFLevJ/BVtg+6+xfAn4HHzKxze64lHZcSgpSTy4CRQJ9UBz3iAeBVIiWK5ONfA08Bp5MiISSd+yVw\nJZHEMiDXQM3swqTeSAk/ae65EvgYiO8NtRuwOOBtOwObARvl4VrSASkhSNmINhjfD/wmy6kTgZFm\ntkWKYxcC/+Xu/0w+YGaXmNmeZtbFzLoR6dlUDyxtQ6xXJvVGSvjJ8NG7gYvNrJeZ/QeRBHhXqhPN\n7EgzG2hmnaJVadcB86KlhZyuJQJKCFJ+fkfq6qAW7r4QeAkYk+LYcnd/Od1HgTuBz4k0Qv8IONjd\n49/o5ye97f+xLV8ig8uA94EPibSF/MHdn44djN5z7+hmH+BpItVAC4Fm4Iig1xJJZlogR0REQCUE\nERGJUkIQERFACUFERKKUEEREBIj0Wy4bm266qW+zzTZhhyEiUlbmzp37ubunG+XfoqwSwjbbbMOc\nOXPCDkNEpKyY2YdBzlOVkYiIAEoIIiISpYQgIiKAEoKIiEQpIYiICKCEICIiUUoIIiICKCGIiJS0\nW265hdmzZxflXmU1ME1EpKN46623+O53vwvAtttuy/vvv1/weyohiIiUEHenU6fEypu33367KPdW\nlZGISIm44447EpLBEUccgbvTo0ePotxfJQQRkZA1NzdTVVWVsK++vp6ePXsWNQ6VEEREQjR+/PiE\nZHDqqafi7kVPBqASgohIKL755htqampa7evatWtIEYVYQjCzrc3sBTN7x8wWm9k5YcUiIlJMI0eO\nTEgGV155Je4eajKAcEsI64Dz3P0tM+sBzDWzZ939nRBjEhEpmE8++YQtt9wyYV9TU1OrXkVhCS0K\nd//Y3d+K/v4VsAToE1Y8IiKFZGYJyeCuu+5K2cU0TCXRhmBm2wCDgddTHDsFOAWgX79+RY1LRKS9\nnnvuOX70ox8l7HP3kKLJLPTUZGYbAg8Do9z9y+Tj7j7Z3Ye4+5DevbMuCSoiUjLMLCEZ/OY3vynZ\nZAAhlxDMrJpIMpjq7o+EGYuISL7ccsstnHnmmQn7SjkRxISWEMzMgCnAEne/Lqw4RKRjmzavjkkz\nl7K8voGtamsYM3wghw9ue3Nm5NG23qWXXsoVV1zR3jCLIswqo2HAr4Afmtnb0Z+fhBiPiHQw0+bV\nMe6RhdTVN+BAXX0Do+5/m8G/e4Zp8+pyutbpp5/eKhm4e9kkAwixhODuLwOW9UQRkQKZNHMpDY1N\nrfav/LqRcY8sBAhUWkhOBE8//TTDhw/PT5BFFHqjsohIWJbXN6Q91tDYxKSZS1MemzavjmETZ9Gp\numvKUkE5JgMokW6nIiJh2Kq2hroMSSFVwpg2r46xDy9g6YTEGu6r736C7+wymGETZ7W7PSLf7RpB\nKSGISIc1ZvhAxj2yMGW1EUQSRrIj9ujbal//C6Zz7wfVrP3b+mvV1TfkVO0UE2vXaO912kIJQUQ6\nrNgD9vLHF1Pf0JhwrKa6ijHDB7Zsr169utW6BFudejvVtVsAtPo8rK92yuVBnqpdo6GxifMemJ8Q\ncyEoIYhIh3b44D4cPrhPxmqa5HYCiJQKgsjUTpHL+U3uBS8pKCGISMULUicfSwzx55/1p6eomzwy\n4bx7XnyH8c9+2Oot3gxSjT1LVe2USaZ2jbaUOHKhhCAiZS3bwz7XOvnY+e+OP6jVsUff+ojDB/eh\ne4+NuOKJxaz8en01UapkkFztFES2do1cSxy5UEIQkbIV5GGfrk4+3Zv2RTffy7tTzk/Y12/MY1in\nqoRrT5q5NCEhxFSZ0eze5t5BsfPPe2A+TSmyTK4ljlwoIYhI2QrysE/3Rl1X39Cqi2i6HkSprp2p\nrr9PO7uKxj6XXFJoS4kjF0oIIlK20j2U4/enq5M3aNm/5Jl7OWLclITj6RqNY9fOVNcfX1IB2jSm\nIL6EU6zxCEoIIlK20j2U46tVUtXJGxCrjPnw6kMSPtt7y75sctJtWccmZKvrb2hs4vLHF7N2XXOb\nxxTEN3QXg6auEJGyNWb4QGqqqxL2JVerHD64D1cdOYg+tTUY0Ke2Bge+eO5/WiWDbS6YzmfLl3HV\nkYPotUF1q/vFXzv+uunUNzSmrdIqRSohiEjZClqtkvymnTyuYKOhR9Fr3xNa3v7jxybE9ybq2jnx\nHTp23rCJszJOgZGskD2F2kMlBBHJWWxytwFjn2TYxFk5TxWdzzhyqWPfc889WyWD/hdMp9e+J6Rt\nsP2msbnl9/qGxpTTY6crqaQqZUBhewq1h0oIIpKTMOfaaU8cyYng3Ak38nqnnTImk6DTY6crqUDr\nnkLVVcaatesYMPbJok5cF4QSgojkJNd+/WHH0a9fP5YtW5ZwXtDlLINMj51t2otYrMvrG6jdoJrV\n36xrmfcorGSajqqMRCQnQbp6lkIcTU1NmFlCMli4cGFOaxtnq9pZXt+QctW1cY8sbKlSOnxwH2aP\n/SH/mHgwG3TpTGNz4v1LqZFZCUFEcpLuIVnsevFMcZgZnTsnVoC4O7vssktO90jVNpB8r0wllZhY\nm0u6hudSaWRWQhCRnATp6hlWHF2aGnhl3P4J+z777LOcSgXxYl1La6pTPyr327F31pJKfAkinVJp\nZFYbgojkJIwRtEHi+GfSmAII3lYA6XssxeYtSvVAf+HdFVkHx6VrmI4xKHoyTUcJQURyVuwRtJni\n2L1XIwMGDEjY39DQQLdu3VJ+JtWDH8jYYylTKeD6X+yecc6hbNVBTmk0KIMSgojkWTHXA061cE2m\nUkG6rqrdqjtl7LGUqRSQrcSUbd3mTCOdi00JQUTyJl9jFLIllZkzZ3LggQcmfKa5uTllgoiXrgE4\n29oDqeYtSp7GIt33yzTnURhtL5moUVlE8iZIj5tssnXjNLOEZLD33nvj7lmTAeTemyd+Kovk+ZCu\nOnJQ4Anq4uc8qorGmcs1ikUlBBHJm3yMUUiXVM69eDxHzPhTwv5cew+lq76pralOmJUUUk+Sl7wS\nW/J6Cuke7qXS5pJNqCUEM7vDzD4zs0VhxiEi+ZGPMQqpkseHVx/CP+OSwVFHHdWmrqTpusxefujO\nOZUAspViylXYJYS7gJuAu0OOQ0TyIFtdexDxb/ErHp/E10teTDje1jEFkL3LbNC3+FKZviPfQk0I\n7v6SmW0TZgwikptMDb6ZJnnLVr0Su25dfQMGrcYVHP+bC/nzDRPaHX8+qm9KZfqOfAu7hJCVmZ0C\nnAKRSapEJDxBehGlqmtP9xkgIQk4rVcwA3j0rY9K6s07yEpt5ajkexm5+2R3H+LuQ3r37h12OCId\nWlt6EaX7zBVPLE6Y0qHZvVUy6HfUONy9pJIBlM70HflW8iUEESkdmapK0lUlpRuUFVuFDFKXCmKL\n3E+bV5e2eimsqTNKZfqOfFNCEJHA0lWV9KypTlstVGVGU5qG4ObGtSy77mcJ+7Y84b/psvm2Ldvn\nPTAfWP8QLpUFesqlK2kuwu52eh/wKjDQzD4ys5PCjEdEMktXVWJG2qqkdMngw6sPaZUM+l8wPSEZ\nADS5J3TpzMfgN0kt1ITg7se6+5buXu3ufd19SpjxiEhm6Ubs1sdV/8RbXt/Qaq6edV+uaFVF1PfM\n/22pIkol/oFfqT18SoGqjEQkJ6mqStJNDx2rW49V8WRqK8gm9sCv1B4+paDkexmJSOnL1Ovm8MF9\nOPE761olgwdf/wc7XvxU4HvEHvi59vCJTTExYOyTDJs4q+xHExeSSggi0m6Zet1kmqK6c3U15z0w\nP207Q0zyzKLp7pWsVBqgy4USgoi0WabunzfffDO2x1kJ5ydPOxE7N3m6i+oqo3uXzqxqaEz5wA/a\nw6dSp5goFCUEEWmTTG/fR+zRN+Hcvn37smzZspTXyXef/vgkla7coQbo1AIlBDPbBdgJaFmTzt01\nIZ1IB5bq7fsffx7DEeMTJy8OMhldvvr0JyepdNQAnVrWhGBmlwH7EkkIM4CDgJfRDKUiHVryW3Zy\no3G3/rvT97gJrUYaF3KUcbYF7aEyppgolCAlhKOA3YB57n6imW0O3FPYsESk1MW6f2bqStrY5Fzx\nxOKijTLOVBVkUDFTTBRKkITQ4O7NZrbOzDYCPgO2LnBcIlJkqd7cAa54YnHLvEO1NdVcfujOHD64\nD2OGD2zVVtB90AFs+pNRCfvi5ywqdCNvujEKfWprmD32h+2+fqULkhDmmFktcBswF1hNZLoJESkx\nba2OSfXmPubB+TQDTc3r2wDqGxoZ8+D8VokAgg0wK/Qo43ws0NORZU0I7n5G9Nc/mdnTwEbuvqCw\nYYlUlmLMztme6phUb+6Nza0bg92bee+qQxP2DTj0bJr/Y3jaa9fWVLf8XuhRxpU6C2mxBGlUft7d\n9wdw938m7xORzIo1OKo91TFB3tBTtRW4e8aePdWdjMsP3blluxhv8JU4C2mxpJ26wsy6mdnGwKZm\n1svMNo7+bAPob1tKUilOU1Cs2TnbUh0T+/vK1DG0ee2aVslgt3Nub+lOGj/hHUSmu4ZIvf2kn+/W\nakBZLovZS3FlKiGcCowCtiLSdhAbf/4lcFOB4xLJWalOU1Cs2TnTVcc4kfWMk6tOLp62kKmv/Stt\nMugE/CNFqWD7cU9y+c93S9iXy1t5WG/wYS+qUw7SJgR3vwG4wczOdvcbixiTSJuU6jQFudSbt+eh\ntd+OvdM+4JOT47R5dRmTQbfVH7P05pEJ+7Ye9QDdN9yQq47ctewepKX6slBqgjQq36iRylIOSnWe\n/KD15u15aE2bV8fDc+syVv3EJ8dJM5emPTfTuIKNu3ct2AO02APWSuFlodRopLJUjFKdJz9dzxeA\n3a94hvqGSD/9TgbJHXuCPrSCjNCFSJKZNq8uZZL8+u+vs+KR3yfs63f+45itb2osVHINa8Ba2C8L\npUYjlaVilHIf9OR682nz6hjz4PyErp0penkCwR5auTzYxj2ykNoNqgMvch+vUMk1rAFrYb8slJog\nC+Q0uHszoJHKUtLKqQfLpJlLU/bzTyXIQyuXB1tDYxPukWT5+ZPXt0oG7s6jb32U0yI07VWMAWvF\n/D7lSiOVpaKUSx/0oA+6oA+tVKWj6iqjsSl10lnV0Mg/04wrgOIP8NKAtdJgQaambTk5MgYhtJHK\nQ4YM8Tlz5oRxa5G8GjZxVsoHIET68Te7B3poxTfE1m5QjTsJi8qkWut42R9/QfPaNQn7Hn3ro7w/\nLHNpJE41uK2muqpkS3jlxszmuvuQbOelLSGY2R6Zjrn7W20NTqSjGzN8YKs2BIi81U86ardAD8Hk\nh+jKrxupqa7i+l/snvD5+HOSq4f22GMPLrn98bw36ObaSKw3+NKQtoRgZi9Ef+0GDAHmExmctisw\nx91/UJQI46iEIJVk2rw6Ln98cUsvo14bVHPZT3cO/BBMV8pIntlz2ry6lJPRxf7fD3qdXBTimtJ2\n7S4huPt+0Qs9Auzh7guj27sAl+cpTpEOq73tHUEbYpOTwbhx47jyyitzvk4hYpPSEqRReWAsGQC4\n+yIz+4983NzMDgRuAKqA2919Yj6uK9IRZGuINbNWx1LVCBSiQVfdPMtTkG6nC8zsdjPbN/pzG9Du\nRmUzqwJuJjLQbSfgWDPbqb3XFeko0nWlHLXfNq2SweOPP552beNCdMlUN8/yFKSEcCJwOnBOdPsl\n4NY83Pt7wHvu/gGAmf0FOAx4Jw/XFql4qRpiXxm3P0ePTzwvW0/CQjToqpG4POXU7TSvNzY7CjjQ\n3U+Obv8KGOruZ6X7jBqVRVJbvnw5ffokPmz//ve/s/322wOa6bOja3ejcqkws1OAUwD69esXcjQi\nbVeoh3K2toKgXUCVNCRIG0Kh1JE4BUbf6L4E7j7Z3Ye4+5DevXsXLTiRfIo9lOvqG3DWP5Tbs4DP\n7NmzWyWDNWvWtKoiCrJATyHik/KTNSGY2c+D7GuDN4EdzGyAmXUBjgEez8N1RUrOFU8szuuqaWbG\nf/7nfybsc3c22GCDVucG6QJarFXdpLQFKSGMC7gvJ+6+DjgLmAksAR5w98Xtva5IqZk2ry5hZtF4\ndfUNOS31ecstt7QqFTQ3N2dsOE7X1TN+v8YNCGSeuuIg4CdAHzP777hDGwHr8nFzd59BZI0FkYqV\n7S076FQRQccVJAsyLXiQcQNqY6h8mUoIy4E5wDdEZjmN/TwODC98aCKVIchbdqbqmWOOOaZVMnD3\nlMlg2rw6hk2cxYCxT7aUPIJMC55t3IDaGDqGTFNXzAfmm9m97t4IYGa9gK3dfWWxAhQpd+nevpOl\nShzJiWDzzTfnk08+Sfn5bL2JMr3NZxs3oCUoO4Yg3U6fNbNDo+fOBT4zs1fc/dzChiYSjnxXjaSq\nskklvnpmww03ZM2axCmq97rqeZZH2xxSxdTeh3ampNGWNgZVMZWfIAmhp7t/aWYnA3e7+2VmFsp6\nCCKFVoi1fePfsjOVFFauWZtyZtL9fvpzPt5tRMtn08VUyIbhbG0MyQ///XbszcNz6wq2RrIURpBe\nRp3NbEvgaKD1IqsiFaSt3S9T1d3HO3xwH2aP/SGtm4XXWzL+J62Sgbuzdq/TAsUUpDdRW2VqY0jV\nvjD1tX+pG2sZCpIQfkeka+j77v6mmW0L/L2wYYkUT/zDPN0bfLaqkaANrqkezu7eauGa++67r6XR\nOOibfyEnlMvUMJ0qiabr+6RurKUta5WRuz8IPBi3/QHws0IGJVIsqZZuTCXTW3Yudfdjhg9k1P1v\nt2wnJwKAbS6YzjHHHJxw7yBTSRd6Qrl0bQy5POQ1/XVpCzJS+Ttm9ryZLYpu72pmFxc+NJHCS/Uw\nT5btLTvXuvuqTkZz4zetksGWJ/43/S+Y3uqhmcubf6xq6vpf7A7Aufe/ndPAt7ZI95BPrh7T9Nel\nL0iV0W1ERiY3Arj7AiLTTIiUvUxvt+n67CcLWnc/bV4d5z0wnw+uOphl1x2VcKz/BdPpstm2KR+a\nQcYRJN8n1zED2dpAMkmXsI77fr/AMUtpCNLLaAN3fyOpP3ReRiqLhC1ddUyvDaqZd+mPA10jyEjg\nafPqGHXbM3x464iEz/Y96x6qutcCkYdmuiqeXJbbzLX7aXt7Vmntg8oRJCF8bmbbEW0niq5j8HFB\noxIpkjHDBzLmofk0NiU2g67+Zh0XT1vIC++uyPqQC/JATLXIff8L1nfay+fi87lWYeVj0Fl714eW\n0hAkIZwJTAZ2NLM64B/AcQWNSqRIDh/ch8sfX0x9Q+Lkc43NztTX/tXSWybbW3O6B+K8efPYY489\nEvb1G/MY1ml9FUs+69anzaujkxlNadZOTkUT20lMkITg7n6AmXUHOrn7V2Y2oNCBiRTLqqRkEJP8\nSM321pw8OOuVcfsnHLfOXel33sOtPpevuvVY1U+qZJAp6QTtxSSVL0ij8sMA7r7G3b+K7nuocCGJ\nFFcuD750b83xDbmrl85ulQx2v2JmymTQa4PqvFW1pOsxVWWWMekUcvyClJdM01/vCOwM9DSzI+MO\nbQR0K3RgIsUSdK4hSJ88Yg/j5K6ktTvsycq/vcG0eXWt2iqqq4zLfrpz+4KPky5ZNbu3a2I76Tgy\nVRkNBA4BaoGfxu3/ChhZyKBE8iXIBGux7fMemJ+yuiUm01vz0r8+wr9n3pywr/8F01v64hfjodue\nqh81Cgtknv76MeAxM/uBu79axJhE8iKX7pSHD+7DuXEjiJNl6hKaPEX1Rj/4Bb32+RWQ+DAu9EM3\nSPdXkUyytiEoGUi5ynWiunRv0rEuockP81GjRrVKBv0vmN6SDIr9MM51AJtIsiC9jETKUq7dKXN5\nw05OBE899RTfbD4o9Hp4Vf1IeyghSMXKtU49SD3/3nvvzcsvv5zwufilLPUwlnKWqZfR6EwfdPfr\n8h+OSP60pU493Rt2c3MzVVWJXTMXLVrEzjvnr5eQSNgylRB6RP8cCOwJPB7d/inwRiGDEsmHfPXs\nqaqqorm5OWFfqgXuRcqdZfuHbWYvAQfHBqWZWQ/gSXffpwjxJRgyZIjPmTOn2LeVDurrr7+me/fu\nCfs+/fRTNttss5AiEmkbM5vr7kOynRekDWFz4Nu47W+j+0QqVnKjMahUIJUvSEK4G3jDzB6Nbh8O\n/LlwIYmE5/PPP6d3794J+xoaGujWTYPzpfIFGYcwATgRWBn9OdHdr2zPTc3s52a22MyazSxrMUak\nGMwsIRkccMABuLuSgXQYQbudbgB86e53mllvMxvg7v9ox30XAUcC/9OOa4jkRV1dHX37Jq5X0Nzc\nnLLaSKSSBVlT+TLgAiLLaAJUA/e056buvsTdUw8XFSkiM0tIBqNGjcLdlQykQwpSQjgCGAy8BeDu\ny6M9jYrCzE4BTgHo169fsW4rFW7hwoXsuuuuCfvUaCwdXZD1EL71yP8psSU0u2c5n+h5z5nZohQ/\nh+USoLtPdvch7j4kubFPpC3MLCEZPPHEE0oGIgQrITxgZv8D1JrZSGAEcHu2D7n7Ae0NTiSfZs2a\nxf77Jy5co0Qgsl7WhODu15jZj4AviYxavtTdny14ZCJ5lNwm8MYbb7DnnnuGFI1IaQrSqHy1uz/r\n7mPc/bfu/qyZXd2em5rZEWb2EfAD4Ekzm9me64mk88Ybb7RKBu5OXeetGDZxFgPGPsmwibOYNq8u\npAhFSkeQNoQfpdh3UHtu6u6Puntfd+/q7pu7+/D2XE8klS5dujB06NCW7ffffx93T1j/2Fm/cI6S\ngnR0aROCmZ1uZguBHc1sQdzPP4CFxQtRJDdPP/00ZkZjYyMAo0ePxt3ZdtttgdwXzhHpKDK1IdwL\nPAVcBYyN2/+Vu39R0KhE2sDdueOOOzj55JNb9n311VdsuOGGCeflunCOSEeRtoTg7qvc/Z/ADcAX\n7v6hu38IrDOzoek+JxKGDz/8kIMOOqglGVx77bW4e6tkAOkXyAmyGL1IJQvShnArsDpue3V0n0jo\nmpubuemmm9h55515+eWXuemmm2hqamL06PTrO40ZPpCa6sTFbrQYvUiwcQjmcZ213b3ZzLT0poRu\n6dKlnHzyybz88sv8+Mc/ZvLkyfTv3z/r5/K1cI5IpQnyYP/AzH7D+lLBGcAHhQtJJLN169Zx7bXX\nctlll1FTU8Odd97Jr3/965zmH9Ji9CKtBakyOg3YC6gDPgKGEp1bSKTY5s+fz9ChQxk7diwHH3ww\nS5Ys4YQTTtBkdCJ5EGSk8mfAMUWIRSrAtHl1BamKWbt2LePHj2fixIlsvPHGPPjggxx11FF5iFhE\nYtImBDM7393/YGY3Ep3YLp67/6agkUnZiQ34ivXxjw34AtqVFF577TVGjBjBkiVLOP7447nuuuvY\nZJNN8hKziKyXqcpoSfTPOcDcFD8iCfI94GvNmjWce+657LXXXqxevZoZM2bw5z//WclApEDSlhDc\n/Ynon1o/WQLJ54CvWbNmMXLkSD744ANOP/10Jk6cyEYbbdTeEEUkg0xVRk+Qoqooxt0PLUhEUra2\nqq2hLsXDP5cBX6tWreK3v/0tt99+O9tvvz0vvvgi++yzTz7DFJE0MlUZXQNcC/wDaABui/6sBt4v\nfGhSbto74OuJJ55gp5124o477mDMmDEsWLBAyUCkiDJVGb0IYGbXuvuQuENPmNmcgkcmZaetA75W\nrFjBOeecw3333cegQYOYNm2a1ioQCUGQgWndzWxbd/8AwMwGAIGW0ZSOJ5cBX+7O/fffz9lnn82q\nVau44oorGDt2LF26dClwlCKSSpCEcC7wVzP7ADCgP3BqQaOSihQ/RmETW02nV+/gzRef4Xvf+x5T\npkxhl112CTtEkQ4tyMC0p81sB2DH6K533X1tYcOSShMbo/D1t418NecJPpx9LzSt44RzL+H2SZdR\nVVWV/SIiUlBZE4KZbQCMBvq7+0gz28HMBrr79MKHJ5Vi0sylfDZnBl88fSMAXbfakU0OGc3fNttO\nyUCkRASpMrqTyEC0H0S364AHASUECWTt2rW8Mm7/hH2b/78/YNZJi9KIlJAgk9tt5+5/ABoB3P1r\nIm0JIllNmDCBbt26tWxvethY+l8wHbPIPz0tSiNSOoKUEL41sxqig9TMbDtAbQiS0apVq6itrU3Y\nN/CiGXwnSjAHAAAP5UlEQVSzrrllW4vSiJSWICWEy4Cnga3NbCrwPHB+QaOSsnbfffclJIPnn38e\nd2fiz3alT20NBvSpreGqIwdpTQKREmJxi6G1PhiZZL4v8DXwfSJVRa+5++fFCS/RkCFDfM4cjYkr\nVd988w1XXHEFEydOBKC2tpaVK1eGHJWImNncpAHGKWUsIUSXzpzh7v929yfdfXpYyUBK28svv8xu\nu+3GxIkTGTFiBF988YWSgUiZCVJl9JaZaR4BSWn16tWcffbZ7LPPPnz77bc888wzTJkyhV69eoUd\nmojkKEij8lDg/5nZP4E1RKqN3N13betNzWwS8FPgWyIT5Z3o7vVtvZ6E45lnnuGUU07hX//6F2ef\nfTYTJkxgww03DDssEWmjIAlheAHu+ywwzt3XmdnVwDjgggLcRwpg5cqVjB49mrvuuouBAwfyf//3\nfwwbNizssESknTKth9ANOA3YHlgITHH3dfm4qbs/E7f5GqDFccvEo48+yhlnnMGKFSu48MILueSS\nSxLGGYhI+crUhvBnYAiRZHAQkbURCmEE8FS6g2Z2ipnNMbM5K1asKFAIks2nn37K0UcfzZFHHskW\nW2zBm2++2WrQmYiUt0xVRju5+yAAM5sCvJHLhc3sOWCLFIcucvfHoudcBKwDpqa7jrtPBiZDpNtp\nLjFI+7k7U6dO5ZxzzmH16tVMmDCBMWPGUF1dHXZoIpJnmRJCY+yXaF1/Thd29wMyHTezE4BDgP09\n02AICc2yZcs47bTTmDFjBnvttRdTpkxhxx13zP5BESlLmRLCbmb2ZfR3A2qi27FeRm1e8dzMDiQy\n2vm/onMjSQlpbm5m8uTJnH/++TQ1NXHDDTdw5plnalZSkQqXaQnNQv7ffxPQFXg2WvJ4zd1PK+D9\nJKD33nuPk08+mRdffJEDDjiAyZMnM2DAgLDDEpEiCNLtNO/cffsw7ivpNTU1cf3113PJJZfQtWtX\npkyZwoknnkiuVYUiUr5CSQhSWhYtWsSIESN48803Oeyww7jlllvYaqutwg5LRIpMCaHMxK9LvFVt\nDWOGD2zzjKHffvstV111FRMmTKC2tpa//OUvHH300SoViHRQSghlJLYucUNjEwB19Q2Me2QhQM5J\n4c0332TEiBEsWrSI4447jj/+8Y9suummeY9ZRMpHkMntJINp8+oYNnEWA8Y+ybCJs5g2r65g95o0\nc2lLMohpaGxi0sylga/x9ddfM2bMGL7//e+zcuVKpk+fzj333KNkICIqIbRHPt/Yg0i3/nDQdYlf\nfPFFTj75ZN577z1OPfVUrr76anr27JnPEEWkjKmE0A75eGPPRbr1h7OtS/zll19y+umns+++++Lu\nzJo1iz/96U9KBiKSQAmhHdr7xp6rMcMHUlOdODwk27rETz31FLvssguTJ09m9OjRLFiwgP32268g\n8YlIeVOVUTtsVVtDXYqHf7Y39raKVUNl62U0bV4dVz7yBu88eiNrFr/A1tsN5JVXXmHo0KEFiUtE\nKoMSQjuMGT4woQ0Bsr+xt9fhg/tkbJ949K2POPK7W7ds99zrWGr2OYaPu/QtWEwiUhlUZdQOhw/u\nw1VHDqJPbQ0G9Kmt4aojBxWkQTmIGTNmJCSDTX4yitq9j2OtVxWsXUNEKodKCO2U7Y29GJqbm1tN\nPLf1uQ/Rqcv6tQoK1a4hIpVDJYQyd/fddyckg20OPp3+F0xPSAZQuHYNEakcKiGUqW+//ZauXbu2\n2vfkos+K3q4hIpVBJYQyNGfOnIRkMHXqVNyd6urqkmvXEJHy0SFKCPmcEC5MjY2NXHPNNVx44YUt\n+5qbm1tNRlcK7RoiUn4qvoQQm16irr4BZ/30EoWcc6gQ3n77bYYOHcqFF17Iz372Mz7++GPcXTOT\nikjeVHxCKPb0Evm2du1aLr74Yvbcc0+WL1/OQw89xEMPPcQWW2wRdmgiUmEqvsqo2NNL5NOrr77K\nSSedxJIlS/j1r3/Nddddx8Ybbxx2WCJSoSq+hNDWCeHCtGbNGkaNGsWwYcNYs2YNTz31FHfddZeS\ngYgUVMUnhLZMCBem559/nkGDBnHDDTdwxhlnsGjRIg488MCwwxKRDqDiE0K5dMOsr69n5MiRHHDA\nAXTu3JmXXnqJm266iR49eoQdmoh0EBXfhgCl3w3z8ccf5/TTT+eTTz7h/PPP5/LLL6empnSrtESk\nMlV8CaGUrVixgmOPPZbDDjuMTTfdlNdff52rr75ayUBEQqGEEAJ357777mOnnXbi4Ycf5ne/+x1v\nvvkmQ4YMCTs0EenAQkkIZvZ7M1tgZm+b2TNmtlUYcYShrq6OQw89lF/+8pdst912zJs3j0suuYQu\nXbqEHZqIdHBhlRAmufuu7r47MB24NKQ4isbdue2229hpp514/vnnue6665g9ezY777xz2KGJiAAh\nNSq7+5dxm90BDyOOYnn//fcZOXIkL7zwAvvttx+33XYb2223XdhhiYgkCK0NwcwmmNky4DgylBDM\n7BQzm2Nmc1asWFG8APOgqamJ66+/nkGDBjF37lwmT57M888/r2QgIiXJ3Avzcm5mzwGpJty5yN0f\niztvHNDN3S/Lds0hQ4b4nDlz8hhl4bzzzjucdNJJvPbaaxxyyCHceuut9O2rdY1FpPjMbK67Z+21\nUrAqI3c/IOCpU4EZQNaEUA4aGxuZOHEi48ePp0ePHkydOpVjjz1Ws5KKSMkLpQ3BzHZw979HNw8D\n3g0jjnybO3cuI0aMYMGCBRxzzDHccMMNbLbZZmGHJSISSFgjlSea2UCgGfgQOC2kOPJi1apVXHnl\nlVx77bVsttlmPPbYYxx66KFhhyUikpOwehn9LIz7FsJVV13VsoLZSSedxDXXXENtbW3IUYmI5K5D\nzGVUCB9//DFbbbV+PN1FF13E+PHjQ4xIRKR9lBDaYPTo0Vx//fUt259++qnaCkSk7Gkuoxy8//77\nmFlLMrjmmmtwdyUDEakIKiEE9Ktf/Yp77rmnZbu+vp6ePXuGGJGISH6phJDFggULMLOWZHDnnXfi\n7koGIlJxVEJIw90ZPnw4zz77LAA9e/bk448/1loFIlKxVEJIYfbs2XTq1KklGUybNo36+nolAxGp\naCohxGlqamKPPfZgwYIFAHznO99h8eLFdO6svyYRqXwqIUTNmDGDzp07tySDv/71ryxdulTJQEQ6\njA7/tFu7di39+vXjs88+A2CfffbhhRdeoFMn5UoR6Vg69FPv3nvvpVu3bi3JYO7cubz44otKBiLS\nIXXIEsLq1avp0aNHy/bRRx/NX/7yF01RLSIdWod7Fb7pppsSksHSpUu5//77lQxEpMPrMCWEzz//\nnN69e7dsn3XWWdx4440hRiQiUlo6REJobGxMSAbLli3TcpYiIkk6REKoqqrijDPOYLPNNuOyyypi\npU4RkbzrEAmhU6dO3HzzzWGHISJS0jpco7KIiKSmhCAiIoASgoiIRCkhiIgIoIQgIiJRSggiIgIo\nIYiISJQSgoiIAGDuHnYMgZnZCuDDsONoo02Bz8MOIo8q6ftU0neByvo+lfRdILzv09/de2c7qawS\nQjkzsznuPiTsOPKlkr5PJX0XqKzvU0nfBUr/+6jKSEREACUEERGJUkIonslhB5BnlfR9Kum7QGV9\nn0r6LlDi30dtCCIiAqiEICIiUUoIIiICKCEUlZn93swWmNnbZvaMmW0VdkxtZWaTzOzd6Pd51Mxq\nw46pPczs52a22MyazaxkuwVmYmYHmtlSM3vPzMaGHU97mNkdZvaZmS0KO5Z8MLOtzewFM3sn+u/s\nnLBjSkUJobgmufuu7r47MB24NOyA2uFZYBd33xX4GzAu5HjaaxFwJPBS2IG0hZlVATcDBwE7Acea\n2U7hRtUudwEHhh1EHq0DznP3nYDvA2eW4n8fJYQicvcv4za7A2Xbou/uz7j7uujma0DfMONpL3df\n4u5Lw46jHb4HvOfuH7j7t8BfgMNCjqnN3P0l4Iuw48gXd//Y3d+K/v4VsAToE25UrXWINZVLiZlN\nAI4HVgH7hRxOvowA7g87iA6uD7AsbvsjYGhIsUgGZrYNMBh4PdxIWlNCyDMzew7YIsWhi9z9MXe/\nCLjIzMYBZwGXFTXAHGT7LtFzLiJSHJ5azNjaIsj3ESkkM9sQeBgYlVRjUBKUEPLM3Q8IeOpUYAYl\nnBCyfRczOwE4BNjfy2BASw7/bcpRHbB13Hbf6D4pEWZWTSQZTHX3R8KOJxW1IRSRme0Qt3kY8G5Y\nsbSXmR0InA8c6u5fhx2P8Cawg5kNMLMuwDHA4yHHJFFmZsAUYIm7Xxd2POlopHIRmdnDwECgmcg0\n3qe5e1m+xZnZe0BX4N/RXa+5+2khhtQuZnYEcCPQG6gH3nb34eFGlRsz+wnwR6AKuMPdJ4QcUpuZ\n2X3AvkSmi/4UuMzdp4QaVDuY2X8C/wcsJPL/P8CF7j4jvKhaU0IQERFAVUYiIhKlhCAiIoASgoiI\nRCkhiIgIoIQgIiJRSghStsxsk+jMsW+b2SdmVhe33SWP99kjOu6iLZ/d2MwCdcc1s/FmNirLOUea\n2Y5tiUUkGyUEKVvu/m933z06e+yfgOtj29EJ3rCI9v4734O2z7y5MZDP8RlHAkoIUhBKCFJxzGz7\n6LzzU4HFwNZmVh93/Bgzuz36++Zm9oiZzTGzN8zs+0nXqiEyTflx0ZLHUWa2oZndFT1/npn9NHru\nIDN7M3reAjPbFpgIDIzum5gi1kvN7G9m9jKwQ9z+06LXmm9mD5pZjZntDfwEuD56vW1SnZf3v1Dp\nMDSXkVSqHYHj3X2OmWX6d/7fwB/c/bXoLJTTgV1iB929wcx+R2Tth1EAZvYH4Gl3P8HMegGvm9mz\nwBnANe5+v5l1BQwYC2wfLcUkMLPvAT8DdgO6AG8Dr0YPP+juf4qeNxE4wd1vNbMZwEPuPi16rNV5\nwK25/mWJgBKCVK733X1OgPMOIPIGH9vuZWY17t6Q4TM/Bg6KW5WsG9APeAW42Mz6A4+4+3tx101l\nH+Dh6L0azOyJuGO7RhNRLdCDSKJKJeh5IlkpIUilWhP3ezORt/WYbnG/G/C9WJtDQAYc7u7vJ+3/\nm5m9ChwMPG1mI4DlOVw33t3AQe6+yMxOJrLKVnvOE8lKbQhS8dy9GVhpZjtEG5iPiDv8HHBmbMPM\nWlXtAF8RefuOmQmcHfeZwdE/t3X399z9BiJv6rum+Gy8l4AjzKybmW1EZCrxmO7AJ9Epk3+ZIZZ0\n54nkTAlBOooLiDzIXyGymljMmcCwaCPwO8DIFJ+dBewWbUA+CrgC6G5mC81sMXB59LxfWmQB9beB\n7wD3uPunwNzouQmNyu7+BvAosAB4Engj7vClRKa0ng28E7f/PuDCWKNyhvNEcqbZTkVEBFAJQURE\nopQQREQEUEIQEZEoJQQREQGUEEREJEoJQUREACUEERGJ+v+zCWsuVeV5eAAAAABJRU5ErkJggg==\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7f1c6b1fec90>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
Expand Down

0 comments on commit 201a4e6

Please sign in to comment.