{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# The Efficient Frontier of Philanthropy\n", "\n", "The concept of an [efficient frontier](https://en.wikipedia.org/wiki/Efficient_frontier) for investment portfolios originated with Harry Markowitz's [Modern Portfolio Theory](https://en.wikipedia.org/wiki/Modern_portfolio_theory). Markowitz's main insight was that you can minimize the risk you take for any level of return by diversifying a portfolio. The end result is a hyperbola, often called the Markowitz Bullet, that demonstrates that greater returns are associated with higher risk. Although this approach to porfolio optimization has [theoretical challenges](http://www.bloombergview.com/articles/2015-08-20/more-financial-risk-equals-higher-returns-not-quite-) and is less popular today, modified versions of MPT still underly many robo-advisors currently on the market ([Wealthfront](http://www.slideshare.net/wealthfront/engineering-your-portfolio-with-etfs/35-Want_us_to_do_this), [WiseBanyan](https://wisebanyan.com/investment-strategy), [Betterment](https://www.betterment.com/portfolio/)). \n", "\n", "![frontier image](https://upload.wikimedia.org/wikipedia/commons/e/e1/Markowitz_frontier.jpg) \n", "\n", "The following code implements this concept in Python, and is based on a [blog post](http://blog.quantopian.com/markowitz-portfolio-optimization-2/) from Quantopian. Their original code returns an optimal portfolio for any given level of risk when given a list of historic returns for a group of assets. I then extend the concept of an efficient frontier to the area charitable giving, and refactor the code to take a list of charitable returns (in DALY/\\$) and a covariance matrix as an input. The final output is a philanthropic efficient frontier. \n", "\n", "The following is a rough sketch of my concept, feel free to critique and contribute!\n", "\n", "### Traditional Investment Example: " ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": false }, "outputs": [], "source": [ "#%matplotlib notebook\n", "%matplotlib inline\n", "\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "import cvxopt as opt\n", "from cvxopt import blas, solvers\n", "import pandas as pd\n", "\n", "import mpld3\n", "from mpld3 import plugins\n", "\n", "#mpld3.enable_notebook()\n", "\n", "\n", "np.random.seed(123)\n", "\n", "# Turn off progress printing \n", "solvers.options['show_progress'] = False\n", "\n", "## NUMBER OF ASSETS\n", "n_assets = 4\n", "\n", "## NUMBER OF OBSERVATIONS\n", "n_obs = 1000 #original 1000\n", "\n", "#results in a n_assets x n_obs vector, with a return for each asset in each observed period\n", "return_vec = np.random.randn(n_assets, n_obs) \n", "\n", "\n", "\n", "def convert_portfolios(portfolios):\n", " ''' Takes in a cvxopt matrix of portfolios, returns a list of portfolios '''\n", " port_list = []\n", " for portfolio in portfolios:\n", " temp = np.array(portfolio).T\n", " port_list.append(temp[0].tolist())\n", " \n", " return port_list\n", "\n", "\n", "def optimal_portfolio(returns):\n", " n = len(returns)\n", " returns = np.asmatrix(returns)\n", " \n", " N = 100\n", " mus = [10**(5.0 * t/N - 1.0) for t in range(N)]\n", " \n", " # Convert to cvxopt matrices\n", " S = opt.matrix(np.cov(returns)) \n", " \n", " #S is the covariance matrix, \n", " # diagonal is the variance of each stock:\n", " #[ 1.00e+00 -3.00e-02 -3.19e-02 -2.72e-02]\n", " #[-3.00e-02 9.19e-01 -1.26e-02 -3.74e-02]\n", " #[-3.19e-02 -1.26e-02 9.64e-01 6.39e-03]\n", " #[-2.72e-02 -3.74e-02 6.39e-03 9.79e-01]\n", " \n", " \n", " pbar = opt.matrix(np.mean(returns, axis=1)) #pbar is the mean return of each stock\n", "\n", " # Create constraint matrices\n", " G = -opt.matrix(np.eye(n)) # negative n x n identity matrix\n", " h = opt.matrix(0.0, (n ,1))\n", " A = opt.matrix(1.0, (1, n))\n", " b = opt.matrix(1.0)\n", " \n", " # Calculate efficient frontier weights using quadratic programming\n", " portfolios = [solvers.qp(mu*S, -pbar, G, h, A, b)['x'] #returns a list of objects\n", " for mu in mus]\n", "\n", " port_list = convert_portfolios(portfolios)\n", " \n", " \n", " ## CALCULATE RISKS AND RETURNS FOR FRONTIER\n", " returns = [blas.dot(pbar, x) for x in portfolios]\n", " risks = [np.sqrt(blas.dot(x, S*x)) for x in portfolios] #np.sqrt returns the stdev, not variance\n", " \n", " ## CALCULATE THE 2ND DEGREE POLYNOMIAL OF THE FRONTIER CURVE\n", " m1 = np.polyfit(returns, risks, 2)\n", " #print m1 # result: [ 159.38531535 -3.32476303 0.4910851 ]\n", " x1 = np.sqrt(m1[2] / m1[0])\n", " # CALCULATE THE OPTIMAL PORTFOLIO\n", " wt = solvers.qp(opt.matrix(x1 * S), -pbar, G, h, A, b)['x'] #Is this the tangency portfolio? X1 = slope from origin? \n", "\n", " return np.asarray(wt), returns, risks, port_list\n", "\n", "\n", "weights, returns, risks, portfolios = optimal_portfolio(return_vec)" ] }, { "cell_type": "code", "execution_count": 7, "metadata": { "collapsed": false, "scrolled": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZgAAAEPCAYAAAB/WNKuAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xt8VPWd//HXJ+F+MaKQAFG5iXirRauI626ZLV0L6hZ1\nKwWLFHpZZKV2/bV9QG1dgm5b6XapuiguihcEG62tVbdWsa1TewW84IWLgCYIAUIAESEJhuTz+2NO\nYAgTGMiczCXv5+Mxj8w58z1nPt8c4MP53o65OyIiIqmWl+4AREQkNynBiIhIKJRgREQkFEowIiIS\nCiUYEREJhRKMiIiEIvQEY2ajzGyNma01s+nNlLnbzNaZ2QozOz9uf7mZvWFmr5vZsrj9M81sk5m9\nFrxGhV0PERE5Nu3CPLmZ5QFzgZHAZmC5mT3t7mviyowGBrn7YDO7GJgHDA8+bgAi7v5BgtPPcfc5\nYcYvIiLHL+w7mGHAOnff4O51QCkwpkmZMcBCAHdfChSYWVHwmR0hRgshXhERSZGwE0wxsDFue1Ow\n70hlKuLKOPCimS03s683OW5a0KT2gJkVpDJoERFpuUzv5L/U3S8ALgduNLO/D/bfCwx096HAVkBN\nZSIiGSbUPhhidyOnxW2fEuxrWubURGXcfUvws8rMniLW5PYnd6+KK38/8GyiLzczLbQmInIc3L3F\n3RBh38EsB043s35m1gEYBzzTpMwzwEQAMxsO7HL3SjPrYmbdgv1dgcuAt4Pt3nHHX9O4PxF3z9nX\nzJkz0x6D6qe6qX6590qVUO9g3L3ezKYBS4glswXuvtrMpsQ+9vnu/pyZXW5m64G9wOTg8CLgqeAu\npB2w2N2XBJ/92MyGEhtlVg5MCbMeIiJy7MJuIsPdnweGNNn3v022pyU4rgwY2sw5J6YyRhERSb1M\n7+SXI4hEIukOIVS5XL9crhuofhJjqWxvyzRm5rlcPxGRMJgZngWd/CIi0kYpwYiISCiUYEREJBRK\nMCIiEgolGBERCYUSjIiIhEIJRkREQqEEIyIioVCCERGRUCjBiIhIKJRgREQkFEowIiISCiUYEREJ\nhRKMiIiEIvQEY2ajzGyNma01s+nNlLnbzNaZ2QozOz9uf7mZvWFmr5vZsrj9PcxsiZm9Y2YvmFlB\n2PUQEZFjE+rzYMwsD1gLjAQ2A8uBce6+Jq7MaGCau19hZhcDd7n78OCz94BPufsHTc47G9jh7j8O\nklYPd5+R4Pv1PBgRyXrl5WXMm3crtbUVdOpUzNSpt9O//4DQvi9bngczDFjn7hvcvQ4oBcY0KTMG\nWAjg7kuBAjMrCj6zZmIcAzwSvH8EuCrVgYuIZILy8jJmzvwnIpHFXH11lEhkMTNn/hPl5WXpDu2o\nwk4wxcDGuO1Nwb4jlamIK+PAi2a23My+Hlem0N0rAdx9K1CY0qhFRDLEvHm3Mm7cu3TuHNvu3BnG\njXuXefNuTW9gSWiX7gCO4lJ332JmvYglmtXu/qcE5ZptByspKTnwPhKJ6FnaIpI1amrK2bXrbweS\nS6POnaG2dnPKvicajRKNRlN2vkZhJ5gK4LS47VOCfU3LnJqojLtvCX5WmdlTxJrc/gRUmlmRu1ea\nWW9gW3MBxCcYEZFMV1u7gW3bfk5V1RPU1paRl1dATQ2HJJmaGujUqW/KvrPpf75nzZqVkvOG3US2\nHDjdzPqZWQdgHPBMkzLPABMBzGw4sCtIHF3MrFuwvytwGfB23DGTgvdfBp4OtRYiIiGqrd3Axo3/\nzauvXsyrr15ITc1aBgz4AZdcspnp01+ktHQQNTWxsjU1UFo6iKlTb09v0EkIdRQZxIYpA3cRS2YL\n3P0OM5sCuLvPD8rMBUYBe4HJ7v6amQ0AniLW/NUOWOzudwTlTwKeIHbnswEY6+67Eny3RpGJSEaq\nrX2fqqon2bbtCWpq1tOr19X06nUtJ574j+TltT+k7MFRZJvp1Klv1owiCz3BpJMSjIhkksakUlX1\nc6qr19Gz51UUFl7LiSd+5rCkkk5KMElQghGRdKut3RiXVNbSs+cYCgvHZlxSiacEkwQlGBFJh9ra\nTUFSeYLq6nfo2fMqevW6lh49RmZsUomnBJMEJRgRaS379lUc6FOprl5Dz55j4pJKh3SHd0yUYJKg\nBCMiYYollV8ESWVVkFTGZmVSiacEkwQlGBE5Xs2t/7Vv3+YDfSp7967k5JM/T2HhWHr0+GxWJ5V4\nSjBJUIIRkePRuP5X4xItNTWwePHJjB/fn4KC94Kkcm2QVDqmO9yUU4JJghKMiByP6dMnEIksPmz2\n/AsvjGDOnBdyMqnES1WCyfS1yEREWk19fTU7djzLjh3PJ1z/yz0v55NLKinBiEib5l7Prl1RKisX\nsX37r+jefRhdu55BTc1fQ13/qy1QE5mItEl79rxBZeUiKisfo0OH3hQVTaCwcBwdO/ZJ2AdTWjqI\nWbNeDHWJlkyhPpgkKMGISLza2o1s2/YzKisXsX//hxQVTaCo6Et07Xr2YWVbe/2vTKIEkwQlGBHZ\nv/9Dqqp+QWXlIvbseYNevf6FoqIJFBT8PbGnuktTSjBJUIIRaZsaGj5m587nqaxcxM6dL9Cjx0iK\niiZw0kmXk5/fKd3hZTwlmCQowYi0He7O7t1/pbJyEVVVP6dLl7MoKppAr15foH37k9IdXlZRgkmC\nEoxI7quuXht01i8mL68DRUXXU1h4HZ079093aFlLCSYJSjAi2S/Rki19+3Zl27bHqax8lNra9ykq\nGk9R0QS6dbsAsxb/u9jmZU2CCZ5oeScHn2g5O0GZu4HRxJ5oOcndV8R9lge8Amxy988H+2YCXwe2\nBcVucffnE5xXCUYkiyUaLrxwYWeuvjqfT3ziKoqKJnDiiSPJy9OUvlTKipn8QXKYC4wENgPLzexp\nd18TV2Y0MMjdB5vZxcB9wPC403wTWAWc0OT0c9x9Tpjxi0h6zZv3/QPJBWKz6SdOrOGll8Yyduyj\n6Q1OjirsMXrDgHXuvsHd64BSYEyTMmOAhQDuvhQoMLMiADM7BbgceCDBuXUfLJKjPv64io0b51BV\n9auES7bs21eVnsDkmISdYIqBjXHbm4J9RypTEVfmp8B3gETtXNPMbIWZPWBmBSmKV0TSxL2eHTue\nZ+XKa1m6dDB79rzBiSf+HTU1h5bTki3ZI2MbLs3sCqDS3VeYWYRD71juBW5zdzez/wTmAF9NdJ6S\nkpID7yORCJFIJKyQReQ41NZuYMuWh9i69UHaty+kT5+vMWTIA7RrV8BNNzW3ZMvt6Q47p0SjUaLR\naMrPG2onv5kNB0rcfVSwPQPw+I5+M7sPeMndHw+21wAjiPW9TAD2A52B7sAv3X1ik+/oBzzr7ucl\n+H518otkoIaGfWzf/gxbtjzARx+9QmHhePr0+Srdu59/WNm2vGRLumTFKDIzywfeIdbJvwVYBox3\n99VxZS4HbnT3K4KEdKe7D29ynhHAt+JGkfV2963B+5uBi9z9ugTfrwQjkkH27l3Jli0LqKxcRNeu\n59Knz1fp2fMa8vM7H/1gaTVZMYrM3evNbBqwhIPDlFeb2ZTYxz7f3Z8zs8vNbD2xYcqTkzj1j81s\nKNAAlANTQqqCiLTQ/v17qKp6nC1bHqC2dgO9e0/i/PP/Qpcup6c7NAmZJlqKSIskmgjZr19/du9e\nypYtD7B9+y8oKPg0ffp8jZNOGq05K1kgK5rI0k0JRiRciSZCLlp0Mtde24PevZ0+fb5G795fpmPH\nPukOVY6BEkwSlGBEwtXcs+uXLBnJnDkvatmWLJWqBKOHIYjIcdm//yM+/HBZwomQDQ0NSi6SufNg\nRCQzVVe/Q0XFXCorFwPdqKlBz66XhHQHIyJH5V7P9u3P8sYbl/H6658mP7+ACy98gxkz/kBp6aAD\ns+0bJ0JOnaqJkKI+GBE5grq6nWzZ8iCbN99L+/Y9KS7+Br16XXvIUyE1ETL3qJM/CUowIsdnz543\nqKiYS1XVk5x88pUUF3+DE04Ylu6wpJVkxURLEckeDQ11bN/+FBUVc6mpeY++fW9g2LA1dOhQlO7Q\nJEspwYi0MU0nRn7lKzfRseMLbN58H507n05x8U307DmGvLz26Q5VspyayETakEQTIx9+2PjXfx3L\nxRffQrduh60ZK22Q5sGIyDGbN+97hz0hctIkZ8mSdkouknJqIhNpA+rra6msXMi2bYmfEFlbuzk9\ngUlO0x2MSA6rq9vFhg0/YunSAWzf/jQ9elyqJ0RKq1GCEclB+/ZV8O6732Hp0kFUV6/ivPOWcN55\nv+amm+ZrYqS0GnXyi+SQvXtXs3Hjf7F9+68oKprIqafeTKdO/Q4po4mRcjSaaJkEJRjJRYmev9Kj\nx2bef//H7N79N4qLp1Fc/G+0b39yukOVLJU1CcbMRgF3cvCJlrMTlLkbGE3siZaT3H1F3Gd5wCvA\nprhHJvcAHgf6EXui5Vh3/zDBeZVgJKckGma8cGEnrr22Jxdf/F16955Efn6XdIcpWS4rhikHyWEu\n8DngHGC8mZ3ZpMxoYJC7Dyb26OP7mpzmm8CqJvtmAL919yHA74HvhhC+SMaZN+/Ww4YZT5xYy9Kl\n/0Bx8b8puUhGCbuTfxiwzt03uHsdUAqMaVJmDLAQwN2XAgVmVgRgZqcAlwMPJDjmkeD9I8BV4YQv\nkjnc69m9+7VmhhlvTU9QIkcQdoIpBjbGbW8K9h2pTEVcmZ8C3wGatnMVunslgLtvBQpTFbBIpnFv\noLKylGXLzmH//ioNM5askbETLc3sCqDS3VeYWQQ4Untgsx0tJSUlB95HIhEikUiKIhQJl3sDVVW/\npLy8hPz8bgwefDe33HI6JSWXHdIHU1o6iFmzNMxYjl80GiUajab8vKF28pvZcKDE3UcF2zMAj+/o\nN7P7gJfc/fFgew0wgljfywRgP9AZ6A780t0nmtlqIOLulWbWOzj+rATfr05+yTruzo4dz1BWNhOz\ndgwYcBsnnTT6wCOINcxYwpYVo8jMLB94BxgJbAGWAePdfXVcmcuBG939iiAh3enuw5ucZwTwrbhR\nZLOBne4+28ymAz3cfUaC71eCkYzVdLjxDTfcxgknrKGs7D9wr2PAgNs4+eTP69n20uqy4nkw7l5v\nZtOAJRwcprzazKbEPvb57v6cmV1uZuuJDVOenMSpZwNPmNlXgA3A2LDqIBKGRMONb7nlSb74xVO5\n9NI76NnzamKDMEWylyZaiqTB9OkTiEQWHzIirKYGotHrmD17cfoCEyFL5sGISGLV1eubGW68JT0B\niYRACUakFdXV7WT9+pvZu/d1DTeWnKcEI9IKGhr2sXHjf7Ns2RAaGmqZMeNPWtVYcp76YERC5O5U\nVT3Be+99l65dz2HgwNl07Xo2oOHGkrmyYphyuinBSDp9+OGfWb/+W7h/zKBBP6FHj8+kOySRpGTF\nMGWRtqi6ej3vvTeDjz5axoABP6Co6EsacixtkhKMSIrU1e2gvPx2KisXceqp3+Kssx4lP7/z0Q8U\nyVFKMCItVF9fS0XFXDZunE2vXmMZNmwVHTpo/VURJRiR4+TubNv2OGVl36Vr1/MYOvSPdO165tEP\nFGkjlGBEjsOuXX/k3Xe/jXs9Q4Y8RI8ekXSHJJJxlGBEjkF19dqgA/9VBg78IYWF49WBL9IMJRiR\nJpqucjx16u307dudDRtuo7LyMU477TucddZideCLHIXmwYjESbTK8aJFJ3HVVfUMHXo9/fr9Bx06\n9Ep3mCKh0kTLJCjByLFqbpXj3/3uSn7yk2fTF5hIK9JqyiIhqK2tSLjKcV3d3vQEJJLFQk8wZjbK\nzNaY2drg6ZOJytxtZuvMbIWZDQ32dTSzpWb2upm9ZWYz48rPNLNNZvZa8BoVdj0k99XVfUBdXYVW\nORZJkVATjMWG18wFPgecA4w3szOblBkNDHL3wcAU4D4Ad98H/KO7nw8MBUab2bC4Q+e4+wXB6/kw\n6yG5zb2BLVseZNmys/jCF4ZRWtpfqxyLpEDYo8iGAevcfQOAmZUCY4A1cWXGAAsB3H2pmRWYWZG7\nV7p7dVCmYxBrfIeKHlQuLfbRR6+ydu2NAJx33q/p3v1TDBx46CrHs2ZplWOR4xF2gikGNsZtbyKW\ndI5UpiLYVxncAb0KDALucfflceWmmdn1wCvAt9z9w1QHL7mrrm4nZWXfo6rqKQYO/CG9e086MJ+l\nf/8BzJ69KM0RimS/jO7kd/eGoInsFOBiMzs7+OheYKC7DwW2AnPSFaNkF/cGNm9+gGXLzgbyGTZs\nNX36fEWTJUVCEPYdTAVwWtz2KcG+pmVOPVIZd99tZi8Bo4BV7l4V9/H9QLPjR0tKSg68j0QiRCKR\n5KOXnLJ79yusW3cjZvmcd95v6N79/HSHJJIRotEo0Wg05edNeh6Mmf0d0J+4pOTuC49yTD7wDjAS\n2AIsA8a7++q4MpcDN7r7FWY2HLjT3YebWU+gzt0/NLPOwAvAHe7+nJn1dvetwfE3Axe5+3UJvl/z\nYIS6uh2899732L79Vwwc+CN69/6y7lhEjqBVHzhmZo8S6wdZAdQHu52gc7457l5vZtOAJcSa4xa4\n+2ozmxL72OcHCeNyM1sP7AUmB4f3AR4J+mHygMfd/bngsx8Hw5kbgHJio89EDhEbHbaAsrLvU1g4\nlmHDVtO+fY90hyXSZiR1B2Nmq4Gzs+12QHcwbdfu3cuD5rD2DB58D927D013SCJZo7Ufmfw20JtY\nM5dIxoo1h93Cjh3PMHDgHRQVXa/mMJE0STbB9ARWmdkyYF/jTnf/fChRiSTh0FWP+3DNNWezf///\nUFj4RS66aDXt25+Y7hBF2rRkm8hGJNrv7n9IeUQppCay3JVo1eOFCzvxve89ybnnXpHu8ESymlZT\nToISTO5qbtXjaPRLmiQp0kKtupqymQ03s+VmtsfMPjazejPb3dIvFzleH330ZsJVj2trN6cnIBE5\nTLK9n3OB8cA6oDPwNeCesIISac7HH1exatV49u0r06rHIhku6eE17r4eyHf3end/iNisepFW4e5U\nVj7G8uWfoGPHU7jllqWUlg7SqsciGSzZTv6Xgc8CDxBb+2sLMMndPxlueC2jPpjcUFu7ibVrb2Df\nvvcZMmQBJ5xwERA/iiy26vHUqVr1WCQVWrWT38z6AZVAB+BmoAC4N7iryVhKMNktNhP/fsrKvk9x\n8Tc47bQZ5OV1SHdYIjmv1UeRBeuBnebu77T0S1uLEkz2qq5ez9q1X6e+vpohQxbQrdu56Q5JpM1o\n7VFk/0xsHbLng+2hZvZMS79cpCn3et5//ye89tpwTj7581xwwV+UXESyVLIz+UuIPSgsCuDuK8xM\njd2SUnv2vM0773yF/PxufOpTS+nceVC6QxKRFkg2wTQumx+/T21PkhINDR+zYcMP2bz5HgYM+CF9\n+nyNJn/WRCQLJZtgVprZdUC+mQ0GbgL+El5YkqsOXT+smAkT/oWamlvp3HkgF164go4di9Mdooik\nSLKjyLoA3wMuC3a9ANzu7vuaPyr91MmfWRKtH/bww/l8+9v/zYUX3qS7FpEM0aqd/MDZwasd0AkY\nAyxv6ZdL2zJv3q0HkgvElnaZNKmeJ59cruQikoOSTTCLgQeBa4Arg9c/J3OgmY0yszVmttbMpjdT\n5m4zW2dmK4InVWJmHc1sqZm9bmZvmdnMuPI9zGyJmb1jZi+YWUGS9ZA0qqnZpPXDRNqQZBNMlbs/\n6+5l7r6h8XW0g4LHHc8FPgecA4w3szOblBkNDHL3wcQefXwfQND89o/ufj4wFBhtZsOCw2YAv3X3\nIcDvge8mWQ9Jk5qa96ipeVvrh4m0IckmmJlm9oCZjTezaxpfSRw3DFgXJKQ6oJRY81q8McBCAHdf\nChSYWVGwXR2U6Uisec7jjnkkeP8IcFWS9ZBW5u5s3foIr712MZMn30Bp6UCtHybSRiQ7imwycCbQ\nHmgI9jnwy6McVwxsjNveRCzpHKlMRbCvMrgDehUYBNzj7o39PoXuXgng7lvNrDDJekgrqqv7gLVr\nb2Dv3pV88pO/o1u38+jb96uHrB82a5bWDxPJVckmmIuC5qhW5e4NwPlmdgLwKzM7291XJSrayqHJ\nUXzwQZQ1aybSs+dVnHnmw+Tnxzpf+vcfoAeCibQRySaYvxzhH/cjqQBOi9s+JdjXtMypRyrj7rvN\n7CVijwhYRezupsjdK82sN7CtuQBKSkoOvI9EIkQikWOsghyLhoaPKSv7DyorFzJkyAJOPnl0ukMS\nkaOIRqNEo9GUnzfZeTCriTVTlQH7AAPc3c87ynH5wDvASGJL/C8Dxrv76rgylwM3uvsVZjYcuNPd\nh5tZTw6uINCZ2NybO9z9OTObDex099nByLQe7j4jwfdrHkwrqq5+h1WrrqNjx74MGbKADh3UcimS\njVI1DybZO5jjeriYu9eb2TRgCbEBBQvcfbWZTYl97PODhHG5ma0H9hLr7wHoAzwS9MPkAY+7+3PB\nZ7OBJ8zsK8AGYOzxxCfHp+ls/BtuuI2OHX/Le+/dwoABt9O37w2a1yIiyS/Xn410B5N6iWbjL1zY\nhfHj+zFy5C/o2vWsdIcoIi3U2jP5RYDEs/EnTqzmz3/+pJKLiBxCCUaOSW1tRTOz8SvTE5CIZCwl\nGDkmHTqcrNn4IpIUJRhJ2u7dy7nkkr/w6KMnaja+iByVOvnlqNydLVvmU1b2fc4443/Zu/f8Q2bj\nT52q2fgiuSRVnfxKMHJE9fU1rFv3b+zevZxzz/0FXbq0+oIOItLKNIpMQldT8x6vv/53NDTs44IL\n/qbkIiLHRAlGEtqx49e89tol9O49mbPOWky7dt3SHZKIZJlkZ/JLDjt0Zn4frryyJ/n5v+Tcc39J\nQcGl6Q5PRLKUEkwbl2hm/r33dmLWrJcpKLgo3eGJSBZTE1kbl3hmfi0LFtyV3sBEJOspwbRxzc/M\n35yegEQkZyjBtHHuezQzX0RCoQTTRrnXs27dv/OZz+zkZz87TTPzRSTlNNGyDdq/fw+rV19Hff1e\nzjnnSSoqdmlmvogcoJn8SVCCOdy+fZt5660r6dbtfM44Yx55eR3SHZKIZJjWfqLlcTOzUcCdHHyi\n5ewEZe4GRhN7ouUkd19hZqcAC4EioAG4393vDsrPBL4ObAtOcYu7Px92XbJN0ydPTpw4nt27p9K3\n71ROO22GnjopIqEK9Q4meNzxWmAksBlYDoxz9zVxZUYD09z9CjO7GLjL3YebWW+gd5BsugGvAmPc\nfU2QYD5y9zlH+f42eweTaH7Lww/n8Z3v3MmFF34j3eGJSAbLlrXIhgHr3H2Du9cBpcCYJmXGELtT\nwd2XAgVmVuTuW919RbB/D7AaKI47Tv/9PoJE81smTWrg5z9fmt7ARKTNCDvBFAMb47Y3cWiSSFSm\nomkZM+sPDAXi/3WcZmYrzOwBMytIVcC5QvNbRCTdMn6pmKB57Engm8GdDMC9wG3u7mb2n8Ac4KuJ\nji8pKTnwPhKJEIlEQo03U3Tq1JeaGg5JMprfIiKJRKNRotFoys8bdh/McKDE3UcF2zMAj+/oN7P7\ngJfc/fFgew0wwt0rzawd8H/Ab9w94dolZtYPeNbdz0vwWZvsg3Gv53e/u44FC55m4sR9B/pgSksH\nMWvWixqCLCJHlC2jyJYDpwdJYAswDhjfpMwzwI3A40FC2uXulcFnDwKrmiYXM+vt7luDzWuAt8Oq\nQLZpaPiY1auvp2fP7dx++yvcf/8dB+a3zJql+S0i0npCnwcTDFO+i4PDlO8wsynE7mTmB2XmAqM4\nOEz5dTO7FHgZeAvw4HWLuz9vZguJ9ck0AOXAlLikFP/dbeoOpr6+hpUrr8Usj7PPfoL8/E7pDklE\nspAmWiYh1xNM/DyXjh0LGTHifQYM6M+ZZy4kL699usMTkSyVLU1kEpJE81wefbQ7P/jBo0ouIpIR\ntNhllko0z+X66z/ivvtK0hqXiEgjJZgspXkuIpLplGCyVMeORXqOi4hkNCWYLFRfX0skspFHH+2u\n57iISMbSKLIs09BQx8qV/0JeXme6dPkB991Xoue4iEhKaZhyEnIlwRwcjryJffvKuOqqQVx22fN6\nlouIhELDlNuIRMORf/azdpx5ZoXuVkQko6kPJsMlGo48btx7zJt3a3oDExE5CiWYDKfhyCKSrZRg\nMlxenmk4sohkJSWYDLZnz1tceukKHnusWMORRSTraBRZhmkcMVZdXUZ19QpuuOGH9Or1+WAUmYYj\ni0j4NEw5CdmWYBKNGNNDwkSktaUqwaiJLIMkHjH2rkaMiUhWUoLJIBoxJiK5JPQEY2ajzGyNma01\ns+nNlLnbzNaZ2QozGxrsO8XMfm9mK83sLTO7Ka58DzNbYmbvmNkLZlYQdj1aQ35+e40YE5GcEWqC\nMbM8YC7wOeAcYLyZndmkzGhgkLsPBqYA9wUf7Qf+n7ufA1wC3Bh37Azgt+4+BPg98N0w6xGm8vIy\npk+fwDe+cSlbt77M/Pk9NWJMRHJC2EvFDAPWufsGADMrBcYAa+LKjAEWArj7UjMrMLMid98KbA32\n7zGz1UBxcOwYYERw/CNAlFjSySqJOvUXLSri6afH0LHjbjp16susWRoxJiLZKewEUwxsjNveRCzp\nHKlMRbCvsnGHmfUHhgJ/C3YVunslgLtvNbPClEbdShJ16k+Y8D7R6D8we/av0huciEgLZfxil2bW\nDXgS+Ka7722mWLNjkUtKSg68j0QiRCKRVIbXIurUF5FMEI1GiUajKT9v2AmmAjgtbvuUYF/TMqcm\nKmNm7Ygll0fd/em4MpVBM1qlmfUGtjUXQHyCyTQdOxZSU8MhSUad+iLS2pr+53vWrFkpOW/Yo8iW\nA6ebWT8z6wCMA55pUuYZYCKAmQ0HdjU2fwEPAqvc/a4Ex0wK3n8ZeJos0tixv3XrS/zgB+0oL4/t\nV6e+iOSS0Gfym9ko4C5iyWyBu99hZlMAd/f5QZm5wChgLzDJ3V83s0uBl4G3iDWBOXCLuz9vZicB\nTxC789kAjHX3XQm+O+Nm8ifq2L/nnm706nUuhYWDtAyMiKSdlopJQiYmmOnTJxCJLD6sWSwa/RKz\nZy9KX2DO8HMJAAALRElEQVQiIgEtFZOl1LEvIm2FEkwraex3Wb/+DRYsgK1bD36mjn0RyUVqImsF\nifpdHnoIrrkGCgq0YrKIZBb1wSQhUxJMc/0uP/lJERde+Fl17ItIRklVgsn4iZa5oLl+l3PPPVsd\n+yKSs5RgQlZeXsbq1eXs3Ant2sHnPge9e6vfRURynxJMiBr7Xm6+ufxA38vDD8Po0fDSS4OYNUsT\nKkUkd6kPJkTN9b389Kf9mT//9+p3EZGMpD6YDFdeXsYrr7zI1q2Ql3ewaaxzZzjrrAFKLiKS85Rg\nQtDYNPbtb287pGns6qtjw5LV9yIibYGayELQXNPYY49BXp7mvIhIZlMTWQbbtm09TzwBDQ2HNo99\n8EERc+cquYhI26AEk2Ll5WVUVa3kxhs5bOTYOed8VslFRNoMNZGlWHPNYz/6UTcWLnxTCUZEMp6a\nyDJUc81jZ511rpKLiLQpSjApdKTmsRNPHJTu8EREWlXoy/Wb2SgzW2Nma81sejNl7jazdWa2wszO\nj9u/wMwqzezNJuVnmtkmM3steI0Kux7JuOOOm+nRYw+lpfDII/DhhzBpEjz4YDc9BllE2pxQ72DM\nLA+YC4wENgPLzexpd18TV2Y0MMjdB5vZxcA8YHjw8UPA/wALE5x+jrvPCTP+Y1FeXsamTS8cdvdy\n9dVqHhORtinsO5hhwDp33+DudUApMKZJmTEECcTdlwIFZlYUbP8J+KCZc7e4AyqVYncvtYfdvfz6\n12oeE5G2KewEUwxsjNveFOw7UpmKBGUSmRY0qT1gZgUtC7NlGu9errsOJk+GsWPhqadiSaayspOa\nx0SkTcrWTv57gdvc3c3sP4E5wFcTFSwpKTnwPhKJEIlEUh7MvHm38oUv1B4yemz06NjdyxlnXKbm\nMRHJaNFolGg0mvLzhjoPxsyGAyXuPirYngG4u8+OK3Mf8JK7Px5srwFGuHtlsN0PeNbdz2vmO5r9\nvLXmwUyePJzq6qVMmnRo/8vmzR149NE1SjAiklVSNQ8m7Cay5cDpZtbPzDoA44BnmpR5BpgIBxLS\nrsbkEjCa9LeYWe+4zWuAt1Md+LFYv34THTpwWP9L+/a9lFxEpM0KNcG4ez0wDVgCrARK3X21mU0x\ns38NyjwHlJnZeuB/gX9rPN7MHgP+ApxhZu+b2eTgox+b2ZtmtgIYAdwcZj2OpLy8jHbtth3Y3r8f\nFi2KJZnBg09NV1giImmnpWJa6IYbrmLnzqeZPPlg89hDD0FdHZxxxpeYPXtRqN8vIpJqqWoiU4Jp\noZEji7j++m384Q8HO/hHjIB77snj5z9fryYyEck62dIHk/Nqamq5//5Y01hjcvnNb6Br1+5KLiLS\npinBtEB5eRl5eR8xcGAsuezfD08+GRuinJfXLd3hiYiklRJMC9xxx80UFBzaBOcem/+iDn4Raeuy\ndaJlRnjzzT/SqROUlx/s4G/fHqqqYOBALQ8jIm2bEkwL7N79Ie3bQ//+sSayhgZ4//3YEGUtDyMi\nbZ0STAs0NNSTn3/4HUx+PurgF5E2TwmmBRoaYNeuWL9Lfj7s2QNmsf0iIm2dOvlbYN++w5NJQ0Ns\nv4hIW6c7mBZoaIDdu2PNYw0N8PHHsWYy3cGIiOgOpkVqamDz5kP3bd4c2y8i0tbpDqYFdu+GT38a\nbr75YCf/T38KL7+c7shERNJPCaYFBg06mFwg9vPmm6GiIr1xiYhkAjWRtcBJJx1MLo06d47tFxFp\n65RgWmDnzsP7W2pqYvtFRNq60BOMmY0yszVmttbMpjdT5m4zW2dmK8zs/Lj9C8ys0szebFK+h5kt\nMbN3zOwFMysIux6JvPturM+lMck09sG8+246ohERySyhJhgzywPmAp8DzgHGm9mZTcqMBga5+2Bg\nCjAv7uOHgmObmgH81t2HAL8HvhtC+EdVW+u8/DJ8+9tw662xny+/HNvfGqLRaKt8T7rkcv1yuW6g\n+klM2Hcww4B17r7B3euAUmBMkzJjgIUA7r4UKDCzomD7T8AHCc47BngkeP8IcFUIsSelttZZudL5\n4x9jP1sruUDu/yHP5frlct1A9ZOYsBNMMbAxbntTsO9IZSoSlGmq0N0rAdx9K1DYwjhFRCTFcqWT\nP3ef+ywikq3cPbQXMBx4Pm57BjC9SZn7gC/Gba8BiuK2+wFvNjlmdWMZoDewupnvd7300ksvvY79\nlYocEPZEy+XA6WbWD9gCjAPGNynzDHAj8LiZDQd2NTZ/BSx4NT1mEjAb+DLwdKIvd/emx4mISCux\n4H/64X2B2SjgLmLNcQvc/Q4zm0IsQ84PyswFRgF7gcnu/lqw/zEgApwMVAIz3f0hMzsJeAI4FdgA\njHX3XaFWREREjknoCUZERNqmrOzkT2byZlDuIjOrM7Nr4vaVm9kbZva6mS1rnYiPzdHqZ2YjzGyX\nmb0WvL6f7LGZoIX1y/rrF5SJBHV428xeOpZj062F9cvo65fEn81vB7G/ZmZvmdl+MzsxmWMzQQvr\nd+zXLsxO/pAGDuQB64l1/rcHVgBnNlPud8D/AdfE7X8P6JHuerSkfsAI4Jnj/d1ka/1y6PoVACuB\n4mC7Z45dv4T1y/Trd6y/f+BKYhO+c+baNVe/47122XgHk8zkTYBvAE8C25rsNzL7zi3Z+iUawJDs\nsenUkvo17s/263cd8At3rwBw9+3HcGy6taR+kNnX71h//+OBnx3nsenQkvrBcVy7TL3QR3LUyZtm\n1he4yt3ncfg/VA68aGbLzezroUZ6fJKZnApwSbB226/N7OxjPDadWlI/yI3rdwZwkpm9FNTj+mM4\nNt1aUj/I7OuX9O/fzDoTG5j0i2M9No1aUj84jmuXq8+DuROIb1+MTzKXuvsWM+tF7Je12mNL0mST\nV4HT3L06WMvtV8T+UueKI9UvF65fO+AC4DNAV+CvZvbX9IaUUgnr5+7ryY3rB/DPwJ88d0evJqrf\nMV+7bLyDqQBOi9s+JdgX70Kg1MzKgC8A95jZ5wHcfUvwswp4ithtYyY5av3cfY+7VwfvfwO0D4Zu\nJ/O7SbeW1C8nrh+x/zm+4O617r4DeBn4ZJLHpltL6pfp1+9Yfv/jOLT5KFeuXaOm9Tu+a5fujqfj\n6KjK52BHVQdiHVVnHaH8QwSd/EAXoFvwvivwZ+CydNfpWOvHoSsdDAPKj+d3k4X1y5XrdybwYlC2\nC/AWcHYOXb/m6pfR1y/Z3z+xQQw7gM7HemwW1++4rl3WNZG5e72ZTQOWcHDy5mprMnkz/pC490XA\nU2bmxG7jF7v7klYJPElJ1u8LZjYVqANqgC8e6di0VKQZLakfOXL93H2Nmb0AvAnUA/PdfRVALly/\n5upnZgPI4Ot3DP+2XEXsDq3maMe2chWOqCX14zj/7mmipYiIhCIb+2BERCQLKMGIiEgolGBERCQU\nSjAiIhIKJRgREQmFEoyIiIRCCUakFZnZN82sUzOffdnM/qe1YxIJixKMSOv6d2KzopujiWmSM7Ju\nJr9ItjCzLsQe7V1MbJmOJ4G+wEtmtt3dR5rZZGAG8AGxme+16YpXJNWUYETCMwqocPcrAczsBGAS\nEHH3D8ysN1ACnA/sBqLAa2mJVCQEaiITCc9bwD+Z2Y/M7O/dfTexR0c0Pj7iYuAld9/p7vuBx9MV\nqEgYdAcjEhJ3X2dmFwCXA7eb2e85vI+luSd3imQ93cGIhMTM+gA17v4Y8BNiD+H6CDghKLIU+LSZ\n9TCz9sC16YlUJBy6gxEJzyeA/zKzBuBjYCpwCfC8mVUEnfyzgL8R6+Rfkb5QRVJPy/WLiEgo1EQm\nIiKhUIIREZFQKMGIiEgolGBERCQUSjAiIhIKJRgREQmFEoyIiIRCCUZERELx/wHvR53EiZb4EQAA\nAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/html": [ "\n", "\n", "\n", "\n", "
\n", "" ], "text/plain": [ "" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# Matplotlib only plot: \n", "plt.ylabel('mean')\n", "plt.xlabel('std')\n", "plt.plot(risks, returns, 'y-o') #risks and returns are just arrays of points along the frontier\n", "plt.show()\n", "\n", "\n", "# Interactive plot using mpld3:\n", "def interactive_plot(risks, returns, portfolios):\n", " fig, ax = plt.subplots()\n", " ax.grid(True, alpha=0.3)\n", "\n", "\n", " labels = []\n", " for i in range(len(risks)):\n", " rounds = map(lambda x: round(x,3), portfolios[i])\n", " label = \" Risk: \" + str(risks[i]) + \"
\" + \" Return: \" + str(returns[i]) + \"
\" + \" Portfolio Weights: \" + str(rounds)\n", " labels.append(str(label))\n", "\n", " points = ax.plot(risks, returns, 'y-o', color='b',\n", " mec='k', ms=15, mew=1, alpha=.6)\n", "\n", " ax.set_xlabel('standard deviation')\n", " ax.set_ylabel('return')\n", " ax.set_title('Efficient Frontier', size=20)\n", "\n", " tooltip = plugins.PointHTMLTooltip(points[0], labels,\n", " voffset=10, hoffset=10)\n", " plugins.connect(fig, tooltip)\n", "\n", " return mpld3.display()\n", "\n", "#Note: if you're in the IPython notebook, mpld3.show() is not the best command\n", "# to use. Consider using mpld3.display(), or mpld3.enable_notebook().\n", "# See more information at http://mpld3.github.io/quickstart.html.\n", " \n", "interactive_plot(risks, returns, portfolios)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Application to Philanthropy\n", "\n", "An interesting next question is, does this same type of risk return relationship apply to philanthropy? Are areas with highly uncertain yet historically high returns like scientific research the best place to put philanthropic money? I outline a framework that might help answer this question below. \n", "\n", "The above code is great if you have a well curated list of historical returns for different stocks, but the same isn't true for the world of philanthropy. Very few charities or philanthropic endeavors have good data on how their interventions impact the world. And even if a charity does have good data on their impact, it is still very hard to compare effectiveness *across* interventions. \n", "\n", "One method that healthcare economists have developed to compare across interventions is the [Disability Adjusted Life Year](https://en.wikipedia.org/wiki/Disability-adjusted_life_year), or DALY. A measure that determines the cost effectiviness of an intervention is the DALY/\\$, which can be though of as the \"return\" to charitable giving. Economists will often perform sensitivity analysis when calculating DALY/\\$ values, which result in a variance or standard deviation for their estimates. The DALY/\\$ metric, combined with the variance from the sensitivity analyses give us all the information we need to create an efficient frontier that combines different interventions. \n", "\n", "Next, I refactor the above code to take covariance and mean return matrices as inputs instead of historical returns. This would allow anyone who has done a cost effectiveness analysis along with a sensitivity analysis to input their values as lists, and output an efficient frontier. " ] }, { "cell_type": "code", "execution_count": 8, "metadata": { "collapsed": false }, "outputs": [], "source": [ "def covmean_portfolio(covariances, mean_returns):\n", " ''' Returns an optimal portfolio given a covariance matrix and matrix of mean returns '''\n", " n = len(mean_returns)\n", " \n", " N = 100\n", " mus = [10**(5.0 * t/N - 1.0) for t in range(N)]\n", "\n", "\n", " S = opt.matrix(covariances) # convert array to matrix \n", "\n", " pbar = opt.matrix(mean_returns) # convert array to matrix\n", "\n", " # Create constraint matrices\n", " G = -opt.matrix(np.eye(n)) # negative n x n identity matrix\n", " h = opt.matrix(0.0, (n ,1))\n", " A = opt.matrix(1.0, (1, n))\n", " b = opt.matrix(1.0)\n", " \n", " # Calculate efficient frontier weights using quadratic programming\n", " portfolios = [solvers.qp(mu*S, -pbar, G, h, A, b)['x'] \n", " for mu in mus]\n", " \n", " #Convert portfolios to list for output\n", " port_list = convert_portfolios(portfolios)\n", " \n", " ## CALCULATE RISKS AND RETURNS FOR FRONTIER\n", " frontier_returns = [blas.dot(pbar, x) for x in portfolios] #differ from those input, renamed\n", " risks = [np.sqrt(blas.dot(x, S*x)) for x in portfolios] #np.sqrt returns the stdev, not variance\n", " \n", " ## CALCULATE THE 2ND DEGREE POLYNOMIAL OF THE FRONTIER CURVE\n", " m1 = np.polyfit(frontier_returns, risks, 2)\n", " x1 = np.sqrt(m1[2] / m1[0])\n", " \n", " # CALCULATE THE OPTIMAL PORTFOLIO\n", " wt = solvers.qp(opt.matrix(x1 * S), -pbar, G, h, A, b)['x'] \n", "\n", " return np.asarray(wt), frontier_returns, risks, port_list\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Basic Application\n", "\n", "Next, I create a simple list of variances and covariances, assuming there are no covariances between charities. The covariances is a list of lists, with each sub-list representing a column in the covariance matrix. The diagonal of this matrix represents the variances for individual interventions. The mean returns list just contains the estimate for the mean impact in DALY/\\$1000. In the future, it might be possible to incorporate positive externalities, or flow-through effects in this matrix by including them as covariances between interventions (although covariances aren't causative or directional). In this case, I assumed each intervention has a variance of 1, and chose a few mean returns. " ] }, { "cell_type": "code", "execution_count": 9, "metadata": { "collapsed": false, "scrolled": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYgAAAEPCAYAAABY9lNGAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xt8VeWV//HPCgkh3JNwFdAAggXRolK11f5IqcUibb0M\n49iqiDNt0dFfHaqWmbZWaK0OHVurdWrrbxwV1AJt1TqjVnEkbbUdBRGwgmiQIAQhQMI1J5DL+v1x\nDiGEHTkkOWefy/f9euXVk3N29lnuJmfxPOt51jZ3R0REpLWcsAMQEZHUpAQhIiKBlCBERCSQEoSI\niARSghARkUBKECIiEijhCcLMKsxslZm9aWavt3HMfWb2npmtNLPxiY5JRESOLTcJ79EElLp7TdCL\nZjYFGOnuo8zsHOAXwLlJiEtERD5CMqaY7BjvczEwH8DdXwP6mNnAJMQlIiIfIRkJwoElZrbMzL4W\n8PoQYFOL7ytjz4mISIiSMcV0nrt/aGb9iSaKte7+ShLeV0REOiDhCcLdP4z973Yzewo4G2iZICqB\nYS2+Hxp77ghmpqZRIiLt4O7Wnp9L6BSTmXU3s56xxz2AycBfWx32DDA9dsy5wC533xZ0PnfXlzu3\n33576DGkypeuha5FNl6L+vpdVFU9ybp11/GXv4zg1VdPYO3aGWzd+gQHDlQxdiw89xwsXdqxz/BE\njyAGAk/F/vWfCzzu7i+a2UzA3f1Bd3/OzC4ys3JgP3BtgmMSEUkr7o3s2bOMmpoXqa5+gf37V9O7\n96coKrqQceNuoEePUzE7PEioqIB77oFZszr2vglNEO6+AThqX4O7/7LV9zcmMg4RkXRTV7eR6uoX\nqal5kZqa/yE/fyhFRRdSUjKHPn3Op0uXgjZ/dv9+p0cPo/Koyfrjk4witXSy0tLSsENIGboWh+la\nHJaO16KhYR+7dpU1jxIaGmooLJxMcfEXOfnk+8jPH3xc59u/P1q2bTmyOF7mnh61XzPzdIlVRORY\n3JvYt+/N5lHC3r3L6dXrbIqKJlNYOJmePT+OWcfLxGaGt7NIrQQhItJJKio28MADt1FXV0m3bkO4\n/vofUFIyvPn1Awcqqa5eEps2WkJeXj8KCy+kqGgyfftOpEuXHp0ekxKEiEjIKio2cPvtn+OKK9ZT\nUACRCCxcOJxZs75Hz55vUVPzIgcObKGw8ILYKOFzdOt2YsLjUoIQEQnZ7NlXUVr6OAUtaseRCDz1\nVH++/e1vUFQ0mV69zsKsS1Lj6kiCUJFaRKQDGhsj7Nr1MjU1S49IDgAFBVBQMI6Sku+GE1wHKUGI\niBynurqN7Nz5LDt3Psvu3X+iZ88zyM8fQCSy5agRRLduJ4QXaAdpiklE5BiamurZs+fPsaTwHPX1\nVRQVTaG4eCqFhZ8jL6+wjRrESObOXXJEoTrZVIMQEelkBw9WUV39e3bufJaamiV06zaC4uKLKC6e\nSq9enwhcgnp4FdMWunU74ahVTGFQghAR6aBD+xIOTR3V1q6jsPCzFBdPpahoynFvVEsVShAiIu3Q\n0LCHmpol7Nz5LNXVz9OlSx+Ki6dSXHwRffp8mpycrmGH2GFKECIirQRtWjvppBIikXebRwl7975O\n796fiiWFqRQUjAw77E6nBCEi0kJQwXjBgt5Mm9aHQYOamkcJfft+ltzcnmGHm1BKECIiLdx66zQm\nTfrtUUtOX3ppCnff/WyHGtilG22UE5GsF4lUsGPH0+zY8TTbt78SuGmtoaEuq5JDRylBiEhacnf2\n71/dnBQOHKikuPiLDBt2CwMGDCISWZRRm9bCoCkmEUkb7o3s3v1qc1IA6NfvUvr1u4Q+fT7V3Oco\nVTethUE1CBHJWI2NEWpqXmLHjqfZufO/yM8f0pwUevQ4rc0po1TctBYGJQgRySj19TXs3PksO3Y8\nTU3NEnr1OpN+/S6huPhiCgpKwg4vrShBiEhaCdqjMGhQHjt3/o7t259i797X6dt3UiwpfIGuXfuF\nHXLaUoIQkbQRVB+YPz+fyy7LZ9y4S+jX7xKKiiYn5O5q2UgJQkTSxs03f4ELLnj2qBVGZWVfZt68\nJ8ILLENpH4SIpLTa2nVUVS2mqmoR1dXvBe5RqKvbGk5w0qaj+9WKiHSCSGQ9GzfeybJl41m58jPU\n1+/glFN+yYAB04hEWh+rPQqpSFNMItJpIpEKtm9fTFXVYg4c2ET//n/DgAF/R58+52uPQkhUgxCR\nhApadXTow7yubhPbt/+aqqrF1NWtp1+/S2NJYSI5OcGz2NqjkDxKECKSMMH/4j+J666bTn7+/1Bb\n+w79+l3CgAGX07fvJHJy8sIOWVpQghCRhJk9+ypKSx8/atXRM88M5447fkZh4ecy4sY6mUqrmEQk\nIRobI+zZ82bgqqOuXUsoLp4aTmCSFEoQInIEd2f37lfYtm0+27f/lqamfCIR1Bk1CylBiAgQXZa6\ndet8tm1bQE5OAYMGXcOECasZOrS+jVVHPwg7ZEkw1SBEslh9/S62b1/M1q3ziUTeZcCALzNo0HR6\n9jzziC6pWnWUvlSkFpE2tV6iOnPm9+jd+122bZtPdfULFBVNZuDA6RQVfV4rkDKQEoSIBApaovrI\nIzlcffV4zjzz6/Tvfzl5eYVhhykJ1JEEoVYbIhns5z+f3ZwcIFponjGjiT/9aQwnnDBTyUE+korU\nIhnG3dm793W2bPklVVVPttEYb0s4wUlaScoIwsxyzGyFmT0T8NpEM9sVe32FmX03GTGJZJqGht1U\nVv6c5cvHs2bNV+je/RT6979YjfGk3ZI1grgJWAP0buP1P7r7l5IUi0haCuqHdNJJJezdu4wtW37J\njh1PUlh4ASNH/pjCwkmY5XDDDRu4/fZVWqIq7ZLwIrWZDQUeBn4IfLN1IjCzicAt7v7FY5xHRWrJ\nWkHF5scf78fll/dnwIADDB78NQYPvpauXQcG/qyWqGavlF7FZGa/Jpoc+gA3t5EgfgtsBiqBW919\nTcB5lCAka7XVD+nFFyfxk58swUzrTSRYyvZiMrOpwDZ3X2lmpUBQkG8AJ7p7rZlNAZ4GRgedb86c\nOc2PS0tLKS0t7eyQRVKOeyN7964OLDY3NbmSgxyhrKyMsrKyTjlXQkcQZnYncBXQABQAvYAn3X36\nR/zMBuAsd69u9bxGEJJVGhr2sXXrw2zefC8PPVTDZZdVB9zH+UrmzXssvCAl5aX0FFPzG0WnkoKm\nmAa6+7bY47OBxe5eEvDzShCScYIKz4MG5bB588/YuvVhCgsnMXToP1FdPZg5cybrLmxy3NIuQZjZ\nTMDd/UEzuwG4HqgHIsAsd38t4OeVICSjBBWeFyzowaWXduGMM77GkCE3UlBQcsTxKjbL8UqLBNFR\nShCSadoqPC9dejk/+tGi8AKTjKJWGyJppqmpgT17VgQWng8c2B5OUCKtqNWGSIIE1RdOPHEI27Yt\nYOPGu2hs3Ksb8UhK0xSTSAK0tbHtsstyGTnydE466Tvs2jWsjRvxqPAsnUc1CJEU01Z9YcmSC/nJ\nT37f/JwKz5JoKbtRTiRbRSKbAusLjY0Hj3iupGS49jFIylKRWqQTNTU18OGHD7F//3J1UZW0pwQh\n0gncm6iqWsyyZaeybdtj3HTTfBYuHNmcJA7VF66/Xl1UJX2oBiHSAe5OdfXv2bDhO0AOI0bcRWHh\nBZiZ6guSElSkFgnBrl2vsGHDt6mv387w4XfQr99lmLXr71AkYVSkFkmivXtXsmHDd9i//21KSuYw\ncOBV5OToT0kyj36rReJUW/suGzZ8j927/8CJJ36bceOeJCcnP+ywRBJGCULkGOrqNrFx4/fZvv0p\nhg37Jqec8h/k5vYMOyyRhFOCEGnDwYPb+eCDu9i69VFOOOHrnHPOu+TlFYUdlkjSKEFIVgvqlzR0\naDGbNv2Yysr7GTDgCj7xib+Snz847FBFkk6rmCRrBfVLeuyxYi691Dn11KmUlMyhoGBE2GGKdIiW\nuYq0Q1v9kl56aSo//vF/hxeYSCfS/SBE2qGurjKwX1JDQ204AYmkGCUIyUr79r1FXd076pck8hGU\nICSrHDy4jXXrZrJq1QXMmDGThQtHqF+SSBtUg5Cs0NhYx+bNP2XTprsZNGg6J510G3l5heqXJBlP\nRWqRNrg727f/mvffn03PnuMZMeJHdO8+KuywRJJGvZhEAuzZ8zrl5bNoaqrllFP+k8LCz4Qdkkha\nUYKQjFNXt4n33/8Xdu1ayvDhdzBo0HTMuoQdlkjaUZFaMkZDwz42bLiN5cvHU1AwgrPPXsfgwdcq\nOYi0k0YQkvbcG9m6dT4bNnyXvn0/w4QJK+nWbVjYYYmkPSUISWs1NWWsXz+LnJzujBv3FL17nx12\nSCIZQwlC0lJt7Xu8//632LdvJSNGzKN//7/V3dxEOpkShKS01t1Wv/rVWzCbz9at8znxxFsZM+ZX\ndOnSLewwRTKS9kFIygrqtvrIIzl8/euX8+lP/5SuXQeGHaJIylOzPslIDzxwW3NygGgjvRkzmnjx\nxS5KDiJJoAQhKau2dn1gt9W6ui3hBCSSZZQgJOU0NR1g48Y72b9/hbqtioRICUJSys6dv2fZstPY\ns+c1vvWtJSxcOFLdVkVCoiK1pIS6uo2Ul89i375VjBp1H8XFUwHUbVWkg9TNVdJWY2MdmzbdzebN\nP2Xo0JsYNuxWLVsV6UQp383VzHKA5cBmd/9SwOv3AVOA/cAMd1+ZjLgkXDt3Pk95+Tfo0WMcZ521\nnIKCkrBDEpEWkrVR7iZgDdC79QtmNgUY6e6jzOwc4BfAuUmKS0IQiVSwfv0s9u//KyeffB/FxVPC\nDklEAiS8SG1mQ4GLgP9o45CLgfkA7v4a0MfMtMg9AzU21lFRcQdvvDGBnj3PYsKEt5QcRFJYMkYQ\n9wC3An3aeH0IsKnF95Wx57YlOC5Jop07n+e99/4vPXuepukkkTSR0ARhZlOBbe6+0sxKAXVTyzKR\nSAXl5f9Ebe3bjBr1M40YRNJIokcQ5wFfMrOLgAKgl5nNd/fpLY6pBFo27x8ae+4oc+bMaX5cWlpK\naWlpZ8crnSS6Ounf2Lz5XoYNm8XYsQu1OkkkCcrKyigrK+uUcyVtmauZTQRubr2KKZY8bnD3qWZ2\nLvBTdz+qSK1lrqmrdcfVK66YxIEDd9Kz5+mcfPI9dOt2UtghimStlF/m2pqZzQTc3R909+fM7CIz\nKye6zPXaMGKS9gnquDpv3iJmz/5/jBs3I+zwRKQDtFFOOmT27KsoLX38iKZ6kQiUlV3JvHmPhReY\niABq9y0h2rdvrTquimQoJQhpl4aGfZSXzyISWaOOqyIZSglCjlu04+o46uur+ed/fkUdV0UylGoQ\nEreDB3ewfv0sdu9+hdGjf0lR0WRAHVdFUpm6uUpCuTtVVU9QXn4zAwdeyfDh36dLlx5hhyUicUi7\nZa6SPurqNvLuu9dx4MAWTjvtv+jd+xNhhyQiSaIahARyb2Tz5ntZvvws+vT5NGedtVzJQSTLxD2C\nMLNPASUtf8bd5ycgJgnZvn1/Zd26r5KTk8+ZZ75K9+6nhB2SiIQgrgRhZguAkcBKoDH2tBNr0y2Z\nobGxjg8+uJMtWx5g+PAfMnjwV4ne60lEslG8I4gJwFhViTNH6/5JV175Jerqbqd79zFMmLCK/Hzt\nYxDJdvEmiL8Cg4APExiLJElQ/6Qf/nAh3/rWzxg37vqwwxORFBHv/EE/YI2ZvWBmzxz6SmRgkjgP\nPHBbc3KAaGuMGTMaWbz41XADE5GUEu8IYk4ig5DkikQ2qn+SiBxTXAnC3f+Q6EAkOWpqXmb//hVE\nIhzVgVX9k0SkpbimmMzsXDNbZmb7zOygmTWa2Z5EByedp7ExQnn5LNaunc6NN96v/kkickxxtdow\ns+XAFcCvia5omg6Mdvd/SWx4R8SgRVTttHfvG6xdezU9epzG6NE/Jy+vWP2TRLJEwnsxmdlyd59g\nZqvd/fTYc2+6+xntedP2UII4fk1NDXzwwZ1UVv47J598LwMHXhF2SCKSZMnoxVRrZl2BlWb2I6LL\nXbWDKoXV1q5j7drp5Ob2ZcKEFeTnDwk7JBFJM/F+yF8dO/ZGoveNHgb8TaKCkvZzb2Lz5vtZseI8\nBg26htNP/72Sg4i0S9ztvs2sADjR3dclNqQ2319TTMdQV7eJdev+noaGvYwZM5/u3UeHHZKIhCzh\nU0xm9kXgbqArMNzMxgPfd/cvtedNpWNat8m47rrv0737Xygvn8XQoTcxbNhscnLUyV1EOibeIvUb\nwCSg7FBh2szecvfTEhxfyxg0giC4TcaCBT24/PLBfPazi+jV68ywQxSRFNKREUS8NYh6d9/d6jl9\nWocgqE3G1Vfv589/nqDkICKdKt55iLfN7CtAFzMbBXwD+HPiwpK21NVVBrbJOHBgWzgBiUjGincE\n8X+BU4EDwBPAbuCmRAUlbcvL6928A/oQtckQkUSIN0GMjX3lAt2Ai4FliQpKgm3f/iTnnvsnHn+8\nWG0yRCTh4i1SrwNuIXpfiKZDz7v7xsSFdlQMWVukbmys4/33b2XnzmcZO3YR1dX91CZDROKSjFYb\nr7j7+e15g86SrQmitracNWsup1u3EZxyyn+Ql9c37JBEJI0kI0F8Fvgy8D9E6xAAuPuT7XnT9sjG\nBFFVtYj33ruRkpI5nHDCP2LWrv+PRSSLJaMX07XAx4A8Dk8xOZC0BJFNDrXmrql5idNPf0HLV0Uk\nFPEmiE+4+ykJjSRLtd4VPWPGDPbtu5nu3T/GhAkryM3tHXaIIpKl4k0Qfzazse6+JqHRZJmgXdFz\n5vyKb35zLmPHfkdTSiISqniXuZ5LtNX3OjNbbWZvmdnqRAaWDYJ2Rc+Y0cSTT76j5CAioYt3BPH5\nhEaRpdraFV1XtyWcgEREWogrQSRzv0M26dIlj0iEI5KEdkWLSKrQXeFCsmXLg3zqU8t54onB2hUt\nIikp7hsGtevkZvnAH4neRyIX+I27z211zETgd8D7saeedPc7As6VEfsgmpoOUl5+E7t2/YFx435H\nVVWudkWLSMIkfKNcR5hZd3evNbMuwKvAN9z99RavTwRuPtbNhzIhQRw8WMXbb08jN7cvY8Y8piWs\nIpJwybgfRLu5e23sYT7RUUTQp3zGL9nZu3cFb7xxNn37TmTcuKeVHEQk5SX8vpRmlgO8AYwE/t3d\ng7rAftLMVgKVwK3pvt+i9ea3adM+wcGDdzBq1AMMGDAt7PBEROKS8Cmm5jcy6w08DdzYMgGYWU+g\nKTYNNQW4191HB/x8WkwxBW1+e/TRXG677WlOPXVq2OGJSJZJRi+mDnP3PWa2lOieijUtnt/X4vHz\nZvZzMyty9+rW55gzZ07z49LSUkpLSxMac3sEbX675poG5s//FfPmKUGISGKVlZVRVlbWKedK9Cqm\nfsTuZ21mBcALwL+6+3Mtjhno7ttij88GFrt7ScC50mIEcdNNn+HSS8uOev6ppz7Dvfe+nPyARCSr\npfIIYjDwaKwOkQMscvfnzGwm4O7+IDDNzK4H6oEI8HcJjimhcnN7aPObiGSEpNUgOiodRhC7dv2J\nl166hKefzuXKK6uaaxALF45k7twl2t8gIkmX0vsgOkuqJ4iqqt/w3nv/yJgxj7Nnz8na/CYiKUEJ\nImSbN9/LBx/8iNNO+2969Toj7HBERJqlcg0iI7Xc51Bfv5XJkw9y4YWvUlBQEnZoIiKdRgniOAXt\nc1i4sITx452SkrCjExHpPOrmepyC9jlccUUFDzxwW7iBiYh0MiWI4xSJbNRNfkQkKyhBHIf6+hoi\nkTXN9284RPscRCQTKUHE6eDBHaxaNYkrr7yEhQtH6iY/IpLxtMw1DgcPVrFq1QUUF09l+PA72bix\nQvscRCQtaB9EJ2u5jLVr10LOO281H//4VZSU3I5Zxt+6QkQyiPZBdKKgZayPPVbE+PHXKDmISFZR\nDaKVoGWsV11VrWWsIpJ1lCBaqaur1DJWERGUII6Snz9Ay1hFRFCCOEJjY4SJEytYsKC3lrGKSNbL\n+lVMh1YsRSKbOXCgnGnTzuTkk+/hF7+4XctYRSTtaZlrOwU33hvB3LkvKSGISEboSILI6imm4MZ7\n72vFkogIWZ4gtGJJRKRtWZ0gzFwrlkRE2pC1CWLfvrc4//zV/OpXQ7RiSUQkQFYVqQ+tWKqt3UBt\n7Squu+6H9O//JTXeE5GMpVVMcQhesTSSuXOXKCGISMbSKqY4BK9YWq8VSyIibciaBKEVSyIixydr\nEkRuboFWLImIHIeMrkEcbqOxkR07VlBbW8DMmTtVgxCRrKEidYDgG/+cSK9eZ5Cfv0crlkQkKyhB\nBJg9+ypKSx8/ou4QiUBZ2ZXMm/dYAiIUEUk9WsUUQEVpEZGOybgEUVGxgdmzr2L9+jU89BBs3Xr4\nNRWlRUTil1FTTEF1h4cfhssugz59VJQWkeyjGkRMW3WHu+8eyIQJF6goLSJZpyMJIrezgwlTW3WH\ncePGqjAtInKcMiZBVFRsYO3aCqqrITcXLrwQBg1S3UFEpL0yIkEcqj3MmlXRXHt45BGYMgWWLh3J\n3Llq3y0icrwSWoMws3zgj0BXosnoN+4+N+C4+4ApwH5ghruvDDimzRpEW7WHe+4p4cEHX1bdQUSy\nVsrWINz9gJl9xt1rzawL8KqZPe/urx86xsymACPdfZSZnQP8Ajg33veoqNjA8uVL2LoVcnIOTy0V\nFMCYMcOVHERE2inhU0zuXht7mB97v9bDgIuB+bFjXzOzPmY20N23Hevch6aWbrml6oippUsvjS5r\nVe1BRKT9Ep4gzCwHeAMYCfy7uy9rdcgQYFOL7ytjzx0zQQTd42HGDHjiCcjJUe1BRKQjkjGCaALO\nMLPewNNmNtbd17TnXHPmzGl+XFpa2uay1pqagdx/vzbEiUj2KSsro6ysrFPOldSNcmZ2G7Df3X/S\n4rlfAEvdfVHs+3eAia2nmFoXqSsqNvD1r09i4MCKo5a1qiGfiEhUyhapzawfUO/uu82sAPgc8K+t\nDnsGuAFYZGbnAruOVX/QslYRkcRL9DLX04BHiTYFzAEWufsPzWwm4O7+YOy4+4HPE13meq27rwg4\nV/MI4rrrLmH//t+Rk3N45VKfPlrWKiLSWsqOINz9LeDMgOd/2er7G+M9Z0XFBjZvfoEbbuColUta\n1ioi0nnSrt33Aw/cxg031B21cunZZ7WsVUSkM6Vdq42qqnIWL4ampiM3xm3b1o0f/EC1BxGRzpJW\nCaKiYgPbt7991PTSlCkwevRkTS+JiHSitLofxMyZFwcWp++6qyfz569WghARaSVli9Sdre3i9Dgl\nBxGRTpZWCWLatLoj6g9TpkSL0yecMDLs0EREMk5aTTFdfnl0xVLLEcSWLV1ZsOAdjSBERAJ0ZIop\nrZa5du0KCxfCo4/C7t3RZJGX11/JQUQkAdIqQRzS0ACPPRZNEqNGDQs7HBGRjJRWNYivfOXw9NLD\nD8OiRTB6tOoPIiKJkFY1iOnTj1ze+r3v5fDrX5driklEpA1ZU4OAI6eX+vUrVHIQEUmQtEsQAO7R\n6aWRI88LOxQRkYyVVjWIltavh3vu+WnYYYiIZKy0HEEA5OTkanpJRCSB0moE8eGH0RVMeXlhRyIi\nkvnSagQxZAj8wz9EN8zV1zeFHY6ISEZLqwRRXg7z58MVV0BuWo19RETST1p9zO7bB2aweDEUFqZV\nbhMRSTtp9Snb0AAHDsC6dZCTMzDscEREMlpajSCammDXrugo4vvffyLscEREMlpatdoYOxYqKiA/\nH6qr0yNuEZEwZU2rjbvvhvPOg9rasCMREcl8aTXFVFAAs2ZBZWXYkYiIZL60GkFANEkUFYUdhYhI\n5ku7BBGJQHV12FGIiGS+tEoQkQjccw/U1PQOOxQRkYyXdquYamp6s2XL7rDDERFJCx1ZxZRWCSJd\nYhURSRVZs8xVRESSRwlCREQCKUGIiEggJQgREQmkBCEiIoGUIEREJFBCE4SZDTWzl83sbTN7y8y+\nEXDMRDPbZWYrYl/fTWRMIiISn0SPIBqAb7r7qcAngRvM7GMBx/3R3c+Mfd2R4JjSXllZWdghpAxd\ni8N0LQ7TtegcCU0Q7r7V3VfGHu8D1gJDAg5t1yaObKVf/sN0LQ7TtThM16JzJK0GYWYlwHjgtYCX\nP2lmK83sWTMbm6yYRESkbUm5H4SZ9QR+A9wUG0m09AZworvXmtkU4GlgdDLiEhGRtiW8F5OZ5QL/\nDTzv7vfGcfwG4Cx3r271vBoxiYi0Q3t7MSVjBPGfwJq2koOZDXT3bbHHZxNNWkfd8aG9/4EiItI+\nCU0QZnYecCXwlpm9CTjwbeAkwN39QWCamV0P1AMR4O8SGZOIiMQnbdp9i4hIcqXcTmoz+7yZvWNm\n75rZ7IDXs2Zj3bGuReyYUjN708z+amZLkx1jssTxe3FL7DqsiG3KbDCzvmHEmmhxXIveZvZMbGXg\nW2Y2I4QwkyKOa9HXzJ40s1Vm9r+ZukrSzB4ys21mtvojjrnPzN6L/V6Mj+vE7p4yX0QTVjnRKag8\nYCXwsVbHTASeCTvWFLkWfYC3gSGx7/uFHXdY16LV8V8AXgo77hB/L/4FuOvQ7wSwE8gNO/aQrsWP\ngNtij0/J4N+L84luI1jdxutTgGdjj88B/jee86baCOJs4D133+ju9cBC4OKA47KhYB3PtfgK8Ft3\nrwRw9x1JjjFZ4v29OOTLwK+SElnyxXMtHOgVe9wL2OnuDUmMMVniuRZjgZcB3H0dUGJm/ZMbZuK5\n+ytAzUcccjEwP3bsa0AfMxt4rPOmWoIYAmxq8f1mgndeZ8PGuniuxWigyMyWmtkyM7s6adElV7y/\nF5hZAfB54LdJiCsM8VyL+4GxZrYFWAXclKTYki2ea7EKuAyaV0meCAxNSnSppfW1qqSNv6GWkrJR\nrpNpY91hucCZwCSgB/AXM/uLu5eHG1aovgi84u67wg4kRBcCb7r7JDMbCSwxs9P96E2q2eBfgXvN\nbAXwFvAm0BhuSOkj1UYQlUQz/CFDY881c/d97l4be/w8kGdmRckLMWmOeS2I/ovpBXevc/edwB+B\njycpvmT6l58PAAACgUlEQVSK51occgWZO70E8V2La4EnAdx9PbABCGqSme7i+bzY6+5/79FGoNcA\nA4D3kxhjqqgEhrX4/qP+hpqlWoJYBpxsZieZWVeif+zPtDyg5bzZR22sywDHvBbA74DzzayLmXUn\nWnxam+Q4kyGea4GZ9SG6iOF3SY4vmeK5FhuBC6D572U0mfmhGM/nRR8zy4s9/hrwhwweSRlt12ef\nAaYDmNm5wC6PbVD+KCk1xeTujWZ2I/Ai0eT1kLuvNbOZZNnGuniuhbu/Y2YvAKuJDpsfdPc1IYad\nEHH+XgBcQnREFQkr1kSL81rcATzSYsnjtzLxH1FxXosxwKNm1kR0xd8/hBdx4pjZE0ApUGxmHwC3\nA105/FnxnJldZGblwH6io8xjnze27ElEROQIqTbFJCIiKUIJQkREAilBiIhIICUIEREJpAQhIiKB\nlCBERCSQEoRIO5nZTWbWrY3XrjGznyU7JpHOpAQh0n7/BHT/iNe1yUjSWkrtpBZJVbFWJouJdsDs\nAvwGOAFYamY73P2zZnYt8M9E2y6vBurCilekMyhBiMTn80Clu38BondtA2YApe5eY2aDgDnAGcAe\noAxYEUqkIp1EU0wi8XkL+JyZ3WVm57v7Ho5sjnYOsNTdq2M351kUVqAinUUjCJE4uPt7ZnYmcBHw\nAzN7maNrDNlwp0PJIhpBiMTBzAYDEXd/Arib6I2a9gK9Y4e8BvwfMyuMtZf+23AiFek8GkGIxOc0\n4N9ibaMPAtcDnwR+b2aVsSL1XOB/iRapV4YXqkjnULtvEREJpCkmEREJpAQhIiKBlCBERCSQEoSI\niARSghARkUBKECIiEkgJQkREAilBiIhIoP8PFQA3RdgPFmcAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/html": [ "\n", "\n", "\n", "\n", "
\n", "" ], "text/plain": [ "" ] }, "execution_count": 9, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# Sub-lists represent matrix columns, diagonal is the variance of each intervention\n", "covariances = [[1,0,0,0],[0,1,0,0],[0,0,1,0],[0,0,0,1]] \n", "\n", "mean_returns = [1.2,5.0,3.0,1.5]\n", "\n", "weights, returns, risks, portfolios = covmean_portfolio(covariances, mean_returns)\n", "\n", "# Matplotlib only plot: \n", "plt.ylabel('mean')\n", "plt.xlabel('std')\n", "plt.plot(risks, returns, 'y-o') #risks and returns are just arrays of points along the frontier\n", "plt.show()\n", "\n", "# Interactive mpld3 plot:\n", "interactive_plot(risks, returns, portfolios)\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Model Shortcomings\n", "\n", "\n", "- Markowitz's version of Modern Portfolio Theory (MPT) assumes past performance is indicative of the future returns. Similarly, if the user of this model inputs an average return to an intervention, this does not take into account the declining marginal returns to doing good in the world (e.g. as more money is poured into scientific research, it might become harder to make transformative discoveries). A person could avoid this problem by only using marginal cost effectiveness estimates, and regularly updating them or using a decay function based on society's investments.\n", "\n", "- MPT is usually limited to publicly traded assets with a long history of returns. It is hard to include returns from private equity or real estate in this type of analysis. High risk philanthropic endeavors might be more analogous to venture capital investments, so MPT might not be the best model to use. \n", "\n", "- Using the standard deviation assumes a normal distribution. There are versions of MPT that use other distributions that might be useful here. \n", "\n", "- There is very little reason to be risk averse when it comes to philanthropy because the returns don't accrue to us. Someone else benefits, so why not take the riskiest option with the highest expected value? \n", "\n", "A more complete list of criticisms to Modern Portfolio Theory is available [here](https://en.wikipedia.org/wiki/Modern_portfolio_theory#Criticisms)." ] }, { "cell_type": "markdown", "metadata": { "collapsed": true }, "source": [ "# Conclusion\n", "\n", "This approach might be a useful way to come up with an optimal portfolio for charitable giving. This allows us to not only combine DALY/\\$ measures, but also account for the uncertainty involved in imperfect estimates. The next step for this analysis is to find a real world dataset that incorporates both DALY/\\$ estimates along with their variances. I think it's likely this dataset whould show that there are returns to risk taking in charitable giving. \n", "\n", "The DALY is an imperfect measure, as there are some things that are very hard to quantify, but it is still a worthwhile exercise to try to come up with back of the envelope calculations to partially inform the decision making process. Note that this approach isn't necessarily tied to the DALY/\\$, as any cross-intervention measure of social return could be used. In some cases back of the envelope calculations could actually harm our judgement by giving us a false sense of certainty. One way to guard against this risk is to only use this as a part of the decision making process, rather than relying solely on this type of analysis. \n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Next Steps\n", "\n", "### Real Datasets\n", "\n", "Although the above example shows a theoretical relationship between risk and return, it should be possible to use a few existing datasets to find if this relationship holds in reality. \n", "\n", "Here are some potential data sets to test this theory:\n", "\n", "- Disease Control Priorities Project ([DCP2](http://dcp-3.org/dcp2/chapters)): This project took a very comprehensive look at the global disease burden, and provided DALY/\\$ [estimates](http://dcp-3.org/sites/default/files/dcp2/DCP02.pdf) for a broad variety of diseases in the developing world. Unfortunately, the estimates don't include sensitivity analyses, and mostly just include a range of point estimates. I could use the spread of the point estimates from their dataset, but the downside of this is that the interventions with very few estimates would automatically have the lowest variance. \n", "\n", "- The Copenhagen Consensus: [The Copenhagen Consensus](http://www.copenhagenconsensus.com/) performs cost benefit analyses for a broad variety of health and policy interventions. They convert everything into monetary cost benefit ratios terms by putting a value on a human life rather than use the DALY/\\$. It would be possible to plot the standard deviation of their estimates against the cost benefit ratio to create an efficient frontier. However, they do not seem to publish a dataset including both cost benefit ratios and sensitivity analyses. Their approach relies heavily on expert judgmement as well, which might not be quantifiable. Converting everything into monetary terms also has the downside of valuing present lives over future lives because future cash flows are discounted. \n", "\n", "- [AidGrade](http://www.aidgrade.org/): This organization performs meta-analyses that combine the treatment effects of most of the known randomized control trials in the development economics literature. This [presentation](http://evavivalt.com/2015/09/effective-altruism-talk) by Eva Vivalt outlines the initial results, with treatment effects in standard deviations. In this case, you could plot standard deviation in treatment effect against the mean treatment effect, as she does in the presentation. This might give you an idea of the curve shape, although these values should be converted to a cost effectiveness figure like DALY/\\$ that can be compared across studies. \n", "\n", "- Crowdsourcing: People could submit and discuss cost effectiveness estimates and critique them publicly. GitHub might be a good venue for these open source cost effectiveness conversations. \n", "\n", "### Black-Litterman Model\n", "\n", "The [Black-Litterman model](https://en.wikipedia.org/wiki/Black%E2%80%93Litterman_model) generates the expected returns of different assets based on the market's current equilibrium asset allocation rather than using their historical returns. This concept could potentially be applied to society's entire allocation of money via markets, governments and philanthropy. This approach would incorporate the entirety of society's wisdom into a philanthropic portfolio. This is kind of an off-the-wall idea and I don't think it makes sense, but I might write a post about it in the future. \n" ] } ], "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.12" } }, "nbformat": 4, "nbformat_minor": 0 }