Skip to content

Commit

Permalink
Browse files Browse the repository at this point in the history
added 4.7.5
  • Loading branch information
jet08013 committed Apr 27, 2018
1 parent 756240b commit 770e7f4
Showing 1 changed file with 166 additions and 0 deletions.
166 changes: 166 additions & 0 deletions BDA 4.7.5.ipynb
@@ -0,0 +1,166 @@
{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Problem 4.7.5.\n",
"\n",
"Approximate mean and variance.\n",
"1. Suppose x and y are independent normally distributed random variables, where x~N(4,1) and y~N(3,2). What are the mean and standard deviations of y/x? Compute this using simulation.\n",
"\n",
"2. Do the same computation without simulation.\n",
"\n",
"3. What assumptions do you need for part (2)?"
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"import matplotlib.pyplot as plt\n",
"import numpy as np\n",
"from scipy.stats import norm"
]
},
{
"cell_type": "code",
"execution_count": 48,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"mean: 0.8079667283774233 sd: 0.6574088802979349\n"
]
}
],
"source": [
"x=norm(4.0,1.0)\n",
"y=norm(3.0,2.0)\n",
"\n",
"samples_x=x.rvs(10000)\n",
"samples_y=y.rvs(10000)\n",
"z=samples_y/samples_x\n",
"print('mean:',np.mean(z),'sd:',np.sqrt(np.var(z)))"
]
},
{
"cell_type": "code",
"execution_count": 57,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAX0AAAD8CAYAAACb4nSYAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzt3XuYXFWZ7/Hvr5skQEAIJiBCQgcNl3ANNAFkJujhYkAn0VGPiccRBE/OzHCZGceZE2RGMzBqvMxB58goyOQIoxIRBaJGMdy8kpAGAiGBkBBCCIkkGAgCIaHT7/lj74ZKpS+7u6v2rur6fZ6nnq7al1pvde16a9Xaa6+liMDMzBpDU9EBmJlZfpz0zcwaiJO+mVkDcdI3M2sgTvpmZg3ESd/MrIE46ZuZNRAnfTOzBuKkb2bWQHYrOoByI0eOjJaWlqLDsEHs/vvvfy4iRuVdro9tq6asx3XNJf2Wlhba2tqKDsMGMUlPFVGuj22rpqzHtZt3zLogabKkFZJWSZrZw3YflBSSWvOMz6y/nPTNykhqBq4GzgHGA9Mlje9iu72BS4FF+UZo1n9O+ma7mgisiojVEbEdmAtM7WK7K4EvAa/mGZzZQGRq05c0Gfga0AxcFxGzy9ZfBbwrfbgnsH9E7Juu2wEsTdetjYgplQjcsmuZ+dPX76+Z/Z4CI6kbBwFPlzxeB5xcuoGkCcDoiPiJpE/lGZwlfFz3T69Jv+Sn7lkkB/9iSfMiYnnnNhHxdyXbXwJMKHmKrRFxfOVCNqs6dbHs9YknJDUBVwHn9/pE0gxgBsCYMWMqFJ5Z/2Vp3sn6U7fTdODGSgRnVpB1wOiSxwcD60se7w0cDdwjaQ1wCjCvq5O5EXFtRLRGROuoUbn3EjXbRZbmnV5/6naSdAgwFrirZPHuktqAdmB2RNzaz1jN8rIYGCdpLPAMMA34SOfKiNgCjOx8LOke4FMR4f6YOSht1rG+y5L0e/ypW2YacHNE7ChZNiYi1ks6FLhL0tKIeGKnAvwT2GpIRLRLuhi4neQ81pyIWCbpCqAtIuYVG2FjcHKvjixJv7efuqWmAReVLoiI9enf1WmNaALwRNk21wLXArS2tnrSXitcRMwH5pct+0w3274zj5jMKiFLm/7rP3UlDSVJ7LvUdCQdDowA7i1ZNkLSsPT+SOA0YHn5vmZmlo9ea/p9+Kk7HZgbEaU19SOBayR1kHzBzC7t9WNmZvnK1E8/y0/diJjVxX6/A44ZQHxmZr1yn/3sfEWumVkDqblRNs3MBsK1/p65pm9m1kBc028w5X2fXRMyayxO+oOUL2wxs6446ZvZoOX2/V25Td/MrIE46ZuZNRAnfTOzBuI2fTOrGe6AUH2u6ZuZNRDX9BvI4VrLpKaH2cowbt1xGi+xZ9EhmVnOnPQbxFFaw9yhV7K3tgLw3uaFnLf9fxcclVnCzTr5cfNOAxjFC3x76Gy2MJxJ267ib7f/Nac0PcrsId8qOjQzy5lr+g3gE7v9lP34I5O3f5G1cQBr4wAObV/PpbvdChseggOPKzpEM8uJa/qD3Jt4if/RfCc/7jiVlXHw68uva38Pf4w94DdfLTA6M8ubk/4g97HmBeylV/lm+5Sdlr/IcL674wxYfitsXl1QdGaWNyf9QS2Yttvd/HLHsTwWY3ZZO6f9HEDwwA35h2ZmhXDSH8SO0xMcrOe4bcc7uly/kREwdhIsvw12mtrYzAYrJ/1B7NzmRWyPZu7oOLH7jY56X9K88+wj+QVmZoVx0h+0gvc0L+I3HcfwIsO73eqEHwxjR4j/+/Wv5BibNaqWmT99/WbFyJT0JU2WtELSKkkzu1h/vqRNkpakt0+UrDtP0sr0dl4lg7fuHavVHKznmN9xco/bbeZNLOwYz7lNi3KKzCzhL4Bi9Jr0JTUDVwPnAOOB6ZLGd7Hp9yPi+PR2XbrvfsBngZOBicBnJY2oWPTWrdObHqIjxJ07JvS67S86Wnlb0wbY/GQOkZlZkbJcnDURWBURqwEkzQWmAssz7PtuYEFEbE73XQBMBm7sX7jWk9Ia0/eHPsKyOITneVOv+/2m4+jkzup7YL+xVYrOzGpBluadg4CnSx6vS5eV+4CkhyXdLGl0H/e1ChrOVk7QSn7TcUym7Z+It7Ih9oPVd1c5MrPiuCkpkSXpq4tl5f37fgy0RMSxwB3A9X3YF0kzJLVJatu0aVOGkKwnJzc9yhDt4NcZkz6I3+w4Glb/Ejp2VDU2MytWlqS/Dhhd8vhgYH3pBhHxh4jYlj78FnBi1n3T/a+NiNaIaB01alTW2K0bf9q0lFdjCPd3HJZ5n193HAOvvgAbllQxMjMrWpY2/cXAOEljgWeAacBHSjeQdGBEbEgfTgEeTe/fDny+5OTt2cBlA47aenRa0yPc13EE2xiaeZ/flrbrH9RDv36zOlfaxLNm9nsKjKQYvdb0I6IduJgkgT8K3BQRyyRdIalzQJdLJS2T9BBwKXB+uu9m4EqSL47FwBWdJ3WtOkbwIoc1PcPCjq46WHXvD+wDo46Ap35XpcjMrBZkGlo5IuYD88uWfabk/mV0U4OPiDnAnAHEaH1wUtMKABZ1HNH3ncecCktvTtr1m5orHJmZ1QJfkTvITGx6jG0xhKVxaN93PuQ02P5H+P3SygdmZjXBSX+QOalpBUvibWxnSN93PuTU5O/aeysblJnVDCf9QWQ4WzlaT/avaQdgn4NhnzFu1zcbxJz0B5ETmlbSrGBxf5M+wCHvSGr6HmrZbFBy0h9EWpseZ0eIBzrG9f9JDjkVXt7k2bTMBilPjD6InKDHeSzG8DJ79P9JRqejcj69CN78tsoEVockTQa+BjQD10XE7LL1nwQ+AbQDm4ALIuKp3AOtA7U89EEj9tl3TX+w6NjBhKZVfboKt0sjD4dh+yRJv0FlHFn2QaA1HXrkZuBL+UZp1j9O+oPFxuXspVdpG2jSb2qC0SfB0/dVJq769PrIshGxHegcWfZ1EXF3RLySPlxIMsSIWc1z0h8s0pr5AzHApA9JE8/GR2HrCwN/rvrU19FhLwR+1tUKDyZotcZJf7B4+j6ejX1ZFyMH/lyjTwYC1rUN/LnqU6bRYQEkfRRoBb7c1XoPJmi1xidyB4u1C9NeO13lq2w6T2rtyass36MJnl4I486sUIB1JdPosJLOBC4HTi8ZZdasprmmPxj88ffwwlMDb89PvcLucMDRjdyu//rIspKGkowsO690A0kTgGuAKRGxsYAYzfrFNf3BIE3OD1Qo6QNJE89DN8KOdmhurMMkItoldY4s2wzM6RxZFmiLiHkkzTl7AT+QBLA2IqZ0+6QNqJa7ajayxvo0D1ZPL4LmYTwSFZzfdvTJsPhbsHEZHHhc5Z63TmQYWbYh272s/rl5ZzB4ehEcdAKvVfI7fPTE9LkbtonHbFBy0q93r70K65e8kaQrZd8xsPeBDX2Rltlg5KRf7zYsgY7X3hg+oVKk5IvESd9sUHHSr3NfuObbAJz47S2Vf/LRJ8MLa+HFDb1va2Z1wSdy69xJTY/xRMeByRy3lTbmlOTv2nvh6D+v/POb1ZBGGXzNNf161tFBa9PjLO44vDrP/5ZjYcienknLbBBx0q9nmx5jX708sElTetI8BA4+CZ5y0jcbLDIlfUmTJa2QtErSzC7Wf1LSckkPS7pT0iEl63ZIWpLe5pXvawOwNpnWcHFUqaYPyUxazz7SyIOvmQ0qvSb9CowtvjUijk9vvmKxktYu5NnYl7Wxf/XKGHMqyeBri6tXhpnlJktN32OL16IIeOretGmn/4Os9ergVmjazZOlW0NpmfnTQTuMRJbeO12NLd5Tp/DyscV3l9RGMq3c7Ii4tc9R2q6eXwMvrmNRR5VHAxg6HN46AZ76bXXLsUFhsCbKwSRL0u/P2OKnlyweExHrJR0K3CVpaUQ8UbbfDGAGwJgxYzIF3vCe/BUAv+s4qipPv1P3tbMmwW+/Btv+CMP2rkp5ZpaPLM07fR1bfErp2OIRsT79uxq4B5hQvq8nmuiHJ38Fex3AE/HW6pc1dhJ0tMPahdUvy8yqKkvS7/fY4pJGSBqW3h8JnAYsr1TwDSsiSfpjJ1HV9vxOo0+G5qHw5C+rX5aZVVWvST8i2oHOscUfBW7qHFtcUmdvnNKxxUu7Zh4JtEl6CLibpE3fSX+gnnscXt6YJv0cDNkjSfxpk5KZ1a9MwzD0d2zxiPgdcMxAArQurE5r3GMnAcvyKXPsJLj78/DKZthzv3zKNCvYYByawVfk1qMn7oQRLcktL4e+EwhYfXd+ZZpZxXnAtXrz2takpn/Cx/It96ATYY8RsHIBHP2BfMu2muZumvXFNf16s+a30L4Vxp2db7lNzfD2M5Ok39GRb9lmVjFO+vVm5S9gtz2g5bT8yx73bnjlOVj/YP5lm1lFOOnXkwhYeXtyUnXIHvmX//YzACUxmFldctKvJxsfTYZfOCznpp1Oe+6XTKH42PzetzWzmuQTufVk2S2gJlp/uCfP/bCgk2fjp8Ltn4bnVsHItxcTgxWuEU/eDpbum67p14sIWPYjaPkTnqvG1Ig96BxxsGXmT2H8+5KFy27JNQYzqwzX9OvFs4/AH1bBqRcl10UXZZ+DkjH2l/0ITv+HAgOxvDVi7X4wck2/XjzyQ1AzHFkD89Ac9eewcXlyjsHM6oqTfj3Y8Ros+R6MOwuGjyw6Gjjq/dA0BB64oehIzKyP3LxTD1bMh5eehdYLio4ksdcoOPLPYMl34b/9Mwzds+iIzHJVzyd1XdOvB21zYJ/RyRWxteKkC+HVLT6h2wAG89SBjcg1/Vq38TFYfQ+86/JkKISCvfHhD9YcfBjcdw0c/xFQDuP6m9mAuaZf6375RRi6F7ReWHQkZQTvuAQ2PASP/7zoYMwsIyf9Wrbx0aT5ZOIMGP7moqPZ1XHTk+Gd7/lCch2BmdU8J/1aFQF3/AsMHZ7UqGtR8xCY9I9Jbd9t+9ag6u2ch9v0a9Xy2+Dxn8FZV8Ce+9XuQXXsh5N2/Z/PhLe9Kxlz38xqlpN+LXr5OZj/D3DgcXDKRUVH063OL6Lx+jDzd/9n+Pmn4f3fKDgqq4SarWTUsHrpxunmnVrTvg2+/1HY9iJM+To01/738vJogT/9e3joe7DomqLDMbMe1H5GaSTt2+BHM2DtvfDBOXDgsUVHlN07Z8Kzy5Jmnj3fDMd8sOiIzApTy7V+J/1a8dJGuPkCWPNrOPtf624e2pZP/5w9+ADXD13FxB9eCC88Baf9bU1cW2DdczNO48mU9CVNBr4GNAPXRcTssvXDgBuAE4E/AB+OiDXpusuAC4EdwKUR4WmXSr22FR78Dtx1ZXL//dfCcR+uyw/jVnbnL7ZfxleGfJM/u/MKHlzwPT7/2kf4wec/WXcXbw3kmDcr1dVnucjaf69JX1IzcDVwFrAOWCxpXkQsL9nsQuD5iHi7pGnAF4EPSxoPTAOOAt4K3CHpsIjYUekXUjciYMvTyTyzq+6ER38MWzfDIafBe78Kow4rOsIB2cZQLnntEu7YcQL/NOQ7/GDYFfDNW+DI90LLn8BbjoHd850PoK8GcsznH2129ViRsMrLUtOfCKyKiNUAkuYCU4HSD8BUYFZ6/2bg65KULp8bEduAJyWtSp/v3m5Le/GZpBdIXUgvSIpI7kdAdEDsgI522NEO7Vth+8vw6ovw8qbk9bW/muw3dO9k5MyTPgGHvIOWy+YDK4t6MRUkbuv4E27fdhJ/3vwbPrDhV0z4/WyalPy/NsdePBsjeD725iX2YCvDmHrioclJazWnTUJKfx2o5FdCbr8W+n3MR/gqNetdd1/AefwCyJL0DwKeLnm8Dji5u20iol3SFuDN6fKFZfseVF6ApBnAjPThSzrnCysyRd+1kcBzA9g/x3JfBL6d3vIuuyJ6LfcL6W1nLwLrq1puLw7pZf1Ajvmd4trl2Jb6e2wX9R6Xcxw7q2gc+uKA4ujtuAayJf2uqlfltZnutsmyLxFxLXBthlh6JaktIlor8Vz1UG6RZQ/icgdyzO+8oELHdpHHl+OomzhasmybpZ/+OmB0yeOD2bWa9vo2knYD9gE2Z9zXrNYM5Jg3q2lZkv5iYJyksZKGkpyYnVe2zTzgvPT+B4G70rbNecA0ScMkjQXGAfdVJnSzqhnIMW9W03pt3knbKy8GbifpvjYnIpZJugJoi4h5wH8C/5WeqN1M8iEh3e4mkhNg7cBFOfTcqUgzUR2VW2TZg7LcgRzzVVTk8VXKceys7uKQKydmZo3DY++YmTUQJ30zswYyaJO+pE9JCkkjcyzzy5Iek/SwpFsk7Vvl8iZLWiFplaSZ1SyrrNzRku6W9KikZZL+Jq+y0/KbJT0o6Sd5lluEot7jshgKfb+7iKcm3n9J+0q6Of3MPyrp1ILi+Lv0fXlE0o2Sdu9p+0GZ9CWNJrmEfm3ORS8Ajo6IY4HHgcuqVVDJUAHnAOOB6emwF3loB/4+Io4ETgEuyrFsgL8BHs2xvEIU/B6XKvr9Llcr7//XgJ9HxBHAcRQQk6SDgEuB1og4mqTjQY+dCgZl0geuAv6RLi6WqaaI+EVEtKcPF5L0766W14cKiIjtQOdQAVUXERsi4oH0/h9JDvZdrrSuBkkHA+8BrsujvIIV9h6XKvL9Llcr77+kNwGTSHpxERHbI+KFgsLZDdgjvV5kT3q5FmrQJX1JU4BnIuKhgkO5APhZFZ+/q6ECcv8gSmoBJgCLciryqyRf6B05lVekmniPSxXwfperlff/UGAT8P/SpqbrJA3PO4iIeAb4CkmrxgZgS0T8oqd96jLpS7ojbb8qv00FLgc+U1DZndtcTvKT+LvVioOMwwBUk6S9gB8CfxsRL+ZQ3nuBjRFxf7XLqhGFv8el8n6/uyi/lt7/3YATgG9ExATgZSD3cy6SRpD8+htLMpLxcEkf7WmfupxEJSLO7Gq5pGNIXvxDySCfHAw8IGliRPy+mmWXxHAe8F7gjCpfoVnoEBeShpAkgO9GxI9yKvY0YIqkc4HdgTdJ+k5E9HiQ17GaGcakoPe7XC29/+uAdRHR+YvnZgpI+sCZwJMRsQlA0o+AdwDf6XaPiBi0N2ANMDLH8iaTXH08KoeydgNWk3zJDQUeAo7K6XWKZAKRrxb43r4T+ElR5ef0Ggt7j2vt/a7F9x/4NXB4en8W8OUCYjgZWEbSli/geuCSnvapy5p+Dfs6MAxYkP7SWBgRf1mNgqKboQKqUVYXTgP+AlgqaUm67NMRMT+n8htCwe9xKb/fXbsE+G46PtNq4ON5BxARiyTdDDxA0qT8IL0MydDrMAyS5pA0V2yMpEtQ+XqRdF06F3gFOD/SM/1pU8c/pZv+a0Rc36dXZGZmFZXlRO63SZotunMOyeiZ40gmi/gGgKT9gM+S/PyYCHw2PelgZmYF6TXpR8Sv6Hmc8KnADZFYCOwr6UDg3cCCiNgcEc+TXLjU05eHmZlVWSW6bHbXl7jm+hibmTW6SpzIHdBUibDzPKLDhw8/8YgjjqhAWGZdu//++5+LiFF5lzty5MhoaWnJu1hrEFmP60ok/e76Eq8j6VZVuvyerp4gSuYRbW1tjba2tgqEZdY1SU8VUW5LSws+tq1ash7XlWjemQd8TIlTSC4D3kDSzexsSSPSE7hnp8vMzKwgvdb0Jd1IUmMfKWkdSY+cIQAR8U1gPkl3zVUkXTY/nq7bLOlKkvlGAa6ICE8cbWZWoCxz5E7vZX0AF3Wzbg4wp3+hmZlZpdXlgGtmNkjN2ie5WdU46ZuZNRAnfTOzBuKkbw3lggsuADhO0iOdyyTNkvSMpCXp7dySdZel89OukPTukuWFz11r1h9O+tZQzj//fICVXay6KiKOT2/zAdJ5YKcBR5EMIfIf6aTctTJ3rVmfeWhlayiTJk2CZAjaLBWeqcDciNgGPClpFcnggZDOXQsgqXPu2uWVj9isslzTN0tcLOlhSXNKRoMd8LhSkmZIapPUtmnTpmrEbdYnTvpmyXDgbwOOJ5lc+t/S5QMeVyoiro2I1ohoHTUq9+F+zHbh5h1reBHxbOd9Sd8CfpI+7GmO2pqYu9asr1zTt4aXzv/Q6f1AZ8+eecA0ScMkjSWZKOg+kqFFxkkam06VNy3d1qzmuaZvDWX69OkAR5DM9Nk5ltQ7JR1P0kSzBvhfABGxTNJNJCdo24GLImIHyc61MHetWZ856VtDufHGG5k7d+7DEdFasvg/u9s+Ij4HfK6L5fNJBhs0qytu3jGz2lA65o7H36kaJ30zswbipG9m1kCc9M3MGkimpN/b4FKSrioZrOpxSS+UrNtRss7d2szMCpRlusTOwaXOIrlYZbGkeRHx+jgjEfF3JdtfAkwoeYqtEXF85UI2M7P+ylLTn0g6uFREbAc6B5fqznTgxkoEZ2ZmlZUl6fdlcKlDgLHAXSWLd08HnFoo6X39jtTMzAYsy8VZmQeXIrkc/ebOqxZTYyJivaRDgbskLY2IJ3YqQJoBzAAYM2ZMhpDMzKw/stT0exp0qtw0ypp2ImJ9+nc1cA87t/d3buORCM3McpAl6WcaXErS4cAI4N6SZSMkDUvvjwROwxNNmJkVptfmnYho72pwKUlXAG0R0fkFMJ1klqHSpp8jgWskdZB8wcwu7fVjZmb5yjTgWleDS0XEZ8oez+piv98BxwwgPjMzqyBfkWtm1kCc9M3MGoiTvplZA3HSN7PaNGsfj6tfBU76ZmYNxEnfzKyBOOmbmTUQJ31rKBdccAHAcZIe6VwmaT9JCyStTP+OSJdL0r+n80g8LOmEkn3OS7dfKem8/F/JIOO2+9w46VtDOf/88wFWli2eCdwZEeOAO9PHAOcA49LbDOAbkHxJAJ8FTiYZevyznV8UZrXOSd8ayqRJkwDayxZPBa5P718PvK9k+Q2RWAjsK+lA4N3AgojYHBHPAwuAyVUP3qwCnPTN4ICI2ACQ/t0/Xd7dXBKZ55gwqzVO+mbd624uicxzTEiakU4i1LZp06aKBmfWH076ZvBs2mxD+ndjury7uSQyzzHhuSKs1jjpmyXzQ3T2wDkPuK1k+cfSXjynAFvS5p/bgbPT+SJGAGeny8xqXqahlc0Gi+nTpwMcQdIjcx1JL5zZwE2SLgTWAh9KN58PnAusAl4BPg4QEZslXUkywRDAFRGxObcXYTYATvrWUG688Ubmzp37cES0lq06o3zbdEKgi7p6noiYA8ypQohmVZWpeUfSZEkr0otUZnax/nxJmyQtSW+fKFnni1jMzGpErzV9Sc3A1cBZJCewFkua18W0h9+PiIvL9u28iKWVpHfD/em+z1ckejMz65MsNf2JwKqIWB0R24G5JBetZOGLWMzMakiWpJ/1QpQPpOOT3Cypsztbpn3dl9nMLB9Zkn6WC1F+DLRExLHAHbxxSXumi1jcl9nMLB9Zkn6vF6JExB8iYlv68FvAiVn3NTOz/GRJ+ouBcZLGShoKTCO5aOV1nVczpqYAj6b3fRGLmVkN6bX3TkS0S7qYJFk3A3MiYpmkK4C2iJgHXCppCsnohZuB89N9fRGLmVkNyXRxVkTMJ7k6sXTZZ0ruXwZc1s2+vojFzPpv1j4wa0vRUQwaHnvHzKyBOOmbmTUQJ30zswbipG9m1kCc9M3MGoiTvplZA/F4+mZWnFn7FB1Bw3FN38ysgTjpm5k1ECd9M7MG4qRvZtZAnPTNUpLWSFqazvPcli7bT9KCdI7nBelosSjx7+m80Q9LOqHY6M2ycdI329m7IuL4iGhNH88E7oyIccCd6WOAc4Bx6W0G8I3cIzXrByd9s55N5Y2Z4K4H3ley/IZILAT2LZtXwqwmOembvSGAX0i6X9KMdNkBEbEBIP27f7o869zRVgnuz18xmZK+pMmSVqTtlzO7WP9JScvTts07JR1Ssm5H2ka6RNK88n3NashpEXECSdPNRZIm9bBtpvmfJc2Q1CapbdOmTZWK06zfek36kpqBq0k+COOB6ZLGl232INCaTox+M/ClknVb0zbS4yNiSoXiNqu4iFif/t0I3AJMBJ7tbLZJ/25MN880/3NEXBsRrRHROmrUqGqGb5ZJlpr+RGBVRKyOiO3AXJL2zNdFxN0R8Ur6cCHJB8CsbkgaLmnvzvsk8zk/QjIf9HnpZucBt6X35wEfS3vxnAJs6WwGMqtlWcbe6art8uQetr8Q+FnJ493T7m/twOyIuLXPUZpV3wHALZIg+Vx8LyJ+LmkxcJOkC4G1wIfS7ecD5wKrgFeAj+cfslnfZUn6mdouASR9FGgFTi9ZPCYi1ks6FLhL0tKIeKJsvxkk3d4YM2ZMpsDNKikiVgPHdbH8D8AZXSwP4KIcQjOrqCzNO5naLiWdCVwOTImIbZ3LS9pJVwP3ABPK93W7p5lZPrIk/cXAOEljJQ0FppG0Z75O0gTgGpKEv7Fk+QhJw9L7I4HTgOWVCt7MzPqm1+adiGiXdDFwO9AMzImIZZKuANoiYh7wZWAv4Adpm+jatKfOkcA1kjpIvmBmR4STvplZQTJNohIR80lOXJUu+0zJ/TO72e93wDEDCdDMzCrHV+SamTUQJ30zK4aHViiEk76Z1YdZ+/iLogKc9M3MGoiTvplZA3HSNzNrIE76ZmYNxEnfzOqLT+YOiJO+mVkDcdI3M2sgTvpmZg3ESd/MrIE46ZtZvipxZa1P5vabk76ZWQNx0jczayBO+mZmDSRT0pc0WdIKSaskzexi/TBJ30/XL5LUUrLusnT5CknvrlzoZsXq7XNhXahkW3znuQG37/dJr0lfUjNwNXAOMB6YLml82WYXAs9HxNuBq4AvpvuOJ5lT9yhgMvAf6fOZ1bWMnwuzmpOlpj8RWBURqyNiOzAXmFq2zVTg+vT+zcAZSibLnQrMjYhtEfEksCp9PrN6l+VzYXlxrT+zLHPkHgQ8XfJ4HXByd9ukE6lvAd6cLl9Ytu9B/Y7WrHZk+VxYEUm4tMxZW954PGtL/rHUoCxJX10si4zbZNkXSTOAGenDbZIeyRBXNYwEnmugcossu8jXfHgFnqM/x/ZLklak94t8/Z1qIQaoZhz/oq7v5xlD3wwkjkOybJQl6a8DRpc8PhhY38026yTtBuwDbM64LxFxLXAtgKS2iGjNEnxuYyE+AAAFKUlEQVSlFVW2X3P+ZVfgafp8bJfHUNTrr6UYaiWOWoghrziytOkvBsZJGitpKMmJ2Xll28wDzkvvfxC4KyIiXT4t7d0zFhgH3FeZ0M0KleVzYVZzeq3pp230FwO3A83AnIhYJukKoC0i5gH/CfyXpFUkNfxp6b7LJN0ELAfagYsiYkeVXotZbrr7XBQcllmvsjTvEBHzgfllyz5Tcv9V4EPd7Ps54HN9iGmXn8I5Kqpsv+Y6LLurz0XeMQxQLcQAtRFHLcQAOcShpBXGzMwagYdhMDNrIIUl/YEM7ZBD2Z+UtFzSw5LulJSpK9RAyy3Z7oOSQlJFzuJnKVfSf09f8zJJ36tEuVnKljRG0t2SHkz/3+dWqNw5kjZ21/1XiX9P43pY0gmVKLeXmD6U/n87yt/booYrkTRL0jOSlqS3ivz/M5ZdE8NYSFojaWn6+ivRsytLmbscn5L2k7RA0sr074iqFB4Rud9ITnw9ARwKDAUeAsaXbfPXwDfT+9OA7+dY9ruAPdP7f1WJsrOUm263N/ArkovaWnN6veOAB4ER6eP9c/xfXwv8VXp/PLCmQmVPAk4AHulm/bnAz0j6258CLKr0cd5FmUeSXCNwT+l7m77uh4BhwNj0f9Zc7XjSsmcBn8qjrL4eGznGsgYYmXOZuxyfwJeAmen9mcAXq1F2UTX9gQztUPWyI+LuiHglfbiQpA921ctNXUny5r9agTKzlvs/gasj4nmAiNiYY9kBvCm9vw9d9HXvj4j4FUlPsu5MBW6IxEJgX0kHVqLsHmJ6NCJWdLGqEYcraehhLLo5Pktz3vXA+6pRdlFJv6tL2MuHZ9hpaAegc2iHPMoudSFJjbDq5UqaAIyOiJ9UoLzM5QKHAYdJ+q2khZIm51j2LOCjktaR9IS5pEJl96avx8FgjuXitIlrTtWaFHZV9GsuFcAvJN2fXkFdlAMiYgNA+nf/ahSSqctmFQxkaIc8yk42lD4KtAKnV7tcSU0kI5SeX4GyMpeb2o2kieedJL9qfi3p6Ih4IYeypwPfjoh/k3QqyfUeR0dExwDLrkRsfX9S6Q7gLV2sujwibsszliwxAd8g+XUZ6d9/Ay6oVNk9hdXFsqK6Ep4WEesl7Q8skPRYWhMflIpK+gMZ2iGPspF0JsmH4vSI2JZDuXsDRwP3pK1YbwHmSZoSEQM5uZT1f70wIl4DnlQyPsw4kqtOByJL2ReSDLtNRNwraXeS8Ucq1cQ0kNj6LCLOrJVYOmWNSdK3gEr+yuxJVV9zX0TE+vTvRkm3kDQ9FZH0n5V0YERsSJsaq/IZKKp5ZyBDO1S97LSZ5RpgSgXbt3ssNyK2RMTIiGiJiBaScwkDTfi9lpu6leTkNZJGkjT3rB5guVnLXguckZZ9JLA7sKkCZfdmHvCxtBfPKcCWzp/WBShsuJKy8xjvB/Ia7LAmhrGQNFzS3p33gbPJ739QrjTnnQd098twYPI8Y1129vpc4HGSM/iXp8uuIEl0kHz4f0ByUus+4NAcy74DeBZYkt7m5VFu2bb3UIHeOxlfr4D/QzJcxlJgWo7/6/HAb0l6bywBzq5QuTcCG4DXSGqVFwJ/CfxlyWu+Oo1raaX+173E9P40lm3p8XV7ybrL01hWAOdUO5aScv8rff0PkySdA3Mse5djI+8bSe+hh9Lbsrzi6Ob4fDNwJ7Ay/btfNcr2FblmZg3EV+SamTUQJ30zswbipG9m1kCc9M3MGoiTvplZA3HSNzNrIE76ZmYNxEnfzKyB/H/GxpPhCZ8T7QAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 432x288 with 4 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"h=norm(.807,.657)\n",
"fig,ax=plt.subplots(2,2)\n",
"ax[0,0].set_xlim(-5,5)\n",
"j=ax[0,0].hist(z,bins=100,density=True)\n",
"j=ax[0,1].hist(samples_x,bins=50,density=True)\n",
"j=ax[1,1].hist(samples_y,bins=50,density=True)\n",
"ax[0,0].plot(np.linspace(-5,5,100),h.pdf(np.linspace(-5,5,100)))\n",
"\n",
"b=plt.hist(z,bins=100)"
]
},
{
"cell_type": "code",
"execution_count": 58,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"(array([1.000e+00, 0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00,\n",
" 0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00,\n",
" 0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00,\n",
" 0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00,\n",
" 0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00,\n",
" 0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00,\n",
" 0.000e+00, 1.000e+00, 1.000e+00, 6.000e+00, 7.000e+00, 1.800e+01,\n",
" 5.800e+01, 1.330e+02, 3.230e+02, 6.530e+02, 1.059e+03, 1.401e+03,\n",
" 1.615e+03, 1.410e+03, 1.124e+03, 7.740e+02, 5.150e+02, 3.100e+02,\n",
" 2.150e+02, 1.320e+02, 7.000e+01, 5.800e+01, 2.700e+01, 2.600e+01,\n",
" 2.200e+01, 7.000e+00, 5.000e+00, 7.000e+00, 4.000e+00, 2.000e+00,\n",
" 3.000e+00, 2.000e+00, 1.000e+00, 0.000e+00, 1.000e+00, 1.000e+00,\n",
" 0.000e+00, 0.000e+00, 2.000e+00, 0.000e+00, 0.000e+00, 1.000e+00,\n",
" 0.000e+00, 1.000e+00, 1.000e+00, 0.000e+00, 0.000e+00, 0.000e+00,\n",
" 0.000e+00, 0.000e+00, 1.000e+00, 0.000e+00, 1.000e+00, 0.000e+00,\n",
" 0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00,\n",
" 0.000e+00, 0.000e+00, 0.000e+00, 1.000e+00]),\n",
" array([-9.34383882, -9.13726736, -8.9306959 , -8.72412443, -8.51755297,\n",
" -8.31098151, -8.10441005, -7.89783858, -7.69126712, -7.48469566,\n",
" -7.2781242 , -7.07155274, -6.86498127, -6.65840981, -6.45183835,\n",
" -6.24526689, -6.03869543, -5.83212396, -5.6255525 , -5.41898104,\n",
" -5.21240958, -5.00583812, -4.79926665, -4.59269519, -4.38612373,\n",
" -4.17955227, -3.97298081, -3.76640934, -3.55983788, -3.35326642,\n",
" -3.14669496, -2.9401235 , -2.73355203, -2.52698057, -2.32040911,\n",
" -2.11383765, -1.90726619, -1.70069472, -1.49412326, -1.2875518 ,\n",
" -1.08098034, -0.87440888, -0.66783741, -0.46126595, -0.25469449,\n",
" -0.04812303, 0.15844844, 0.3650199 , 0.57159136, 0.77816282,\n",
" 0.98473428, 1.19130575, 1.39787721, 1.60444867, 1.81102013,\n",
" 2.01759159, 2.22416306, 2.43073452, 2.63730598, 2.84387744,\n",
" 3.0504489 , 3.25702037, 3.46359183, 3.67016329, 3.87673475,\n",
" 4.08330621, 4.28987768, 4.49644914, 4.7030206 , 4.90959206,\n",
" 5.11616352, 5.32273499, 5.52930645, 5.73587791, 5.94244937,\n",
" 6.14902083, 6.3555923 , 6.56216376, 6.76873522, 6.97530668,\n",
" 7.18187815, 7.38844961, 7.59502107, 7.80159253, 8.00816399,\n",
" 8.21473546, 8.42130692, 8.62787838, 8.83444984, 9.0410213 ,\n",
" 9.24759277, 9.45416423, 9.66073569, 9.86730715, 10.07387861,\n",
" 10.28045008, 10.48702154, 10.693593 , 10.90016446, 11.10673592,\n",
" 11.31330739]),\n",
" <a list of 100 Patch objects>)"
]
},
"execution_count": 58,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"b\n"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.6.4"
}
},
"nbformat": 4,
"nbformat_minor": 2
}

0 comments on commit 770e7f4

Please sign in to comment.