{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# PCA Analysis" ] }, { "cell_type": "code", "execution_count": 212, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "from sklearn.decomposition import PCA\n", "from scipy.linalg import svd\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "code", "execution_count": 213, "metadata": {}, "outputs": [], "source": [ "country_correlation = np.load('Data/country_correlation.npy')\n", "country_names = np.load('Data/country_names.npy')\n", "data = np.load('Data/country_capability_dict.npy').item()\n", "number_of_features = 50000\n", "\n", "X = []\n", "\n", "for country in data: \n", " X.append(data[country])" ] }, { "cell_type": "code", "execution_count": 214, "metadata": {}, "outputs": [], "source": [ "X = np.array(X)\n", "X = X[:, 0:number_of_features]\n", "N = X.shape[0]\n", "M = X.shape[1]" ] }, { "cell_type": "code", "execution_count": 215, "metadata": { "collapsed": true }, "outputs": [], "source": [ "Y = X - np.ones((N,1))*X.mean(0)" ] }, { "cell_type": "code", "execution_count": 216, "metadata": {}, "outputs": [], "source": [ "U,S,V = svd(Y,full_matrices=False)\n", "\n", "rho = (S*S) / (S*S).sum() " ] }, { "cell_type": "code", "execution_count": 217, "metadata": { "scrolled": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA4IAAAG5CAYAAADbK4omAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAIABJREFUeJzs3XmYXGWZ9/HvnQSysIYQhG6ygCAj\ngqK0wYQlKIMiIriggK2CMgZ3HcZBmKgwjHHEfUN9UVDEuCCKxhHFGRFERSCg6ADqxJiVLQkhIE2A\nJPf7x6k2lU5Vpzrp01Xd/f1cV11V5zmnqu6q7pD8eLbITCRJkiRJw8eIZhcgSZIkSRpYBkFJkiRJ\nGmYMgpIkSZI0zBgEJUmSJGmYMQhKkiRJ0jBjEJQkSZKkYcYgKEktIiImR8TfImJks2spW0RcHxH/\n1OC1d0bE0SXUcHRELOvlfEbEfv39vo2IiB9HxOnb+BpHRsSf+qmehn9ekqTBwSAoSVshIn4SERfW\naD8pIu6LiFF9fc3MXJKZO2bm+v6pcmjIzGdk5vXNrmMgZeaLM/PybXyNGzPzgP6qSZuLiAsi4uvN\nrkOStoZBUJK2zuXAayMierS/Dpibmev68mJbExw19ETBv5slSaXzLxtJ2jrfByYAR3Y3RMR44ATg\na5Xjl0TEbyPi4YhYGhEXVF07tTL08MyIWAJcV9U2qnLNGyLi7oh4JCIWRsRZVc8/OiKWRcS/RMQD\nEXFvRLyh6vzYiPh4RCyOiDUR8cuIGFs597yI+HVEPBQRd/Q27DIi2iLiuxGxIiL+GhHvrDp3TUR8\nvOr4WxFxWeXxGRHxq4j4XOX9/xgRx9R5j6dGxHURsSoiVkbE3IjYter8ooj4x8rjCyLiyoj4WuV7\nuTMiOhqsd2xEfDUiVkfEXcBz633uKsdXvvuVEfHRiBgREdtHxIMRcXDVa+8REV0RMbHG5+v1u6gM\nu5wTEb8CuoB9q4diVp7/y4j4WKX2v0bEi6uev1tEfCUi7qmc/36lfZOhr5Xv8byIuKty3VciYkzl\n3PiI+K/K97a68njvBr4fImJkRPxbRPyl8jO5LSImVc7NiIhbK5/71oiY0eNzf7Dyu/i3iPhhREyo\n/Pwfrlw/ter6jIh39vx5VM6NiIj3VX7fH6j8fuxSOdf95+r0iFhSee7sqtcdERHnVupfVfn92m1L\nz42I44B/A06p1H9H1c9rYeW7+GtEdDbyPUrSQDMIStJWyMzHgCuB11c1vxr4Y2beUTl+tHJ+V+Al\nwFsi4mU9Xmom8HTgRTXe5gGKYLkz8AbgkxHxnKrzewK7AO3AmcDFUYRRgI8BhwIzgN2Ac4ANEdEO\n/Aj4YKX9PcB36wSYEcAPgTsq73EM8O6I6K71jcDrIuIFlX/sTgPeVfUShwF/AXYHzge+1/0P7J5v\nBfwn0Fb5LiYBF9S4rtuJwLcovtd5wOcarPd84KmV24uARubgvRzoAJ4DnAS8MTOfqLz/a6uuOw34\nWWauqPM6W/ouXgfMAnYCFtd5/p8qz/8IcGnE33ujrwDGAc8A9gA+2cvn6aT47E8Fnga8r9I+AvgK\nMAWYDDxG5XttwNkUn/94it/VNwJdlc/3I+AzFP/T5BPAjyJiQtVzT6X47O2Vmm6q1LEbcDfFd1Vt\ns59Hpf2Myu35wL7AjjXqPwI4gOL34gMR8fRK+zuAl1H8WWwDVgMXb+m5mfkT4EPAtytDup8VETtU\nPu+LM3Mnij9/v6v/1UlSE2WmN2/evHnbihvFPw4fAsZUjn8F/HMv138K+GTl8VQggX2rzne3jarz\n/O8D76o8PpriH+ujqs4/ADyP4h/1jwHPqvEa7wWu6NF2LXB6jWsPA5b0aDsP+ErV8SuBpcBK4Iiq\n9jOAe4CoarsFeF3l8fXAP9X5nC8Dflt1vAj4x8rjC4D/qTp3IPBYI/UCC4Hjqs7NApb18vPKHte/\nlSLs/f29uj8fMB94dZ3XaeS7uLDHc/7+/VSev6Dq3LhKbXsCewEbgPE13vfo6s9X+R7fXHV8PPCX\nOjUfAqyuVU+Na/8EnFSj/XXALT3abgLOqHrN2VXnPg78uOr4pcDvGvx5/Ax4a9W5A4AngVFs/HO1\nd4/v/9TK47uBY6rO7dWH514AfL3q3A4U/014JTC23u+WN2/evLXCzTkpkrSVMvOXEbESeFlE3ErR\nI/aK7vMRcRjwYeAgYHtgNPCdHi+ztN7rV4b/nU/RczOCIgD8oeqSVbnpXMQuip6Q3YExFD1QPU0B\nXhURL61q2w74eZ1r2yLioaq2kcCNVcc/BD4L/Ckzf9nj+cszM6uOF1P0uGwiIp4CfJpimO1OFJ91\ndY16ut1X9bgLGBPFcNot1dvGpt93rZ63nnpe3waQmTdHRBdwdETcC+xH0TtZz5a+i7q/BxV//8yZ\n2VXpDNyRoufswczs7fuqVvPzRMQ4ip7E44DuXuWdImJkbnnxoknU/l1rY/PveDFF71+3+6seP1bj\neMdG6q/xXospgtxTqtp6/t50v/YU4OqI2FB1fn2Dz91EZj4aEadQ9LRfWhnu+y+Z+cda10tSMzk0\nVJK2zdcohn++Frg2M6v/IfsNinAwKTN3Ab5IMQyyWlJDRIwGvksxxPMpmbkrcE2N59eyElhLMdSu\np6UUPYK7Vt12yMwP17n2rz2u3Skzj6+6Zg5Fj8peEXFaj+e3Vw1fhGLI4T013udDFN/DwZm5M8V3\n2cjn7Gu991KElup6tqTn9dX1X16p9XXAVZm5tpfX2dJ3UfP3oAFLgd2iak7lFtT7PP9C0Yt2WOVn\ncFSlvZGfw1Jq/67dQxGyqk0GljdYay316u/5XpOBdWwaLOtZSjGUs/r3ZkxmNlLnZj+3zLw2M4+l\n6Fn8I/ClBl5HkgacQVCSts3XgH8E3kQRDKrtRNFbszYipgGv6cPrdvcgrgDWVXoHX9jIEzNzA3AZ\n8IkoFk8ZGRHTK+Hy68BLI+JFlfYxlUVFai0McgvwSES8N4qFVkZGxEER8VyAiDiKYu7i6ynm2322\nMgex2x7AOyNiu4h4FcX8v2tqvM9OwN+ANZXn/2sjn7Ov9VLM6TwvioVR9qaYG7Yl/1q5fhLF/Mdv\nV537OsWctddSWSCoF41+F32SmfcCPwY+X6lzu8rPpZ63RcTelfl7s9n4eXai6IF7qHKu59y83nwZ\n+I+I2D8Kz6zMA7wGeFpEvCYiRlV6yg4E/quPH7NavZ/HN4F/joh9ImJHNs7da2T13i8CcyJiCkBE\nTIyIkxqs535gamxctOYpUWwhswPwOMXv9YbeXkCSmsUgKEnbIDMXAb+mmBvUc2jgW4ELI+IR4AMU\nQaTR130EeGflOaspQmRvQw97eg/FMNJbgQeBi4ARmbmUYpGNf6MImUspgtdmfx9UhgSeQDFf7K8U\nPY1fBnaJiJ0pws/bM3N5Zt4IXAp8parn62Zg/8rz5gAnZ+aqGrX+O8XiH2soFhf5Xh8+Z0P1Vr3P\n4sq5n1IssrIlPwBuo1jw40cUn7H7/ZYCt1P0Ct1Y89kbNfpdbI3XUcxp+yPFPNF393LtNyg++0KK\n4ZwfrLR/Chhbqe83wE/68P6foPg9/SnwMMV3NLby+U6g6G1cRbFg0QmZubIPr91TvZ/HZRQ/z19Q\n/HzX0ljQh2JY8jzgp5U/q7+hmAPaiO6h3qsi4naKP0dnU/RQPkixAM1bGnwtSRpQ3ZPcJUnqNxFx\nBsXiIkc0u5YyRbFdxj2Z+b5erjmDFvguImJRpY7/aWYdWysiEtg/Mxc0uxZJGgpcLEaSpK0QxR53\nrwCe3dxKJEnqO4eGSpLURxHxH8D/Ah/NzL82ux5JkvrKoaGSJEmSNMzYIyhJkiRJw0ypcwQj4jiK\n1bhGAl/uuU9VZSnzrwGHUqwodkpmLoqI7ShWentOpcavZeZ/9vZeu+++e06dOrX/P4QkSZIkDQK3\n3Xbbysyc2Mi1pQXBiBgJXAwcCywDbo2IeZl5V9VlZwKrM3O/iDiVYnnzU4BXAaMz8+CIGAfcFRHf\nrCzTXtPUqVOZP39+WR9HkiRJklpaRCxu9Noyh4ZOAxZk5sLMfAL4FsXeVdVOYuMGzFcBx1T2n0pg\nh4gYRbGv0RMUexNJkiRJkrZRmUGwnWKj4m7LKm01r8nMdRSbCU+gCIWPAvcCS4CPZeaDPd8gImZF\nxPyImL9ixYr+/wSSJEmSNAS16mIx04D1QBuwD/AvEbFvz4sy85LM7MjMjokTGxoKK0mSJEnDXplB\ncDkwqep470pbzWsqw0B3oVg05jXATzLzycx8APgV0FFirZIkSZI0bJQZBG8F9o+IfSJie+BUYF6P\na+YBp1cenwxcl8XGhkuAFwBExA7A84A/llirJEmSJA0bpQXBypy/twPXAncDV2bmnRFxYUScWLns\nUmBCRCwAzgbOrbRfDOwYEXdSBMqvZObvy6pVkiRJkoaTKDrgBr+Ojo50+whJkiRJw1VE3JaZDU2p\na9XFYiRJkiRJJTEISpIkSdIwYxCUJEmSpGHGIChJkiRJw4xBUJIkSZKGGYOgJEmSJA0zBsGyzJ0L\nU6fCiBHF/dy5za5IkiRJkgAY1ewChqS5c2HWLOjqKo4XLy6OATo7m1eXJEmSJGGPYDlmz94YArt1\ndRXtkiRJktRkBsEyLFnSt3ZJkiRJGkAGwTJMnty3dkmSJEkaQAbBMsyZA+PGbdo2blzRLkmSJElN\nZhAsQ2cnXHIJTJmyse2ii1woRpIkSVJLMAiWpbMTFi2CW28tjidMaGo5kiRJktTNIFi2Qw6BnXaC\nX/yi2ZVIkiRJEmAQLN+oUXD44XDDDc2uRJIkSZIAg+DAmDkT7r4bHnig2ZVIkiRJkkFwQMycWdzf\neGNz65AkSZIkDIID49BDYexY5wlKkiRJagkGwYGw/fYwY4bzBCVJkiS1BIPgQDnqKPj972H16mZX\nIkmSJGmYMwgOlJkzIRN++ctmVyJJkiRpmDMIDpRp04ohos4TlCRJktRkBsGBMnYsHHaY8wQlSZIk\nNZ1BcCAddRTcfjs88kizK5EkSZI0jBkEB9LMmbB+Pfz6182uRJIkSdIwZhAcSNOnw8iRDg+VJEmS\n1FQGwYG0447Q0eGCMZIkSZKayiA40I46Cm65Bbq6ml2JJEmSpGHKIDjQZs6EJ5+Em29udiWSJEmS\nhimD4EA7/HCIcJ6gJEmSpKYxCA60XXeFQw5xnqAkSZKkpjEINsNRR8FNN8Hjjze7EkmSJEnDkEGw\nGWbOhLVrYf78ZlciSZIkaRgyCDbDkUcW984TlCRJktQEBsFm2H13eMYzDIKSJEmSmqLUIBgRx0XE\nnyJiQUScW+P86Ij4duX8zRExtdLeGRG/q7ptiIhDyqx1wB11FPzqV7BuXbMrkSRJkjTMlBYEI2Ik\ncDHwYuBA4LSIOLDHZWcCqzNzP+CTwEUAmTk3Mw/JzEOA1wF/zczflVVrU8ycCY8+Crff3uxKJEmS\nJA0zZfYITgMWZObCzHwC+BZwUo9rTgIurzy+CjgmIqLHNadVnju0HHVUce82EpIkSZIGWJlBsB1Y\nWnW8rNJW85rMXAesASb0uOYU4Ju13iAiZkXE/IiYv2LFin4pesDstRfsv7/zBCVJkiQNuJZeLCYi\nDgO6MvN/a53PzEsysyMzOyZOnDjA1fWDmTPhxhth/fpmVyJJkiRpGCkzCC4HJlUd711pq3lNRIwC\ndgFWVZ0/lTq9gUPCUUfBmjXwhz80uxJJkiRJw0iZQfBWYP+I2CcitqcIdfN6XDMPOL3y+GTgusxM\ngIgYAbyaoTg/sNvMmcW98wQlSZIkDaDSgmBlzt/bgWuBu4ErM/POiLgwIk6sXHYpMCEiFgBnA9Vb\nTBwFLM3MhWXV2HSTJ8OUKc4TlCRJkjSgRpX54pl5DXBNj7YPVD1eC7yqznOvB55XZn0tYeZMuOYa\nyITNFkyVJEmSpP7X0ovFDAtHHQUrV8Lddze7EkmSJEnDhEGw2brnCTo8VJIkSdIAMQg221OfCrvs\nAu95D4wYAVOnwty5za5KkiRJ0hBW6hxBNeAb34C//W3jXoKLF8OsWcXjzs7m1SVJkiRpyLJHsNlm\nz958Q/murqJdkiRJkkpgEGy2JUv61i5JkiRJ28gg2GyTJ/etXZIkSZK2kUGw2ebMgXHjNm0bN65o\nlyRJkqQSGASbrbMTLrkERo4sjqdMKY5dKEaSJElSSVw1tBV0dsJnPlNsI/HTnza7GkmSJElDnD2C\nraK9HZYvb3YVkiRJkoYBg2CraGszCEqSJEkaEAbBVtHeDmvWwKOPNrsSSZIkSUOcQbBVtLcX9/fc\n09w6JEmSJA15BsFW0dZW3Ds8VJIkSVLJDIKtwh5BSZIkSQPEINgquoOgPYKSJEmSSmYQbBU77ww7\n7mgQlCRJklQ6g2AraWtzaKgkSZKk0hkEW4mbykuSJEkaAAbBVmIQlCRJkjQADIKtpHtoaGazK5Ek\nSZI0hBkEW0l7Ozz5JKxc2exKJEmSJA1hBsFW4hYSkiRJkgaAQbCVtLUV964cKkmSJKlEBsFWYo+g\nJEmSpAFgEGwle+1V3BsEJUmSJJXIINhKttsO9tjDIChJkiSpVAbBVtPe7hxBSZIkSaUyCLYaN5WX\nJEmSVDKDYKtpazMISpIkSSqVQbDVtLcXG8o//nizK5EkSZI0RBkEW033FhL33tvcOiRJkiQNWQbB\nVuNegpIkSZJKZhBsNW1txb0rh0qSJEkqiUGw1dgjKEmSJKlkpQbBiDguIv4UEQsi4twa50dHxLcr\n52+OiKlV554ZETdFxJ0R8YeIGFNmrS1jt91g9GiDoCRJkqTSlBYEI2IkcDHwYuBA4LSIOLDHZWcC\nqzNzP+CTwEWV544Cvg68OTOfARwNPFlWrS0lohge6tBQSZIkSSUps0dwGrAgMxdm5hPAt4CTelxz\nEnB55fFVwDEREcALgd9n5h0AmbkqM9eXWGtrcVN5SZIkSSUqMwi2A0urjpdV2mpek5nrgDXABOBp\nQEbEtRFxe0ScU+sNImJWRMyPiPkrVqzo9w/QNAZBSZIkSSVq1cViRgFHAJ2V+5dHxDE9L8rMSzKz\nIzM7Jk6cONA1lqd7aGhmsyuRJEmSNASVGQSXA5OqjveutNW8pjIvcBdgFUXv4S8yc2VmdgHXAM8p\nsdbW0t4OXV2wZk2zK5EkSZI0BJUZBG8F9o+IfSJie+BUYF6Pa+YBp1cenwxcl5kJXAscHBHjKgFx\nJnBXibW2FreQkCRJklSi0oJgZc7f2ylC3d3AlZl5Z0RcGBEnVi67FJgQEQuAs4FzK89dDXyCIkz+\nDrg9M39UVq0tx03lJUmSJJVoVJkvnpnXUAzrrG77QNXjtcCr6jz36xRbSAw/9ghKkiRJKlGrLhYz\nvHX3CBoEJUmSJJXAINiKxo6F8eMNgpIkSZJKYRBsVe3tzhGUJEmSVAqDYKtyU3lJkiRJJTEItiqD\noCRJkqSSGARbVVsb3H8/rFvX7EokSZIkDTEGwVbV3g4bNhRhUJIkSZL6kUGwVbmXoCRJkqSSGARb\nVfdegq4cKkmSJKmfGQRblT2CkiRJkkpiEGxVe+wBI0caBCVJkiT1O4NgqxoxAvbay6GhkiRJkvqd\nQbCVuZegJEmSpBIYBFuZQVCSJElSCQyCraytzaGhkiRJkvqdQbCVtbfDmjXw6KPNrkSSJEnSEGIQ\nbGVuISFJkiSpBAbBVtYdBB0eKkmSJKkfGQRbWVtbcW+PoCRJkqR+ZBBsZQ4NlSRJklQCg2Ar22mn\n4mYQlCRJktSPDIKtzi0kJEmSJPUzg2Crc1N5SZIkSf3MINjqDIKSJEmS+plBsNW1tcG998KGDc2u\nRJIkSdIQYRBsde3t8OSTsHJlsyuRJEmSNEQYBFudW0hIkiRJ6mcGwVbXvam8K4dKkiRJ6icGwVZn\nj6AkSZKkfmYQbHV77gkRBkFJkiRJ/cYg2Oq22w722MOhoZIkSZL6jUFwMHAvQUmSJEn9yCA4GBgE\nJUmSJPUjg+Bg0N7u0FBJkiRJ/cYgOBi0tRUbyj/+eLMrkSRJkjQEGAQHg+4tJOwVlCRJktQPSg2C\nEXFcRPwpIhZExLk1zo+OiG9Xzt8cEVMr7VMj4rGI+F3l9sUy62x57iUoSZIkqR+NKuuFI2IkcDFw\nLLAMuDUi5mXmXVWXnQmszsz9IuJU4CLglMq5v2TmIWXVN6i0tRX39ghKkiRJ6gdl9ghOAxZk5sLM\nfAL4FnBSj2tOAi6vPL4KOCYiosSaBid7BCVJkiT1ozKDYDuwtOp4WaWt5jWZuQ5YA0yonNsnIn4b\nETdExJG13iAiZkXE/IiYv2LFiv6tvpWMHw9jxhgEJUmSJPWLVl0s5l5gcmY+Gzgb+EZE7Nzzosy8\nJDM7MrNj4sSJA17kgIkohoc6NFSSJElSPygzCC4HJlUd711pq3lNRIwCdgFWZebjmbkKIDNvA/4C\nPK3EWlufm8pLkiRJ6idlBsFbgf0jYp+I2B44FZjX45p5wOmVxycD12VmRsTEymIzRMS+wP7AwhJr\nbX0GQUmSJEn9pLRVQzNzXUS8HbgWGAlclpl3RsSFwPzMnAdcClwREQuABynCIsBRwIUR8SSwAXhz\nZj5YVq2DQvfQ0MxiqKgkSZIkbaXSgiBAZl4DXNOj7QNVj9cCr6rxvO8C3y2ztkGnvR0eewweeqhY\nPEaSJEmStlKrLhajntxCQpIkSVI/qdsjGBGv6O2Jmfm9/i9HdVVvKn/QQc2tRZIkSdKg1tvQ0JdW\n7vcAZgDXVY6fD/waMAgOJHsEJUmSJPWTukEwM98AEBE/BQ7MzHsrx3sBXx2Q6rRRd4+gQVCSJEnS\nNmpkjuCk7hBYcT8wuaR6VM+YMTBhgpvKS5IkSdpmjawa+rOIuBb4ZuX4FOB/yitJdbW12SMoSZIk\naZttMQhm5tsj4uUUe/sBXJKZV5dblmpyU3lJkiRJ/aDRfQRvBx7JzP+JiHERsVNmPlJmYaqhvR3u\nuKPZVUiSJEka5LY4RzAi3gRcBfy/SlM78P0yi1IdbW1w//2wbl2zK5EkSZI0iDWyWMzbgMOBhwEy\n8/8otpTQQFu2DDZsgO23h6lTYe7cZlckSZIkaRBqJAg+nplPdB9ExCggyytJNc2dC9/4RvE4ExYv\nhlmzDIOSJEmS+qyRIHhDRPwbMDYijgW+A/yw3LK0mdmz4fHHN23r6iraJUmSJKkPGgmC5wIrgD8A\nZwHXAO8rsyjVsGRJ39olSZIkqY5Gto/YAHypclOzTJ5cDAet1S5JkiRJfdDIqqGHR8R/R8SfI2Jh\nRPw1IhYORHGqMmcOjBu3aduIEXD++c2pR5IkSdKg1cjQ0EuBTwBHAM8FOir3GkidnXDJJTBlCkTA\nhAnFCqLf/jasXdvs6iRJkiQNIo0EwTWZ+ePMfCAzV3XfSq9Mm+vshEWLigC4ciVceilcey288pWb\nLyQjSZIkSXVscY4g8POI+CjwPeDvaSMzby+tKjXmjW8sNpc/6yx41avgqquKPQYlSZIkqReNBMHD\nKvcdVW0JvKD/y1GfzZoF69fDW98Kp5wCV14J223X7KokSZIktbBGVg19/kAUom3wlrcUYfAd74DT\nToNvftMwKEmSJKmuukEwIl6bmV+PiLNrnc/MT5RXlvrs7W8vhon+8z/DEUfAfffB0qXF9hJz5hTz\nCyVJkiSJ3nsEd6jc7zQQhagfvPvdcMstRY9gt8WLi+GjYBiUJEmSBEBkZrNr6BcdHR05f/78ZpfR\nfFOn1t54fsqUYsVRSZIkSUNSRNyWmR1bvrKBOYIRMQY4E3gGMKa7PTPfuNUVqjxLlvStXZIkSdKw\n08g+glcAewIvAm4A9gYeKbMobYPJk/vWLkmSJGnYaSQI7peZ7wcezczLgZewcUsJtZo5c2DcuE3b\nxo0r2iVJkiSJxoLgk5X7hyLiIGAXYI/yStI26eyESy6B9vbiePz44tiFYiRJkiRVNBIEL4mI8cD7\ngXnAXcBHSq1K26azE5YtK8Lgi15kCJQkSZK0iUY2lP9y5eENwL7llqN+NWMG3HRTs6uQJEmS1GJ6\n21C+5kby3dxQfhCYPh2+8x245x5oa2t2NZIkSZJaRG9DQ3fawk2tbsaM4t5eQUmSJElV6vYIZua/\nD2QhKsGznw2jR8Ovfw2vfGWzq5EkSZLUIra4WExE7BsRP4yIFRHxQET8ICKcKzgYbL89dHQUQVCS\nJEmSKhpZNfQbwJXAXkAb8B3gm2UWpX40YwbcfjusXdvsSiRJkiS1iEaC4LjMvCIz11VuXwfGlF2Y\n+sn06fDEE0UYlCRJkiQaC4I/johzI2JqREyJiHOAayJit4jYrewCtY2mTy/uXTBGkiRJUkUjQfDV\nwFnAz4HrgbcApwK3AfN7e2JEHBcRf4qIBRFxbo3zoyPi25XzN0fE1B7nJ0fE3yLiPQ19Gm1uzz1h\nn32cJyhJkiTp7xrZUH6frXnhiBgJXAwcCywDbo2IeZl5V9VlZwKrM3O/iDgVuAg4per8J4Afb837\nq8qMGfCzn0EmRDS7GkmSJElN1siqof9RCXXdxztHxFcaeO1pwILMXJiZTwDfAk7qcc1JwOWVx1cB\nx0QUSSUiXgb8FbizgfdSb6ZPh/vug8WLm12JJEmSpBbQyNDQUcAtEfHMiDgWuJViWOiWtANLq46X\nVdpqXpOZ64A1wISI2BF4L9DrXoYRMSsi5kfE/BUrVjRQ0jDVvbG8w0MlSZIk0UAQzMzzgHOAmyl6\n716SmZ8rua4LgE9m5t+2UNslmdmRmR0TJ04suaRB7OCDYYcdXDBGkiRJEtDAHMGIOAr4DHAhcDDw\n2Yg4MzPv2cJTlwOTqo73rrTVumZZRIwCdgFWAYcBJ0fER4BdgQ0RsXYAAujQNGoUTJtmj6AkSZIk\noIEgCHwMeFX3Ii8R8QrgOuAftvC8W4H9I2IfisB3KvCaHtfMA04HbgJOBq7LzASO7L4gIi4A/mYI\n3EYzZsCHPwyPPlr0DkqSJEkathqZIzi9eqXPzPwecPiWnlSZ8/d24FrgbuDKzLwzIi6MiBMrl11K\nMSdwAXA2sNkWE+on06fD+vVw663NrkSSJElSkzXSI7h7RHwIaM/M4yLiQGA6RYjrVWZeA1zTo+0D\nVY/XAq/awmtc0ECN2pLnPa/PtdPsAAAgAElEQVS4v+kmOProppYiSZIkqbka6RH8KkWv3l6V4z8D\n7y6rIJVkwgQ44ADnCUqSJElqKAjunplXAhvg70M+15dalcoxY0bRI5jZ7EokSZIkNVEjQfDRiJgA\nJEBEPI9ivz8NNjNmwKpV8H//1+xKJEmSJDVRI3MEz6ZY3fOpEfErYCLFCp8abKZPL+5//Wt42tOa\nW4skSZKkpmlkQ/nbgZnADOAs4BmZ+fuyC1MJnv502GUXN5aXJEmShrlGegS75wXeWXItKtuIEcXq\noS4YI0mSJA1rjcwR1FAyYwbceSescZqnJEmSNFwZBIeb6dOLVUNvvrnZlUiSJElqki0GwSi8NiI+\nUDmeHBHTyi9NpTjsMIhweKgkSZI0jDXSI/h5YDpwWuX4EeDi0ipSuXbeGQ4+2AVjJEmSpGGskSB4\nWGa+DVgLkJmrge1LrUrlmj4dfvMb2LCh2ZVIkiRJaoJGguCTETGSjRvKTwRMEIPZjBnw8MNw113N\nrkSSJElSEzQSBD8DXA3sERFzgF8CHyq1KpWremN5SZIkScNOIxvKzwXOAf4TuBd4WWZ+p+zCVKL9\n9oPdd3eeoCRJkjRMbXFD+Yh4HnBnZl5cOd45Ig7LTPcfGKwiil5BewQlSZKkYamRoaFfAP5Wdfy3\nSpsGsxkz4M9/hpUrm12JJEmSpAHWSBCMzMzug8zcQAM9iWpxM2YU97/5TXPrkCRJkjTgGgmCCyPi\nnRGxXeX2LmBh2YWpZB0dMGqUw0MlSZKkYaiRIPhmYAawHFgGHAbMKrMoDYBx4+CQQ1wwRpIkSRqG\ntjjEMzMfAE4dgFo00KZPh0svhXXrit5BSZIkScNCI6uGTgTeBEytvj4z31heWRoQGzZAVxdsvz1M\nngxz5kBnZ7OrkiRJklSyRrqBfgDcCPwPsL7ccjRg5s6Fyy4rHmfC4sUwqzLi1zAoSZIkDWlRtSBo\n7QsifpeZhwxQPVuto6Mj58+f3+wyBo+pU4vw19OUKbBo0UBXI0mSJGkbRcRtmdnRyLWNLBbzXxFx\n/DbWpFazZEnf2iVJkiQNGY0EwXdRhMHHIuLhiHgkIh4uuzCVbPLkvrVLkiRJGjK2GAQzc6fMHJGZ\nYzNz58rxzgNRnEo0Z06xhURPRxwx8LVIkiRJGlAN7RkQEeOB/YEx3W2Z+YuyitIA6F4QZvbsYjjo\npEmw557FIjLHHgunn97c+iRJkiSVppHtI/6JYnjo3sDvgOcBNwEvKLc0la6zc9MVQh9/HE44Ac48\nE8aPhxNPbF5tkiRJkkrT6BzB5wKLM/P5wLOBh0qtSs0xejRcfTUceii8+tVw/fXNrkiSJElSCRoJ\ngmszcy1ARIzOzD8CB5Rblppmxx3hmmtg332LHsHbb292RZIkSZL6WSNBcFlE7Ap8H/jviPgBUGMD\nOg0ZEybAT39aDA897jj485+bXZEkSZKkftTIqqEvz8yHMvMC4P3ApcDLyi5MTbb33vDf/108nj69\nOB4xotiIfu7cppYmSZIkadvUDYIRsXPlfrfuG/AH4JfAjgNUn5rpaU+Dd70LHnwQli+HTFi8GGbN\nMgxKkiRJg1hvPYLfqNzfBsyvca/h4Etf2rytq6vYdkKSJEnSoFR3+4jMPCEiApiZmUsGsCa1kiV1\nfvT12iVJkiS1vF7nCGZmAj8aoFrUiiZPrt3e1jawdUiSJEnqN42sGnp7RDx3a148Io6LiD9FxIKI\nOLfG+dER8e3K+ZsjYmqlfVpE/K5yuyMiXr41769+MGcOjBu3eXtXl6uJSpIkSYNUI0HwMOCmiPhL\nRPw+Iv4QEb/f0pMiYiRwMfBi4EDgtIg4sMdlZwKrM3M/4JPARZX2/wU6MvMQ4Djg/0VE3WGsKlFn\nJ1xyCUyZAhHF/Zw5MGoUHHkk/H6LvwqSJEmSWkwj4epFW/na04AFmbkQICK+BZwE3FV1zUnABZXH\nVwGfi4jIzK6qa8YAuZU1qD90dha3aq98JfzjP8LMmfCTn8BhhzWnNkmSJEl91sg+goszczHwGEUg\n675tSTuwtOp4WaWt5jWZuQ5YA0wAiIjDIuJOii0r3lw5v4mImBUR8yNi/ooVKxooSf3mgAPgxhuL\nzeePOQZ+/vNmVyRJkiSpQVsMghFxYkT8H/BX4AZgEfDjkusiM2/OzGcAzwXOi4gxNa65JDM7MrNj\n4sSJZZeknqZOLcLg1Klw/PHwnvcUj914XpIkSWppjcwR/A/gecCfM3Mf4BjgNw08bzkwqep470pb\nzWsqcwB3AVZVX5CZdwN/Aw5q4D010PbaC66/HvbcEz7+8WLDeTeelyRJklpaI0HwycxcBYyIiBGZ\n+XOgo4Hn3QrsHxH7RMT2wKnAvB7XzANOrzw+GbguM7PynFEAETEF+AeKnki1ot13h/XrN29343lJ\nkiSpJTWyWMxDEbEj8AtgbkQ8ADy6pSdl5rqIeDtwLTASuCwz74yIC4H5mTkPuBS4IiIWAA9ShEWA\nI4BzI+JJYAPw1sxc2dcPpwG0bFntdjeelyRJklpOFHvG93JBxA7AWiCATorhm3MrvYQto6OjI+fP\nn9/sMoavqVOL4aA9TZkCixYNdDWSJEnSsBMRt2VmI6M36w8NjYiLI+LwzHw0M9dn5rrMvDwzP9Nq\nIVAtoNbG8xFw/vnNqUeSJElSXb3NEfwz8LGIWBQRH4mIZw9UURqEem48P3FisWjMDTcU95IkSZJa\nRt0gmJmfzszpwEyKlTwvi4g/RsT5EfG0AatQg0dnZzEMdMMGeOABuOACuPxy+OIXm12ZJEmSpCqN\nbih/UWY+GzgNeBlwd+mVafB7//vhJS+Bd70Lbrqp2dVIkiRJqmhkQ/lREfHSiJhLsZH8n4BXlF6Z\nBr8RI+CKK2DSJDj5ZLjvvmZXJEmSJIneF4s5NiIuA5YBbwJ+BDw1M0/NzB8MVIEa5MaPh6uvhtWr\n4dWvhiefbHZFkiRJ0rDXW4/gecCvgadn5omZ+Y3M3OL+gdJmnvlM+PKX4cYb4Zxzml2NJEmSNOzV\n3VA+M18wkIVoiHvNa+CWW+BTn4Jp0+C005pdkSRJkjRsbXGOoNRvPvpROPJIOP10aGsr5hBOnQpz\n5za7MkmSJGlYMQhq4Gy3HZx6ajFP8N57i/0FFy+GWbMMg5IkSdIAMghqYH3kI5u3dXXB7NkDX4sk\nSZI0TBkENbCWLOlbuyRJkqR+ZxDUwJo8uW/tkiRJkvqdQVADa84cGDdu8/Zjjhn4WiRJkqRhyiCo\ngdXZCZdcAlOmQETRE/isZ8FXv1psPC9JkiSpdAZBDbzOTli0CDZsKFYN/dWvNu4teMMNza5OkiRJ\nGvIMgmq+HXaA//ov2HdfOPFEuOOOZlckSZIkDWkGQbWGCRPg2mth553huONg4cJmVyRJkiQNWQZB\ntY5Jk4ow+MQT8KIXwQMPNLsiSZIkaUgyCKq1HHhgMUx0+XI4/ni49FKYOhVGjCju585tdoWSJEnS\noDeq2QVIm5k+Ha66Ck44AWbNKhaVgWJhmVmzisednc2rT5IkSRrk7BFUazr+eNhtt40hsFtXF8ye\n3ZyaJEmSpCHCIKjW9eCDtduXLBnYOiRJkqQhxiCo1jV5ct/aJUmSJDXEIKjWNWcOjBu3efuRR0Lm\nwNcjSZIkDREGQbWuzk645BKYMgUiYO+94TnPga9/HU45BR5+uNkVSpIkSYOSQVCtrbMTFi0qFo1Z\nuhTmz4eLLoLvfQ+e+1z43/9tdoWSJEnSoGMQ1OASAeecAz/7WdEjOG0aXHFFsb+g+w1KkiRJDTEI\nanCaORNuv73oFXz96+GMM4p9BjM37jdYKwwaGCVJkiQ3lNcgttdeRc/ghAmbzxfs6oKzz4YDDoBd\ndy1u11wDb3lLcQ7coF6SJEnDVuQQWX2xo6Mj58+f3+wy1AwjRmzbKqJTphTzECVJkqRBLCJuy8yO\nRq61R1CD3+TJRe9eT095SrHq6Jo18NBD8M531n7+4sWwdi2MGVNunZIkSVKLcI6gBr9a+w2OGwcf\n/ziceCK87nXwjncUPX/1TJoE//ZvsGRJcexcQkmSJA1hBkENfj33G5wypTjuOe+vXmA87zw44ohi\nW4p99ikWoDnzzMYWn5EkSZIGIYOghobq/QYXLaq9+Eu9wPihD8HVV8PChfCv/1qsRvr445s+t6sL\nzj1387mI9hxKkiRpEHKxGKmn3hafGTcO9t23uD3xRLFq6ZNPbnq+Vm+kJEmSVLK+LBZTao9gRBwX\nEX+KiAURcW6N86Mj4tuV8zdHxNRK+7ERcVtE/KFy/4Iy65Q2MXly7fbx4+Gss4oQuHAhXHvtpiEQ\nip7Dd71r41zDbvYcSpIkqYWUtmpoRIwELgaOBZYBt0bEvMy8q+qyM4HVmblfRJwKXAScAqwEXpqZ\n90TEQcC1QHtZtUqbmDOnmBPYvd8gFD19n/3spj19I+r8f5RVq4php/vuC89/Pmy/PXz1q/DYY8V5\n9y+UJElSk5XZIzgNWJCZCzPzCeBbwEk9rjkJuLzy+CrgmIiIzPxtZt5Tab8TGBsRo0usVdqo0cVn\n6vUc7rUXfOpTcPDB8N3vwhe+sDEEduvqgtmzy6lfkiRJ2oIyg2A7sLTqeBmb9+r9/ZrMXAesASb0\nuOaVwO2Z2WP1DoiIWRExPyLmr1ixot8KlxpafKbeKqQf/WgxPPT734eVK4swWcuSJfXnIkqSJEkl\naulVQyPiGRTDRc+qdT4zL8nMjszsmDhx4sAWJzXSczhyZP2ew0w47DC48kpYt25gapYkSZIoNwgu\nByZVHe9daat5TUSMAnYBVlWO9wauBl6fmX8psU5p621tz+HYsfCGN8BDD8Epp8D++8NnPgOXXeai\nMpIkSSpdaYvFALcC+0fEPhSB71TgNT2umQecDtwEnAxcl5kZEbsCPwLOzcxflVijVL7ucDh7djEc\ndPLkIhx2dsL69fDDH8LHPlYMJ63mojKSJEkqSan7CEbE8cCngJHAZZk5JyIuBOZn5ryIGANcATwb\neBA4NTMXRsT7gPOA/6t6uRdm5gP13st9BDXo7bkn3H9/7fbly+uvUipJkiTRt30E3VBeahW9bWQ/\nYQIceywcdxy88IXFyqRz59buZZQkSdKw1DIbykvqg3qLyuy+O7zkJfDzn8MZZ0BbW3HtGWcUw0cz\nNw4jrTWn0M3sJUmS1EOZcwQl9UW9jew/9amip2/DBvj97+Haa+H88zdfabSrC848E372s2IF0ylT\nYMEC+MQn3MxekiRJm3BoqNRKGh3u2dsw0rY2uPfe3vconDSpeI+teW9JkiS1JOcISkPd1KlF715P\nU6YU21g8/jgsXQpPe1r9QHj00XDMMcVtwQJ485s3743suS+iJEmSWpZzBKWhrtbehOPGFe0Ao0fD\nfvvVn3e4886wZg28//0wYwacfvqmIRCK49mz+792SZIkNZ1BUBqMOjuL3ropUyCiuK/Ve1cvMH7+\n83D77bBiBVx5Zf1ew8WLix7Gai4+I0mSNOg5NFQa6hqZ+1dvqGm3ffcthpCOHg2XXrpx8RlwCKkk\nSVKLcI6gpL6ZO7f2iqXnnw9jxxYrkV5/fTGctJbuuYmSJElqGucISuqbekNNzzkH3vEO+P73YeXK\n4lwtixfDu98NV18Nq1YVbQ4hlSRJaln2CEpqXL0hpGPGFCGxe8jo3nvDffdtutehQ0glSZJKZY+g\npHLUW3zmy1+G1avhxhvhgx8seg9rbXh/7rkDV6skSZLqMghKalxvq5WOHg1HHFEsTPP447Wfv2wZ\nHHts8ZwVK4o2h5BKkiQNOIeGSup/9YaQ7rwz7LFHsYH9yJFwwAHF4yee2HiNQ0glSZK2ikNDJTVX\nb/sX/vnP8NvfwnvfWzyuDoFQDCE955zN9za051CSJKnf2CMoqRyN7F84YkT9zezHj4dDD4WOjiIc\nfulL7l8oSZLUC3sEJTVfZ2ext+CGDcV9rcA2eXLt5+62G5x8Mjz4IHzsY/CZz2waAqEIh+95z6Z7\nH3az91CSJKlXo5pdgKRhbM6c2hvZf+YzG4Pj2rVFW62ew/vugx13hKc+FQ4+GA46CB56qOg9XLu2\nuGbx4uI9wN5DSZKkCnsEJTVPb6uQdhszpn7P4e67wwc+AM96Ftx1VxEsP/vZjSGwW1cXnH32xs3u\nu9lzKEmShinnCEpqfXPn1u457BkaH3sMdtih/rxDgPb2IjiOGgU/+YkrlkqSpCHDOYKShpZGeg4B\nxo6t33u4xx7w0Y/C858PS5fCvHm1VyydPbuczyBJktRC7BGUNLQ02ntYb8XSCFi/vriXJEkaROwR\nlDR8Ndp7WK/nMBOOOgp+/vPya5UkSWoSg6CkoaeRrSvqbXp/+umwcCG84AXFMNIbb3RRGUmSNOQY\nBCUNT/V6Dr/6VfjLX+DTn4a77y56B1//+mIbisyN21HUCoONBkaDpSRJajLnCEpSPV1dxSqjDz20\n+bkdd4TzzoO2tuL229/Cv//7phvf15qb2OgcRkmSpD7qyxxBg6Ak9abeojKNGjUK9t23uB81quhl\nfPLJza+bMqUYxipJkrSVXCxGkvpLvUVlpkwpevUWLoRf/rL+KqPr1sFzngMHHlgEwlohEIohp3fe\n2T81S5IkbYFBUJJ6U29RmTlzin0L99kHDj+898D4zW/Cd74DV19dHNdz0EEwbRp8/vOwerVzCSVJ\nUmkMgpLUm0a3o+gtMDZy3Re+AJ/8JDz+OLztbTBxYrGCaX8uUiNJklThHEFJ6i9z58Ls2bBkSdFD\nOGdO7QVgersus1h45uij4ZFHNn/uuHFw1lnwlKcUtzvvhM99Dtau3fSaeovPNFqjJEkadFwsRpIG\nu94WqRk3btNVR2sZO7boQZw8GSZNKu5vvRXe+97GViw1MEqSNOj0JQiOKrsYSdJWmDy5GA7aU/fq\noo8+CvffD/vtVzswPvYYfPnLxXW96eqCd78bnv70YjGbXXfdfIuL7mGpYBiUJGmIcI6gJLWiLc05\n3GGHIrj1tkjNI4/Agw/CHXfAD39Y/71WroRDD4Xx42HCBHjjGzfvcezqKnoIe3J+oiRJg5JBUJJa\nUX8sUhNRhLtnPhNOOKH+iqV77QXf/S589KPw6lfDE0/Uvm7xYrjgArjmGlixYmPPYSML2kiSpJbi\nHEFJGuz6skhN9ZBPqD1HcOrU2sNSt9uu2Bex+++NkSNh/frNr+sevro1NUqSpK3mhvKSNJx0dhbB\na8OG4r5ewNrWXsavfAUefhhuuKHoPawVAqEIkdOmFb2L55wDZ5wBZ57pVhiSJLWQUnsEI+I44NPA\nSODLmfnhHudHA18DDgVWAadk5qKImABcBTwX+Gpmvn1L72WPoCT1o0Z68Or1HO64I8yYUYTSxYuL\nvRFr2Wkn+NCH4IAD4B/+oQiYZ53lqqaSJG2lltg+IiJGAn8GjgWWAbcCp2XmXVXXvBV4Zma+OSJO\nBV6emadExA7As4GDgIMMgpLUghoZarphA4waVX8rjGoRta/bc0+4/XbYY49iOGqjQ1y7azQwSpKG\niVYZGjoNWJCZCzPzCeBbwEk9rjkJuLzy+CrgmIiIzHw0M38JrEWS1JoaGWo6YkT9lU0nT4bly+G6\n6+ALX6gfFu+7D9raYPRo2Hvv+quannPOpttl9HUxG4elSpKGkTKDYDuwtOp4WaWt5jWZuQ5YA0xo\n9A0iYlZEzI+I+StWrNjGciVJfdbI/MR6cw4/9KEi4D3/+fDmN9df1XTiRPj85+Hcc+HYY+uvanrP\nPcWw1PHj4eCD4Z/+qXZgPPfczUNnX0KjgVGSNAQM6sViMvOSzOzIzI6JEyc2uxxJUi3bukjNJz8J\nb3kLfPCDxYI19QLjhAnwn/8Jr30t7LcfrK0zqGTZMhgzBvbZBw4/HF71quL1a4XG887btM3AKEka\nIkaV+NrLgUlVx3tX2mpdsywiRgG7UCwaI0kaSjo7tzw3r/v8lub0zZlTe47gpz/d2DYY48cXz1++\nvOhF/MMf4JFHate0dGnx2rvvXtzuvnvzgNnVBWefDU9/ehFGJ0yAH/xg0xq7A2P155QkqYnK7BG8\nFdg/IvaJiO2BU4F5Pa6ZB5xeeXwycF0OlY0NJUl918hQ023tYfzsZ+HDH4YrroCf/Qz++Mf6vYy7\n7gpvexsccwy0t9fvZXzgATj00CJ87rRT0StZq4fxPe+BNWs2bbfnUJLUBGVvH3E88CmK7SMuy8w5\nEXEhMD8z50XEGOAKihVCHwROzcyFlecuAnYGtgceAl5YveJoT64aKknaTKOrhja6Emm9XsanPAW+\n+EVYtaq4vfe9vde1115FD+LIkXD99fDkk72/b18+iyRp2GqJ7SMGmkFQkrRNGgla2xoYJ04segXv\nvru43XJL7dVSx44t5i3utx889alw553wvve5ZYYkqVcGQUmSytKfgXHEiPrbZowdC4891nstO+0E\n73gH7LxzcfvDH+Cyy+Dxx3t/X0nSkNQq+whKkjT09Oc8xnp7LE6ZUuyJuHw53HBD/VoeeQQuuqjY\nEuOtby32Y6wOgVCE0Te/udiC46abNu612Je5ic5jlKQhxx5BSZKaZVuHmk6ZAn/9a9Fz+PDDxb6M\nW/p7fcQI2HNPuP9+WL9+Y/vYscU8x9e/futqlCQ1nT2CkiQNBtu6AuqcOcXzxo0rwl29HsbJk4sg\n+f3vw/vfD6tXbxoCoQiTp59erJS6337wvOfBCSfAWWfVXgF19uzN36fRnkN7GCWp6QyCkiQ100Bs\nmfGhDxVh8KST4IIL6m+DAUUYnDatmHN4zz0bh5L2tHgxvOIVRSD8+tfhgx+EN72paM/cuHdiz5DX\n3cO4peu6rzVYSlIpHBoqSdJQ0shiNr0NNV206P+3d+9BclV1Ase/Pwh5GmNeREJCEgziBjZGSaG4\navlYAZUCXHEF45oIJWUVumjtikgVPrBQ1HJ9IEaUt8ozCqTWVRGiJYUQCRoMT8VAAgQJBBIeiYGQ\n3/5x7pjOpHsyw2Sme6a/n6quvn3uubdPz6npyS/nnN/pXt0RI8r977tvx9HFWqNGwYIFJbAcPRq+\n9rUyItnZ5Mlw662l/qhRcMUV3ZuS6tRVSfoHs4ZKkqTGehI87azuc8+VYPDAAxuvTxw7tiS22bKl\n920fMQKOPXZbptSzz4b163esVy+odWsNSYNcTwLBIX3dGEmS1GI6gp/uBEU7qzt0KMyatW0dYmcd\nAVlmmZL6ylfCQw/tWG/8+DK9dOPGMh31s5+t3/ZNm+D662HDhhJcNgo+V62Ct78d9t23PNasgfPO\n2zYttmNKau1n7GDAKKkNOCIoSZJ6r7ujjLsiU2rHSN/WraXegw/uWG/kSJg9G1auhLVrG7d7+PCy\n1nH8eJgwodS/7LIy0tlV+2o/T3eCRoNLSf3ArKGSJKl/dTehza7IlNpht93gy1+uX+/73y/7Jj76\naBk5jKjf7r//HZYuhUsugc99Di6+ePsgEErQOn8+vO518J73wEknlXaceCKccML2iW8+8pFyj9r/\naDdBjqQW5IigJElqTbtytK07I4xbtpSpro3+bXTYYWWK6Zo1sG5d122PKKONI0aUaaz1EupMmAA/\n+xnMmFGOL73UBDmSesVkMZIkSbV25ZRUKCOJI0c2DhpPP72sZ9y0Cc45Z+ftGzWqjEQ+//yO58aM\ngVNOKaN/HaOg9RLk1Fun6ZRUqa0YCEqSJHXWnaCoJ6Nt3Q0aG9Xbay/43vfg/vvL41vfepEfrMas\nWTBzZnmsWweXXw6bN+/8sxgwSoOCawQlSZI6mzevBGhbt5bnRllSu7OGEbq3jrGrel/7Ghx5JJx8\nMnzzm+W96tlnnzIC2ZFRderU+vVGjy5ZWVeuhIULy1rF2iAQyj1OOgmuvBLuvLOMQvZkDSO4PlEa\nJAwEJUmSanUnYOyo1x8Jcr70JRg2rKw3HDmycYKchQvh6qthxYoSMDZKkLNhA7z//WXvx5EjYcGC\n7UdAobz+1KfKHpFr15ZANLNvEt9IagqnhkqSJLWK/kiQM3UqLF4Md98Nd91V9m/sjj32KElvtm7d\n8dzo0XDqqTB2bHksX16munbs2wi9n5bq9FVpp1wjKEmS1O56myBn4kT4+tfhqae2Pc46q3dtGjYM\njj66rI+cPLlMY73wwp2vY+zJ2k0DRrUxA0FJkiT1b4Kce+6BJ58sjwMPbJxRdeZMeOSRMn21kd12\nK+196UtL1tRly0oG1s723BOWLIFJk2DcOLjsMgNGtTWTxUiSJKl/E+QMH15G+mbNKoFVPdOmwV/+\nAs88U0YYG61j3LoV3vSmssfikCH1g0Ao6xcPPLCMXg4bBh/6UP31jh//eAn8fv5zWLq0jHR+5CO7\nfq2j6yI1gDgiKEmSpO7blaOMvd2CY889y1rEtWvh0UdLYp3eGDkSjj++jDBOmlRGOb/73e6vdezu\naKTUR5waKkmSpObalQFjbwPLKVPghhvgiSfK493vbtzuMWNKZtWuRMDLX14yuXY8VqzYcbsOKMl5\nVq/evswEOeojPQkEh/R1YyRJktSG5s3bedDScX5nwU536515Zv2A8ayzyh6LHaZN63okcvPmMso4\nbVr9tY6ZcMQRZcpqx6NeEAjw4INlXeQBB5THhg1wwQXbRhk7pqXWfk7YMfhtVK+jrgGjesgRQUmS\nJA0ezZi62lXdMWPg0EPhzjvhz3+GLVvqt3vIkBIoDh9eRhhvuWX7Kakd9twTrr++rMccPx4uvdQE\nOfoHk8VIkiSpPe3KBDldJcfprFHdc86BK68sgeCzzzZOkLNlSwkmx4wpba8XBEIZqZw9e1uCnPnz\n6yfIOfnksl/kTTeVPSMXLiwB465OkKMByxFBSZIkqZGejKJ1p25vE+RMmgRnnw1r1pRtOL7ylRf5\nwSrjx8N118H++8OoUe7ZOMCZLEaSJElqRf2VIGfyZLj2Wli3riTI+cAHdt626dPhb3+rPxo5eXKZ\nrjpiRJm+evXV8NGPmnraZn8AAA0lSURBVCW1xTg1VJIkSWpF3Z2W2tvpq1/9KsydC4cdBscdV66v\nZ6+94Kqr4Iwz4JBDGk9JXbOmjPpNnAijRzfes/GTnyzBZGdON205jghKkiRJA1l/JMgZP75kX920\nqQSLp5zSdZumTCmB6Ny58PTT8O1vl2u7em/1mlNDJUmSJG2vPwLGSZPg05+GZcvK489/btyevfcu\n22vUJtBx3WGvGAhKkiRJenF2ZcC4fj2MG1d/P0Yo5+bMKY+NG+Gii7afnmqimh4xEJQkSZLUt7ob\njDUaPRw7Fo45BpYvhxUrGq9PHD0aPvvZkrBm8mS47TY4/XSnmtZhIChJkiSpNXRn9HDLFhg6tPHI\nYXd03oKjDZk1VJIkSVJr6E4G1CFDyqhiPdOmwVNPwT33wJIljd9n9epd2+5BzkBQkiRJUt+aN6+M\n1m3dWp7rTeFstBXGmWeW6aH77w9vfWvjrTAAvvhFePbZXdnyQctAUJIkSVLz9WbvxOHD4aCDylrC\nmTPh3HPLdFP3L2zINYKSJEmSBpZGiWp+97uyx+FNN8HLXw5PPgmbN2+7bpAnlWmZNYIRcXhE3BsR\n90XEqXXOD4uIK6rzSyNies25z1Tl90bEYX3ZTkmSJEkDSKOppm94A9x4I1xzDTz++PZBIJSENaed\ntuP92nDksM8CwYjYHTgHeCcwCzguImZ1qnYC8GRmzgS+AXylunYWcCxwAHA48N3qfpIkSZLUWAQc\ndRS88EL986tXw+zZpc4nPwkf+hCccELZ4iKzPJ94Yv1gcBAFjEP68N4HA/dl5kqAiLgcOAq4q6bO\nUcDnq+NFwHciIqryyzNzM3B/RNxX3e/mPmyvJEmSpMFin33q7184ejTMmAErV8L112+/rUWHjRvh\nwx+GSy6BCRNg4kR4+GFYvBiee67U6QgYYUBONe3LqaF7Aw/WvH6oKqtbJzO3ABuA8d28VpIkSZLq\na5SFdOFCuPbason9M8+UEcR6nn8e1q+Hm2+G88+HRYu2BYEdNm4saxUHoAGdNTQiToyIZRGx7LHH\nHmt2cyRJkiS1iu5kIY3oev/CpUvLyOHTTzcOGAfo/oV9GQg+DEyteT2lKqtbJyKGAGOAdd28lsz8\nfmbOzcy5EydO3IVNlyRJkjTg9Xb/wlqNAsZG5S2uLwPBW4H9ImJGRAylJH9Z3KnOYmB+dXwMsCTL\nfhaLgWOrrKIzgP2A3/dhWyVJkiS1o97sX1gvYBwg+ixZTGZuiYiPAb8EdgcuyMw7I+IMYFlmLgbO\nB35YJYN5ghIsUtW7kpJYZgtwUmY2SPsjSZIkSb0wb97OE750nK+3f+EA5IbykiRJkjQItMyG8pIk\nSZKk1mMgKEmSJEltxkBQkiRJktqMgaAkSZIktRkDQUmSJElqMwaCkiRJktRmDAQlSZIkqc0YCEqS\nJElSmzEQlCRJkqQ2YyAoSZIkSW3GQFCSJEmS2oyBoCRJkiS1mcjMZrdhl4iIx4BV/fiWE4DH+/H9\n1D32S+uxT1qT/dJ67JPWY5+0Jvul9dgnrWNaZk7sTsVBEwj2t4hYlplzm90Obc9+aT32SWuyX1qP\nfdJ67JPWZL+0HvtkYHJqqCRJkiS1GQNBSZIkSWozBoIv3veb3QDVZb+0HvukNdkvrcc+aT32SWuy\nX1qPfTIAuUZQkiRJktqMI4KSJEmS1GYMBCVJkiSpzRgIvggRcXhE3BsR90XEqc1uTzuKiAsiYm1E\n3FFTNi4ifhURf6mexzazje0oIqZGxK8j4q6IuDMiTq7K7ZsmiYjhEfH7iLi96pMvVOUzImJp9T12\nRUQMbXZb201E7B4Rf4yI/61e2ydNFhEPRMSKiFgeEcuqMr+/migiXhYRiyLinoi4OyIOsU+aKyL2\nr35HOh5PRcQn7JeBx0CwhyJid+Ac4J3ALOC4iJjV3Fa1pYuAwzuVnQrckJn7ATdUr9W/tgD/lZmz\ngNcDJ1W/H/ZN82wG3paZrwbmAIdHxOuBrwDfyMyZwJPACU1sY7s6Gbi75rV90hremplzavZE8/ur\nub4F/CIzXwW8mvI7Y580UWbeW/2OzAEOAjYCV2O/DDgGgj13MHBfZq7MzOeAy4GjmtymtpOZvwWe\n6FR8FHBxdXwxcHS/Nkpk5iOZ+Yfq+GnKH+y9sW+aJotnqpd7VI8E3gYsqsrtk34WEVOAdwPnVa8D\n+6RV+f3VJBExBngzcD5AZj6XmeuxT1rJ24G/ZuYq7JcBx0Cw5/YGHqx5/VBVpuablJmPVMd/AyY1\nszHtLiKmA68BlmLfNFU1BXE5sBb4FfBXYH1mbqmq+D3W/74JnAJsrV6Pxz5pBQlcFxG3RcSJVZnf\nX80zA3gMuLCaRn1eRIzCPmklxwKXVcf2ywBjIKhBKcu+KO6N0iQR8RLgJ8AnMvOp2nP2Tf/LzBeq\nKTxTKLMaXtXkJrW1iDgCWJuZtzW7LdrBGzPztZTlHydFxJtrT/r91e+GAK8FFmbma4Bn6TTd0D5p\nnmod85HAVZ3P2S8Dg4Fgzz0MTK15PaUqU/M9GhF7AVTPa5vcnrYUEXtQgsAfZ+ZPq2L7pgVUU6p+\nDRwCvCwihlSn/B7rX/8CHBkRD1CWF7yNsg7KPmmyzHy4el5LWfN0MH5/NdNDwEOZubR6vYgSGNon\nreGdwB8y89Hqtf0ywBgI9tytwH5VdrehlCHxxU1uk4rFwPzqeD5wbRPb0paqdU7nA3dn5v/UnLJv\nmiQiJkbEy6rjEcA7KGs3fw0cU1WzT/pRZn4mM6dk5nTK35AlmTkP+6SpImJURIzuOAYOBe7A76+m\nycy/AQ9GxP5V0duBu7BPWsVxbJsWCvbLgBNl5FY9ERHvoqzv2B24IDPPbHKT2k5EXAa8BZgAPAp8\nDrgGuBLYB1gF/Htmdk4ooz4UEW8EbgRWsG3t02mUdYL2TRNExGzKov3dKf/5d2VmnhER+1JGo8YB\nfwQ+mJmbm9fS9hQRbwH+OzOPsE+aq/r5X129HAJcmplnRsR4/P5qmoiYQ0mqNBRYCXyY6rsM+6Rp\nqv8sWQ3sm5kbqjJ/VwYYA0FJkiRJajNODZUkSZKkNmMgKEmSJEltxkBQkiRJktqMgaAkSZIktRkD\nQUmSJElqMwaCkqSmi4gXImJ5RNwREVdFxMgG9f6vY1/EHt5/ckQs6kX7HoiICS/2+oEiIhZExORm\nt0OS1PcMBCVJrWBTZs7JzAOB54CP1p6MYrfMfFdmru/pzTNzTWYes/OabW8BYCAoSW3AQFCS1Gpu\nBGZGxPSIuDciLgHuAKZ2jMxV5+6OiB9ExJ0RcV1EjACIiJkRcX1E3B4Rf4iIV1T176jOL4iIayPi\nNxHxl4j4XMcbR8Q1EXFbdc8Td9bQiDi8eo/bI+KGqmxcdZ8/RcQtETG7Kv98RFwcETdGxKqI+LeI\n+GpErIiIX0TEHlW9B2rKfx8RM6vy6RGxpLrvDRGxT1V+UUR8OyJ+FxErI+KYmvZ9KiJura75Qs19\ndvjZVdfNBX5cjc6O2AV9KUlqUQaCkqSWERFDgHcCK6qi/YDvZuYBmbmqU/X9gHMy8wBgPfDeqvzH\nVfmrgTcAj9R5q4Or+rOB90XE3Kr8+Mw8iBIQ/WdEjO+irROBHwDvrd7rfdWpLwB/zMzZwGnAJTWX\nvQJ4G3Ak8CPg15n5z8Am4N019TZU5d8BvlmVnQ1cXN33x8C3a+rvBbwROAI4q2rfodXP6GBgDnBQ\nRLy50c8uMxcBy4B51ejspkafXZI08BkISpJawYiIWE4JRFYD51flqzLzlgbX3J+Zy6vj24DpETEa\n2DszrwbIzL9n5sY61/4qM9dVwc5PKUEUlODvduAWYColYGrk9cBvM/P+6r2eqMrfCPywKlsCjI+I\nl1bnfp6Zz1MC3d2BX1TlK4DpNfe+rOb5kOr4EODS6viHNW0GuCYzt2bmXcCkquzQ6vFH4A/Aq2o+\nzw4/uy4+pyRpEBrS7AZIkkS1RrC2ICIAnu3ims01xy8APZnKmJ1fR8RbgH8FDsnMjRHxG2B4D+7Z\nHZsBMnNrRDyfmR3t2Mr2f5OzwXGX961EzfOXM/Pc2ooRMZ3e/ewkSYOAI4KSpEEjM58GHoqIowEi\nYliDDKTvqNbyjQCOBm4CxgBPVkHgqygjfl25BXhzRMyo3mtcVX4jMK8qewvweGY+1cOP8v6a55ur\n498Bx1bH86r36covgeMj4iVVW/aOiD13cs3TwOgetlWSNAA5IihJGmz+Azg3Is4Anqes3dvaqc7v\ngZ8AU4AfZeayiFgBfDQi7gbupQR6DWXmY1VCmZ9GxG7AWuAdwOeBCyLiT8BGYP6L+Axjq+s3A8dV\nZR8HLoyITwGPAR/eSfuui4h/Am6uRlefAT5IGQFs5CLgexGxiTIy6jpBSRqkYtusFEmSBr+IWADM\nzcyPNbst9UTEA5T2Pd7stkiSBi+nhkqSJElSm3FEUJIkSZLajCOCkiRJktRmDAQlSZIkqc0YCEqS\nJElSmzEQlCRJkqQ2YyAoSZIkSW3m/wFQBFBOD29GuAAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.figure(figsize=(15, 7))\n", "plt.plot(range(1,len(rho)+1),rho,'ro-')\n", "plt.title('Variance explained by principal components');\n", "plt.xlabel('Principal component');\n", "plt.ylabel('Variance explained');\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Here plotted is the vcariance explaned of every principal component. From 0(0%) to 1(100%). " ] }, { "cell_type": "code", "execution_count": 218, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAgAAAAHwCAYAAADQAtd+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAIABJREFUeJzt3X+cXHV97/HXJyGAEVw00GCB7HIB\nW4EYlBXrvdSi8VZCpQgtFF0VULtNH/JQ7KN6oest0nar1+o1tlf0bisP0K6CXkWh/qrEIlq0uLFI\nAFHB7sYgP0KEiI1gIN/7x5yNs5vdndndmTlnznk9H495ZOecMzPfzMzueZ/vz0gpIUmSqmVJ3gWQ\nJEmdZwCQJKmCDACSJFWQAUCSpAoyAEiSVEEGAEmSKsgAIKmjIuKUiNg6x/4PRcT/nOnYiLgjIk6Z\n47FfiIjzWlpgqaQMANI8RMSrImIsIn4WEfdlJ5yTO/C6KSKObvFzPi0iNkTEluz/c092/+BWvs58\npZTWp5T+cpZ9x6WUbgSIiHdExD9O278upXRVO8sXETdGxGPZe/ZQRHw6Ip5Zt/+kiPh8RDwSET+J\niFsi4oJpz3FkROyOiA+2s6zSXAwAUpMi4k+ADcBfAyuBVcDlwBl5lgsgIvaZ5/H7AhuB44BTgacB\nLwS2Aye1vIDlc2FK6QDgWcBBwPsAIuKFwFeArwJHAyuAPwbWTXv8a4GHgT+IiP06VWhpipSSN2/e\nGtyAHuBnwNlzHLMftYDw4+y2Adgv23c+8PVpxyfg6OznK4EPAJ8DHgX+DTgq23dTdux/ZmX4A+AU\nYCvwP4D7gY8CtwOn1z3/MuAh4LkzlPUNwAPAAXP8fy4G7snKcydwZt2+84F/Bf4PsAO4C1hbt/8C\n4LvZY38I/FHdvsmy/1lWvnFgoG7/lcBf1R9bt28ceCm10PILYFf2nnwn238j8Ia641+XleNh4EtA\nb7Y9qJ20HwR+CmwGjm/yuzD9Nd4I3J79/HXgAw0eH9n7+sfZZ/D7eX+/vVXzZg2A1JwXAvsD185x\nzBDwG8AJwBpqV9Jvn8drnAtcBjwduBsYBkgpvSjbvyaldEBK6Zrs/qHAM4BeYBD4CPDquuc7Dbgv\npfTvM7zWS4EvppR+Nkd57gF+k1r4uQz4x/qqbuAF2TEHA5cCn46IZ2T7HgReTq1m4QLgfRHxvLrH\nHpo97jDgPGAkIn5tjrJMkVL6IrWamGuy92TN9GMi4gxqIeMs4BDga8DHs92/DbyI2hV8D3AOtdqP\nyWae25opR9Zc8nvAv0fEcmrfk//X4GEnA4cDVwOfoPb/lzrOACA1ZwXwUErpiTmOGQD+IqX0YEpp\nG7WT5mvm8RrXppRuyV5jlFqQmMtu4NKU0uMppZ8D/wicFhFPy/a/hlrNwExWAPfN9eQppU+mlH6c\nUtqdhY4fMLV54EFgQ0ppV7b/e8DvZI/9XErpnlTzVeCfqYWJev8zK/tXqdV8nNPg/ztf64F3ppS+\nm72nfw2cEBG91GoODgR+HYjsmPuysn8spfScBs/9txHxCPAdau/jn1ALbkto8L5SO+F/IaX0MPAx\n4NSI+JWF/RelhTMASM3ZDhzcoK39V4GJuvsT2bZm3V/3807ggAbHb0spPTZ5J6X0Y2rV8r8XEQdR\na3ceneWx24FnzrIPgIh4bUTcmnVmewQ4ntpV+6R7U0r1q4nt+f9GxLqI+GbWCe4RarUR9Y99OKX0\nnzM9toV6gffXlf8n1KrfD0spfYVa88UHgAcjYqQuODXjTSmlg1JKh6WUBrLA9zC1UDbr+xoRTwHO\nJvtcUkrfALYAr1rA/09aFAOA1JxvAI8Dr5jjmB9TO+lMWpVtg1r7/fLJHRFxaAvKNNNSnldRawY4\nG/hGSuneWR57A/CyiHjqTDuzq+S/By4EVqSUDqLWxyDqDjssIurvrwJ+nHVq+xTwHmBl9tjPT3vs\n06e9dv171axGS5n+iFrfg4Pqbk9JKd0MkFL625TSicCx1JoC3jrP159amJR2Uvue/N4ch51JrVnk\n8oi4PyLu55fNIFJHGQCkJqSUdgB/DnwgIl4REcsjYll2pfvu7LCPA2+PiEOytuE/p1YtD7Wq4uMi\n4oSI2B94xzyL8ADwX5o47jPA84A3U+sTMJuPUjtBfioifj0ilkTEioj4s4g4DXgqtRPsNoBsGNvx\n057jV4A3Ze/D2cCzqZ3o96XWIXIb8ERErKPW5j7dZRGxb0T8JrX+Ap9s4v9X7wGgLyJm+zv2IeCS\niDgu+z/0ZOUkIp4fES+IiGXUwtlj1K7eF+ttwPkR8daIWJG91pqIuDrbfx5wBbCaWhPPCcB/A9ZE\nxOoWvL7UNAOA1KSU0nuptfW+ndrJ7UfUrpA/kx3yV8AYcBu1XuXfzraRUvo+8BfUrrx/QK23+Hy8\nA7gqq86eta086wvwKeBI4NNzHPc4tY6AdwFfptYT/hZq1fT/llK6E3gvtSvaB6idsP512tP8G3AM\ntZ78w9R6s29PKT0KvIlaB7eHqVVvXzftsfdn+35MrTp8fUrprrnfgr1MBobtEfHtGf6P1wL/C7g6\nIn5KrQZjcjje06jVcDxMrflhO/A3ABExEBF3zLMsk695M/CS7PbDiPgJMAJ8PiIOA9ZS6zdxf91t\nE/BFrAVQh8XUJjxJ3S4i/hx4Vkrp1Q0PXvhrnE9tKFzbJ0GS1B7zmjxEUrFlw/Bez/xGH0iqIJsA\npJKIiD+k1izxhZTSTXmXR1Kx2QQgSVIFWQMgSVIFGQAkSaqgUncCPPjgg1NfX1/exZAkqWM2bdr0\nUErpkEbHlToA9PX1MTY2lncxJEnqmIiYaHyUTQCSJFWSAUCSpAoyAEiSVEGl7gMgSeouu3btYuvW\nrTz22GOND664/fffn8MPP5xly5Yt6PEGAElSYWzdupUDDzyQvr4+pq42rXopJbZv387WrVs58sgj\nF/QcNgFIkgrjscceY8WKFZ78G4gIVqxYsaiaEgOAJKlQPPk3Z7HvkwFAkqRp7r//fs4991yOOuoo\nTjzxRE477TS+//3vt+z5b7zxRm6++eaWPd9CGAAkSV1rdPMofRv6WHLZEvo29DG6eXTRz5lS4swz\nz+SUU07hnnvuYdOmTbzzne/kgQceaEGJa+YKAE888UTLXmcuBgBJUlca3TzK4PWDTOyYIJGY2DHB\n4PWDiw4B//Iv/8KyZctYv379nm1r1qzh5JNP5q1vfSvHH388q1ev5pprrgFqJ/OXv/zle4698MIL\nufLKK4HajLSXXnopz3ve81i9ejV33XUX4+PjfOhDH+J973sfJ5xwAl/72tc4//zzWb9+PS94wQt4\n29vexjHHHMO2bdsA2L17N0cfffSe+63iKABJUlca2jjEzl07p2zbuWsnQxuHGFg9sODnvf322znx\nxBP32v7pT3+aW2+9le985zs89NBDPP/5z+dFL3pRw+c7+OCD+fa3v83ll1/Oe97zHv7hH/6B9evX\nc8ABB/Cnf/qnAHz4wx9m69at3HzzzSxdupSenh5GR0e56KKLuOGGG1izZg2HHNJwev95sQZAktSV\ntuzYMq/ti/X1r3+dV77ylSxdupSVK1fyW7/1W3zrW99q+LizzjoLgBNPPJHx8fFZjzv77LNZunQp\nAK973ev4yEc+AsAVV1zBBRdcsPj/wDQGAElSV1rVs2pe25t13HHHsWnTpqaP32effdi9e/ee+9OH\n5u23334ALF26dM72/ac+9al7fj7iiCNYuXIlX/nKV7jllltYt25d0+VplgFAktSVhtcOs3zZ8inb\nli9bzvDa4UU970te8hIef/xxRkZG9my77bbbOOigg7jmmmt48skn2bZtGzfddBMnnXQSvb293Hnn\nnTz++OM88sgjbNy4seFrHHjggTz66KNzHvOGN7yBV7/61VNqBlrJACBJ6koDqwcYOX2E3p5egqC3\np5eR00cW1f4PtfH11157LTfccANHHXUUxx13HJdccgmvetWreM5znsOaNWt4yUtewrvf/W4OPfRQ\njjjiCM455xyOP/54zjnnHJ773Oc2fI3TTz+da6+9dk8nwJn87u/+Lj/72c/aUv0PECmltjxxEfT3\n96exsbG8i5Gb0c2jDG0cYsuOLazqWcXw2uFF/2JIUjt997vf5dnPfnbexSiEsbEx3vKWt8waEGDm\n9ysiNqWU+hs9v6MASmpyeMxkD9nJ4TGAIUCSCu5d73oXH/zgBxkdXfy8BrOxCaCk5hoeI0kqtosv\nvpiJiQlOPvnktr2GAaCkOj08RpLUXQwAJdWu4TGS1G5l7pvWSot9nwwAJTW8dph9l+47Zdu+S/dd\n9PAYSWqn/fffn+3btxsCGkgpsX37dvbff/8FP4edAEts+i+Qv1CSiu7www9n69atLZ/3voz2339/\nDj/88AU/3mGAJdW3oY+JHRN7be/t6WX8ovHOF0iS1BHNDgO0CaCk7AQoSZqLAaCk7AQoSZqLAaCk\n2jVHtiSpHAwAJdWuObIlSeVgJ0BJkkrEToCSJGlWBgBJkirIACBJUgUZACRJqiADgPYY3TxK34Y+\nlly2hL4NfYxubt861JKkfLkWgIDayX/w+kF27toJwMSOCQavHwRw6KAklZA1AAJgaOPQnpP/pJ27\ndjK0cSinEkmS2skAIMC1AySpagwAAlw7QJKqxgAgwLUDJKlqDAACXDtAkqrGtQAkSSoR1wKQJEmz\nMgBIklRBBgBJkirIACBJUgUZACRJqiADgCRJFWQAkCSpggwAkiRVkAFAkqQKMgBIklRBBgAxunmU\nvg19LLlsCX0b+hjdPJp3kSRJbbZP3gVQvkY3jzJ4/SA7d+0EYGLHBIPXDwK4EJAklZg1ABU3tHFo\nz8l/0s5dOxnaOJRTiSRJnWAAqLgtO7bMa7skqRwMABW3qmfVvLZLksrBAFBxw2uHWb5s+ZRty5ct\nZ3jtcE4lkiR1ggGg4gZWDzBy+gi9Pb0EQW9PLyOnj9gBUJJKLlJKeZehbfr7+9PY2FjexZAkqWMi\nYlNKqb/RcdYASJJUQQYANcXJgiSpXJwISA05WZAklY81AGrIyYIkqXwMAGrIyYIkqXwMAGrIyYIk\nqXwMAGrIyYIkqXwMAGrIyYIkqXxynQgoIq4AXg48mFI6Ptv2DOAaoA8YB85JKT0cEQG8HzgN2Amc\nn1L69lzP70RAkqSq6ZaJgK4ETp227WJgY0rpGGBjdh9gHXBMdhsEPtihMkqSVDq5BoCU0k3AT6Zt\nPgO4Kvv5KuAVdds/kmq+CRwUEc/sTEklSSqXvGsAZrIypXRf9vP9wMrs58OAH9UdtzXbJkmS5qmI\nAWCPVOugMK9OChExGBFjETG2bdu2NpVMkqTuVsQA8MBk1X7274PZ9nuBI+qOOzzbNkVKaSSl1J9S\n6j/kkEPaXlhJkrpREQPAdcB52c/nAZ+t2/7aqPkNYEddU4EkSZqHXBcDioiPA6cAB0fEVuBS4F3A\nJyLi9cAEcE52+OepDQG8m9owwAs6XmBJkkoi1wCQUnrlLLvWznBsAt7Y3hJJklQNRWwCkCRJbWYA\nkCSpggwAkiRVkAFAkqQKMgBIklRBBgBJkirIACBJUgUZACRJqiADgCRJFWQAkCSpggwAkiRVkAFA\nkqQKMgBIklRBBgBJkirIACBJUgUZACRJqiADgCRJFWQAkCSpggwAkiRVkAFAkqQKMgBIklRBBgBJ\nkirIACBJUgUZACRJqiADgCRJFWQAkCSpggwAkiRVkAFAkqQKMgBIklRBBgBJkirIACBJUgUZACRJ\nqiADgCRJFWQAkCSpggwAkiRVkAFAkqQKMgBIklRBBgBJkirIACBJUgUZACRJqiADgCRJFWQAkCSp\nggwAkiRVkAFAkqQKMgBIklRBBgBJkirIACBJUgUZACRJqiADgCRJFWQAkCSpggwAkiRVkAFAkqQK\nMgBIklRBBgBJkirIACBJUgUZACRJqiADgCRJFWQAkCSpggwAkiRVkAFAkqQKMgBIklRBBgBJkirI\nACCpK4xuHqVvQx9LLltC34Y+RjeP5l0kqavtk3cBJKmR0c2jDF4/yM5dOwGY2DHB4PWDAAysHsiz\naFLXsgZAUuENbRzac/KftHPXToY2DuVUIqn7GQAkFd6WHVvmtV1SYwYASYW3qmfVvLZLaqywASAi\nxiNic0TcGhFj2bZnRMSXI+IH2b9Pz7ucktpveO0wy5ctn7Jt+bLlDK8dzqlEUvcrbADIvDildEJK\nqT+7fzGwMaV0DLAxuy+p5AZWDzBy+gi9Pb0EQW9PLyOnj9gBUFqESCnlXYYZRcQ40J9Seqhu2/eA\nU1JK90XEM4EbU0q/Nttz9Pf3p7GxsfYXVl1hdPMoQxuH2LJjC6t6VjG8dtgTiKTSiYhNdRfOsyry\nMMAE/HNEJOD/ppRGgJUppfuy/fcDK3MrnbqKw8gkaaoiNwGcnFJ6HrAOeGNEvKh+Z6pVXexVfRER\ngxExFhFj27Zt61BRVXQOI5OkqQobAFJK92b/PghcC5wEPJBV/ZP9++AMjxtJKfWnlPoPOeSQThZZ\nBeYwMkmaqpABICKeGhEHTv4M/DZwO3AdcF522HnAZ/MpobqNw8gkaapCBgBqbftfj4jvALcAn0sp\nfRF4F/DfI+IHwEuz+1JDDiOTpKkK2QkwpfRDYM0M27cDaztfInW7yY5+jgKQpJrCDgNsBYcBSpKq\nptlhgEVtApAWzGVjJamxQjYBSAvleH9Jao41ACoVx/tLUnMMACoVx/tLUnMMACoVx/tLUnMMACoV\nx/tLUnMMACoVl42VpOY4D4AKy+V7JWn+nAdAXW1yON/EjgkSac9wvrKO6XfuAkmdZgBQIVVpOF/V\nwo6kYjAAqJCqNJyvSmFHUnEYAFRIVRrOV6WwI6k4DAAqpCoN56tS2JFUHAYAFVKVhvNVKexIKg6H\nAUoF4JBHSa3S7DBAA4AkSSXiPACSJGlWBgB1LSfPkaSF2yfvAkgLMTl5zuT4+cnJcwDbziWpCdYA\nqCvNZ/IcawokaW/WAKgrNTt5jjUFkjQzawDUlZqdPMdpdpUna59UZAYAdaVmJ89xmt1qKdIJ10We\nVHQGAHWlZmcKXOw0u0U6oUxX5LLloWgnXGufVHROBKRSm94HAGo1Bc1MK7yYx7Zbkcs2m3bPdti3\noY+JHRN7be/t6WX8ovGWvU6zlly2hMTef1+DYPeluzteHlWHEwFJLG5NgSJfwXXbKIhOXJ0XrbnH\nRZ5UdAYAld7A6gHGLxpn96W7Gb9ovOmrzvmcUDp9kp3vKIi8q8U7EaaKdsJ1kScVnQFAmkWzJ5Q8\nTrLdNgqiE1fnRTvhVmlFS3UnA4A0i2ZPKHmcZLttFEQnrs6LeMJdaO2T1AlOBCTNYvKPdaOOa3mc\nZJst26qeVTN2jOt0tfjw2uEZOy22+up8YPWAJ1mpSQYAVVYzvdIbnVBGN4+yJJbwZHpyr33tPsk2\nc7Lr1Im3kWYDi6TOcRigKmmmYXRBsL5/PZf/zuULfo5JRRqS1+7hd5KKpdlhgAYAVdJsY8aD4KNn\nfbSpE+Rsz7E0lnLVmVfteQ5PwJI6yXkAVGmNhuXN1j6fSE133pvtOXan3VNO/kUYhtdqRZhbQNLi\nGABUOs2cdOdqn2+2814zPdtbNUKgSCfcsoYaqWoMACqdZk66w2uHCWLGx880zn+mk28zQ/FaMUKg\naCfcoswtIGlxDAAqnWZOugOrB1jfv36vEDD9BD66eZQLPnPBlJPvBZ+5gNHNo02NO3/GU54xY1lm\n2z6Top1wizK3gKTFMQCodJqddOby37mcj5710TlP4G/+wpvZtXvXlMft2r2LN3/hzUBnJnop2gm3\naFPuSloYA4BKZz5TwjY6gW//+fYZX2O27dP95Oc/mdf2mRTthFu0KXeLrkj9N6R6BgCVTpGmhG3F\nybtoJ9wivb9FV7T+G1I95wGQ5nDwuw+e8Wp/xVNW8NDbHmr4+JkmC1rIJEHOJdCdZpsrorenl/GL\nxjtfIFVCs/MAOBWwNIf3r3s/r/vs6/jFk7/Ys23fpfvy/nXvb+rxrZoC1znuu1PR+m9I9QwA0hxa\ncQL35F1dRVmMSZqJfQCkBlzSVQtVtP4bRWZnyc4zAEjqKt10orDDZHPsLJkPOwFKLWJHvfZrVadK\nFYudJVvLxYCkFpvrytMrmM4o2qyIag07S+bDACA1odEJ3hPTws2nSt8TRTkVbbKrqjAASE1odIIv\n2ompW9rJ51tz4ominKraWTLv31MDgNSERif4Ip2Yuqk5Yr41J1U9UZRdFTtLFuH31AAgNaHRCb5I\nJ6Zuao6Yb81JFU8UVVG14bZF+D11IiCpCcNrh2fsfT55gm/VjH+tULTmiLksZKIcJ1ZSGRTh99Qa\nAKkJzVx5FuUKpkjNEY0UqeZE6qQi/J4aAKQmFeUE30g3nVSt0ldVFeH31ImApBJyUiKVSVm/z+36\nfzU7EZABQJJUWM7+OH8tmwkwIp4WEUfNsP05Cy2cJEnNKEJv+bKaMwBExDnAXcCnIuKOiHh+3e4r\n21kwSZKK0Fu+rBrVAPwZcGJK6QTgAuCjEXFmti/aWjJJUiHkOWNdEXrLl1WjALA0pXQfQErpFuDF\nwNsj4k1AeTsPSJKA/GesK0Jv+bJqFAAerW//z8LAKcAZwHFtLJckqQDyboN3qGj7NJoJ8I+ZVtWf\nUno0Ik4FzmlbqSRJhVCENnhnf2yPRjUA/wmsnGH7ScA3W18cSRLkv1LcJNvgy6tRANgA/HSG7T/N\n9kmSWizvdvd6tsGXV6MAsDKltHn6xmxbX1tKJGmKolwJqnPybnevZxt8eTXqA3DQHPue0sqCSNrb\n9FnQJq8EAf8AL1KRp5ctQrt7Pdvgy6lRDcBYRPzh9I0R8QZgU3uKJGlSka4Ey6RIVewzsd1dndAo\nAFwEXBARN0bEe7PbV4HXA29uf/H2FhGnRsT3IuLuiLg4jzJInVK0K8GyKHqwst1dnTBnAEgpPZBS\n+q/AZcB4drsspfTClNL97S/eVBGxFPgAsA44FnhlRBzb6XKovWzz/iWvBNujUbDK+ztou7s6Yc4+\nABGxP7AeOBrYDHw4pfREJwo2i5OAu1NKP8zKdzW1SYnuzLFMaiHbvKcaXjs840poXgkuzqqeVUzs\nmJhxe1G+g7a7q90aNQFcBfRTO/mvA97T9hLN7TDgR3X3t2bb9oiIwYgYi4ixbdu2dbRwWryiV812\nmleC7TFXFbvfQVVFo1EAx6aUVgNExIeBW9pfpMVJKY0AIwD9/f2uV9BlbPPem1eCrTf5fs40CuA1\nn37NjI+p8ndQ5dQoAOya/CGl9ERE7gsA3gscUXf/8GybSmKuqlmplWYLVn4HVRWNmgDWRMRPs9uj\nwHMmf46ImWYIbLdvAcdExJERsS9wLnBdDuVQm7Sq93PenbjUveyBr6qYswYgpbS0UwVpRlYLcSHw\nJWApcEVK6Y6ci6UWmqtqtllF6cSl7tSK76DUDSKl8jaT9/f3p7GxsbyLoQ7r29A3YxVub08v4xeN\nd75AkvZS5JkYu11EbEop9Tc6rlETgNR17EgoFVsnZ2K0OXB2BgCVjpPndAf/MFdXp4Za5jXlc7d8\ntw0AKh07cRVf0efiV3t1qpYujzkduum7bQBQ6Th5TmN5X6E42U61daqWLo/mwG76bjeaB0DqSk6e\nM7sijJKwn0a1dWqK6zzmdOim77Y1AFLFFOEKxX4a1dapWro8mgO76bttAJAqpghXKPbT0MDqAcYv\nGmf3pbsZv2i8LbVPeTQHdtN32yYAqWKKMNWtk+2oUzrdHNhN320nApIqZnofAKhdodhRUioHJwKS\nNCNHSUgCawAkSSoVawAkSdKsDACS1KXyntBJ3c1RAJLUhYowoZO6mzUAktSFijChk7qbAUBdz2pQ\nVVERJnRSdzMAqKt108pbUit105Sz3aJqFxMGAHU1q0FVtT/ak7ppytluUMWLCQOAuprVoNVWxT/a\nk5zQqbWqeDFhAFBXsxq02qr4R7teJxbUqYpOXEwUrbbKAKCuZjVotVkDpFZp98VEEWurDADqalaD\nVps1QGqVdl9MFLG2yomA1PU6vdynimN47fCMKxtaA6T5avcyvkWsrTIASOpa3bT2uoqvnRcTq3pW\nMbFjYsbteTEASOpq1gCpGxSxtso+AJIktVkR+ytFSim3F2+3/v7+NDY2lncxJEnqmIjYlFLqb3Sc\nNQCSJFWQAUCSpAoyAEiSVEEGAEmSKsgAIElSBRkAJEmqIAOAJCl3RVsprwqcCVCSlKvJlfImZ8mb\nXCkPcJbHNrIGQJKUqyKulFcFBgBJUq6KuFJeFRgAJEm5mm1FvDxXyqsCA4CklrNDl+ZjeO0wy5ct\nn7It75XyqsAAUCD+0VQZTHbomtgxQSLt6dDl91mzKeJKeVXgaoAFMb0XLNQSsL8E6jZ9G/qY2DGx\n1/benl7GLxrvfIGkinE1wC5jL1iVhR26pO5gACgI/2iqLOzQpVawSbT9DAAF4R9NlYUdurRY9iPp\nDANAQfhHU2Vhhy4tlk2ineFUwAUx+cdxaOMQW3ZsYVXPKobXDvtHU11pYPWA310tmE2inWEAKBD/\naBbP6OZRQ5nUYat6Vs04ksQm0dayCUCahe2QUj5sEu0MA4A0C9shpXzYj6QzbAKQZmE7pJQfm0Tb\nzxoAaRYOzZRUZgYAaRa2Q0oqMwOANItOt0M685mkTnIxIKkAXAxKUqu4GJDURRxxIKnTDABSATji\nQFKnGQCkAnDEgaROMwBIBeCIA0mdZgCQCsCZzyR1mqMAJEkqEUcBSJKkWRkAJEmqIAOAJGkvzkxZ\nfq4GKEmaYvrMlBM7Jhi8fhDAjqklYg2AJGkKZ6asBgOAJGkKZ6asBgOAJGkKZ6ashsIFgIh4R0Tc\nGxG3ZrfT6vZdEhF3R8T3IuJleZZTksrKmSmroaidAN+XUnpP/YaIOBY4FzgO+FXghoh4VkrpyTwK\nKEllNdnRb2jjEFt2bGFVzyqG1w7bAbBkihoAZnIGcHVK6XHgPyLibuAk4Bv5FkvNGt086h8UqUsM\nrB7w97PkCtcEkLkwIm6LiCsi4unZtsOAH9UdszXbpi4wOaxoYscEibRnWJFjiyUpH7kEgIi4ISJu\nn+F2BvBB4CjgBOA+4L3zfO6+fmrDAAALq0lEQVTBiBiLiLFt27a1ofRaCIcVSVKx5NIEkFJ6aTPH\nRcTfA/+U3b0XOKJu9+HZtunPPQKMQG0xoMWVVK3isCJJKpbCNQFExDPr7p4J3J79fB1wbkTsFxFH\nAscAt3S6fFoYhxW1l9O2SpqvwgUA4N0RsTkibgNeDLwFIKV0B/AJ4E7gi8AbHQHQPRxW1D72r5C0\nEJFSeWvJ+/v709jYWN7FUMZRAO3Rt6GPiR0Te23v7ell/KLxzhdIUq4iYlNKqb/Rcd00DFBdzmFF\n7WH/CkkLUcQmAHU526M7y/4VkhbCAKCWsj268+xfIWkhDABqKcf7d97A6gFGTh+ht6eXIOjt6WXk\n9BGbWyTNyT4Aainbo/Nh/wpJ82UNgFrK9mhJ6g4GALWU7dHS3uwYqyIyALRY1X/RbY+WprJjrIrK\niYBaaPIXvb4T3PJlyz0BShXmRE3qtGYnArIGoIXsAS9pOjvGqqgMAC3kL7qk6ewYq6IyALSQv+iS\nprNjrIrKANBC/qKrm1S9w2qn2DFWRWUnwBZzxTt1AzusSuXVbCdAA4DUpRYTNu2ZLpWXywFLJTb9\nCn5ybDnQVAiww6ok+wBIXWixQ07tsCrJACB1ocVewdthVZIBQOpCi72Ct2e6JPsASF1oeO3wjL34\n53MF7xLCUrVZAyB1Ia/gJS2WwwAlSSoRFwOSJEmzMgBIklRBBgBJkirIACBJUgUZACSpwFy1Ue3i\nPACSVFCLXfNBmos1AJJUUItd80GaiwFAkgrKVRvVTgYASSooV21UOxkAJKmgXLVR7WQAkKSCcs0H\ntZNrAUiSVCKuBSBJkmZlAJAkqYIMAJIkVZABQJKkCjIASJJUQQYAFZaLoEhS+7gYkArJRVAkqb2s\nAVAhuQiKJLWXAUCF5CIoktReBgAVkougSFJ7GQBUSC6CIkntZQBQIbkIiiS1l4sBSZJUIi4GJEmS\nZmUAkCSpggwAkiRVkAFAkqQKMgBIklRBBgBJkirIACBJUgUZACRJqiADgCRJFWQAkCSpggwAkiRV\nkAFAkqQKMgBIklRBBgBJkirIANBio5tH6dvQx5LLltC3oY/RzaN5F0mSpL3sk3cBymR08yiD1w+y\nc9dOACZ2TDB4/SAAA6sH8iyaJElTWAPQQkMbh/ac/Cft3LWToY1DOZVIkqSZGQBaaMuOLfPaLklS\nXgwALbSqZ9W8tkuSlBcDQAsNrx1m+bLle23/2S9+ZmdASVKh5BIAIuLsiLgjInZHRP+0fZdExN0R\n8b2IeFnd9lOzbXdHxMWdL3VjA6sHGDl9hBVPWTFl+/afb2fw+kFDgCSpMPKqAbgdOAu4qX5jRBwL\nnAscB5wKXB4RSyNiKfABYB1wLPDK7NjCGVg9wAH7HrDXdjsDSpKKJJdhgCml7wJExPRdZwBXp5Qe\nB/4jIu4GTsr23Z1S+mH2uKuzY+/sTInnx86AkqSiK1ofgMOAH9Xd35ptm217IdkZUJJUdG0LABFx\nQ0TcPsPtjHa9Zva6gxExFhFj27Zta+dLzWqmzoDLly1neO1wLuWRJGm6tjUBpJReuoCH3QscUXf/\n8Gwbc2yf/rojwAhAf39/WkAZFm1y1r+hjUNs2bGFVT2rGF477GyAkqTCKNpUwNcBH4uI/w38KnAM\ncAsQwDERcSS1E/+5wKtyK2UTBlYPeMKXJBVWLgEgIs4E/g44BPhcRNyaUnpZSumOiPgEtc59TwBv\nTCk9mT3mQuBLwFLgipTSHXmUXZKkMoiUcqkl74j+/v40NjaWdzEkSeqYiNiUUupvdFzRRgFIkqQO\nMABIklRBBgBJkirIACBJUgUZACRJqiADgCRJFWQAaNLo5lH6NvSx5LIl9G3oc2lfSVJXK9pMgIU0\nunmUwesH2blrJwATOyYYvH4QwNn+JEldyRqAJgxtHNpz8p+0c9dOhjYO5VQiSZIWxwDQhC07tsxr\nuyRJRWcAaMKqnlXz2i5JUtEZAJowvHaY5cuWT9m2fNlyhtcO51QiSZIWxwDQhIHVA4ycPkJvTy9B\n0NvTy8jpI3YAlCR1LVcDlCSpRFwNUJIkzcoAIElSBRkAJEmqIAOAJEkVZABoM9cQkCQVkWsBtJFr\nCEiSisoagDZyDQFJUlEZANrINQQkSUVlAGgj1xCQJBWVAaCNXENAklRUBoA2cg0BSVJRuRaAJEkl\n4loAkiRpVgYASZIqyAAgSVIFGQAkSaogA0AbOP+/JKnoXAugxZz/X5LUDawBaDHn/5ckdQMDQIs5\n/78kqRsYAFrM+f8lSd3AANBizv8vSeoGBoAWc/5/SVI3cC0ASZJKxLUAJEnSrAwAkiRVkAFAkqQK\nMgBIklRBBgBJkirIACBJUgUZACRJqiADgCRJFWQAkCSpggwAkiRVkAFAkqQKMgBIklRBBgBJkirI\nACBJUgUZACRJqqBIKeVdhraJiG3ARHb3YOChHIsjP4Mi8DPIn59B/sr+GfSmlA5pdFCpA0C9iBhL\nKfXnXY4q8zPIn59B/vwM8udnUGMTgCRJFWQAkCSpgqoUAEbyLoD8DArAzyB/fgb58zOgQn0AJEnS\nL1WpBkCSJGVKGQAi4uyIuCMidkdE/7R9l0TE3RHxvYh4Wd32U7Ntd0fExZ0vdXlFxDsi4t6IuDW7\nnVa3b8bPQ63l9zsfETEeEZuz7/1Ytu0ZEfHliPhB9u/T8y5n2UTEFRHxYETcXrdtxvc9av42+924\nLSKel1/JO6uUAQC4HTgLuKl+Y0QcC5wLHAecClweEUsjYinwAWAdcCzwyuxYtc77UkonZLfPw+yf\nR56FLCO/37l7cfa9n7wYuRjYmFI6BtiY3VdrXUntb0q92d73dcAx2W0Q+GCHypi7UgaAlNJ3U0rf\nm2HXGcDVKaXHU0r/AdwNnJTd7k4p/TCl9Avg6uxYtddsn4day+93sZwBXJX9fBXwihzLUkoppZuA\nn0zbPNv7fgbwkVTzTeCgiHhmZ0qar1IGgDkcBvyo7v7WbNts29U6F2bVa1fUVXn6vneG73N+EvDP\nEbEpIgazbStTSvdlP98PrMynaJUz2/te2d+PffIuwEJFxA3AoTPsGkopfbbT5am6uT4PalVqf0nt\nj+FfAu8FXte50km5OTmldG9E/Arw5Yi4q35nSilFhEOxOsz3vaZrA0BK6aULeNi9wBF19w/PtjHH\ndjWh2c8jIv4e+Kfs7lyfh1rH9zknKaV7s38fjIhrqTXHPBARz0wp3ZdVNT+YayGrY7b3vbK/H1Vr\nArgOODci9ouII6l1+rgF+BZwTEQcGRH7UuuYdl2O5SyVae1pZ1LrpAmzfx5qLb/fOYiIp0bEgZM/\nA79N7bt/HXBedth5gDWWnTHb+34d8NpsNMBvADvqmgpKrWtrAOYSEWcCfwccAnwuIm5NKb0spXRH\nRHwCuBN4AnhjSunJ7DEXAl8ClgJXpJTuyKn4ZfTuiDiBWhPAOPBHAHN9HmqdlNITfr9zsRK4NiKg\n9rf2YymlL0bEt4BPRMTrqa1Wek6OZSyliPg4cApwcERsBS4F3sXM7/vngdOodULeCVzQ8QLnxJkA\nJUmqoKo1AUiSJAwAkiRVkgFAkqQKMgBIklRBBgBJkirIACBp3iLiyWyFu9sj4pMRsTzbfmhEXB0R\n92TT334+Ip6V7ftiRDwSEf8097NL6gQDgKSF+Hm2wt3xwC+A9VEb8H4tcGNK6aiU0onAJfxyzvW/\nAV6TT3ElTWcAkLRYXwOOBl4M7EopfWhyR0rpOymlr2U/bwQezaeIkqYzAEhasIjYh9p66puB44FN\n+ZZIUrMMAJIW4ikRcSswBmwBPpxzeSTNUynXApDUdj9PKZ1QvyEi7gB+P6fySJonawAktcpXgP0i\nYnByQ0Q8JyJ+M8cySZqFAUBSS6TaymJnAi/NhgHeAbwTuB8gIr4GfBJYGxFbI+Jl+ZVWkqsBSpJU\nQdYASJJUQQYASZIqyAAgSVIFGQAkSaogA4AkSRVkAJAkqYIMAJIkVZABQJKkCvr/F0W1CpjAaFMA\nAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "V = V.T\n", "Z = np.matmul(Y, V)\n", "\n", "\n", "# Indices of the principal components to be plotted\n", "i = 0\n", "j = 1\n", "\n", "# Plot PCA of the data\n", "plt.figure(figsize=(8, 8))\n", "plt.title('Country Capabilities: PCA')\n", "plt.plot(Z[:,i], Z[:,j], 'go', label= 'Country')\n", "plt.legend()\n", "plt.xlabel('PC{0}'.format(i+1))\n", "plt.ylabel('PC{0}'.format(j+1))\n", "\n", "# Output result to screen\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "In the following graph the data is projected onto the principal directions selected. " ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true }, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 2", "language": "python", "name": "python2" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 2 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython2", "version": "2.7.13" } }, "nbformat": 4, "nbformat_minor": 2 }