4-plot-histogram.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "c300bafe-1ddc-4aa6-a36a-246e43f8fe5f",
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "import xlrd\n",
    "from scipy.stats import norm\n",
    "import matplotlib.mlab as mlab\n",
    "from math import sqrt\n",
    "from scipy import stats\n",
    "from numpy import mean\n",
    "from scipy.stats import sem\n",
    "from scipy.stats import t\n",
    "from scipy.stats import f\n",
    "from scipy.stats import bartlett\n",
    "from scipy.stats import levene\n",
    "from t_test_equSampleSize import *\n",
    "import csv\n",
    "\n",
    "def readcsv(filename):\n",
    "    with open(filename, 'r') as dest_f:\n",
    "        data_iter = csv.reader(dest_f, delimiter=',', quotechar='\"')\n",
    "        data = [data for data in data_iter]\n",
    "    data_array = np.asarray(data)\n",
    "    return data_array\n",
    "\n",
    "fig, ax = plt.subplots(figsize=(20, 16))\n",
    "data = readcsv(file1)\n",
    "output = file1[0:-29] + '.png'\n",
    "FigLabel = 'Antarctica'\n",
    "Title = ''\n",
    "font_size = 70\n",
    "(a, b) = np.shape(data)\n",
    "data = data[1:a,1]\n",
    "data[data == ''] = np.nan\n",
    "data = data.astype(np.float)\n",
    "data = data[np.logical_not(np.isnan(data))]\n",
    "x = data\n",
    "N, bins, patches = plt.hist(x, bins=500, alpha=0.5, color = \"gray\",normed='True')\n",
    "for i in range(0,132):\n",
    "    patches[i].set_facecolor('royalblue')\n",
    "\n",
    "plt.setp(ax, xticks=[-111,-90,-60,-30],\n",
    "                         xticklabels=['-111','-90', '-60', '-30'])\n",
    "plt.xlim(-111, -30)\n",
    "plt.ylim(0, 0.07)\n",
    "yticks = ax.yaxis.get_major_ticks()\n",
    "yticks[0].label1.set_visible(False)\n",
    "plt.yticks(np.arange(0, 0.07, step=0.03))\n",
    "plt.gca().set(title='Frequency Histogram', ylabel='Frequency')\n",
    "box = ax.get_position()\n",
    "ax.set_position([box.x0+0.01, box.y0 + box.height * 0.2,\n",
    "                 box.width, box.height * 0.9])\n",
    "legend1 = ax.legend(loc='upper left', bbox_to_anchor=(0, 1),\n",
    "                    fancybox=True, shadow=False, ncol=5, title=FigLabel, edgecolor='k', fontsize=font_size,frameon=False,markerscale=0)\n",
    "ax.add_artist(legend1)\n",
    "legend1.get_title().set_fontsize(font_size)\n",
    "plt.xticks(list(plt.xticks()[0]) + [-111])\n",
    "plt.yticks(fontsize=font_size,rotation=90)\n",
    "plt.xticks(fontsize=font_size)\n",
    "ax.set_xlabel('Daily minimum temperature (°C)', fontsize=font_size + 2,labelpad=20, fontweight='bold')\n",
    "ax.set_ylabel('Frequency (%)', fontsize=font_size + 2,labelpad=20, fontweight='bold')\n",
    "plt.title(Title, fontsize=font_size + 2)\n",
    "ax.spines['top'].set_linewidth(0)\n",
    "ax.spines['bottom'].set_linewidth(3)\n",
    "ax.spines['left'].set_linewidth(3)\n",
    "ax.tick_params(direction='out', length=16, width=3)\n",
    "ax.spines['right'].set_linewidth(0)\n",
    "plt.savefig(output, dpi=100, transparent=True)\n",
    "plt.show()\n",
    "\n"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "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.10.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}