Earthquakes_v2.ipynb

{"nbformat":4,"nbformat_minor":0,"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"},"colab":{"name":"Earthquakes_v2.ipynb","provenance":[]}},"cells":[{"cell_type":"code","metadata":{"id":"WoKaxKn3XtJ8","colab_type":"code","colab":{}},"source":["import math\n","from datetime import *\n","import pandas as pd\n","import csv\n","import matplotlib.pyplot as plt"],"execution_count":0,"outputs":[]},{"cell_type":"code","metadata":{"id":"_UdbolqlXtKC","colab_type":"code","colab":{}},"source":["PI = 3.14159\n","R = 6378.1 # radius of the earth in km\n","SECONDS_IN_A_DAY = 24 * 60 * 60\n","ORIGINAL_CSV_FILE = \"./data.csv\""],"execution_count":0,"outputs":[]},{"cell_type":"code","metadata":{"id":"dH_9qt1sXtKE","colab_type":"code","colab":{}},"source":["class Earthquake:\n","    \n","    \"\"\"This class defines the concept of an earthquake for all subsequent uses in this project.\n","       An earthquake has three main attributes:\n","       timestamp : a datetime object that represents the time of occurence\n","       coordinates : a dict containing the latitude and longitude values \n","       magnitude : a floating point value, the magnitude of the earthquake\"\"\"\n","    \n","    def __init__(self, timestamp, coordinates, magnitude):\n","        \n","        self.timestamp = timestamp\n","        self.coordinates = coordinates\n","        self.magnitude = magnitude\n","        \n","    def __str__(self):\n","        \n","        return \"Timestamp \\t\\t Lat. \\t\\t Long. \\t\\t Magnitude \\n{} \\t {} \\t {} \\t{}\".format(\n","                self.timestamp, \n","                self.coordinates[\"latitude\"], self.coordinates[\"longitude\"], \n","                self.magnitude)"],"execution_count":0,"outputs":[]},{"cell_type":"code","metadata":{"id":"L9vW_nSDXtKH","colab_type":"code","colab":{}},"source":["def degrees_to_radians(deg):\n","    return (math.pi / 180) * deg\n","\n","def spatial_distance(earthquake1, earthquake2):\n","    \n","    \"\"\"calculates and returns the spatial distance in kms between two given earthquakes,\n","       using haversine formula\"\"\"\n","\n","    coordinates1 = earthquake1.coordinates\n","    coordinates2 = earthquake2.coordinates\n","    \n","    latitude1 = degrees_to_radians(coordinates1[\"latitude\"])\n","    latitude2 = degrees_to_radians(coordinates2[\"latitude\"])\n","\n","    longitude1 = degrees_to_radians(coordinates1[\"longitude\"])\n","    longitude2 = degrees_to_radians(coordinates2[\"longitude\"])\n","\n","    delta_lat = latitude2 - latitude1\n","    delta_long = longitude2 - longitude1\n","\n","    a = math.cos(latitude1) * math.cos(latitude2) * math.pow(math.sin(delta_long/2), 2) + math.pow(math.sin(delta_lat/2), 2)\n","\n","    c = 2 * math.atan2(math.sqrt(a), math.sqrt(1-a))\n","\n","    d = R * c\n","\n","    return d"],"execution_count":0,"outputs":[]},{"cell_type":"code","metadata":{"id":"TBMtYXj0XtKJ","colab_type":"code","colab":{}},"source":["def absolute_temporal_distance(earthquake1, earthquake2):\n","    \n","    \"\"\"returns the absolute time difference in days between two given earthquakes\"\"\"\n","    \n","    timestamp1 = earthquake1.timestamp\n","    timestamp2 = earthquake2.timestamp\n","    \n","    duration = timestamp2 - timestamp1\n","    duration_in_seconds = duration.total_seconds()\n","    duration_in_days = duration_in_seconds / (SECONDS_IN_A_DAY)\n","    return abs(duration_in_days)"],"execution_count":0,"outputs":[]},{"cell_type":"code","metadata":{"id":"f7q7uBGEXtKM","colab_type":"code","colab":{},"outputId":"7948180d-61ec-4e94-8d20-4e9c197be0f4"},"source":["#EXAMPLE usage of the class Earthquake and the spatio-temporal distance functions\n","\n","timestamp_example1 = datetime(1999, 4, 29, 20, 0, 1)\n","coordinates_example1 = {\"latitude\" : 24.323, \"longitude\" : 73.422}\n","magnitude_example1 = 9.5\n","earthquake_example1 = Earthquake(timestamp_example1, coordinates_example1, magnitude_example1)\n","\n","timestamp_example2 = datetime(1999, 4, 28, 20, 2, 1)\n","coordinates_example2 = {\"latitude\" : 34.323, \"longitude\" : 93.422}\n","magnitude_example2 = 6.9\n","earthquake_example2 = Earthquake(timestamp_example2, coordinates_example2, magnitude_example2)\n","\n","print(earthquake_example1)\n","print(\"\\n\")\n","print(earthquake_example2)\n","print(\"\\n\")\n","print(\"Spatial distance :\", spatial_distance(earthquake_example1, earthquake_example2), \"km\")\n","print(\"Temporal distance :\", absolute_temporal_distance(earthquake_example1, earthquake_example2), \"days\")"],"execution_count":0,"outputs":[{"output_type":"stream","text":["Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1999-04-29 20:00:01 \t 24.323 \t 73.422 \t9.5\n","\n","\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1999-04-28 20:02:01 \t 34.323 \t 93.422 \t6.9\n","\n","\n","Spatial distance : 2231.409274019451 km\n","Temporal distance : 0.9986111111111111 days\n"],"name":"stdout"}]},{"cell_type":"code","metadata":{"id":"aD3OKtrPXtKP","colab_type":"code","colab":{}},"source":["def timestamp_from_string(date_, time_):\n","    \n","    \"\"\"Takes as input two parameters:\n","    date_ : yyyy-mm-dd\n","    time_ : hh:mm:ss\n","    At present, it ignores the floating point part of the seconds value\"\"\" \n","    \n","    tokens = date_.split(\"-\")\n","    year = int(tokens[0])\n","    month = int(tokens[1])\n","    day = int(tokens[2])\n","    \n","    tokens = time_.split(\":\")\n","    hours = int(tokens[0])\n","    minutes = int(tokens[1])\n","    seconds = int(float(tokens[2]))\n","    \n","    return datetime(year, month, day, hours, minutes, seconds)"],"execution_count":0,"outputs":[]},{"cell_type":"code","metadata":{"id":"-9nXBulEXtKS","colab_type":"code","colab":{},"outputId":"b0dd1cbd-9bb5-446f-b985-550a9145280e"},"source":["#EXAMPLE USAGE OF timestamp_from_string(date_, time_)\n","\n","timestamp_example = timestamp_from_string(\"1999-04-29\", \"20:02:34.4\")\n","print(\"Year: {} \\nMonth: {} \\nDay: {} \\nHours: {} \\nMinutes: {} \\nSeconds: {}\".format(\n","       timestamp_example.year, timestamp_example.month, timestamp_example.day,\n","       timestamp_example.hour, timestamp_example.minute, timestamp_example.second))"],"execution_count":0,"outputs":[{"output_type":"stream","text":["Year: 1999 \n","Month: 4 \n","Day: 29 \n","Hours: 20 \n","Minutes: 2 \n","Seconds: 34\n"],"name":"stdout"}]},{"cell_type":"code","metadata":{"id":"lIfDvZBrXtKU","colab_type":"code","colab":{}},"source":["def extract_from_csv(file):\n","    \n","    earthquakes = []\n","    with open(file) as f:\n","        data = csv.reader(f)\n","\n","        for row in data:\n","            if row[0] == None:\n","                return earthquakes\n","#             print(row)\n","            if row[0] == \"EVENT ID\": #Skipping the first row\n","                continue\n","            datestring = row[2]\n","            timestring = row[3]\n","            latitude = row[4]\n","            longitude = row[5]\n","            magnitude = float(row[10])\n","            \n","            timestamp = timestamp_from_string(datestring, timestring)\n","            coordinates = {\"latitude\" : float(latitude), \"longitude\" : float(longitude)}\n","            \n","            earthquake = Earthquake(timestamp, coordinates, magnitude)\n","            \n","            earthquakes.append(earthquake)\n","        return earthquakes"],"execution_count":0,"outputs":[]},{"cell_type":"code","metadata":{"id":"q5WLEtk5XtKX","colab_type":"code","colab":{}},"source":["\"\"\"Defining a spatio-temporal window for declustering\"\"\"\n","def spatial_window(earthquake):\n","    return math.exp(-1.024 + 0.804 * earthquake.magnitude)\n","\n","def temporal_window(earthquake):\n","    return math.exp(-2.870 + 1.235 * earthquake.magnitude)"],"execution_count":0,"outputs":[]},{"cell_type":"code","metadata":{"id":"xB25u41GXtKZ","colab_type":"code","colab":{}},"source":["def is_dependent(shock, mainshock, spatial_radius, temporal_radius):\n","        return(spatial_distance(shock, mainshock) <= spatial_radius and\n","           absolute_temporal_distance(shock, mainshock) <= temporal_radius)    "],"execution_count":0,"outputs":[]},{"cell_type":"code","metadata":{"id":"Vm2bwQvfXtKb","colab_type":"code","colab":{}},"source":["def decluster_earthquakes(earthquakes):\n","    \n","    declustered_earthquakes = []\n","    \n","    while True:\n","        \n","        earthquakes.sort(key=lambda x : x.magnitude, reverse=True)\n","        \n","        mainshock = earthquakes[0]\n","        \n","        if(len(earthquakes) == 1):\n","            declustered_earthquakes.append(mainshock)\n","            return declustered_earthquakes\n","        \n","        spatial_radius = spatial_window(mainshock)\n","        temporal_radius = temporal_window(mainshock)\n","        \n","        for earthquake in earthquakes[1:]:\n","            \n","            if(is_dependent(earthquake, mainshock, spatial_radius, temporal_radius)):\n","               earthquakes.remove(earthquake)\n","                \n","        declustered_earthquakes.append(mainshock)\n","        earthquakes.remove(mainshock)\n","        \n","        if(earthquakes == []):\n","            return declustered_earthquakes"],"execution_count":0,"outputs":[]},{"cell_type":"code","metadata":{"id":"CSFJ1XlHXtKe","colab_type":"code","colab":{}},"source":["def find_interevent_times(earthquakes):\n","    \n","    \"\"\"returns a list of interevent times in days\"\"\"\n","    \n","    intervent_times = []\n","    earthquakes.sort(key=lambda x : x.timestamp)\n","    l = len(earthquakes)\n","    \n","    print(len(earthquakes))\n","    for earthquake in earthquakes:\n","        print(earthquake)\n","\n","    for prev_event, next_event in zip(earthquakes[:l-1], earthquakes[1:]):\n","        diff = absolute_temporal_distance(prev_event, next_event)\n","        intervent_times.append(diff)\n","    return(intervent_times)\n"],"execution_count":0,"outputs":[]},{"cell_type":"code","metadata":{"id":"w4kog9sxXtKg","colab_type":"code","colab":{},"outputId":"46dcfce3-a731-4576-e571-4f63712bfc3d"},"source":["earthquakes = extract_from_csv(ORIGINAL_CSV_FILE) # FIRST STEP IN WORKFLOW\n","print(\"Events found: {}\".format(len(earthquakes)))\n","for earthquake in earthquakes :\n","    print(earthquake)"],"execution_count":0,"outputs":[{"output_type":"stream","text":["Events found: 86\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1905-04-04 00:49:59 \t 32.6359 \t 76.7878 \t7.8\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1908-08-20 09:53:24 \t 31.2276 \t 87.8353 \t6.9\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1911-10-14 23:23:59 \t 30.7593 \t 80.2794 \t6.5\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1914-10-09 02:39:14 \t 32.6041 \t 76.1809 \t6.2\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1915-12-03 02:39:24 \t 27.6997 \t 91.6401 \t6.5\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1916-08-28 06:39:41 \t 29.73 \t 80.7454 \t7.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1918-07-08 10:22:08 \t 24.599 \t 90.8528 \t7.2\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1923-09-09 22:03:46 \t 24.9372 \t 90.7207 \t6.9\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1924-10-08 20:32:55 \t 30.8815 \t 89.6514 \t6.5\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1930-07-02 21:03:43 \t 25.9288 \t 90.1765 \t7.2\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1934-01-15 08:43:25 \t 26.8847 \t 86.5885 \t8.2\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1934-12-15 01:57:42 \t 31.2496 \t 89.162 \t7.3\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1935-01-03 01:50:11 \t 30.7366 \t 88.3183 \t6.5\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1936-05-27 06:19:15 \t 28.3969 \t 83.311 \t6.8\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1941-01-21 12:41:46 \t 27.1659 \t 91.8596 \t6.6\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1943-10-23 17:23:21 \t 26.6373 \t 93.8518 \t7.2\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1944-10-17 18:36:57 \t 31.2656 \t 83.1143 \t6.8\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1945-06-04 12:08:59 \t 30.218 \t 80.0822 \t6.4\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1945-06-22 18:00:53 \t 32.5089 \t 76.2474 \t6.4\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1947-02-10 04:02:03 \t 31.5368 \t 85.2782 \t6.6\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1947-07-29 13:43:25 \t 28.576 \t 93.628 \t7.8\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1950-08-15 21:42:22 \t 25.3428 \t 92.9443 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1950-08-17 01:54:17 \t 28.4265 \t 94.8443 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1950-08-26 06:33:12 \t 27.2325 \t 94.9646 \t6.1\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1950-09-30 07:29:00 \t 28.8732 \t 94.4095 \t6.3\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1950-10-08 04:50:18 \t 28.5497 \t 94.4008 \t6.3\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1951-04-14 23:40:56 \t 28.2801 \t 93.861 \t6.5\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1951-11-17 04:46:02 \t 30.9084 \t 91.4983 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1951-11-18 09:35:54 \t 31.0559 \t 91.2612 \t7.7\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1951-12-26 10:07:00 \t 31.2439 \t 90.7841 \t6.3\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1952-08-17 16:02:14 \t 30.6482 \t 91.601 \t7.5\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1953-10-08 19:11:04 \t 32.2583 \t 82.8921 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1953-10-11 17:08:09 \t 32.1163 \t 82.7862 \t6.3\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1953-12-03 14:54:11 \t 31.1514 \t 85.7426 \t6.5\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1954-02-23 06:40:35 \t 27.6723 \t 91.6049 \t6.1\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1955-01-28 17:02:42 \t 33.1359 \t 82.5786 \t6.5\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1955-06-27 10:14:12 \t 32.2943 \t 78.5057 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1957-04-14 07:11:57 \t 30.5206 \t 84.3482 \t6.5\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1957-07-01 19:30:24 \t 24.3104 \t 93.8943 \t6.1\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1958-10-28 10:46:33 \t 30.4718 \t 84.5533 \t6.2\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1964-09-26 00:46:01 \t 29.8553 \t 80.452 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1964-10-21 23:09:18 \t 28.0827 \t 93.8315 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1965-01-12 13:32:25 \t 27.3829 \t 87.8286 \t6.1\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1966-03-06 02:15:55 \t 31.373 \t 80.5448 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1966-06-27 10:41:08 \t 29.4815 \t 80.8388 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1966-06-27 10:59:18 \t 29.7075 \t 80.9233 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1966-12-16 20:52:16 \t 29.5431 \t 80.8488 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1967-03-14 06:58:05 \t 28.4311 \t 94.2144 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1967-09-15 10:32:44 \t 27.3343 \t 91.877 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1971-05-22 20:03:30 \t 32.304 \t 92.1004 \t6.2\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1973-09-08 07:25:42 \t 33.1915 \t 86.8693 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1975-01-19 08:02:01 \t 32.3511 \t 78.487 \t6.4\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1975-05-05 05:18:47 \t 33.0852 \t 92.8762 \t6.3\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1977-11-18 05:20:10 \t 32.645 \t 88.5119 \t6.5\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1978-04-04 00:40:28 \t 32.7793 \t 82.2939 \t6.2\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1980-02-22 03:02:45 \t 30.5052 \t 88.6477 \t6.3\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1980-07-29 14:58:42 \t 29.58 \t 81.0481 \t6.1\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1980-11-19 19:00:51 \t 27.3596 \t 88.7832 \t6.1\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1982-01-23 17:37:29 \t 31.6787 \t 82.2471 \t6.1\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1986-06-20 17:12:46 \t 31.1279 \t 86.7831 \t6.1\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1987-05-18 01:53:53 \t 25.2727 \t 94.2182 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1988-08-20 23:09:10 \t 26.7162 \t 86.6247 \t6.5\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1991-10-19 21:23:16 \t 30.756 \t 78.7907 \t6.6\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1992-07-30 08:24:49 \t 29.5266 \t 90.1712 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1993-03-20 14:52:01 \t 28.9776 \t 87.4743 \t6.1\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1998-09-03 18:15:52 \t 27.8631 \t 86.9496 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1999-03-28 19:05:12 \t 30.511 \t 79.421 \t6.3\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1999-03-28 19:36:09 \t 30.31 \t 79.364 \t6.1\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2004-03-27 18:45:27 \t 33.9433 \t 89.2131 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2004-03-27 18:47:27 \t 33.9886 \t 89.1815 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2004-07-11 23:08:41 \t 30.719 \t 83.6663 \t6.2\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2005-04-07 20:04:40 \t 30.5171 \t 83.6553 \t6.1\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2008-01-09 08:26:48 \t 32.4042 \t 85.2551 \t6.2\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2008-08-25 13:22:01 \t 31.0612 \t 83.6521 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2008-10-06 08:30:44 \t 29.8445 \t 90.3787 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2009-09-21 08:53:06 \t 27.3693 \t 91.4597 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2011-02-04 13:53:46 \t 24.6094 \t 94.7347 \t6.4\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2011-09-18 12:40:49 \t 27.8039 \t 88.1536 \t6.5\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2015-04-25 06:11:26 \t 28.1302 \t 84.7168 \t6.9\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2015-04-25 06:15:23 \t 27.6342 \t 85.4762 \t6.1\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2015-04-25 06:45:21 \t 28.1603 \t 84.8433 \t6.5\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2015-04-26 07:09:09 \t 27.7365 \t 85.9788 \t6.6\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2015-05-12 07:05:18 \t 27.8014 \t 86.126 \t6.7\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2015-05-12 07:36:53 \t 27.5248 \t 86.1418 \t6.1\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2016-01-03 23:05:20 \t 24.8308 \t 93.6047 \t6.6\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2016-10-17 07:14:49 \t 32.8558 \t 94.8768 \t6.1\n"],"name":"stdout"}]},{"cell_type":"code","metadata":{"id":"XafRW9KOXtKj","colab_type":"code","colab":{},"outputId":"130c6c1a-4df6-41ca-f91a-2b3ae765f9ec"},"source":["declustered_earthquakes = decluster_earthquakes(earthquakes) #SECOND STEP IN WORKFLOW\n","print(\"i.i.d. events found: {}\".format(len(declustered_earthquakes)))\n","for earthquake in declustered_earthquakes:\n","    print(earthquake)"],"execution_count":0,"outputs":[{"output_type":"stream","text":["i.i.d. events found: 73\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1934-01-15 08:43:25 \t 26.8847 \t 86.5885 \t8.2\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1905-04-04 00:49:59 \t 32.6359 \t 76.7878 \t7.8\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1947-07-29 13:43:25 \t 28.576 \t 93.628 \t7.8\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1951-11-18 09:35:54 \t 31.0559 \t 91.2612 \t7.7\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1934-12-15 01:57:42 \t 31.2496 \t 89.162 \t7.3\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1918-07-08 10:22:08 \t 24.599 \t 90.8528 \t7.2\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1930-07-02 21:03:43 \t 25.9288 \t 90.1765 \t7.2\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1943-10-23 17:23:21 \t 26.6373 \t 93.8518 \t7.2\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1916-08-28 06:39:41 \t 29.73 \t 80.7454 \t7.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1908-08-20 09:53:24 \t 31.2276 \t 87.8353 \t6.9\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1923-09-09 22:03:46 \t 24.9372 \t 90.7207 \t6.9\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2015-04-25 06:11:26 \t 28.1302 \t 84.7168 \t6.9\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1936-05-27 06:19:15 \t 28.3969 \t 83.311 \t6.8\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1944-10-17 18:36:57 \t 31.2656 \t 83.1143 \t6.8\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2015-05-12 07:05:18 \t 27.8014 \t 86.126 \t6.7\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1941-01-21 12:41:46 \t 27.1659 \t 91.8596 \t6.6\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1947-02-10 04:02:03 \t 31.5368 \t 85.2782 \t6.6\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1991-10-19 21:23:16 \t 30.756 \t 78.7907 \t6.6\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2016-01-03 23:05:20 \t 24.8308 \t 93.6047 \t6.6\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1911-10-14 23:23:59 \t 30.7593 \t 80.2794 \t6.5\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1915-12-03 02:39:24 \t 27.6997 \t 91.6401 \t6.5\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1924-10-08 20:32:55 \t 30.8815 \t 89.6514 \t6.5\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1951-04-14 23:40:56 \t 28.2801 \t 93.861 \t6.5\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1953-12-03 14:54:11 \t 31.1514 \t 85.7426 \t6.5\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1955-01-28 17:02:42 \t 33.1359 \t 82.5786 \t6.5\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1957-04-14 07:11:57 \t 30.5206 \t 84.3482 \t6.5\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1977-11-18 05:20:10 \t 32.645 \t 88.5119 \t6.5\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1988-08-20 23:09:10 \t 26.7162 \t 86.6247 \t6.5\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2011-09-18 12:40:49 \t 27.8039 \t 88.1536 \t6.5\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1945-06-04 12:08:59 \t 30.218 \t 80.0822 \t6.4\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1945-06-22 18:00:53 \t 32.5089 \t 76.2474 \t6.4\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1975-01-19 08:02:01 \t 32.3511 \t 78.487 \t6.4\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2011-02-04 13:53:46 \t 24.6094 \t 94.7347 \t6.4\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1950-09-30 07:29:00 \t 28.8732 \t 94.4095 \t6.3\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1953-10-11 17:08:09 \t 32.1163 \t 82.7862 \t6.3\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1975-05-05 05:18:47 \t 33.0852 \t 92.8762 \t6.3\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1980-02-22 03:02:45 \t 30.5052 \t 88.6477 \t6.3\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1999-03-28 19:05:12 \t 30.511 \t 79.421 \t6.3\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1914-10-09 02:39:14 \t 32.6041 \t 76.1809 \t6.2\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1958-10-28 10:46:33 \t 30.4718 \t 84.5533 \t6.2\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1971-05-22 20:03:30 \t 32.304 \t 92.1004 \t6.2\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1978-04-04 00:40:28 \t 32.7793 \t 82.2939 \t6.2\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2004-07-11 23:08:41 \t 30.719 \t 83.6663 \t6.2\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2008-01-09 08:26:48 \t 32.4042 \t 85.2551 \t6.2\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1950-08-26 06:33:12 \t 27.2325 \t 94.9646 \t6.1\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1954-02-23 06:40:35 \t 27.6723 \t 91.6049 \t6.1\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1957-07-01 19:30:24 \t 24.3104 \t 93.8943 \t6.1\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1965-01-12 13:32:25 \t 27.3829 \t 87.8286 \t6.1\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1980-07-29 14:58:42 \t 29.58 \t 81.0481 \t6.1\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1980-11-19 19:00:51 \t 27.3596 \t 88.7832 \t6.1\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1982-01-23 17:37:29 \t 31.6787 \t 82.2471 \t6.1\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1986-06-20 17:12:46 \t 31.1279 \t 86.7831 \t6.1\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1993-03-20 14:52:01 \t 28.9776 \t 87.4743 \t6.1\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2005-04-07 20:04:40 \t 30.5171 \t 83.6553 \t6.1\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2016-10-17 07:14:49 \t 32.8558 \t 94.8768 \t6.1\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1950-08-15 21:42:22 \t 25.3428 \t 92.9443 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1950-08-17 01:54:17 \t 28.4265 \t 94.8443 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1955-06-27 10:14:12 \t 32.2943 \t 78.5057 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1964-09-26 00:46:01 \t 29.8553 \t 80.452 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1964-10-21 23:09:18 \t 28.0827 \t 93.8315 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1966-03-06 02:15:55 \t 31.373 \t 80.5448 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1966-06-27 10:41:08 \t 29.4815 \t 80.8388 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1966-12-16 20:52:16 \t 29.5431 \t 80.8488 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1967-03-14 06:58:05 \t 28.4311 \t 94.2144 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1967-09-15 10:32:44 \t 27.3343 \t 91.877 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1973-09-08 07:25:42 \t 33.1915 \t 86.8693 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1987-05-18 01:53:53 \t 25.2727 \t 94.2182 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1992-07-30 08:24:49 \t 29.5266 \t 90.1712 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1998-09-03 18:15:52 \t 27.8631 \t 86.9496 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2004-03-27 18:45:27 \t 33.9433 \t 89.2131 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2008-08-25 13:22:01 \t 31.0612 \t 83.6521 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2008-10-06 08:30:44 \t 29.8445 \t 90.3787 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2009-09-21 08:53:06 \t 27.3693 \t 91.4597 \t6.0\n"],"name":"stdout"}]},{"cell_type":"code","metadata":{"id":"O6mC8Ga4XtKl","colab_type":"code","colab":{},"outputId":"58f4bb6c-fd1e-4dba-891d-3d83dbc90f1d"},"source":["#THIRD STEP IN WORKFLOW: Plotting the data\n","fig = plt.figure()\n","fig.set_size_inches(7, 5)\n","plt.plot([earthquake.timestamp for earthquake in declustered_earthquakes], \n","        [earthquake.magnitude for earthquake in declustered_earthquakes], '.')\n","plt.xlabel(\"Time\")\n","plt.ylabel(\"Magnitude\")"],"execution_count":0,"outputs":[{"output_type":"stream","text":["/home/ankit/.local/lib/python3.6/site-packages/pandas/plotting/_converter.py:129: FutureWarning: Using an implicitly registered datetime converter for a matplotlib plotting method. The converter was registered by pandas on import. Future versions of pandas will require you to explicitly register matplotlib converters.\n","\n","To register the converters:\n","\t>>> from pandas.plotting import register_matplotlib_converters\n","\t>>> register_matplotlib_converters()\n","  warnings.warn(msg, FutureWarning)\n"],"name":"stderr"},{"output_type":"execute_result","data":{"text/plain":["Text(0, 0.5, 'Magnitude')"]},"metadata":{"tags":[]},"execution_count":80},{"output_type":"display_data","data":{"image/png":"iVBORw0KGgoAAAANSUhEUgAAAb8AAAE9CAYAAACba2yhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAcEUlEQVR4nO3dfZRcdZ3n8fc3aaKAIjHEAUESsqvuCDMo6YXGnfFhGB+ijrorexbEJ5R1dMZR19WRcTwedM/ujq7HWRlQhoPjwjHgHFFGdgZ8mmFWjppoNyoEEAlxEoGMhtgLgg9Jk+/+UTemqE26qzp1q27V7/06p0+q7r1V/f3dX9Ef7q1b34rMRJKkkiwZdgGSJA2a4SdJKo7hJ0kqjuEnSSqO4SdJKo7hJ0kqzsSwC+jVUUcdlatXrx52GZKkBpmZmbkvM1d2u/3Ihd/q1auZnp4edhmSpAaJiK29bO9pT0lScQw/SVJxDD9JUnEMP0lScQw/SVJxDD9JUnEMP0lScQw/SVJxDD/VbmbrLBffsJmZrbPDLkWSgBHs8KLRMrN1lnMu28CuuT0sm1jC+vOmWLtq+bDLklQ4j/xUqw1bdrJrbg97EnbP7WHDlp3DLkmSDD/Va2rNCpZNLGFpwCETS5has2LYJUmSpz1Vr7WrlrP+vCk2bNnJ1JoVnvKU1AiGn2q3dtVyQ09So3jaU5JUHMNPklQcw0+SVBzDT5JUHMNPklQcw0+SVBzDT5JUHMNPklQcw0+SVBzDT5JUHMNPklQcw0+SVBzDT5JUHMNPklQcw0+SVJxawy8i/lNE3BoRmyLiqoh4dMf6R0XEX0fE5ojYGBGr66xHkiSoMfwi4ljgrcBkZp4ELAXO6tjsDcBsZv5L4M+BD9ZVjyRJe9V92nMCODQiJoDDgHs71r8MuLy6fTVwRkREzTVJkgpXW/hl5j3Ah4FtwHbg/sz8UsdmxwI/rLafA+4HVtRVkyRJUO9pz+W0juxOAJ4IHB4Rr1rkc70xIqYjYnrHjh39LFOSVKA6T3v+LvCDzNyRmbuBzwHP7NjmHuBJANWp0ccBOzufKDMvzczJzJxcuXJljSVLkkpQZ/htA6Yi4rDqfbwzgNs7trkWeG11+0zgHzIza6xJkqRa3/PbSOsilpuAW6rfdWlEfCAiXlpt9glgRURsBt4BnF9XPZIk7RWjdqA1OTmZ09PTwy5DktQgETGTmZPdbm+HF0lScQw/SVJxDD9JUnEMP0lScQw/SVJxDD9JUnEMP0lScQw/SVJxDD9JUnEMP0lScQw/SVJxDD9JUnEMP0lScQw/SVJxDD9JUnGKDL+ZrbNcfMNmZrbODrsUNYyvDakME8MuYNBmts5yzmUb2DW3h2UTS1h/3hRrVy0fdllqAF8bUjmKO/LbsGUnu+b2sCdh99weNmzZOeyS1BC+NqRyFBd+U2tWsGxiCUsDDplYwtSaFcMuSQ3ha0MqR2TmsGvoyeTkZE5PTx/Uc8xsnWXDlp1MrVnhaS09gq8NaTRFxExmTna7fXHv+QGsXbXcP2zaL18bUhmKO+0pSZLhJ0kqjuEnSSqO4SdJKo7hJ0kqjuEnSSqO4ScdJPuBSqOnyM/5Sf1iP1BpNHnkJx0E+4FKo8nwkw6C/UCl0eRpT+kgrF21nPXnTdkPVBoxhp90kOwHKo0eT3tKkopj+EmSimP4SZKKY/hJkopj+EmSimP4SZKKY/hJkopj+EmSimP4SZKKY/hJkopj+EmSimP4SZKKY/hJkopj+EmSilNb+EXEUyPiO20/D0TE2zu2eU5E3N+2zfvqqkeSpL1q+z6/zLwDeDpARCwF7gGu2c+mN2bmS+qqQ5KkToM67XkGcFdmbh3Q75Mk6YAGFX5nAVcdYN3pEfHdiLg+Ik4cUD2SpILVHn4RsQx4KfCZ/ay+CViVmScDfwH8zQGe440RMR0R0zt27KivWElSEQZx5LcOuCkzf9S5IjMfyMwHq9vXAYdExFH72e7SzJzMzMmVK1fWX7EkaawNIvzO5gCnPCPi6IiI6vapVT07B1CTJKlgtV3tCRARhwPPA36/bdmbADLzEuBM4M0RMQf8HDgrM7POmiRJqjX8MvMhYEXHskvabl8EXFRnDZIkdbLDiySpOIafJKk4hp8kqTiGnySpOIafJKk4hp8kqTiGnySpOIafJKk4hp8kqTiGnySpOIafJKk4hp8kqTiGnySpOIafJKk4hp8kqTiGnwZmZussF9+wmZmts8MuRVLhuvoy24gI4BxgTWZ+ICKOB47OzG/WWp3GxszWWc65bAO75vawbGIJ68+bYu2q5cMuS1Khuj3y+xhwOnB2df+nwMW1VKSxtGHLTnbN7WFPwu65PWzYsnPYJUkqWFdHfsBpmXlKRHwbIDNnI2JZjXVpzEytWcGyiSXsntvDIRNLmFqzYtglSSpYt+G3OyKWAgkQESuBPbVVpbGzdtVy1p83xYYtO5las8JTnpKGqtvwuxC4BnhCRPxX4EzgvbVVpbG0dtVyQ09SI3QVfpm5PiJmgDOAAF6embfXWpkkSTWZN/wi4vFtd38MXNW+LjN/UldhkiTVZaEjvxla7/MFcDwwW90+EtgGnFBrdZIk1WDejzpk5gmZuQb4CvB7mXlUZq4AXgJ8aRAFSpLUb91+zm8qM6/beyczrweeWU9JkiTVq9urPe+NiPcCn6runwPcW09JkiTVq9sjv7OBlbQ+7nAN8AT2dXtRTZraC7Opde3PKNUqaXC6/ajDT4C31VyL2jS1F2ZT69qfUapV0mB1deQXETdExD90/tRdXMma2guzqXXtzyjVKmmwun3P751ttx8NvAKY63852qupvTCbWtf+jFKtkgYrMnNxD4z4Zmae2ud6FjQ5OZnT09OD/rVDMbN1tpG9MJta1/6MUq2SFi8iZjJzsuvtuwm/jk4vS4C1wIWZ+dTeSzw4JYWfJKk7vYZft6c92zu9zAE/AN7Qe3mSJA1ft+H365n5i/YFEfGoGuqRJKl23X7O7+v7WfaNfhYiSdKgLPStDkcDxwKHRsQzaJ32BDgCOKzm2iRJqsVCpz1fALwOOA74SNvynwLvqakmSZJqNW/4ZeblwOUR8YrM/OyAapIkqVYLnfZ8VWZ+ClgdEe/oXJ+ZH9nPwyRJarSFTnseXv37mLoLkSRpUBY67fmX1b/vH0w5kiTVr6vP+UXESuA/AqvbH5OZr6+nLEmS6tPth9w/D9wIfAV4uL5yJEmqX7fhd1hmvrvWSiRJGpBuO7z8bUS8qNZKJEkakG7D7220AvDnEfFARPw0Ih6oszBJkurSVfhl5mMzc0lmHpqZR1T3j5jvMRHx1Ij4TtvPAxHx9o5tIiIujIjNEXFzRJxyMIORJKkb3V7tub9Quh/Ympn7/Ub3zLwDeHr1+KXAPcA1HZutA55c/ZwGfLz6V5Kk2nR7wcvHgFOAW6r7vwFsAh4XEW/OzC8t8PgzgLsyc2vH8pcBV2TrG3U3RMSREXFMZm7vsq6x4LeNS9JgdRt+9wJvyMxbASLiacAHgD8GPgcsFH5nAVftZ/mxwA/b7t9dLSsm/Ga2znLOZRvYNbeHZRNLWH/elAEoSTXr9oKXp+wNPoDMvA34V5m5ZaEHRsQy4KXAZxZXIkTEGyNiOiKmd+zYsdinaaQNW3aya24PexJ2z+1hw5adwy5JksZet+F3a0R8PCKeXf18DLit+jb33Qs8dh1wU2b+aD/r7gGe1Hb/uGrZI2TmpZk5mZmTK1eu7LLk0TC1ZgXLJpawNOCQiSVMrVkx7JIkaex1e9rzdcAfAHuv1vwa8E5awffcBR57Nvs/5QlwLfCWiPg0rQtd7i/t/b61q5az/rwp3/OTpAGK1rUmNT15xOHANmBNZt5fLXsTQGZeEhEBXAS8EPgZcG5mTs/3nJOTkzk9Pe8mkqTCRMRMZk52u323H3V4MvDfgacBj967PDPXzPe4zHwIWNGx7JK22wn8YbfFSpLUD92+5/dJWp/Bm6N1mvMK4FN1FSVJUp26Db9DM/PvaZ0m3ZqZFwAvrq8sSZLq0+0FL7+MiCXAnRHxFlpXZPrt7pKkkdRLY+vDgLcCa4FXA6+tqyhJkurU1ZFfZn6ruvkgcG595TTPqLceG/X6JY2Hpv0tmjf8IuLa+dZn5kv7W06zjHrrsVGvX9J4aOLfooWO/E6n1XvzKmAjELVX1CD7az027AnrxajXL2k8NPFv0ULhdzTwPFpdWl4J/B1wVXufz3G2t/XY7rk9I9l6bNTrlzQemvi3qOsOL1Ufz7OB/wG8PzMvqrOwAxl0h5emnafu1ajXL2k81P23qNcOLwuGXxV6L6YVfKtp9eP8q8z8/xpQD4LtzSRJnfra3iwirgBOAq6jdbS36SDrkyRp6BZ6z+9VwEO0Puf31lYfaqB14Utm5hE11iZJUi3mDb/M7PZD8JIkjQzDTZJUHMNPklQcw0+SVBzDT0M1s3WWi2/YzMzW2bH6XZKarduvNJL6bpD9/prYW1DS8Hjkp6HZX7+/cfhdkprP8NPQ7O33tzSovd/fIH+XpObrurdnU9jebLwMsveofU6l8dXX9mZS3dauWj6wIBrk75LUbJ72lCQVx/CTJBXH8JMkFcfwkyQVx/CTJBXH8JMk1a5p7QX9qIMkqVZNbC/okZ8kqVZNbC9o+EmSatXE9oKe9pQk1WrtquWsP2+qUe0FDT9JUu2a1l7Q056SpOIYfpKk4hh+kqTiGH6SpOIYfpKk4hh+Gln9bJfUtNZLGg2+bhanCfvNjzpoJPWzXVITWy+p+XzdLE5T9ptHfhpJ/WyX1MTWS2o+XzeL05T9ZvhpJPWzXVITWy+p+XzdLE5T9ltk5lB+8WJNTk7m9PT0sMtQA8xsne1bu6R+PpfK4etmcerYbxExk5mTXW9v+EmSRl2v4edpT0lScQw/SVJxDD9JUnFqDb+IODIiro6I70XE7RFxesf650TE/RHxnernfXXWI0kS1P8h948CX8jMMyNiGXDYfra5MTNfUnMdkiT9Sm1HfhHxOOBZwCcAMnNXZv7fun7fsC2mXc+gW/w0oaXQOOlmf3ZuM99jnJ/m7IOm1NF0o7yf6jzyOwHYAXwyIk4GZoC3ZeZDHdudHhHfBe4F3pmZt9ZYUy0W065n0C1+mtJSaFx0sz87t3nfS07kA397634f4/w0Zx80pY6mG/X9VOd7fhPAKcDHM/MZwEPA+R3b3ASsysyTgb8A/mZ/TxQRb4yI6YiY3rFjR40lL85i2vUMusVPU1oKjYtu9mfnNtdv2n7Axzg/zdkHTamj6UZ9P9UZfncDd2fmxur+1bTC8Fcy84HMfLC6fR1wSEQc1flEmXlpZk5m5uTKlStrLHlxFtOuZ9AtfprSUmhcdLM/O7dZd9IxB3yM89OcfdCUOppu1PdTrR1eIuJG4LzMvCMiLgAOz8x3ta0/GvhRZmZEnEorIFflPEU1tcPLYtr1DLo1kq2Y+qub/dm5zXyPcX6asw+aUkfTNWk/Naq9WUQ8HbgMWAZsAc4F/gNAZl4SEW8B3gzMAT8H3pGZX5/vOZsafpKk4WlU+NXB8JMkdbK3pyRJCzD8JEnFMfwkScUx/CRJxTH8JEnFMfxUnCs3buPVn9jIlRu3DbuUooxyH0g90jjMZd3f6iA1ypUbt/Gea24B4MY77wPglacdP8ySijDqfSC1z7jMpUd+Ksr1m7bPe1/1GPU+kNpnXObS8FNR1p10zLz3VY9R7wOpfcZlLu3wouJcuXEb12/azrqTjvGU5wA1qQ+kDk4T59L2ZpKk4tjeTJKkBRh+kqTiGH6SpOIYfpKk4hh+kqTiGH7SInTT3mkcWkD128HskxL3Z4ljHhTbm0k96qa907i0gOqng9knJe7PEsc8SB75ST3qpr3TuLSA6qeD2Scl7s8SxzxIhp/Uo27aO41LC6h+Oph9UuL+LHHMg2SHF2kRumnv1MQWUMN2MPukxP1Z4pgXy/ZmkqTi2N5MkqQFGH6SpOIYfpKk4hh+kqTiGH6SpOIYfgNkqyL1oumvl7rqa/q4R00v+7OkfW97swGxVZF60fTXS131NX3co6aX/VnavvfIb0BsVaReNP31Uld9TR/3qOllf5a27w2/AbFVkXrR9NdLXfU1fdyjppf9Wdq+t8PLANmqSL1o+uulrvqaPu5R08v+HOV9b3szSVJxbG8mSdICDD9JUnEMP0lScQw/SVJxDD9JUnHs8CL1qNfLwbvdvu7LzEf5MvZ2M1tn+exNdxPAvzvlOIBFjau0j2o0sa5h1mT4ST3otQVUt9vX3VpqXFpXzWyd5exLv8Guh1sf0frr6R+yJIK5h3sbV2nt2ZpY17Br8rSn1INeW0B1u33draXGpXXVhi072f3wvs8mzz2c7F7EuEprz9bEuoZdk+En9aDXFlDdbl93a6lxaV01tWYFhyyNX92fWBocsohxldaerYl1DbsmO7xIPfI9v+HyPb/FaWJd/azJ9maSpOLY3kySpAUYfpKk4hh+kqTiGH6SpOLUGn4RcWREXB0R34uI2yPi9I71EREXRsTmiLg5Ik6psx5JkqD+Di8fBb6QmWdGxDLgsI7164AnVz+nAR+v/pUa68qN27h+03bWnXQMrzzt+J4fP4xLzjs/HrB21fJfLbvvp78EYOVjH8WJT3wcsz/b1Zfamnhp/YH0o9Z+jrdfz7XY5+l8jY/SXHartvCLiMcBzwJeB5CZu4BdHZu9DLgiW5+32FAdKR6Tmdvrqks6GFdu3MZ7rrkFgBvvvA+gpwAcRkunzpZgn5m5mwt+70QuuHbTr5a1WxIcdG3Dbl3Vi37U2s/x9uu5Fvs8na/xbTsf4n99459GYi57UedpzxOAHcAnI+LbEXFZRBzesc2xwA/b7t9dLXuEiHhjRExHxPSOHTvqq1hawPWbts97fyHDaOnU2RJs99wert+0/RHL2vWjtmG3rupFP2rt53j79VyLfZ7O1/QXbv3nkZnLXtQZfhPAKcDHM/MZwEPA+Yt5osy8NDMnM3Ny5cqV/axR6sm6k46Z9/5ChtHSqbMl2CETS1h30jGPWNZuSR9qG3brql70o9Z+jrdfz7XY5+l8Tb/wxKNHZi57UVuHl4g4GtiQmaur+78NnJ+ZL27b5i+Bf8zMq6r7dwDPme+0px1eNGy+59f97xyV94l8z++RRvE9v0a1N4uIG4HzMvOOiLgAODwz39W2/sXAW4AX0brQ5cLMPHW+5zT8JEmdeg2/uq/2/CNgfXWl5xbg3Ih4E0BmXgJcRyv4NgM/A86tuR5JkuoNv8z8DtCZxJe0rU/gD+usQZKkTnZ4kSQVx/CTJBXH8JMkFcfwkyQVx/CTJBXH8JMkFcfwkyQVp9YOL3WIiB20+oTeN+xahugoHH+p4y957OD4Sx7/QmNflZldN38eufADiIjpXtrYjBvHX+74Sx47OP6Sx9/vsXvaU5JUHMNPklScUQ2/S4ddwJA5/nKVPHZw/CWPv69jH8n3/CRJOhijeuQnSdKiNSb8IuKvIuLHEbGpbdnJEfGNiLglIv53RBzRtu5PImJzRNwRES9oW/7CatnmiDh/0ONYjF7GHhHPi4iZavlMRPxO22PWVss3R8SFERHDGE+vep37av3xEfFgRLyzbdnIzT0s6rX/m9W6W6v1j66Wj9z89/jaPyQiLq+W3x4Rf9L2mFGd+ydFxA0RcVs1n2+rlj8+Ir4cEXdW/y6vlkc1t5sj4uaIOKXtuV5bbX9nRLx2WGPq1iLGfk415lsi4usRcXLbc/U+/5nZiB/gWcApwKa2Zd8Cnl3dfj3wX6rbTwO+CzwKOAG4C1ha/dwFrAGWVds8bdhj6/PYnwE8sbp9EnBP22O+CUwBAVwPrBv22Po9/rb1VwOfAd5Z3R/JuV/E/E8ANwMnV/dXAEtHdf57HPsrgU9Xtw8D/glYPeJzfwxwSnX7scD3q79vHwLOr5afD3ywuv2iam6jmuuN1fLH0/rC8McDy6vby4c9vj6P/Zl7xwSsaxv7oua/MUd+mflV4Ccdi58CfLW6/WXgFdXtl9H6j+CXmfkDWt8Ef2r1szkzt2TmLuDT1baN1svYM/PbmXlvtfxW4NCIeFREHAMckZkbsvWKuAJ4ef3VH7we556IeDnwA1rj32sk5x56Hv/zgZsz87vVY3dm5sOjOv89jj2BwyNiAjgU2AU8wGjP/fbMvKm6/VPgduBYWvVfXm12Ofvm8mXAFdmyATiymvsXAF/OzJ9k5iyt/fbCAQ6lZ72OPTO/Xo0NYANwXHV7UfPfmPA7gFvZN4h/Dzypun0s8MO27e6ulh1o+Sg60NjbvQK4KTN/SWucd7etG+WxwwHGHxGPAd4NvL9j+3Gaezjw/D8FyIj4YkTcFBF/XC0fp/k/0NivptXdaTuwDfhwZv6EMZn7iFhN68zORuDXMnN7teqfgV+rbo/l374ux97uDbSOgGGRY296+L0e+IOImKF1WLxryPUM0rxjj4gTgQ8Cvz+E2gbhQOO/APjzzHxwWIUNyIHGPwH8FnBO9e+/jYgzhlNibQ409lOBh4En0nq74z9HxJrhlNhf1f/UfRZ4e2Y+0L6uOpIf28vyex17RDyXVvi9+2B+78TBPLhumfk9Wqd5iIinAC+uVt3DI4+EjquWMc/ykTLP2ImI44BrgNdk5l3V4nvYdxoARnjsMO/4TwPOjIgPAUcCeyLiF8AMYzL3MO/47wa+mpn3Veuuo/We2acYk/mfZ+yvBL6QmbuBH0fE14BJWv/XP7JzHxGH0Prjvz4zP1ct/lFEHJOZ26vTmj+ulh/ob989wHM6lv9jnXX3Q49jJyJ+E7iM1vvZO6vF8+XBATX6yC8inlD9uwR4L3BJtepa4Kzqva4TgCfTerP/W8CTI+KEiFgGnFVtO3IONPaIOBL4O1pvCH9t7/bVaYIHImKqusrvNcDnB154nxxo/Jn525m5OjNXA/8T+G+ZeRFjNPcw72v/i8BvRMRh1XtfzwZuG6f5n2fs24DfqdYdTuuCj+8xwnNfzdUngNsz8yNtq64F9l6x+Vr2zeW1wGuqqz6ngPuruf8i8PyIWF5dHfn8allj9Tr2iDge+Bzw6sz8ftv2i5v/YV/x03blz1W0zuXvpvV/t28A3kbrCqDvA39G9aH8avs/pXWFzx20XdVG62qo71fr/nTY4+r32Gn9MXgI+E7bzxOqdZPApmrsF7Xvryb/9Dr3bY+7gOpqz1Gd+8WMH3gVrffFNgEfals+cvPf42v/MbSu8L0VuA141xjM/W/ROq13c9t/zy+idRXv3wN3Al8BHl9tH8DF1ThvASbbnuv1tC7+2wycO+yx1TD2y4DZtm2nD2b+7fAiSSpOo097SpJUB8NPklQcw0+SVBzDT5JUHMNPklScRn/IXSpJROy9xBvgaFrdTHZU93+Wmc8cSmHSGPKjDlIDRcQFwIOZ+eFh1yKNI097SiMgIh6s/n1ORPyfiPh8RGyJiD+rvufsm9X3nP2LaruVEfHZiPhW9fNvhjsCqVkMP2n0nAy8Cfh14NXAUzLzVFodMP6o2uajtBqA/2ta3/5x2TAKlZrK9/yk0fOtrL7yJSLuAr5ULb8FeG51+3eBp8W+L3M/IiIek+P/bRhSVww/afT8su32nrb7e9j33/QSYCozfzHIwqRR4WlPaTx9iX2nQImIpw+xFqlxDD9pPL0VmIyImyPiNlrvEUqq+FEHSVJxPPKTJBXH8JMkFcfwkyQVx/CTJBXH8JMkFcfwkyQVx/CTJBXH8JMkFef/AcFbRaBHORXeAAAAAElFTkSuQmCC\n","text/plain":["<Figure size 504x360 with 1 Axes>"]},"metadata":{"tags":[],"needs_background":"light"}}]},{"cell_type":"code","metadata":{"collapsed":true,"id":"f3nmrmBUXtKo","colab_type":"code","colab":{},"outputId":"201b65ef-da36-4c2c-fbb1-a03be7c38229"},"source":["interevent_times = find_interevent_times(declustered_earthquakes) #STEP 4 in WORKFLOW"],"execution_count":0,"outputs":[{"output_type":"stream","text":["73\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1905-04-04 00:49:59 \t 32.6359 \t 76.7878 \t7.8\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1908-08-20 09:53:24 \t 31.2276 \t 87.8353 \t6.9\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1911-10-14 23:23:59 \t 30.7593 \t 80.2794 \t6.5\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1914-10-09 02:39:14 \t 32.6041 \t 76.1809 \t6.2\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1915-12-03 02:39:24 \t 27.6997 \t 91.6401 \t6.5\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1916-08-28 06:39:41 \t 29.73 \t 80.7454 \t7.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1918-07-08 10:22:08 \t 24.599 \t 90.8528 \t7.2\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1923-09-09 22:03:46 \t 24.9372 \t 90.7207 \t6.9\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1924-10-08 20:32:55 \t 30.8815 \t 89.6514 \t6.5\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1930-07-02 21:03:43 \t 25.9288 \t 90.1765 \t7.2\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1934-01-15 08:43:25 \t 26.8847 \t 86.5885 \t8.2\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1934-12-15 01:57:42 \t 31.2496 \t 89.162 \t7.3\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1936-05-27 06:19:15 \t 28.3969 \t 83.311 \t6.8\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1941-01-21 12:41:46 \t 27.1659 \t 91.8596 \t6.6\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1943-10-23 17:23:21 \t 26.6373 \t 93.8518 \t7.2\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1944-10-17 18:36:57 \t 31.2656 \t 83.1143 \t6.8\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1945-06-04 12:08:59 \t 30.218 \t 80.0822 \t6.4\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1945-06-22 18:00:53 \t 32.5089 \t 76.2474 \t6.4\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1947-02-10 04:02:03 \t 31.5368 \t 85.2782 \t6.6\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1947-07-29 13:43:25 \t 28.576 \t 93.628 \t7.8\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1950-08-15 21:42:22 \t 25.3428 \t 92.9443 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1950-08-17 01:54:17 \t 28.4265 \t 94.8443 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1950-08-26 06:33:12 \t 27.2325 \t 94.9646 \t6.1\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1950-09-30 07:29:00 \t 28.8732 \t 94.4095 \t6.3\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1951-04-14 23:40:56 \t 28.2801 \t 93.861 \t6.5\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1951-11-18 09:35:54 \t 31.0559 \t 91.2612 \t7.7\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1953-10-11 17:08:09 \t 32.1163 \t 82.7862 \t6.3\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1953-12-03 14:54:11 \t 31.1514 \t 85.7426 \t6.5\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1954-02-23 06:40:35 \t 27.6723 \t 91.6049 \t6.1\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1955-01-28 17:02:42 \t 33.1359 \t 82.5786 \t6.5\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1955-06-27 10:14:12 \t 32.2943 \t 78.5057 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1957-04-14 07:11:57 \t 30.5206 \t 84.3482 \t6.5\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1957-07-01 19:30:24 \t 24.3104 \t 93.8943 \t6.1\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1958-10-28 10:46:33 \t 30.4718 \t 84.5533 \t6.2\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1964-09-26 00:46:01 \t 29.8553 \t 80.452 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1964-10-21 23:09:18 \t 28.0827 \t 93.8315 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1965-01-12 13:32:25 \t 27.3829 \t 87.8286 \t6.1\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1966-03-06 02:15:55 \t 31.373 \t 80.5448 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1966-06-27 10:41:08 \t 29.4815 \t 80.8388 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1966-12-16 20:52:16 \t 29.5431 \t 80.8488 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1967-03-14 06:58:05 \t 28.4311 \t 94.2144 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1967-09-15 10:32:44 \t 27.3343 \t 91.877 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1971-05-22 20:03:30 \t 32.304 \t 92.1004 \t6.2\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1973-09-08 07:25:42 \t 33.1915 \t 86.8693 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1975-01-19 08:02:01 \t 32.3511 \t 78.487 \t6.4\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1975-05-05 05:18:47 \t 33.0852 \t 92.8762 \t6.3\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1977-11-18 05:20:10 \t 32.645 \t 88.5119 \t6.5\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1978-04-04 00:40:28 \t 32.7793 \t 82.2939 \t6.2\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1980-02-22 03:02:45 \t 30.5052 \t 88.6477 \t6.3\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1980-07-29 14:58:42 \t 29.58 \t 81.0481 \t6.1\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1980-11-19 19:00:51 \t 27.3596 \t 88.7832 \t6.1\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1982-01-23 17:37:29 \t 31.6787 \t 82.2471 \t6.1\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1986-06-20 17:12:46 \t 31.1279 \t 86.7831 \t6.1\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1987-05-18 01:53:53 \t 25.2727 \t 94.2182 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1988-08-20 23:09:10 \t 26.7162 \t 86.6247 \t6.5\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1991-10-19 21:23:16 \t 30.756 \t 78.7907 \t6.6\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1992-07-30 08:24:49 \t 29.5266 \t 90.1712 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1993-03-20 14:52:01 \t 28.9776 \t 87.4743 \t6.1\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1998-09-03 18:15:52 \t 27.8631 \t 86.9496 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","1999-03-28 19:05:12 \t 30.511 \t 79.421 \t6.3\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2004-03-27 18:45:27 \t 33.9433 \t 89.2131 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2004-07-11 23:08:41 \t 30.719 \t 83.6663 \t6.2\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2005-04-07 20:04:40 \t 30.5171 \t 83.6553 \t6.1\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2008-01-09 08:26:48 \t 32.4042 \t 85.2551 \t6.2\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2008-08-25 13:22:01 \t 31.0612 \t 83.6521 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2008-10-06 08:30:44 \t 29.8445 \t 90.3787 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2009-09-21 08:53:06 \t 27.3693 \t 91.4597 \t6.0\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2011-02-04 13:53:46 \t 24.6094 \t 94.7347 \t6.4\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2011-09-18 12:40:49 \t 27.8039 \t 88.1536 \t6.5\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2015-04-25 06:11:26 \t 28.1302 \t 84.7168 \t6.9\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2015-05-12 07:05:18 \t 27.8014 \t 86.126 \t6.7\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2016-01-03 23:05:20 \t 24.8308 \t 93.6047 \t6.6\n","Timestamp \t\t Lat. \t\t Long. \t\t Magnitude \n","2016-10-17 07:14:49 \t 32.8558 \t 94.8768 \t6.1\n"],"name":"stdout"}]},{"cell_type":"code","metadata":{"id":"8V03tubxXtKq","colab_type":"code","colab":{},"outputId":"12ff3e94-f075-4424-b28a-0024ab7e7de2"},"source":["print(len(interevent_times))\n","print(interevent_times)"],"execution_count":0,"outputs":[{"output_type":"stream","text":["72\n","[1234.3773726851853, 1150.5629050925927, 1090.1355902777777, 420.00011574074074, 269.16686342592595, 679.1544791666666, 1889.4872453703704, 394.9369097222222, 2093.021388888889, 1292.4859027777777, 333.71825231481483, 529.1816319444445, 1700.265636574074, 1005.1955439814815, 360.0511111111111, 229.7305787037037, 18.244375, 597.4174768518518, 169.40372685185184, 1113.3326041666667, 1.1749421296296296, 9.19369212962963, 35.03875, 196.6749537037037, 217.4131712962963, 693.3140625, 52.90696759259259, 81.65722222222222, 339.43202546296294, 149.71631944444445, 656.8734375, 78.5128125, 483.63621527777775, 2159.582962962963, 25.93283564814815, 82.59938657407407, 417.53020833333335, 113.3508449074074, 172.42439814814816, 87.42070601851852, 185.1490625, 1345.3963657407407, 839.47375, 498.0252199074074, 105.88664351851853, 928.0009606481482, 136.8057638888889, 689.0988078703704, 158.4971875, 113.16815972222223, 429.94210648148146, 1608.9828356481482, 331.36188657407405, 460.88561342592595, 1154.9264583333334, 284.45940972222223, 233.26888888888888, 1993.1415625, 206.03425925925927, 1825.9862847222223, 106.18280092592593, 269.87221064814815, 1006.5153703703704, 229.20501157407406, 41.797719907407405, 350.0155324074074, 501.20879629629627, 225.94934027777776, 1314.7295949074073, 17.037407407407407, 236.66668981481482, 287.3399189814815]\n"],"name":"stdout"}]},{"cell_type":"code","metadata":{"id":"nMqhh_jjXtKs","colab_type":"code","colab":{"base_uri":"https://localhost:8080/","height":34},"outputId":"f6a5e88a-a1f4-4bb7-cc28-e259dd7a2fe4","executionInfo":{"status":"ok","timestamp":1574348462503,"user_tz":-330,"elapsed":2597,"user":{"displayName":"SUBHAM KUMAR DASH","photoUrl":"https://lh3.googleusercontent.com/a-/AAuE7mCoM7QRiLbW5V1SqI-lIlCj9rythg-c1ahg8C_LbQ=s64","userId":"11090621726530319638"}}},"source":["1234/31\n"],"execution_count":2,"outputs":[{"output_type":"execute_result","data":{"text/plain":["39.806451612903224"]},"metadata":{"tags":[]},"execution_count":2}]},{"cell_type":"code","metadata":{"id":"DZGBG1WAemw6","colab_type":"code","colab":{}},"source":[""],"execution_count":0,"outputs":[]}]}