{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "This notebook briefly walks through fitting the four linear Bayesian regression calibrations for global core top $\\delta^{18}\\mathrm{O}_{\\mathrm{c}}$ and SSTs." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Notebook settings for graphics and load some libraries." ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "%matplotlib inline" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "import pandas as pd\n", "import pymc3 as pm" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Set paths to our core top data. These are also given as supplemental information in the paper." ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "coretop_path = '../data/parsed/coretops.csv'\n", "coretop_grid_path = '../data/parsed/coretops_grid.csv'" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Read coretop data and parse" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "This should be straightforward." ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "coretops_raw = pd.read_csv(coretop_path)\n", "coretops_grid = pd.read_csv(coretop_grid_path)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "coretops_raw['foramtype'] = coretops_raw['species'].astype('category')\n", "coretops_grid['foramtype'] = coretops_grid['species'].astype('category')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Take a quick look at the data to see what we're dealing with... (and make sure everything is there)" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[G. bulloides, G. ruber, N. incompta, N. pachyderma, T. sacculifer]\n", "Categories (5, object): [G. bulloides, G. ruber, N. incompta, N. pachyderma, T. sacculifer]" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" } ], "source": [ "coretops_raw['foramtype'].unique()" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "G. ruber 1002\n", "G. bulloides 635\n", "T. sacculifer 442\n", "N. pachyderma 425\n", "N. incompta 132\n", "Name: foramtype, dtype: int64" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" } ], "source": [ "coretops_raw['foramtype'].value_counts()" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
| \n", " | species | \n", "gridlat | \n", "gridlon | \n", "d18oc | \n", "d18osw | \n", "t_annual | \n", "t_seasonal | \n", "m1 | \n", "m2 | \n", "m3 | \n", "m4 | \n", "m5 | \n", "m6 | \n", "m7 | \n", "m8 | \n", "m9 | \n", "m10 | \n", "m11 | \n", "m12 | \n", "foramtype | \n", "
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 0 | \n", "G. bulloides | \n", "-57.5 | \n", "-53.5 | \n", "2.240000 | \n", "-0.297468 | \n", "2.100175 | \n", "3.783593 | \n", "1 | \n", "1 | \n", "1 | \n", "0 | \n", "0 | \n", "0 | \n", "0 | \n", "0 | \n", "0 | \n", "0 | \n", "0 | \n", "0 | \n", "G. bulloides | \n", "
| 1 | \n", "G. bulloides | \n", "-56.5 | \n", "-115.5 | \n", "2.280000 | \n", "-0.292819 | \n", "5.167257 | \n", "5.167257 | \n", "1 | \n", "1 | \n", "1 | \n", "1 | \n", "1 | \n", "1 | \n", "1 | \n", "1 | \n", "1 | \n", "1 | \n", "1 | \n", "1 | \n", "G. bulloides | \n", "
| 2 | \n", "G. bulloides | \n", "-56.5 | \n", "-56.5 | \n", "2.245000 | \n", "-0.281994 | \n", "3.576477 | \n", "4.666713 | \n", "1 | \n", "1 | \n", "1 | \n", "1 | \n", "1 | \n", "0 | \n", "0 | \n", "0 | \n", "0 | \n", "0 | \n", "0 | \n", "1 | \n", "G. bulloides | \n", "
| 3 | \n", "G. bulloides | \n", "-55.5 | \n", "-58.5 | \n", "2.146000 | \n", "-0.272228 | \n", "4.832495 | \n", "4.983686 | \n", "1 | \n", "1 | \n", "1 | \n", "1 | \n", "1 | \n", "1 | \n", "1 | \n", "1 | \n", "0 | \n", "1 | \n", "1 | \n", "1 | \n", "G. bulloides | \n", "
| 4 | \n", "G. bulloides | \n", "-55.5 | \n", "-57.5 | \n", "2.348571 | \n", "-0.266258 | \n", "4.671475 | \n", "5.106034 | \n", "1 | \n", "1 | \n", "1 | \n", "1 | \n", "1 | \n", "1 | \n", "0 | \n", "1 | \n", "0 | \n", "0 | \n", "1 | \n", "1 | \n", "G. bulloides | \n", "
| \n", " | species | \n", "gridlat | \n", "gridlon | \n", "d18oc | \n", "d18osw | \n", "t_annual | \n", "t_seasonal | \n", "m1 | \n", "m2 | \n", "m3 | \n", "... | \n", "m5 | \n", "m6 | \n", "m7 | \n", "m8 | \n", "m9 | \n", "m10 | \n", "m11 | \n", "m12 | \n", "foramtype | \n", "temp | \n", "
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 0 | \n", "G. bulloides | \n", "-57.5 | \n", "-53.5 | \n", "2.240000 | \n", "-0.297468 | \n", "2.100175 | \n", "3.783593 | \n", "1 | \n", "1 | \n", "1 | \n", "... | \n", "0 | \n", "0 | \n", "0 | \n", "0 | \n", "0 | \n", "0 | \n", "0 | \n", "0 | \n", "G. bulloides | \n", "2.100175 | \n", "
| 1 | \n", "G. bulloides | \n", "-56.5 | \n", "-115.5 | \n", "2.280000 | \n", "-0.292819 | \n", "5.167257 | \n", "5.167257 | \n", "1 | \n", "1 | \n", "1 | \n", "... | \n", "1 | \n", "1 | \n", "1 | \n", "1 | \n", "1 | \n", "1 | \n", "1 | \n", "1 | \n", "G. bulloides | \n", "5.167257 | \n", "
| 2 | \n", "G. bulloides | \n", "-56.5 | \n", "-56.5 | \n", "2.245000 | \n", "-0.281994 | \n", "3.576477 | \n", "4.666713 | \n", "1 | \n", "1 | \n", "1 | \n", "... | \n", "1 | \n", "0 | \n", "0 | \n", "0 | \n", "0 | \n", "0 | \n", "0 | \n", "1 | \n", "G. bulloides | \n", "3.576477 | \n", "
| 3 | \n", "G. bulloides | \n", "-55.5 | \n", "-58.5 | \n", "2.146000 | \n", "-0.272228 | \n", "4.832495 | \n", "4.983686 | \n", "1 | \n", "1 | \n", "1 | \n", "... | \n", "1 | \n", "1 | \n", "1 | \n", "1 | \n", "0 | \n", "1 | \n", "1 | \n", "1 | \n", "G. bulloides | \n", "4.832495 | \n", "
| 4 | \n", "G. bulloides | \n", "-55.5 | \n", "-57.5 | \n", "2.348571 | \n", "-0.266258 | \n", "4.671475 | \n", "5.106034 | \n", "1 | \n", "1 | \n", "1 | \n", "... | \n", "1 | \n", "1 | \n", "0 | \n", "1 | \n", "0 | \n", "0 | \n", "1 | \n", "1 | \n", "G. bulloides | \n", "4.671475 | \n", "
5 rows × 21 columns
\n", "| \n", " | species | \n", "gridlat | \n", "gridlon | \n", "d18oc | \n", "d18osw | \n", "t_annual | \n", "t_seasonal | \n", "m1 | \n", "m2 | \n", "m3 | \n", "... | \n", "m5 | \n", "m6 | \n", "m7 | \n", "m8 | \n", "m9 | \n", "m10 | \n", "m11 | \n", "m12 | \n", "foramtype | \n", "temp | \n", "
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 0 | \n", "G. bulloides | \n", "-57.5 | \n", "-53.5 | \n", "2.240000 | \n", "-0.297468 | \n", "2.100175 | \n", "3.783593 | \n", "1 | \n", "1 | \n", "1 | \n", "... | \n", "0 | \n", "0 | \n", "0 | \n", "0 | \n", "0 | \n", "0 | \n", "0 | \n", "0 | \n", "G. bulloides | \n", "3.783593 | \n", "
| 1 | \n", "G. bulloides | \n", "-56.5 | \n", "-115.5 | \n", "2.280000 | \n", "-0.292819 | \n", "5.167257 | \n", "5.167257 | \n", "1 | \n", "1 | \n", "1 | \n", "... | \n", "1 | \n", "1 | \n", "1 | \n", "1 | \n", "1 | \n", "1 | \n", "1 | \n", "1 | \n", "G. bulloides | \n", "5.167257 | \n", "
| 2 | \n", "G. bulloides | \n", "-56.5 | \n", "-56.5 | \n", "2.245000 | \n", "-0.281994 | \n", "3.576477 | \n", "4.666713 | \n", "1 | \n", "1 | \n", "1 | \n", "... | \n", "1 | \n", "0 | \n", "0 | \n", "0 | \n", "0 | \n", "0 | \n", "0 | \n", "1 | \n", "G. bulloides | \n", "4.666713 | \n", "
| 3 | \n", "G. bulloides | \n", "-55.5 | \n", "-58.5 | \n", "2.146000 | \n", "-0.272228 | \n", "4.832495 | \n", "4.983686 | \n", "1 | \n", "1 | \n", "1 | \n", "... | \n", "1 | \n", "1 | \n", "1 | \n", "1 | \n", "0 | \n", "1 | \n", "1 | \n", "1 | \n", "G. bulloides | \n", "4.983686 | \n", "
| 4 | \n", "G. bulloides | \n", "-55.5 | \n", "-57.5 | \n", "2.348571 | \n", "-0.266258 | \n", "4.671475 | \n", "5.106034 | \n", "1 | \n", "1 | \n", "1 | \n", "... | \n", "1 | \n", "1 | \n", "0 | \n", "1 | \n", "0 | \n", "0 | \n", "1 | \n", "1 | \n", "G. bulloides | \n", "5.106034 | \n", "
5 rows × 21 columns
\n", "