diff --git a/wk7/wk7.ipynb b/wk7/wk7.ipynb index 81f9806..c94124c 100644 --- a/wk7/wk7.ipynb +++ b/wk7/wk7.ipynb @@ -583,7 +583,7 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", - "version": "3.6.5" + "version": "3.5.2" } }, "nbformat": 4, diff --git a/wk8/wk8.ipynb b/wk8/wk8.ipynb index ee766c0..68a5eae 100644 --- a/wk8/wk8.ipynb +++ b/wk8/wk8.ipynb @@ -9,10 +9,494 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 10, + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "Using matplotlib backend: TkAgg\n", + "Populating the interactive namespace from numpy and matplotlib\n" + ] + } + ], + "source": [ + "%pylab\n", + "%matplotlib inline\n", + "import pandas as pd\n", + "from pandas import ExcelWriter\n", + "from pandas import ExcelFile\n", + "import numpy as np\n", + "import matplotlib.pyplot as plt\n", + "from scipy import stats" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "We start the analysis by importing the data. Note, we are only importing the data provided in Table 1 of the Appendix of Dong et al.1 as we are only analysing data from the references days and Day 62.\n", + "\n", + "**\\*\\* Double Check this with Murray \\*\\***\n", + "\n", + "\\*\\* Now we need to read the excel file in a nice way given that the spreadsheet has double headings in some places" + ] + }, + { + "cell_type": "code", + "execution_count": 11, "metadata": {}, "outputs": [], - "source": [] + "source": [ + "data = pd.read_excel('Dong_etal_2018_data.xlsx', sheetname='A1')" + ] + }, + { + "cell_type": "code", + "execution_count": 14, + "metadata": {}, + "outputs": [ + { + "name": "stderr", + "output_type": "stream", + "text": [ + "/usr/local/lib/python3.5/dist-packages/matplotlib/colors.py:494: RuntimeWarning: invalid value encountered in less\n", + " cbook._putmask(xa, xa < 0.0, -1)\n" + ] + }, + { + "data": { + "text/html": [ + "\n", + " \n", + "\n", + "
\n", + " \n", + " \n", + " \n", + " \n", + " | \n", + " Distance from event center (km)\n", + " \n", + " \n", + " \n", + " | \n", + " Unnamed: 1\n", + " \n", + " \n", + " \n", + " | \n", + " Unnamed: 5\n", + " \n", + " \n", + " \n", + " | \n", + " Unnamed: 9\n", + " \n", + " \n", + " \n", + " | \n", + " Unnamed: 13\n", + " \n", + " \n", + " |
---|---|---|---|---|---|
\n", + " Distance from event center (km)\n", + " \n", + " \n", + " \n", + " | \n", + " 1\n", + " \n", + " \n", + " \n", + " | \n", + " nan\n", + " \n", + " \n", + " \n", + " | \n", + " nan\n", + " \n", + " \n", + " \n", + " | \n", + " nan\n", + " \n", + " \n", + " \n", + " | \n", + " nan\n", + " \n", + " \n", + " |
\n", + " Unnamed: 1\n", + " \n", + " \n", + " \n", + " | \n", + " nan\n", + " \n", + " \n", + " \n", + " | \n", + " nan\n", + " \n", + " \n", + " \n", + " | \n", + " nan\n", + " \n", + " \n", + " \n", + " | \n", + " nan\n", + " \n", + " \n", + " \n", + " | \n", + " nan\n", + " \n", + " \n", + " |
\n", + " Unnamed: 5\n", + " \n", + " \n", + " \n", + " | \n", + " nan\n", + " \n", + " \n", + " \n", + " | \n", + " nan\n", + " \n", + " \n", + " \n", + " | \n", + " nan\n", + " \n", + " \n", + " \n", + " | \n", + " nan\n", + " \n", + " \n", + " \n", + " | \n", + " nan\n", + " \n", + " \n", + " |
\n", + " Unnamed: 9\n", + " \n", + " \n", + " \n", + " | \n", + " nan\n", + " \n", + " \n", + " \n", + " | \n", + " nan\n", + " \n", + " \n", + " \n", + " | \n", + " nan\n", + " \n", + " \n", + " \n", + " | \n", + " nan\n", + " \n", + " \n", + " \n", + " | \n", + " nan\n", + " \n", + " \n", + " |
\n", + " Unnamed: 13\n", + " \n", + " \n", + " \n", + " | \n", + " nan\n", + " \n", + " \n", + " \n", + " | \n", + " nan\n", + " \n", + " \n", + " \n", + " | \n", + " nan\n", + " \n", + " \n", + " \n", + " | \n", + " nan\n", + " \n", + " \n", + " \n", + " | \n", + " nan\n", + " \n", + " \n", + " |