Skip to content

Instantly share code, notes, and snippets.

@cdw
Last active March 18, 2020 21:08
Show Gist options
  • Save cdw/242f509d3db64f9ba982370b007db8f5 to your computer and use it in GitHub Desktop.
Save cdw/242f509d3db64f9ba982370b007db8f5 to your computer and use it in GitHub Desktop.
Bar plot with error
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"import matplotlib.pyplot as plt\n",
"import numpy as np\n",
"import pandas as pd"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<div>\n",
"<style scoped>\n",
" .dataframe tbody tr th:only-of-type {\n",
" vertical-align: middle;\n",
" }\n",
"\n",
" .dataframe tbody tr th {\n",
" vertical-align: top;\n",
" }\n",
"\n",
" .dataframe thead th {\n",
" text-align: right;\n",
" }\n",
"</style>\n",
"<table border=\"1\" class=\"dataframe\">\n",
" <thead>\n",
" <tr style=\"text-align: right;\">\n",
" <th></th>\n",
" <th>cell_age</th>\n",
" <th>cell_mean</th>\n",
" <th>ci_lo</th>\n",
" <th>ci_hi</th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>0</th>\n",
" <td>18</td>\n",
" <td>0.427711</td>\n",
" <td>0.395802</td>\n",
" <td>0.458543</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1</th>\n",
" <td>32</td>\n",
" <td>0.448410</td>\n",
" <td>0.412520</td>\n",
" <td>0.484413</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"</div>"
],
"text/plain": [
" cell_age cell_mean ci_lo ci_hi\n",
"0 18 0.427711 0.395802 0.458543\n",
"1 32 0.448410 0.412520 0.484413"
]
},
"execution_count": 2,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"dat = dict(\n",
" cell_age= (18,32),\n",
" cell_mean = (0.4277114502555463,0.44841020643768326),\n",
" ci_lo = (0.3958021553533243, 0.412519854795916), \n",
" ci_hi = (0.4585425031102684, 0.4844134134481588)\n",
")\n",
"df = pd.DataFrame(dat)\n",
"df"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [],
"source": [
"err = np.array(((df.cell_mean - df.ci_lo).values,\n",
" (df.cell_mean - df.ci_hi).values)).T"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAIgAAAGoCAYAAACdYa3QAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAUDElEQVR4nO3df5AfdX3H8eeLSAQMwpScCEk0QVOdaBXkiOg4IShq0DEpBcfgj4KgUccIirZGW0FiLahVi05UomBTBSKIliPGBkWjtBXMJfzQRCLXFMiZtgR/oFAjBN79Y/fkm8s379tcbu82+b4eM9+53c/u7fd9k1d2P/tbEYHZruw31gVYszkglnJALOWAWMoBsZQDYqlaAyJpjqSNkvokLWoz/UxJWyXdVn7eUmc9tvueUNeCJY0DlgAvB/qBNZJ6ImLDoFm/FhELqy534sSJMXXq1JEr1ABYu3bt/RHRNbi9toAAM4G+iNgEIGk5MA8YHJDdMnXqVHp7e0egPGsl6Z527XVuYiYBm1vG+8u2wU6VdIekr0uaUmM9Ngx1BkRt2gYf178emBoRzwO+CyxruyBpgaReSb1bt24d4TItU2dA+oHWNcJkYEvrDBHxy4j4Qzn6ReDYdguKiKUR0R0R3V1dO20mrUZ1BmQNMF3SNEnjgflAT+sMko5oGZ0L/KzGemwYauukRsR2SQuBVcA44PKIWC9pMdAbET3AOZLmAtuBXwFn1lWPDY/2ttP93d3d4b2YkSdpbUR0D273kVRLOSCWckAs5YBYygGxlANiKQfEUg6IpRwQSzkglnJA9iKzZ89m9uzZo/qdDoilHBBLOSCWckAs5YBYygGxlANiKQfEUg6IpRwQSzkglnJAGJtzHHsLB8RSDoilHBBLOSCWckAs5YBYygGxlANiKQfEUg6Ipep8Tqq1+sEIPBXpN78buWWdsNPDhNryGsRSDoilHBBLOSCWckAs5YBYygGxlANiKQfEUg6IpRwQSzkglnJALOWAWMoBsdRefz3Ixbfev8fLuPfBR0ZsWYuOmbjHy2gSr0Es5YBYygGxlANiKQfEUg6IpRwQSzkglnJALOWAWGqvP9TeSVZfcumof6fXIJZyQCzlgFjKAbGUA2KpWgMiaY6kjZL6JC1K5jtNUkiq9lQTGzW1BUTSOGAJcDIwAzhd0ow28x0MnAPcUlctNnx1rkFmAn0RsSkiHgaWA/PazPcR4OPAthprsWGqMyCTgM0t4/1l2x9JOgaYEhEraqzD9kCdAVGbtvjjRGk/4NPAe4dckLRAUq+k3q1bt45giTaUOgPSD0xpGZ8MbGkZPxh4LrBa0t3A8UBPu45qRCyNiO6I6O7q6qqxZBuszoCsAaZLmiZpPDAf6BmYGBEPRMTEiJgaEVOBm4G5ETECz3i0kVLbybqI2C5pIbAKGAdcHhHrJS0GeiOiJ1/C6FnwxevGuoTGqvVsbkSsBFYOajt/F/POrrMWGx4fSbWUA2IpB8RSDoilHBBLOSCWckAs5YBYygGxlANiKQfEUg6IpRwQSzkglnJALOWAWMoBsZQDYikHxFIOiKUcEEs5IJZyQCzlgFjKAbGUA2IpB8RSDoilHBBLOSCWckAs5YBYygGxlANiKQfEUg6IpRwQSzkglnJALOWAWMoBsZQDYikHxFIOiKUcEEs5IJZyQCzlgFjKAbGUA2IpB8RSDoilHBBLDfnWS0ldwFuBqa3zR8RZ9ZVlTVHltajXATcB3wUerbcca5oqATkoIt5feyXWSFX6ICskvar2SqyRqgTkXIqQbJP0u/Lz27oLs2YYchMTEQePRiHWTFX6IEiaC8wqR1dHxIr6SrImGXITI+liis3MhvJzbtlmHaDKGuRVwNER8RiApGXArcCiOguzZqh6JPXQluFD6ijEmqnKGuQi4FZJ3wdE0Rf5QK1VWWNU2Yu5StJq4DiKgLw/Iv6n7sKsGXa5iZH07PLnC4AjgH5gM3Bk2WYdIFuDnAcsAD7ZZloALx1q4ZLmAJcA44AvRcTFg6a/HXgnxTmeB4EFEbGhWuk2GnYZkIhYUA6eHBHbWqdJOmCoBUsaBywBXk6x9lkjqWdQAK6MiC+U888FPgXM2b0/wepUZS/mPyq2DTYT6IuITRHxMLAcmNc6Q0S0HrJ/EsWayRpkl2sQSU8FJgEHSjqGooMK8GTgoArLnkTRZxnQD7ywzfe8k2JzNp4Kmy0bXVkf5JXAmcBkilX/gN8BH6ywbLVp22kNERFLgCWSXg/8LXDGTguSFlD0h3ja055W4attpGR9kGXAMkmnRsS1w1h2PzClZXwysCWZfznw+V3UshRYCtDd3e3N0CiqchzkWkmvBp4DHNDSvniIX10DTJc0DfgFMB94fesMkqZHxF3l6KuBu7BGqXJN6hco+hwnAl8CTgN+PNTvRcR2SQuBVRS7uZdHxHpJi4HeiOgBFko6CXgE+DVtNi82tqocan9xRDxP0h0RcaGkTwLfqLLwiFgJrBzUdn7L8Lm7Va2Nuiq7ub8vf/6fpCMp/rdPq68ka5Iqa5AVkg4FPgGso9gT+VKtVVljVOmkfqQcvFbSCuCAiHig3rKsKbIDZX+RTCMiKvVDbO+WrUFek0wLKnZUbe+WHSh782gWYs1U5aLlwyVdJunb5fgMSWfXX5o1QZXd3H+iONh1ZDn+c+DddRVkzVIlIBMj4mrgMSiOkOKbuDtGlYA8JOkwyjOxko4HvJvbIaocKDsP6AGeIenfgS6K8zHWAdKASNqP4gzuCcCzKK7x2BgRj4xCbdYAaUAi4jFJn4yIFwHrR6kma5AqfZAbJJ0qqd0VYraPq9oHeRKwXdI2is1MRMSTa63MGmGoPoiA50TEvaNUjzVMuomJiAC+OUq1WANV6YPcLOm42iuxRqrSBzkReJuke4CHeLwP8rxaK7NGqBKQk2uvwhpryE1MRNxD8QCZ15SfQ8s26wBVTvefC1wBPKX8fFXSu+ouzJqhyibmbOCFEfEQgKSPAT8CPltnYdYMVfZixI6n9x+l/X23tg+qsgb5MnCLpIHjIX8OXFZfSdYkVW57+FT5jLKXUKw53hwRt9ZdmDVDdtvDcRRXk307ItZR3DSFpLmS9ouItaNVpI2drA/yCeBnbdo3lNOsA2QBOSwi7h7cGBF9wGG1VWSNkgXkwGTak0a6EGumLCDflfTRwRcKSboQ+F69ZVlTZHsx76W4i79P0m1l2/OBXuAtdRdmzZDdevkQcLqkoygePwWwPiI2jUpl1ghVjoNsAhyKDuUXK1vKAbFU1XfWjQMOb53fFzJ3hiqPwXwXcAHwv5Q3cFPcp+tLDjtAlTXIucCzIuKXdRdjzVOlD7IZ383fsaqsQTYBqyV9C/jDQGNEfGrXv2L7iioBubf8jC8/1kGqHCi7cDQKsWaqshfTBfw1O7/twS//6QBVOqlXAHdSPJ/9QuBuild9WAeoEpDDIuIy4JGI+EFEnAUcX3Nd1hBVOqkDj5v67/LFQlso3h5lHaBKQP5O0iEU14d8luKlhu+ptSprjCp7MSvKwQco7vS3DlJlL2Ya8C5gKjuerJtbX1nWFFU2Mf9CcSfd9Tx+ss46RJWAbIuIz9ReiTVSlYBcIukC4AZ2PBezrraqrDGqBOTPgDdRvDa99XoQH0ntAFUCcgpwVEQ8XHcx1jxVjqTeTvEIKutAVdYghwN3SlrDjn0Q7+Z2gCoBuaD2KqyxhnoU9zjgQxFx0ijVYw0z1KO4H6V4Jfsho1SPNUylA2XATyR9h+JJywBExDm1VWWNUSUg3yo/1oGqnM1dNhqFWDNVOZs7HbgImMGO16QeVWNd1hBVDpR9Gfg8sJ3iepB/Br5SZ1HWHFUCcmBE3AgoIu6JiA/j8zAdo0pAtpWvR71L0kJJp1A81H9IkuZI2iipT9KiNtPPk7RB0h2SbpT09N2s32pWJSDvBg4CzgGOBd4InDHUL5UH2ZZQvG9mBsXjrGYMmu1WoLt8OdHXgY9XL91GQ5W9mDUAkiIi3rwby54J9A0800zScmAexYN4B5b9/Zb5b6YInzVIlffFvEjSBsqnLkt6vqTPVVj2JIonAwzoL9t25Wzg2xWWa6OoyibmH4FXAr8EiIjbgVkVfq/dK0Oi7YzSG4FudvGIb0kLJPVK6t26dWuFr7aRUukZZRGxeVDTo21n3FE/MKVlfDLFTVc7kHQS8DfA3Ij4w+Dp5fcvjYjuiOju6uqqUrKNkEoPkJH0YiAkjZf0Pto/5H+wNcB0SdMkjQfmAz2tM0g6BriUIhz37WbtNgqqBOTtwDsp+g+/AI4ux1MRsR1YCKyiCNTVEbFe0mJJAxcbfQKYAFwj6TZJPbtYnI2RKnsx9wNvGM7CI2IlsHJQ2/ktw77OpOGq7MUcJel6SVsl3SfpuvLx3NYBqmxirgSuBo4AjgSuAa6qsyhrjkpvvYyIr0TE9vLzVXaxu2r7nioXDH2/PI+ynCIYrwO+JelPACLiVzXWZ2OsSkBeV/5826D2sygC4/7IPqzKXsy00SjEmmmXfRBJx0l6asv4X5Z7MJ8Z2LzYvi/rpF4KPAwgaRZwMcXVZA8AS+svzZog28SMa+mAvg5YGhHXAte2vMPO9nHZGmScpIEAvYwd33RZ6T0ztvfL/qGvAn4g6X7g98BNAJKeid/+0DGyt15+VNKNFEdQb4iIgYNj+1E81M46QLqpiIib27T9vL5yrGn8UkNLOSCWckAs5YBYygGxlANiKQfEUg6IpRwQSzkglnJALOWAWMoBsZQDYikHxFIOiKUcEEs5IJZyQCzlgFjKAbGUA2IpB8RSDoilHBBLOSCWckAs5YBYygGxlANiKQfEUg6IpRwQSzkglnJALOWAWMoBsZQDYikHxFIOiKUcEEs5IJZyQCzlgFjKAbGUA2IpB8RSDoilHBBLOSCWckAs5YBYygGxVK0BkTRH0kZJfZIWtZk+S9I6SdslnVZnLTY8tQVE0jhgCXAyMAM4XdKMQbPdC5wJXFlXHbZn6nzF+kygLyI2AUhaDswDNgzMEBF3l9Meq7EO2wN1bmImAZtbxvvLNtuL1BkQtWmLNm1DL0haIKlXUu/WrVv3sCzbHXUGpB+Y0jI+GdgynAVFxNKI6I6I7q6urhEpzqqpMyBrgOmSpkkaD8wHemr8PqtBbQGJiO3AQmAV8DPg6ohYL2mxpLkAko6T1A+8FrhU0vq66rHhqXMvhohYCawc1HZ+y/Aaik2PNZSPpFrKAbGUA2IpB8RSDoilHBBLOSCWckAs5YBYygGxlANiKQfEUg6IpRwQSzkglnJALOWAWMoBsZQDYikHxFIOiKUcEEs5IJZyQCzlgFjKAbGUA2IpB8RSDoilHBBLOSCWckAs5YBYygGxlANiKQfEUg6IpRwQSzkglnJALOWAWMoBsZQDYikHxFIOiKUcEEs5IJZyQCzlgFjKAbGUA2IpB8RSDoilHBBLOSCWckAs5YBYygGxlANiKQfEUg6IpRwQSzkglnJALOWAWMoBsVStAZE0R9JGSX2SFrWZ/kRJXyun3yJpap312O6rLSCSxgFLgJOBGcDpkmYMmu1s4NcR8Uzg08DH6qrHhqfONchMoC8iNkXEw8ByYN6geeYBy8rhrwMvk6Qaa7LdVGdAJgGbW8b7y7a280TEduAB4LAaa7Ld9IQal91uTRDDmAdJC4AF5eiDkjbuYW3tTATu39OFfGAEChnCiNTZxtPbNdYZkH5gSsv4ZGDLLubpl/QE4BDgV4MXFBFLgaU11QmApN6I6K7zO0bCaNdZ5yZmDTBd0jRJ44H5QM+geXqAM8rh04DvRcROaxAbO7WtQSJiu6SFwCpgHHB5RKyXtBjojYge4DLgK5L6KNYc8+uqx4ZH/g9bkLSg3JQ12mjX6YBYyofaLdWRAZF0uaT7JP20pe1oSTdLuk1Sr6SZY1zjAZJ+LOl2SeslXVi2X1Gevvhp+XfsX2shEdFxH2AW8ALgpy1tNwAnl8OvAlaPcY0CJpTD+wO3AMeXtan8XAW8o846OnINEhE/ZOfjLQE8uRw+hJ2P2YyqKDxYju5ffiIiVpbTAvgxxfGl2tR5oGxv825glaR/oNj0vniM6xk44bkWeCawJCJuaZm2P/Am4Nw6a+jINcguvAN4T0RMAd5DcYxmTEXEoxFxNMVaYqak57ZM/hzww4i4qc4aHJDHnQF8oxy+huJsdCNExG+A1cAcAEkXAF3AeXV/twPyuC3ACeXwS4G7xrAWJHVJOrQcPhA4CbhT0luAVwKnR8RjddfRkX0QSVcBs4GJkvqBC4C3ApeUJw238fjZ47FyBLCs7IfsB1wdESskbQfuAX5UXjrzjYhYXFcRPpJqKW9iLOWAWMoBsZQDYikHxFIOyCCSnippuaT/lLRB0kpJfzrE7zxY/pzaeoZ4X+CAtCjvyfkmxZncZ0TEDOCDwOFjW9nYcUB2dCLwSER8YaAhIm4bON8h6a8krZF0x8D1GVVImiDpRknrJP1E0ryWaR+SdKek70i6StL7yvZnSPpXSWsl3STp2SP4d1bWkUdSE8+lOHu6E0mvAKZTnKMR0CNpVnnpwFC2AadExG8lTQRultQDHAucChxD8W+xruX7lwJvj4i7JL2Q4uTcS4f/pw2PA1LdK8rPreX4BIrAVAmIgL+XNAt4jOKOwsOBlwDXRcTvASRdX/6cQHG5wTUtd6I+cWT+jN3jgOxoPcX9Oe0IuCgiLh3Gct9Acfb12Ih4RNLdwAG0v7MQik3/b8pT/WPKfZAdfQ94oqS3DjRIOk7SCRT395xV/u9G0iRJT6m43EOA+8pwnMjjtzn+G/Ca8vrTCcCrASLit8B/SXpt+V2S9PyR+AN3lwPSoryM7xTg5eVu7nrgw8CWiLgBuJLiLOpPKJ5GcHDFRV8BdEvqpVib3Fl+3xqKuwtvp7gWpZfiBnbK+c6WdDvFmm3wkxFGhc/mjjFJEyLiQUkHUfRnFkTEurGua4D7IGNvaflgnQOAZU0KB3gNYkNwH8RSDoilHBBLOSCWckAs5YBY6v8BGhuFiQwLGUYAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 144x432 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"fig, ax = plt.subplots(1, 1, figsize=(2, 6))\n",
"ax.bar(df.cell_age, \n",
" df.cell_mean, \n",
" width = 12, \n",
" yerr=err,\n",
" color=(\"skyblue\", \"pink\"))\n",
"ax.set(xticks=df.cell_age,\n",
" xlabel=\"Cell age\", \n",
" ylabel=\"Spearman Correlation\", \n",
" )\n",
"plt.tight_layout()"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"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.7.6"
}
},
"nbformat": 4,
"nbformat_minor": 4
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment