From 037fed135d9608f23ced80e37a8555b44898f5cb Mon Sep 17 00:00:00 2001 From: Morten Hjorth-Jensen Date: Sat, 11 Sep 2021 23:20:29 +0200 Subject: [PATCH] change of week35 --- doc/pub/week35/ipynb/week35.ipynb | 325 ++++++++++++++++-------------- 1 file changed, 178 insertions(+), 147 deletions(-) diff --git a/doc/pub/week35/ipynb/week35.ipynb b/doc/pub/week35/ipynb/week35.ipynb index 36440b523..9adee6d3b 100644 --- a/doc/pub/week35/ipynb/week35.ipynb +++ b/doc/pub/week35/ipynb/week35.ipynb @@ -402,10 +402,7 @@ { "cell_type": "code", "execution_count": null, - "metadata": { - "collapsed": false, - "editable": true - }, + "metadata": {}, "outputs": [], "source": [ "%matplotlib inline\n", @@ -956,10 +953,7 @@ { "cell_type": "code", "execution_count": null, - "metadata": { - "collapsed": false, - "editable": true - }, + "metadata": {}, "outputs": [], "source": [ "# matrix inversion to find beta\n", @@ -978,10 +972,7 @@ { "cell_type": "code", "execution_count": null, - "metadata": { - "collapsed": false, - "editable": true - }, + "metadata": {}, "outputs": [], "source": [ "fit = np.linalg.lstsq(X, Energies, rcond =None)[0]\n", @@ -998,10 +989,7 @@ { "cell_type": "code", "execution_count": null, - "metadata": { - "collapsed": false, - "editable": true - }, + "metadata": {}, "outputs": [], "source": [ "Masses['Eapprox'] = ytilde\n", @@ -1031,10 +1019,7 @@ { "cell_type": "code", "execution_count": null, - "metadata": { - "collapsed": false, - "editable": true - }, + "metadata": {}, "outputs": [], "source": [ "def R2(y_data, y_model):\n", @@ -1051,10 +1036,7 @@ { "cell_type": "code", "execution_count": null, - "metadata": { - "collapsed": false, - "editable": true - }, + "metadata": {}, "outputs": [], "source": [ "print(R2(Energies,ytilde))" @@ -1070,10 +1052,7 @@ { "cell_type": "code", "execution_count": null, - "metadata": { - "collapsed": false, - "editable": true - }, + "metadata": {}, "outputs": [], "source": [ "def MSE(y_data,y_model):\n", @@ -1093,10 +1072,7 @@ { "cell_type": "code", "execution_count": null, - "metadata": { - "collapsed": false, - "editable": true - }, + "metadata": {}, "outputs": [], "source": [ "def RelativeError(y_data,y_model):\n", @@ -1131,10 +1107,7 @@ { "cell_type": "code", "execution_count": null, - "metadata": { - "collapsed": false, - "editable": true - }, + "metadata": {}, "outputs": [], "source": [ "import os\n", @@ -1187,10 +1160,7 @@ { "cell_type": "code", "execution_count": null, - "metadata": { - "collapsed": false, - "editable": true - }, + "metadata": {}, "outputs": [], "source": [ "# equivalently in numpy\n", @@ -1262,10 +1232,7 @@ { "cell_type": "code", "execution_count": null, - "metadata": { - "collapsed": false, - "editable": true - }, + "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", @@ -1285,10 +1252,7 @@ { "cell_type": "code", "execution_count": null, - "metadata": { - "collapsed": false, - "editable": true - }, + "metadata": {}, "outputs": [], "source": [ "from sklearn.datasets import load_boston\n", @@ -1310,10 +1274,7 @@ { "cell_type": "code", "execution_count": null, - "metadata": { - "collapsed": false, - "editable": true - }, + "metadata": {}, "outputs": [], "source": [ "boston = pd.DataFrame(boston_dataset.data, columns=boston_dataset.feature_names)\n", @@ -1331,10 +1292,7 @@ { "cell_type": "code", "execution_count": null, - "metadata": { - "collapsed": false, - "editable": true - }, + "metadata": {}, "outputs": [], "source": [ "# check for missing values in all the columns\n", @@ -1351,10 +1309,7 @@ { "cell_type": "code", "execution_count": null, - "metadata": { - "collapsed": false, - "editable": true - }, + "metadata": {}, "outputs": [], "source": [ "# set the size of the figure\n", @@ -1375,10 +1330,7 @@ { "cell_type": "code", "execution_count": null, - "metadata": { - "collapsed": false, - "editable": true - }, + "metadata": {}, "outputs": [], "source": [ "# compute the pair wise correlation for all columns \n", @@ -1398,10 +1350,7 @@ { "cell_type": "code", "execution_count": null, - "metadata": { - "collapsed": false, - "editable": true - }, + "metadata": {}, "outputs": [], "source": [ "plt.figure(figsize=(20, 5))\n", @@ -1429,10 +1378,7 @@ { "cell_type": "code", "execution_count": null, - "metadata": { - "collapsed": false, - "editable": true - }, + "metadata": {}, "outputs": [], "source": [ "X = pd.DataFrame(np.c_[boston['LSTAT'], boston['RM']], columns = ['LSTAT','RM'])\n", @@ -1449,10 +1395,7 @@ { "cell_type": "code", "execution_count": null, - "metadata": { - "collapsed": false, - "editable": true - }, + "metadata": {}, "outputs": [], "source": [ "from sklearn.model_selection import train_test_split\n", @@ -1476,10 +1419,7 @@ { "cell_type": "code", "execution_count": null, - "metadata": { - "collapsed": false, - "editable": true - }, + "metadata": {}, "outputs": [], "source": [ "from sklearn.linear_model import LinearRegression\n", @@ -1518,10 +1458,7 @@ { "cell_type": "code", "execution_count": null, - "metadata": { - "collapsed": false, - "editable": true - }, + "metadata": {}, "outputs": [], "source": [ "# plotting the y_test vs y_pred\n", @@ -1649,10 +1586,7 @@ { "cell_type": "code", "execution_count": null, - "metadata": { - "collapsed": false, - "editable": true - }, + "metadata": {}, "outputs": [], "source": [ "import sklearn.linear_model as skl\n", @@ -1722,10 +1656,7 @@ { "cell_type": "code", "execution_count": null, - "metadata": { - "collapsed": false, - "editable": true - }, + "metadata": {}, "outputs": [], "source": [ "np.random.seed()\n", @@ -1748,10 +1679,7 @@ { "cell_type": "code", "execution_count": null, - "metadata": { - "collapsed": false, - "editable": true - }, + "metadata": {}, "outputs": [], "source": [ "import matplotlib.pyplot as plt\n", @@ -1802,10 +1730,7 @@ { "cell_type": "code", "execution_count": null, - "metadata": { - "collapsed": false, - "editable": true - }, + "metadata": {}, "outputs": [], "source": [ "# Common imports\n", @@ -2379,12 +2304,31 @@ }, { "cell_type": "code", - "execution_count": null, - "metadata": { - "collapsed": false, - "editable": true - }, - "outputs": [], + "execution_count": 1, + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "[[ 1. -1.]\n", + " [ 1. -1.]]\n", + "test U\n", + "[[0. 0.]\n", + " [0. 0.]]\n", + "test VT\n", + "[[0. 0.]\n", + " [0. 0.]]\n", + "[[-0.70710678 -0.70710678]\n", + " [-0.70710678 0.70710678]]\n", + "[2.00000000e+00 3.35470445e-17]\n", + "[[-0.70710678 0.70710678]\n", + " [ 0.70710678 0.70710678]]\n", + "[[-3.33066907e-16 4.44089210e-16]\n", + " [ 0.00000000e+00 2.22044605e-16]]\n" + ] + } + ], "source": [ "import numpy as np\n", "# SVD inversion\n", @@ -3182,12 +3126,20 @@ }, { "cell_type": "code", - "execution_count": null, - "metadata": { - "collapsed": false, - "editable": true - }, - "outputs": [], + "execution_count": 2, + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "-0.08839767027751376\n", + "3.8294285924714866\n", + "[[0.92932998 2.63805954]\n", + " [2.63805954 8.61477117]]\n" + ] + } + ], "source": [ "# Importing various packages\n", "import numpy as np\n", @@ -3216,12 +3168,20 @@ }, { "cell_type": "code", - "execution_count": null, - "metadata": { - "collapsed": false, - "editable": true - }, - "outputs": [], + "execution_count": 3, + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "0.0850713352585812\n", + "1.5846946541436007\n", + "[[1. 0.63837291]\n", + " [0.63837291 1. ]]\n" + ] + } + ], "source": [ "import numpy as np\n", "n = 100\n", @@ -3263,12 +3223,40 @@ }, { "cell_type": "code", - "execution_count": null, - "metadata": { - "collapsed": false, - "editable": true - }, - "outputs": [], + "execution_count": 4, + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "[[ 1.09669785 3.23795259]\n", + " [-0.22166894 -1.30593687]\n", + " [ 0.10192631 0.04245426]\n", + " [ 0.82011099 2.55486566]\n", + " [ 0.32105408 0.96706068]\n", + " [ 0.60361795 1.47703672]\n", + " [-1.87875598 -5.24764141]\n", + " [ 0.03658513 0.37688017]\n", + " [-0.57033315 -1.70980756]\n", + " [-0.30923424 -0.39286425]]\n", + " 0 1\n", + "0 1.096698 3.237953\n", + "1 -0.221669 -1.305937\n", + "2 0.101926 0.042454\n", + "3 0.820111 2.554866\n", + "4 0.321054 0.967061\n", + "5 0.603618 1.477037\n", + "6 -1.878756 -5.247641\n", + "7 0.036585 0.376880\n", + "8 -0.570333 -1.709808\n", + "9 -0.309234 -0.392864\n", + " 0 1\n", + "0 1.000000 0.990742\n", + "1 0.990742 1.000000\n" + ] + } + ], "source": [ "import numpy as np\n", "import pandas as pd\n", @@ -3297,12 +3285,49 @@ }, { "cell_type": "code", - "execution_count": null, - "metadata": { - "collapsed": false, - "editable": true - }, - "outputs": [], + "execution_count": 5, + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + " 0 1 2 3 4 5 6 7 \\\n", + "0 0.0 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 \n", + "1 0.0 0.079955 0.083252 0.080480 0.081776 0.083055 0.071202 0.072287 \n", + "2 0.0 0.083252 0.087624 0.083966 0.085810 0.087616 0.074576 0.076060 \n", + "3 0.0 0.080480 0.083966 0.085125 0.086868 0.088629 0.077734 0.079241 \n", + "4 0.0 0.081776 0.085810 0.086868 0.089002 0.091153 0.079679 0.081504 \n", + "5 0.0 0.083055 0.087616 0.088629 0.091153 0.093695 0.081663 0.083808 \n", + "6 0.0 0.071202 0.074576 0.077734 0.079679 0.081663 0.072591 0.074289 \n", + "7 0.0 0.072287 0.076060 0.079241 0.081504 0.083808 0.074289 0.076255 \n", + "8 0.0 0.073485 0.077660 0.080883 0.083467 0.086097 0.076120 0.078358 \n", + "9 0.0 0.074811 0.079394 0.082672 0.085583 0.088544 0.078096 0.080610 \n", + "10 0.0 0.061639 0.064871 0.068789 0.070813 0.072887 0.065314 0.067084 \n", + "11 0.0 0.062699 0.066259 0.070225 0.072518 0.074865 0.066904 0.068904 \n", + "12 0.0 0.063878 0.067775 0.071800 0.074368 0.076991 0.068629 0.070864 \n", + "13 0.0 0.065183 0.069422 0.073519 0.076366 0.079274 0.070494 0.072970 \n", + "14 0.0 0.066615 0.071207 0.075386 0.078520 0.081718 0.072505 0.075226 \n", + "\n", + " 8 9 10 11 12 13 14 \n", + "0 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 \n", + "1 0.073485 0.074811 0.061639 0.062699 0.063878 0.065183 0.066615 \n", + "2 0.077660 0.079394 0.064871 0.066259 0.067775 0.069422 0.071207 \n", + "3 0.080883 0.082672 0.068789 0.070225 0.071800 0.073519 0.075386 \n", + "4 0.083467 0.085583 0.070813 0.072518 0.074368 0.076366 0.078520 \n", + "5 0.086097 0.088544 0.072887 0.074865 0.076991 0.079274 0.081718 \n", + "6 0.076120 0.078096 0.065314 0.066904 0.068629 0.070494 0.072505 \n", + "7 0.078358 0.080610 0.067084 0.068904 0.070864 0.072970 0.075226 \n", + "8 0.080737 0.083272 0.068979 0.071034 0.073233 0.075583 0.078091 \n", + "9 0.083272 0.086096 0.071010 0.073303 0.075747 0.078348 0.081114 \n", + "10 0.068979 0.071010 0.059527 0.061166 0.062930 0.064825 0.066855 \n", + "11 0.071034 0.073303 0.061166 0.063004 0.064972 0.067075 0.069319 \n", + "12 0.073233 0.075747 0.062930 0.064972 0.067148 0.069465 0.071929 \n", + "13 0.075583 0.078348 0.064825 0.067075 0.069465 0.072001 0.074691 \n", + "14 0.078091 0.081114 0.066855 0.069319 0.071929 0.074691 0.077614 \n" + ] + } + ], "source": [ "# Common imports\n", "import numpy as np\n", @@ -4086,10 +4111,7 @@ { "cell_type": "code", "execution_count": null, - "metadata": { - "collapsed": false, - "editable": true - }, + "metadata": {}, "outputs": [], "source": [ "x = np.random.rand(100)\n", @@ -4111,10 +4133,7 @@ { "cell_type": "code", "execution_count": null, - "metadata": { - "collapsed": false, - "editable": true - }, + "metadata": {}, "outputs": [], "source": [ "import os\n", @@ -4309,10 +4328,7 @@ { "cell_type": "code", "execution_count": null, - "metadata": { - "collapsed": false, - "editable": true - }, + "metadata": {}, "outputs": [], "source": [ "from mpl_toolkits.mplot3d import Axes3D\n", @@ -4430,10 +4446,7 @@ { "cell_type": "code", "execution_count": null, - "metadata": { - "collapsed": false, - "editable": true - }, + "metadata": {}, "outputs": [], "source": [ "def FrankeFunction(x,y):\n", @@ -4491,7 +4504,25 @@ ] } ], - "metadata": {}, + "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.8" + } + }, "nbformat": 4, "nbformat_minor": 4 }