{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "%matplotlib inline\n",
    "import io\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "from scipy.cluster import hierarchy\n",
    "from scipy.spatial import distance"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Downloaded 'dissimilarities.npz'.\n"
     ]
    }
   ],
   "source": [
    "%download https://examples.obspy.org/dissimilarities.npz"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "url = \"dissimilarities.npz\"\n",
    "with io.BytesIO(open(url, \"rb\").read()) as fh, np.load(fh) as data:\n",
    "     dissimilarity = data['dissimilarity']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([[ 0.        ,  0.93888207,  0.96150459, ...,  0.96971566,\n",
       "         0.65356486,  0.73541168],\n",
       "       [ 0.93888207,  0.        ,  0.98167743, ...,  0.91176233,\n",
       "         0.60265653,  0.58080099],\n",
       "       [ 0.96150459,  0.98167743,  0.        , ...,  0.85596097,\n",
       "         0.93469123,  0.77658697],\n",
       "       ..., \n",
       "       [ 0.96971566,  0.91176233,  0.85596097, ...,  0.        ,\n",
       "         0.97555222,  0.65943323],\n",
       "       [ 0.65356486,  0.60265653,  0.93469123, ...,  0.97555222,\n",
       "         0.        ,  0.21069725],\n",
       "       [ 0.73541168,  0.58080099,  0.77658697, ...,  0.65943323,\n",
       "         0.21069725,  0.        ]])"
      ]
     },
     "execution_count": 18,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "dissimilarity"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "33"
      ]
     },
     "execution_count": 19,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "len(dissimilarity)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "dissimilarity = distance.squareform(dissimilarity)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([ 0.93888207,  0.96150459,  1.        ,  0.68538174,  0.97862019,\n",
       "        0.95732425,  0.88929751,  0.74416172,  0.8120546 ,  0.63233502,\n",
       "        0.68508224,  0.68878307,  0.82971361,  0.93221387,  0.99861152,\n",
       "        0.8815054 ,  0.88068856,  0.76391142,  1.        ,  0.97721853,\n",
       "        0.88318716,  0.84955382,  0.86369877,  0.98763165,  0.91343057,\n",
       "        1.        ,  0.84156442,  1.        ,  0.86324721,  0.96971566,\n",
       "        0.65356486,  0.73541168,  0.98167743,  0.78085595,  0.90033752,\n",
       "        0.90494312,  0.88419423,  0.96923988,  0.71998968,  0.93512688,\n",
       "        0.84986643,  0.91862085,  0.73623252,  0.90527974,  0.60132527,\n",
       "        1.        ,  0.82492802,  0.83972715,  0.61912734,  0.84393098,\n",
       "        0.98633653,  0.98452806,  0.76259394,  0.98720732,  0.82518932,\n",
       "        0.85128888,  1.        ,  0.92485085,  0.91238983,  0.93038732,\n",
       "        0.91176233,  0.60265653,  0.58080099,  0.91306048,  0.73423293,\n",
       "        0.86543878,  0.67716302,  1.        ,  0.95931128,  0.58181014,\n",
       "        0.71968672,  0.70648874,  0.6685183 ,  0.95320051,  0.98842509,\n",
       "        0.75611159,  0.52758392,  0.57014619,  0.65980726,  0.63492586,\n",
       "        0.79565094,  1.        ,  1.        ,  1.        ,  0.7841689 ,\n",
       "        0.6117646 ,  0.63130788,  1.        ,  0.88230793,  1.        ,\n",
       "        0.85596097,  0.93469123,  0.77658697,  0.85032548,  0.89224028,\n",
       "        0.58790815,  0.80479222,  0.96821465,  0.97745347,  0.60526953,\n",
       "        0.97708644,  0.9784995 ,  0.98046401,  0.76511608,  0.78267888,\n",
       "        0.73015123,  0.79286143,  0.79494161,  0.93774377,  1.        ,\n",
       "        1.        ,  0.99570001,  0.64777174,  0.74406892,  0.81944359,\n",
       "        0.7548781 ,  0.96143823,  0.80918253,  0.92018793,  0.79635791,\n",
       "        0.83017116,  0.8414229 ,  0.56429289,  0.65044684,  0.55269893,\n",
       "        0.44188927,  0.20488142,  0.52200707,  0.18073403,  0.1675657 ,\n",
       "        0.2129899 ,  0.90483297,  0.3831195 ,  0.65499954,  0.67533379,\n",
       "        0.45446101,  0.69960382,  0.84161212,  1.        ,  1.        ,\n",
       "        0.77143381,  0.85191795,  0.68257506,  0.59801396,  0.69865218,\n",
       "        0.69884547,  0.71579616,  0.68030409,  0.41243143,  0.50046058,\n",
       "        0.6485806 ,  0.77454182,  0.72265667,  0.5845938 ,  0.66623526,\n",
       "        0.60942704,  0.59818574,  0.52799929,  0.95400164,  0.81868374,\n",
       "        1.        ,  0.76452384,  0.6706793 ,  0.78791189,  0.91275991,\n",
       "        0.66950683,  0.66052778,  0.70392136,  0.66445939,  0.7875532 ,\n",
       "        0.65700753,  0.64884699,  0.72279969,  0.79475795,  0.61223415,\n",
       "        0.67304787,  0.58596424,  0.77602711,  0.77793322,  0.69631784,\n",
       "        0.15026725,  0.79026875,  0.80325117,  0.65473106,  0.81691666,\n",
       "        0.62044313,  0.55035417,  0.39784863,  0.68867162,  0.40246332,\n",
       "        0.93868431,  0.68824351,  0.60481299,  0.92112504,  0.7648137 ,\n",
       "        0.40721604,  0.23316495,  0.50911266,  0.49754708,  0.5996995 ,\n",
       "        0.48176781,  0.68963577,  0.75404799,  0.62741588,  0.49153694,\n",
       "        0.7032785 ,  0.43661662,  0.40371082,  0.46613514,  0.95871575,\n",
       "        0.63426474,  0.51362222,  0.60535603,  0.52962333,  0.6543014 ,\n",
       "        1.        ,  0.89691283,  0.84816973,  0.84709477,  1.        ,\n",
       "        0.59878968,  0.62989347,  0.83261837,  0.66968924,  0.63845773,\n",
       "        0.86176237,  0.52765558,  0.63713423,  0.36799243,  0.61850343,\n",
       "        0.38251971,  0.39445395,  0.42481065,  0.93068673,  0.54074455,\n",
       "        0.60453251,  0.68092583,  0.49502888,  0.7260657 ,  1.        ,\n",
       "        1.        ,  0.84658622,  0.84728396,  1.        ,  0.68150066,\n",
       "        0.7013547 ,  0.62382352,  0.81169668,  0.83847222,  0.85886975,\n",
       "        0.46391951,  0.55635821,  0.54545511,  0.1301536 ,  0.10193699,\n",
       "        0.12946733,  0.91936604,  0.33514303,  0.9217109 ,  0.56505263,\n",
       "        0.30643902,  0.61977694,  0.81924785,  1.        ,  0.82117565,\n",
       "        1.        ,  1.        ,  0.57035375,  0.55282904,  0.72367949,\n",
       "        0.76617387,  0.74889084,  0.83658191,  0.35461064,  0.49880495,\n",
       "        0.4955816 ,  0.56116595,  0.56758363,  0.7028769 ,  0.4673118 ,\n",
       "        0.9830755 ,  1.        ,  0.4593218 ,  0.97750609,  0.95329368,\n",
       "        0.71414558,  0.69370974,  0.73634846,  0.98689414,  1.        ,\n",
       "        1.        ,  0.74206515,  0.88981931,  0.77797241,  0.7608349 ,\n",
       "        0.46098103,  0.40075586,  0.06797162,  0.11862359,  0.89975063,\n",
       "        0.3378624 ,  0.47588603,  0.56875416,  0.32915934,  0.61759867,\n",
       "        0.8518421 ,  1.        ,  0.84827515,  1.        ,  0.84970493,\n",
       "        0.55410556,  0.77874446,  0.51555066,  0.74160474,  0.76918101,\n",
       "        0.8037141 ,  0.32750667,  0.48368975,  0.07343822,  0.90263383,\n",
       "        0.3882076 ,  0.40510246,  0.52763538,  0.34809961,  0.56664056,\n",
       "        0.85266674,  1.        ,  0.82484363,  1.        ,  0.85062975,\n",
       "        0.51763942,  0.78087713,  0.47654548,  0.93058261,  0.7678651 ,\n",
       "        0.80557119,  0.33992337,  0.50533063,  0.87620205,  0.40608619,\n",
       "        0.43308953,  0.50416051,  0.37417889,  0.58365476,  1.        ,\n",
       "        1.        ,  0.81221815,  1.        ,  1.        ,  0.51204949,\n",
       "        0.50497999,  0.41758213,  0.51313095,  0.75057162,  0.58867394,\n",
       "        0.38080921,  0.53036769,  1.        ,  0.95937965,  0.96950625,\n",
       "        0.92087194,  1.        ,  0.89956495,  0.66176862,  0.60361294,\n",
       "        0.91541492,  0.72233389,  0.99052916,  0.93150791,  1.        ,\n",
       "        0.77277337,  1.        ,  0.83554832,  0.92709572,  0.9619748 ,\n",
       "        0.99260295,  1.        ,  0.52523932,  0.9939295 ,  0.99557776,\n",
       "        0.9894934 ,  0.79671547,  0.73075086,  0.78691815,  1.        ,\n",
       "        0.98442829,  0.50845067,  0.54864295,  0.56270932,  0.66111684,\n",
       "        0.51732841,  0.62035284,  0.4633079 ,  0.52776121,  0.50285288,\n",
       "        0.78911551,  1.        ,  0.77622309,  0.81331657,  1.        ,\n",
       "        0.43488541,  0.48816976,  0.40181325,  0.38796986,  0.36513637,\n",
       "        0.5629852 ,  0.53723105,  0.60884071,  0.941744  ,  0.17758826,\n",
       "        1.        ,  0.82670183,  0.7792973 ,  1.        ,  0.83013155,\n",
       "        0.06945653,  0.23163494,  0.46629163,  0.51609293,  0.46235202,\n",
       "        0.60879421,  0.5894823 ,  0.67732176,  0.91178597,  0.84473057,\n",
       "        0.86712963,  0.84229942,  1.        ,  1.        ,  0.96407039,\n",
       "        0.59460724,  1.        ,  1.        ,  0.99723952,  0.93564962,\n",
       "        0.12139526,  0.273916  ,  1.        ,  1.        ,  0.98164858,\n",
       "        0.80080135,  0.69915668,  0.13621555,  0.33228821,  0.51096018,\n",
       "        0.49964611,  0.49665298,  0.65518214,  0.62679682,  0.68233381,\n",
       "        0.89383896,  0.62470824,  0.23048196,  0.42331155,  1.        ,\n",
       "        0.69563822,  0.794525  ,  0.68649458,  0.82732507,  0.76385348,\n",
       "        0.85207077,  1.        ,  0.89550569,  0.43040175,  0.57946445,\n",
       "        0.83229735,  0.99488546,  0.8160061 ,  0.98378634,  0.87217807,\n",
       "        0.95581388,  1.        ,  1.        ,  0.93143121,  0.94222351,\n",
       "        0.85519252,  0.75398807,  0.74132101,  0.68524077,  0.77859882,\n",
       "        0.64539027,  0.84710621,  0.85474693,  0.45494841,  0.83777608,\n",
       "        1.        ,  0.72506533,  0.9943554 ,  1.        ,  0.96791288,\n",
       "        0.82300733,  1.        ,  0.82623169,  0.7775087 ,  1.        ,\n",
       "        0.69513196,  1.        ,  0.73483582,  0.83641542,  1.        ,\n",
       "        0.25149369,  0.45094747,  0.50467077,  0.44533345,  0.6356088 ,\n",
       "        0.58389733,  0.66432989,  0.45051552,  0.48473151,  0.47771943,\n",
       "        0.59385131,  0.60490077,  0.64042403,  0.05676754,  0.0767261 ,\n",
       "        0.2181077 ,  0.7667143 ,  0.76339372,  0.08073411,  0.1686845 ,\n",
       "        0.78191313,  0.62577551,  0.2925989 ,  0.80794125,  0.63562642,\n",
       "        0.97555222,  0.65943323,  0.21069725])"
      ]
     },
     "execution_count": 21,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "dissimilarity"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "528"
      ]
     },
     "execution_count": 27,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "len(dissimilarity)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "plt.rcParams[\"figure.figsize\"] = (13, 5)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "threshold = 0.3\n",
    "linkage = hierarchy.linkage(dissimilarity, method=\"single\")\n",
    "clusters = hierarchy.fcluster(linkage, threshold, criterion=\"distance\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.text.Text at 0x7f2762b46780>"
      ]
     },
     "execution_count": 30,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAwYAAAFGCAYAAADQAYiRAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XmcJHV5+PHPszsbDhUyEwhGblFANMEgoIb96SgimCio\nMZFDJL7UHzFiws8LjyQMajyIEQ88QDniFcADkaCCRlpd0YBENOguh0ZkdxElO14gurM8vz++NWvv\n7MxuT3fV9PTM5/16zW53dfXTT1d/63iqvlUVmYkkSZKkxW1JvxOQJEmS1H8WBpIkSZIsDCRJkiRZ\nGEiSJEnCwkCSJEkSFgaSJEmSgKF+J9CpiPC6qpIkSVINMjOmDhuYwgDAey5IkiRJvYnYrCYA7Eok\nSZIkCQsDSZIkSVgYSJIkScLCQJIkSRIWBpIkSZKwMJAkSZKEhYEkSZIkLAwkSZIkYWEgSZIkCQsD\nSZIkSVgYSJIkSQKG+p1AE0ZGYHy831ksTMPDsG5dv7OQJElS3SIz+51DRyIiO801Agbkaw0cp60k\nSdJgiwgyM6YOtyuRJEmSJAsDSZIkSQv0HINBMMjnQcRmB57mP8+NkCRJ2jLPMeiTQchxIXF6S5Ik\nFZ5jIEmSJGlGFgaSJEmSLAwkSZIkzUFhEBFHRcSqiLg5Ik6bYZzRiPhmRNwYEVc3nZMkSZKkTTV6\n8nFELAFuBg4H1gLXAcdm5qq2cXYErgGenJlrImKnzLxrmliefKyuOb0lSZKKfp18fChwS2belpnr\ngYuAY6aMczzwicxcAzBdUSBJkiSpWU0XBrsCt7c9X10Na7cvMBIRV0fEdRFxYsM5SZIkSZpiPtzg\nbAg4CHgicD/gaxHxtcy8tb9paa41fdO3pm7M5s3TJEnSQtB0YbAG2KPt+W7VsHargbsy817g3oj4\nMnAgsFlhMDY2tvHx6Ogoo6OjNaerfhofb+Y8gKYLjvHx5u8GbfEhSZK61Wq1aLVaWx2v6ZOPlwI3\nUU4+vgO4FjguM1e2jbM/8C7gKGAb4D+BZ2fmd6fE8uTjBa6pabIQpvVC+A6SJGl+mOnk40aPGGTm\nhog4BbiKcj7DeZm5MiJOLi/nuZm5KiKuBL4NbADOnVoUSJIkSWpWo0cM6uQRg4XPIwYzWwjfQZIk\nzQ/9ulypJEmSpAEwH65KJA2kpk9qnqrpE5zBk5wlSVrMLAykLjV1FaV+moviQ5IkzU92JZIkSZLk\nEQNpvpnrLkpT9fOogV2ZJEnqHwsDaZ5ZiF2UOmVXJkmS+seuRJIkSZIsDCRJkiRZGEiSJEnCwkCS\nJEkSFgaSJEmSsDCQJEmShIWBJEmSJCwMJEmSJGFhIEmSJAkLA0mSJElYGEiSJEnCwkCSJEkSFgaS\nJEmSsDCQJEmShIWBJEmSJCwMJEmSJGFhIEmSJAkLA0mSJElYGEiSJEnCwkCSJEkSMNTvBCTNjZER\nGB/vdxZbF9HvDAQwPAzr1vU7C0nSXLIwkBaJ8XHIrD9uRDNx1V8WaJK0+NiVSJIkSZKFgSRJkiQL\nA0mSJElYGEiSJEnCwkCSJEkSc1AYRMRREbEqIm6OiNOmef3xEfHTiPiv6u/vm85JkiRJ0qYavVxp\nRCwBzgYOB9YC10XEZZm5asqoX87Mo5vMRZIkSdLMmj5icChwS2belpnrgYuAY6YZzytmS5IkSX3U\ndGGwK3B72/PV1bCpHhsRN0TEFRFxQMM5SZIkSZpiPtz5+Hpgj8y8JyKeAnwK2LfPOUmSJEmLStOF\nwRpgj7bnu1XDNsrMX7Y9/mxEvCciRjJz3dRgY2NjGx+Pjo4yOjpad76SJEnSgtJqtWi1WlsdLzKz\nsSQiYilwE+Xk4zuAa4HjMnNl2zi7ZOad1eNDgUsyc69pYmWnuUZAg1+rFoOQ41xrapoYdzDjqr/8\nXSVp4YoIMnOzc3wbPWKQmRsi4hTgKsr5DOdl5sqIOLm8nOcCz4qIFwHrgV8Bz24yJ0mSJEmba/SI\nQZ08YrDwDdoebeM2G1f95e8qSQvXTEcMvPOxJEmSJAsDSZIkSRYGkiRJkrAwkCRJkoSFgSRJkiQs\nDCRJkiRhYSBJkiQJCwNJkiRJWBhIkiRJwsJAkiRJEhYGkiRJkrAwkCRJkoSFgSRJkiQsDCRJkiRh\nYSBJkiQJCwNJkiRJWBhIkiRJwsJAkiRJEhYGkiRJkrAwkCRJkoSFgSRJkiQsDCRJkiRhYSBJkiQJ\nCwNJkiRJWBhIkiRJAob6ncAgGRmB8fH64kXUF2t4GNatqy+eJEmSFhcLg1kYH4fMfmcxvTqLDEmS\nJC0+diWSJEmSZGEwH4yMlD3+vfxBb+8fGenvNJAkSVJ/2ZVoHpgPXZTsiiRJkrS4ecRAkiRJkoWB\nJEmSJAsDSZIkScxBYRARR0XEqoi4OSJO28J4h0TE+oh4ZtM5SZIkSdpUo4VBRCwBzgaOBB4OHBcR\n+88w3puBK5vMR5IkSdL0mj5icChwS2belpnrgYuAY6YZ7yXAx4EfN5yPJEmSpGk0fbnSXYHb256v\nphQLG0XEg4CnZ+YTImKT1ySpTiMj5fLA6oyXMe7c8DCsW9fvLCSpN/PhPgZvB9rPPXBVJKkR8+Ge\nIQvdYi2+xscXZyFlQSQtLE0XBmuAPdqe71YNa3cwcFFEBLAT8JSIWJ+Zn54abGxsbOPj0dFRRkdH\n685XktSDxVR8LdYiqN1iLYjaWRxpELRaLVqt1lbHi2xwCR4RS4GbgMOBO4BrgeMyc+UM418AXJ6Z\nn5zmtew014hmVkyDFtccjDvIcZswSLkOqsU0jRfTd9XMbAcaRBFBZm5W1jd6xCAzN0TEKcBVlBOd\nz8vMlRFxcnk5z536libzkSRJkjS9Ro8Y1MkjBs1ayDkYdzDjNmGQch1Ui2kaL6bvqpnZDjSIZjpi\n4J2PJUmSJFkYSJIkSbIwkCRJkoSFgSRJkiQsDCRJkiQxi8IgIvaMiCdVj7eLiAc0l5YkSZKkudRR\nYRARLwQ+DpxTDdoN+FRTSUmSJEmaW50eMXgxcBjwc4DMvAX4/aaSkiRJkjS3Oi0Mfp2Zv5l8EhFD\neJdiSZIkacHotDD4UkS8BtguIo4APgZc3lxakiRJkuZSZAf38Y6IJcDzgScDAVwJfCA7eXNNIqLj\nj2vq9uSDFtccjDvIcZswSLkOqsU0jRfTd51vRkZgfLzfWcwvw8Owbl2/s9CgiAgyMzYb3mFhcD/g\n3szcUD1fCmyTmffUnunMOVgYNGgh52DcwYw7CCt+V8Sbmw/LkrmymL7rfOO035zTRLMxU2HQaVei\n/wC2a3u+HfCFOhKTpOmMj5eV3Hz+m++FiyRJs9FpYbBtZv5y8kn1ePtmUpIkSZI01zotDO6OiIMm\nn0TEo4BfNZOSJEmSpLk21OF4pwIfi4i1lJOPHwg8u7GsJEmSJM2pjk4+BoiIZcB+1dObMnN9Y1lN\n//mefNyghZyDcY3blEHIca4tpmmymL7rfOO035zTRLMx08nHnR4xADgE2Kt6z0FVwA/WlJ8kSZKk\nPuqoMIiIDwH7ADcAG6rBCVgYSJIkSQtAp0cMDgYOmMsbmkmSJEmaO51elehGygnHkiRJkhagTo8Y\n7AR8NyKuBX49OTAzj24kK0mSJElzqtPCYKzJJCRJkiT1V0eFQWZ+qelEJEmS+mlkBMbH+51F92Kz\ni08OjuFhWLeu31mo06sSPQZ4F/Aw4HeApcDdmblDg7lJkiTNmfFx7wXQL4Nc1CwknZ58fDZwHHAL\nsB3wAuDdTSUlSZIkaW51WhiQmbcCSzNzQ2ZeABzVXFqSJEmS5lKnJx/fExG/A9wQEWcCdzCLokKS\nJEnS/Nbpxv2J1binAHcDuwPPbCopSZIkSXOr08Lg6Zl5b2b+PDPPyMyXAk9tMjFJkiRJc6fTwuCk\naYb9VY15SJIkSeqjLZ5jEBHHAccDe0fEp9te2gHwarOSJEnSArG1k4+voZxovBPwL23DfwF8u6mk\nJEmSJM2tLXYlyszbMrMFPAn4SnUH5DuA3YCObkUREUdFxKqIuDkiTpvm9aMj4lsR8c2IuDYiDpv9\n15AkSZLUi8gObvEXEdcD/wcYBr4KXAf8JjNP2Mr7lgA3A4cDa6v3HZuZq9rG2T4z76ke/yFwSWY+\nbJpY2UmuZdxm7lw4aHHNwbjGbdYg5DjXFtM0WUzfdb5ZzMudhcppP7cigszcbCd/p/cxiMy8JyKe\nD7wnM8+MiBs6eN+hwC2ZeVuVxEXAMcDGwmCyKKjcH7ivw5wkSZIWnZERGB/vdxb1i476ogyO4WFY\nN2Bn5HZcGETEY4ETgOdXw5Z28L5dgdvbnq+mFAtTgz8deBOwM/BnHeYkSZK06IyPu3d9EAxiodPp\n5UpPBV4NXJqZ34mIBwNX15VEZn6q6j70dOANdcWVJEmS1JmOjhhUJx1/qe3594G/7eCta4A92p7v\nVg2b6XNWRMSDI2IkMzc7+DI2Nrbx8ejoKKOjox2kIEmSJC1erVaLVqu11fG2ePJxRLw9M0+NiMuB\nzUbMzKO3GDxiKXAT5eTjO4BrgeMyc2XbOPtk5veqxwcBl2Xm7tPE8uTjBi3kHIxr3KYMQo5zbTFN\nk8X0XeebxbzcgcHJc7Gbz79Ttycff6j6/63dfGhmboiIU4CrKN2WzsvMlRFxcnk5zwX+PCKeC/wG\n+BXwl918liRJkqTudXS5UoCI2BkgM3/SaEYzf75HDBq0kHMwrnGbMgg5zrXFNE0W03edbxbzcgcG\nJ8/Fbj7/TjMdMdjqyccRMRYRd1G6BN0cET+JiH9sIklJkiRJ/bHFwiAiXgocBhySmSOZOQw8Gjgs\nIv7fXCQoSZIkqXlbO/n4m8ARmXnXlOE7A1dl5h83nF/7Z9qVqEELOQfjGrcpg5DjXFtM02Qxfdf5\nZjEvd2Bw8lzs5vPv1G1XomVTiwLYeJ7BsrqSkyRJktRfWysMftPla5IkSZIGyNYuV3pgRPx8muEB\nbNtAPpIkSZL6YIuFQWYunatEJEmSJPXPVi9XKkmSJGnhszCQJEmSZGEgSZIkycJAkiRJEhYGkiRJ\nkrAwkCRJkoSFgSRJkiQsDCRJkiRhYSBJkiQJCwNJkiRJWBhIkiRJwsJAkiRJEhYGkiRJkrAwkCRJ\nkoSFgSRJkiQsDCRJkiRhYSBJkiQJGOp3ApIk1WlkxQrGJyaa/6CT9iJaP2j8Y4aHhli3fHnjnyNJ\nkZn9zqEjEZGd5hoBTXytQYs78pYRxu8d72zkq0+HJ5zR0ajD2w6z7rR1PWQ2vUGbvsYdzLh1GoQc\n59p8mCbRapGjo33NYc6Kk1laKEXGYl7uwODkudjN598pIsjMmDrcIwYL2Pi94+TpHbbI0wHGOho1\nztisHUmS2oxPTPS9OJlOtFr9TkHSPOY5BpIkSZIsDCRJkiRZGEiSJEnCwkCSJEkSFgaSJEmSsDCQ\nJEmSxBwUBhFxVESsioibI+K0aV4/PiK+Vf2tiIg/bDonSZIkSZtq9D4GEbEEOBs4HFgLXBcRl2Xm\nqrbRvg88LjN/FhFHAe8HHtNkXpIkSdKWrFgxwsREhzeKncZJJ51Oq9XZzWNnMjQ0zPLl9d9UdsbP\nazj+ocAtmXkbQERcBBwDbCwMMvPrbeN/Hdi14ZwkSZKkLZqYGGd0tPtbF5d7HI71lEOrNbc3lW26\nK9GuwO1tz1ez5Q3/FwCfbTQjSZIkSZtp+ohBxyLiCcDzgOX9zkWSJElabJouDNYAe7Q9360atomI\n+CPgXOCozJyxM9fY2NjGx6Ojo4yWYzSSJEmSZtBqtWi1Wlsdr+nC4DrgIRGxJ3AHcCxwXPsIEbEH\n8AngxMz83paCtRcGkiRJkrZu6g71M86Y/qToRguDzNwQEacAV1HOZzgvM1dGxMnl5TwX+AdgBHhP\nRASwPjMPbTIvSZIkSZtq/ByDzPwcsN+UYee0PX4h8MKm85AkSZI0M+98LEmSJMnCQJIkSZKFgSRJ\nkiQsDCRJkiRhYSBJkiQJCwNJkiRJWBhIkiRJwsJAkiRJEhYGkiRJkrAwkLTYjIxARC1/pzNWWyxG\nRvo9ZSRJi9xQvxOQpDk1Pg6ZtYQaa/u3ZxH1xJE076wYWcHE+ERt8U5iL1rxg9riAQwND7F83fJa\nY2rwWBhIkiQ1aGJ8gtEcrS1eibRXbfEAWtGqNZ4Gk12JJEmSJFkYSJIkSbIwkCRJkoSFgSRJkiQs\nDCRJkiThVYkkSQNgZMUKxic6v9xjtFodjTc8NMS65V6iUZLAwkCSNADGJybI0dHa43ZaQEjSYmBX\nIkmSJEkWBpIkSZIsDCRJkiRhYSBJkiQJTz6WJGlgzPbqTNPp9YRrr+QkLVwWBpIkDYimrs40G17J\nSVq47EokSZIkycJAkiRJkoWBJEmSJCwMJEmSJOHJx5IkSdKcWbFihImJ8Y7Hb7Wio/GGhoZZvnxd\nt2mVGD29W5IkSVLHJibGGR3N2uN2WkBsiV2JJEmSJFkYSJIkSZqDwiAijoqIVRFxc0ScNs3r+0XE\nNRFxb0S8tOl8JEmSJG2u0XMMImIJcDZwOLAWuC4iLsvMVW2j/S/wEuDpTeYiSZIkaWZNHzE4FLgl\nM2/LzPXARcAx7SNk5l2ZeT0w0XAukiRJkmbQdGGwK3B72/PV1TBJkiRJ84iXK5UkSRpAK0ZWMDFe\nX4eLVrRqiTM0PMTydctriaW51XRhsAbYo+35btWwroyNjW18PDo6yujoaLehJEmSBtrE+ASjOdrv\nNDZTV4Gh+rRaLVqt1lbHa7owuA54SETsCdwBHAsct4Xxt3hnhvbCQJIkSdLWTd2hfsYZZ0w7XqOF\nQWZuiIhTgKso5zOcl5krI+Lk8nKeGxG7AN8AHgDcFxF/BxyQmb9sMjdJkiRJv9X4OQaZ+TlgvynD\nzml7fCewe9N5SJIkSZqZdz6WJEmSZGEgSZIkycJAkiRJEhYGkiRJkrAwkCRJkoSFgSRJkiQsDCSp\nHiMjENH9H/T2/oiSgyRJXWr8PgaStCiMj0Nmf3OYLDAkqQcrRlYwMT7RU4xWtHp6/9DwEMvXLe8p\nhmbPwkCSJEkbTYxPMJqjfc2h18JC3bEwkCRpkRtZsYLxic73EEer1dF4w0NDrFvuXl9pUFgYSJK0\nyI1PTJCjo7XH7bSA0OIw2y5KnR41sNtRfSwMJEmS1LimuijZ7ag+XpVIkiRJkoWBJEmSJAsDSZrf\nZnN/BPCeB5KkrnmOgSTNZ03dH8F7HkiSpvCIgSRJkiQLA0mSJEkWBpIkSZKwMJAkSZKEhYEkSZIk\nLAwkSZIkYWEgSZIkCQsDSZIkSVgYSJIkScLCQJIkSRIWBpIkSZKwMJAkSZKEhYEkSZIkLAwkSZIk\nYWEgSZIkCQsDSZIkSVgYSJIkSWIOCoOIOCoiVkXEzRFx2gzjvDMibomIGyLikU3nJEmSJGlTjRYG\nEbEEOBs4Eng4cFxE7D9lnKcA+2TmQ4GTgfc1mZMkSZKkzTV9xOBQ4JbMvC0z1wMXAcdMGecY4IMA\nmfmfwI4RsUvDeUmSJElq03RhsCtwe9vz1dWwLY2zZppxJEmSJDXIk48lSZIkMdRw/DXAHm3Pd6uG\nTR1n962MA0BEdPzBsxh1VgYu7lgzgRuLO2jT17gDGXfgEh6wuI2l20xY4w5q3MGaLQZwAhu32cDz\nM27ThcF1wEMiYk/gDuBY4Lgp43waeDFwcUQ8BvhpZt45NVBmNvaTS5IkSYtdo4VBZm6IiFOAqyjd\nls7LzJURcXJ5Oc/NzM9ExJ9GxK3A3cDzmsxJkiRJ0uYiM/udgyRJkqQ+8+TjeaC638MgxW2kW9eg\nxdVgGrT5rSnOF8WgTYdBy7dug/b9By3fSYOWdxP5NjkN5vP0HagV2dZExI79zmE2JhtGZt4XEUvr\njp+Z99Uds9JIg86aD1+1Td+sc6MtIkbqijVN7FrnyYhYEhEPrjNmW+xG5rcoHtHEPAHNzRd1x62m\nw3OamM4RsQzYpqZYSyLieRHxO3XEmxL30XXGnE4Ty52G228T+dbezqq4tcWs4j0IWFZXzGni153v\ng6v1T63rzIjYoc54M33MHHxGbeqeL5qKCfUuf5vQ9MnHc6Ka6S4AhiLi88Cq6mZpdcReApwEfDQz\nf11HzDbviIj9M/PJ1fkYS3rduKjy/VfgTuBG4OuZuarXRKu4FwC/BL4cEZ/OzF/VFPcc4OfANcAX\nM3O817iV9ul7X0RELzN6lesHgB0i4pOZ+dE6kqzivj4zX1tHnlPiXgx8DXhbr/GmxH0vsH1E/Dtw\naWb+pqbYAXwe+Ehm3lhHzCpuk/PFOcCPq7jXZ+bNNcX9GLAX8EXgZ3W0i7a4a4EfRsRHMnNtDfEO\nB1ZR2lrPqnbwGcr5abUsy6fEb6o9NNl+m2pnHwf2pP52dmn1d2FNeX4IWAp8JSIuy8zVvcadEr+2\nfCuvAJ4TEcdm5ndrnK7vAZZV2zq3ZOb1dSTbFr+p9Xzt21FNbD80uU3SNr+toYblb1vMWqftQjli\n8FzKicsvARI4MiKeUlPszwNnUWbw7WuKSbU3aSWwe0S8NyKW1nTk4HmURvdeYAPw0oj44x5jAlxO\nuUHdRcA+1LfX5pPAj4DPAg8HhusIWs0sK4Hd2qZvr0cOLqCsmF8FHBYRD60jV8qK6OURcTZsPMJR\nx96aTwC/yMy3RcQzIuKwiLhfDXHPB/4XeAtwBPAnNcSc9CfAFZl5QUS8MCKeGBH71hC3qfni74G7\nKL/hdsCpEbFfDXE/BtxAuSv8WVDb3qszgf8BXg3cR+/z8QcoG+5HA6+NcgW6OhwKfLlquy+PiKdF\nxEE1xYayIl1N2ciqsz0sp5n2+1rgJ9Tfzi4GvkkpkupsZ1+lFC8XRrm4yP4RsXMP8V5H2VB7DvBA\n4A9qyLHdV6g3X4DvA+uAd0fEo2s6cv1KYD3wGuABwAkRUefy9zLKfPFv1Luev5LSvk6oczuKspH9\nI8pOhLq2Hy6lgW2Sypsp7aKu5S80sI26UAqD/6Ic9rqP0lC+DvxRryupaiP95cDDgD8CnhsR9+8x\nVyJix8zcQFn5vw24h7IBAPC7PcTdljJT7wz8IDM/BFwBnBgRu/UQ9wHAP1d7tL8CPB54Q0T8bUTs\n3UPc36dcqeofMvMLwG+A/9ttvHbVkZdPUGaY9unb1UxeLdDPzsxXZeatlGl8RkScUn2PXrwhM7cB\n7omIc6v861g5vx/YLiJuoWxo/Q1wUkR0vaCrDoF+JDNfU+0R/RTwzBpynXQfcEhETO7FPJBSND28\n24BV+/0h8PvUO18sA66ltK9bKSuSHSkL6F7m430o88XrM/MdwJqIeGq38driLqXsAZvIzF8Cj6Ns\nzL8mIh7ZRbztKW3hzMz8MvAfwEOq13otbFdTLnX9ScrK8wDgaRFxWI9xJ/2MsoF5W43tYQdgHDg4\nIj5GTe238kNKO/seZSOojna2HXBJ1c7eCayuo51VLgIeFxEt4BjKcv0FEbFrl/FWUAr77Sn3OnpN\nRJwaEYfXkSzl9x+NiC9RT74A/04pxN8EvCUiHsWm93XqxjXADygFx0cp02J5HdsllfdV6/kV1LSe\nr5xB2cg+kJq2oyqXVNsP/0EN2w9Ruqqd28Q2SeV7wPo6lr8AUbpvngbsT43bqAulMLgFuInSkIco\nM8/9KBXvrEXp2/p+yt6Z/TPzDkql90jKHgsi4sHVBsds455D2YNwPHB/YO/MfBlwY0R8G3h6l/me\nC5xLWWHcBDy72jt8BaXbwAN7iPse4IERsW1EHELZc3U5ZYU96/7rbdPhbcDObXuxr6YcwiQi9uly\n+v5T26AlwF69Tl/YWGh8o/qcAyldEN4G7AIc3E3MNt+rPuOVtBUHEbFHtfLuSmZ+hrIn8B2Z+WrK\nXop9gEf1EHM90IpiiLIBN1zlu383+Va/2xur+F8Drge2y8y/z8yzKPPzrPe6VnEvoBzh2KaKW9d8\n8X7gPGAnYIKycH4YZQdFUpY/3cR9H2Uler+I2K6axjfQxfefJt8LKMvHn0TEpZSNrDdRNjQeM9u4\nmXlPtUKetI5q+VhDYbuWMr9toOT9L5QjHb1suLe3s08C3wX+oipw6lhOnl3ldyNwv5ra75uqp1dS\nNlJeQT3t7A2Z+avM/Fg1Ly+lnnb2gYj4IPAtyvL8C5l5MmXaQPd7+q8FvkPZ2bEPpXfAWrpY/0zJ\n980AmflGqq5rNeULsC3wF5Tp+g+UHVW9Xo59JeWI0YuAp1C2IfalFKFdmTJfXFENexT1rOffVMVd\nkZlrKMubA+l9O+qfqrgXxW97WXyRLrcfqpjnA++kTNNJX6T0Rulqm6Qt9uR8/G/Aj6odHl0vf6uY\nFwIfoWyj/ohyBP9A4MRqnFlP20kLojDIzHspey4fCTyVsvL7MWUPQHSx9+p8ykzxFuBJEXF4VRyM\nAXtHxEcoK+/ZTr/JLhhvBp5IWah/LSJOoCyAgmojc5Y5n0/pznAmcCRlI2gZcHxmTlAWJMtnmWt7\nvm8BngQcmpnXZeZzMvPzVdxuTg48n7IR8SZK15HJDdUfAfdGxHHA6yn9SWfjQjbtkvMj4PtVEdbL\n9KWKN7mxsxL4u8z8BmUP4eHRwyHi6vySyROlTwXWRsR/UgqybbuNW8X7HPDu6vEPKe2kp0ImM9dn\nMQH8N7AyIp4BvJ3SJmbrQuBlEfHu6vlngLsi4lXV8xGgm24kF1Da1GuBP62GbQucUOW+Pd3PFz+h\nzG9PBO6lFEjPpGwA3AV00y3lfEp7ejNlPn5MleeVwNHVcqIb7cuHJwMtym91bWb+D6VLQq97tMnM\nf6XMv6+pIVZSlulrKBtXI5Tl2hERsbSb+ZfftrPJjb6rgd35bXvoZTl5F/DWKtcNwHdqaL8XUro4\nvaNaln2cUoT+Ob21swuBV7QtJ7M6gn0VvbezH1OKuBOAKzPzDdVnfJ8yHWZdgFbvX5eZl1KOAH++\nOr9gR+CoiFjWQ3v4fxHx3uoz3kqZ93rOt4rxU0rXxVHgGcC3gYN6aL9k5k8oOyTuoKx7n0fZeO2l\nu9qFbDpfANxUw3r+Qkr73Ri3asdj9LYddSFlPf+uKuaGangv2w/t3YQfF7/tJrwauLvahuhmm2Qy\n35dGxDuY+g1lAAAH60lEQVQz8+fVUeCL6G35O7luezXwJ1HOpVxLWdft0cO0LTJzwfwBu1L6jn4I\n+Bzw0C5iLAOOaHv+VOCdbc+fR9kwfFiPcY8G3gE8tPqRx6rhwz3GfTJlb8eBlL2Yl1A2LPateToc\nS1kgzWoabykupbvAasqGy35d/Hb7Vv+fSemKAeXQ7YXdTt+tfN6zKH1pZ93OZog3eV+RI4DbZtvG\nOoh/LGVP3j515QvsQFlJfaOb32zK7/bPlGIoKCu6l1P2iFwBHDDLmEuAQ9qef5iyt+ZFlI2Liykb\nQr3OF0dOmS+OBL4M7N5j3Knz2+MoC/pte4w7udzZj9Jv+5xqOsxq+m6h7R4MvGy2eW4h7oMoexfP\npxQK+/cQq72dnV21s4dRNgbqWk4eCbwLeEQ1HT7cTfudku9bgfe0T+tu29mUuGdSuk20v1ZXO3va\nlPZ7IqUf/4N7bA+7UDbQzgK+VFN7OBP4wJTXnltjvh8EzqieD/USb5r4T6Qs0/eqaTq8f8prXa3n\nO2hnXW1HbSkuXW4/sPm64hJKN62/q9rxbbONubV2BuxdteNZL3+3ku/9q2l7UzfTdmPMOhvpfPmj\nVPo79fD+ZdXCd4hyFOJD1fB9KHuEutoQnCbuBdXw51AOPQMs6THuQcD51fDdKZX+79U8HfYF/no2\njbnDuL9LOanq4V3GXdr2+O1UK1PKntyup+8Mn7UNpe/hrDYkOoy9N/CQmmNuCzyfHlakW4j9yl6m\nw5Tf7R3Au6vHv0fZ0BzpMu7kxuqBlKMmh1A2ho+knG+wc5dxp7bfD1fD9wFe0EP7nWm+2J+yIb9j\nTXEvqIY/AfhLetiomOazduz299pK3OFul2MdtrO6lpPty989qjbXbfudujyb3KjYpZqX61pOnls9\n3r3mdta+vvgb4BE1tYUdqjz3rLE9tE+H3WrOd6fpPrOm2I/uZdm7lenwB9XyrNv1/Jbi9rIdNVPc\n7YDTu5kv2HRd8S7Keud1lO5a/9jtvDZDvpPbJYdQjv7tVVO+b6DsVDy41zbhnY+3ouq/dhrlUnwv\nBI7NzJ/VFPfVlGr/xcBfZubPa4p7GqVinMz3pzXGbWI6nAbcTFnZPbuX6RDx20vCRcTrKJdRHKd0\nq+p5+k75rJ4vLzuXmso3yhWfNmx9zC3GaP/dxoA/oxwqPbHX9hvlBK2JLFf9ehnwB5n58l5itsWe\nOl/01H6nibsSOBk4Lmu4ZF4V91WUfvB/TVnu/KLXuINiSjs7g9LFrJZ2VsWcuvzt6XebZr6YzPeE\nXn63GeLeSf3TYXK+mJftbAvT9/gsJ4g28lnzTR/a2Qk1rufHKEdY1wLP6THu1HXFg4DXZumu3rVp\ntkuOoHQFPK6XdtbYuq2XqmKh/7F5V4la9rYat/m41f+NdMnxr5m/pn83Guj6NWjzRRNxB+2vqXY2\nh8uznrp9Dep0mMP2UMv0HbS/OWxn8zpuFfNZlIsG1Hbkfg7yrW3d5hGDDkTEK4FPZQ03lDHunMbd\nm3IY79Y646pZTfxuEbEN5fyj1gC134GKO2iaWj4M2vJs0KZDU1xfFAPYzmqN2+S6ooo/7/O1MOhA\nHV0ljDv3caV287kr1UKIq8LpWzgdNKgWezdhCwNJkiRJC+M+BpIkSZJ6Y2EgSZIkycJAkiRJkoWB\nJEmSJCwMJEmSJGFhIEkLQkRsiIj/iohvVv+/sub4j4+Ix9YZs4PP3DMi/nsuP1OSFrOhficgSarF\n3Zl5UIPxR4FfAl9r8DOm0/U1tb2WviTNjkcMJGlhiM0GRBwZEZe0PX98RFxePX5yRFwTEd+IiIsj\nYvtq+P9ExFhEXB8R34qIfSNiT+CvgVOroxGHTfmc0yPivIi4OiJujYiXVMM32eMfES+LiH+sHl8d\nEW+LiOsi4jsRcXBEfCIiboqI17eFXxYRH46I70bEJRGxbfX+gyKiVb3/sxGxS1vcsyLiWuBv65m0\nkrQ4WBhI0sKw3ZSuRH8BfAE4NCK2q8Z5NvDRiPg94LXA4Zl5MHA98NK2WD/OzEcB7wNenpm3VY/P\nysyDMvOr03z+fsARwKOB0yNiaTV8S3v8f52ZhwDnAJcBLwL+EPiriBhui3t2Zh4A/AL4m4gYAt4F\n/Hn1/guAN7bFXZaZh2bmWVucYpKkTdiVSJIWhnum60oUEZ8DnhYRnwD+DHgFpVvQAcBXIyKAZcA1\nbW+7tPr/euAZHX7+FZk5AfxvRNwJ7NLBez5d/f/fwI2Z+eMq5+8BuwM/A36YmV+vxvsw8BLgSuAR\nwOer/JcAa9viXtxhzpKkNhYGkrSwXQycAowD12Xm3dXG9FWZecIM7/l19f8GOl9P/Lrt8X3V+yaA\npW3Dt53hPfdNeX9u4XOT0m3qxsw8bIZx7u4kYUnSpuxKJEkLw2bnGFS+BBwEvBC4qBr2deCwiNgH\nICK2j4iHbiX+L4AdZpnTncDOETEcEdsAT53l+wH2iIhHV4+PB74C3FTFfQxARAxFxAFdxJYktbEw\nkKSFYdsp5xi8ESAz7wP+HTiq+p/MvAv4K+DfIuJblG5E+1VxZjon4HLgGdOdfDyNrD5nAngdcB2l\n+8/KqeNs6f2VVcCLI+K7wO8C78vM9cCzgLdExA3AN4HHTvNeSdIsRKbLUEmSJGmx84iBJEmSJAsD\nSZIkSRYGkiRJkrAwkCRJkoSFgSRJkiQsDCRJkiRhYSBJkiQJCwNJkiRJwP8Hou8Uxxr1n1sAAAAA\nSUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f27400d2d30>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "hierarchy.dendrogram(linkage, color_threshold=0.3)\n",
    "plt.xlabel(\"Event number\")\n",
    "plt.ylabel(\"Distance\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "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.5.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 1
}