{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# MJO and vectors"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"* Wenchang Yang (wenchang@princeton.edu)\n",
"* Department of Geosciences, Princeton University"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {
"ExecuteTime": {
"end_time": "2018-07-31T18:27:03.484048Z",
"start_time": "2018-07-31T18:27:03.478770Z"
}
},
"outputs": [],
"source": [
"import matplotlib.pyplot as plt\n",
"import pandas as pd\n",
"import seaborn as sns\n",
"\n",
"from scipy import stats"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {
"ExecuteTime": {
"end_time": "2018-07-31T18:27:04.885148Z",
"start_time": "2018-07-31T18:27:04.878113Z"
}
},
"outputs": [],
"source": [
"%matplotlib inline\n",
"plt.rcParams['figure.figsize'] = [10, 7]"
]
},
{
"cell_type": "code",
"execution_count": 62,
"metadata": {
"ExecuteTime": {
"end_time": "2018-07-31T19:25:11.226121Z",
"start_time": "2018-07-31T19:25:11.215573Z"
}
},
"outputs": [],
"source": [
"def nan_mode(x):\n",
" m = stats.mode(x)[0]\n",
" if not m:\n",
" m = np.NaN\n",
" return m"
]
},
{
"cell_type": "code",
"execution_count": 64,
"metadata": {
"ExecuteTime": {
"end_time": "2018-07-31T19:26:38.086457Z",
"start_time": "2018-07-31T19:26:38.081677Z"
}
},
"outputs": [],
"source": [
"dfs = dict()"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {
"ExecuteTime": {
"end_time": "2018-07-31T18:26:38.458048Z",
"start_time": "2018-07-31T18:26:38.325151Z"
}
},
"outputs": [],
"source": [
"df = pd.read_csv('MJO_Abundance.csv', index_col=0, parse_dates=True)\n",
"df['month'] = df.index.month"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"ExecuteTime": {
"end_time": "2018-07-30T21:08:12.355614Z",
"start_time": "2018-07-30T21:08:12.333570Z"
}
},
"outputs": [
{
"data": {
"text/html": [
"\n",
"\n",
"
\n",
" \n",
" \n",
" | \n",
" Date_Number | \n",
" MJO Phase | \n",
" Average specimens | \n",
" month | \n",
"
\n",
" \n",
" | Date | \n",
" | \n",
" | \n",
" | \n",
" | \n",
"
\n",
" \n",
" \n",
" \n",
" | 2011-10-01 | \n",
" 40817 | \n",
" 6 | \n",
" 1.000000 | \n",
" 10 | \n",
"
\n",
" \n",
" | 2011-10-11 | \n",
" 40827 | \n",
" 8 | \n",
" 18.625000 | \n",
" 10 | \n",
"
\n",
" \n",
" | 2011-10-12 | \n",
" 40828 | \n",
" 8 | \n",
" 17.375000 | \n",
" 10 | \n",
"
\n",
" \n",
" | 2011-10-13 | \n",
" 40829 | \n",
" 8 | \n",
" 37.181818 | \n",
" 10 | \n",
"
\n",
" \n",
" | 2011-10-15 | \n",
" 40831 | \n",
" 1 | \n",
" 39.238095 | \n",
" 10 | \n",
"
\n",
" \n",
"
\n",
"
\n",
"\n",
"
\n",
" \n",
" \n",
" | \n",
" MJO Phase | \n",
" Average specimens | \n",
"
\n",
" \n",
" | Date | \n",
" | \n",
" | \n",
"
\n",
" \n",
" \n",
" \n",
" | 2011-10-02 | \n",
" 6.0 | \n",
" 1.000000 | \n",
"
\n",
" \n",
" | 2011-10-09 | \n",
" NaN | \n",
" NaN | \n",
"
\n",
" \n",
" | 2011-10-16 | \n",
" 8.0 | \n",
" 31.683983 | \n",
"
\n",
" \n",
" | 2011-10-23 | \n",
" 1.0 | \n",
" 2.000000 | \n",
"
\n",
" \n",
" | 2011-10-30 | \n",
" 2.0 | \n",
" 7.263889 | \n",
"
\n",
" \n",
"
\n",
"
\n",
"\n",
"
\n",
" \n",
" \n",
" | \n",
" MJO Phase | \n",
" Average specimens | \n",
"
\n",
" \n",
" | Date | \n",
" | \n",
" | \n",
"
\n",
" \n",
" \n",
" \n",
" | 2011-10-02 | \n",
" NaN | \n",
" 1.000000 | \n",
"
\n",
" \n",
" | 2011-10-09 | \n",
" 6.0 | \n",
" NaN | \n",
"
\n",
" \n",
" | 2011-10-16 | \n",
" NaN | \n",
" 31.683983 | \n",
"
\n",
" \n",
" | 2011-10-23 | \n",
" 8.0 | \n",
" 2.000000 | \n",
"
\n",
" \n",
" | 2011-10-30 | \n",
" 1.0 | \n",
" 7.263889 | \n",
"
\n",
" \n",
"
\n",
"