2023-01-12 16:13:39 +02:00
|
|
|
{
|
|
|
|
|
"cells": [
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2023-01-12 21:45:54 +02:00
|
|
|
"execution_count": 1,
|
2023-01-12 20:25:12 +02:00
|
|
|
"id": "41405f3d",
|
2023-01-12 16:13:39 +02:00
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [],
|
|
|
|
|
"source": [
|
|
|
|
|
"from matplotlib import pyplot as plt\n",
|
|
|
|
|
"import numpy as np"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "markdown",
|
2023-01-12 20:25:12 +02:00
|
|
|
"id": "dfc88764",
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"source": [
|
|
|
|
|
"## General conversions \n",
|
|
|
|
|
"\n",
|
|
|
|
|
"Conversions\n",
|
|
|
|
|
"\n",
|
|
|
|
|
"Sunflower Seeds, sunflower seed kernels, dried\n",
|
|
|
|
|
"584kcal/100g\n",
|
|
|
|
|
"https://fdc.nal.usda.gov/fdc-app.html#/food-details/170562/nutrients\n",
|
|
|
|
|
"FDC ID: 170562\n",
|
|
|
|
|
"USDA Unit LB/acre\n",
|
|
|
|
|
" \n",
|
|
|
|
|
"Potatoes, flesh and skin, raw\n",
|
|
|
|
|
"https://fdc.nal.usda.gov/fdc-app.html#/food-details/170026/nutrients\n",
|
|
|
|
|
"77kcal/100g\n",
|
|
|
|
|
"FDC ID: 170026 \n",
|
|
|
|
|
"USDA Unit CWT/acre\n",
|
|
|
|
|
"1 CWT = 100lbs\n",
|
|
|
|
|
" \n",
|
|
|
|
|
"Wheat, hard red winter\n",
|
|
|
|
|
"FDC ID: 168890\n",
|
|
|
|
|
"327kcal/100g\n",
|
|
|
|
|
"https://fdc.nal.usda.gov/fdc-app.html#/food-details/168890/nutrients\n",
|
|
|
|
|
"USDA Unit BU/acre\n",
|
|
|
|
|
"60lb / BU https://www.plainsgrains.org/wp-content/uploads/2018/10/2017-HRWW-Report__HiRes-FINAL.pdf\n",
|
|
|
|
|
"\n",
|
|
|
|
|
"Corn grain, yellow\n",
|
|
|
|
|
"365kcal/100g\n",
|
|
|
|
|
"https://fdc.nal.usda.gov/fdc-app.html#/food-details/170288/nutrients\n",
|
|
|
|
|
"FDC ID: 170288 \n",
|
|
|
|
|
"USDA Unit BU/acre\n",
|
|
|
|
|
"Shelled, 56lbs/BU\n",
|
|
|
|
|
"\n",
|
|
|
|
|
"Soybeans, mature seeds, raw\n",
|
|
|
|
|
"446kcal/100g\n",
|
|
|
|
|
"https://fdc.nal.usda.gov/fdc-app.html#/food-details/174270/nutrients\n",
|
|
|
|
|
"FDC ID: 174270 \n",
|
|
|
|
|
"USDA Unit BU/acre\n",
|
|
|
|
|
"60LB/BU https://soybeanresearchinfo.com/research-highlight/exploring-the-impact-of-genetics-on-soybean-test-weight\n"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "markdown",
|
|
|
|
|
"id": "72029669",
|
2023-01-12 16:13:39 +02:00
|
|
|
"metadata": {},
|
|
|
|
|
"source": [
|
|
|
|
|
"## Corn"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2023-01-12 21:45:54 +02:00
|
|
|
"execution_count": 2,
|
2023-01-12 20:25:12 +02:00
|
|
|
"id": "116f0267",
|
2023-01-12 16:13:39 +02:00
|
|
|
"metadata": {},
|
2023-01-12 21:45:54 +02:00
|
|
|
"outputs": [
|
|
|
|
|
{
|
|
|
|
|
"data": {
|
|
|
|
|
"text/plain": [
|
|
|
|
|
"array([[ nan, 2022. , nan, nan, nan, nan, nan, nan,\n",
|
|
|
|
|
" nan, nan, nan, nan, nan, 0. , nan, nan,\n",
|
|
|
|
|
" nan, nan, 172.3, nan],\n",
|
|
|
|
|
" [ nan, 2021. , nan, nan, nan, nan, nan, nan,\n",
|
|
|
|
|
" nan, nan, nan, nan, nan, 0. , nan, nan,\n",
|
|
|
|
|
" nan, nan, 176.7, nan],\n",
|
|
|
|
|
" [ nan, 2020. , nan, nan, nan, nan, nan, nan,\n",
|
|
|
|
|
" nan, nan, nan, nan, nan, 0. , nan, nan,\n",
|
|
|
|
|
" nan, nan, 171.4, nan]])"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
"execution_count": 2,
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"output_type": "execute_result"
|
|
|
|
|
}
|
|
|
|
|
],
|
2023-01-12 16:13:39 +02:00
|
|
|
"source": [
|
|
|
|
|
"csv=np.genfromtxt ('./data/corn.csv',\n",
|
|
|
|
|
" skip_header=1,\n",
|
|
|
|
|
" delimiter=\",\")\n",
|
|
|
|
|
"csv[0:3]\n",
|
|
|
|
|
"#csv[0:3].astype(float)"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2023-01-12 21:45:54 +02:00
|
|
|
"execution_count": 3,
|
2023-01-12 20:25:12 +02:00
|
|
|
"id": "0db69a97",
|
2023-01-12 16:13:39 +02:00
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [],
|
|
|
|
|
"source": [
|
|
|
|
|
"corn_year=csv[:,1]\n",
|
2023-01-12 20:25:12 +02:00
|
|
|
"corn_bu_per_acre=csv[:,-2]\n",
|
|
|
|
|
"\n",
|
|
|
|
|
"grams_per_lbs=453.592\n",
|
|
|
|
|
"corn_lbs_per_bu=(56.0/1.0)\n",
|
|
|
|
|
"corn_kcal_per_gram=(365/100)\n",
|
|
|
|
|
"corn_kcal_per_acre = corn_bu_per_acre*corn_lbs_per_bu*grams_per_lbs*corn_kcal_per_gram"
|
2023-01-12 16:13:39 +02:00
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2023-01-13 20:54:15 +02:00
|
|
|
"execution_count": 32,
|
2023-01-12 20:25:12 +02:00
|
|
|
"id": "eb375958",
|
2023-01-12 16:13:39 +02:00
|
|
|
"metadata": {},
|
2023-01-12 21:45:54 +02:00
|
|
|
"outputs": [
|
|
|
|
|
{
|
|
|
|
|
"data": {
|
2023-01-13 20:54:15 +02:00
|
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAjsAAAHFCAYAAAAUpjivAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAA9hAAAPYQGoP6dpAACp40lEQVR4nOzdeXiMV/vA8e9kXySybxI7sYQQStHWHpSWt1otql7aaqPaoq1WdaFFi1/RovFqtVpbV1RL7Vvb2JIIRcQWEZFIZJFE9szz+yNmzCSTjey5P9flIs8888w5M8NzO+c+91EpiqIghBBCCFFHGVV3A4QQQgghKpMEO0IIIYSo0yTYEUIIIUSdJsGOEEIIIeo0CXaEEEIIUadJsCOEEEKIOk2CHSGEEELUaRLsCCGEEKJOk2BHCCGEEHWaBDuiRps9ezYqlYqbN28afNzHx4c+ffroHYuOjmby5Mm0bt0aS0tLHBwc6NChAy+++CLR0dFlfu2tW7eiUqlwdHQkOzv7frpRIk0fK8uVK1dQqVSsWbNGeywoKIjZs2eTkpJSaa9b2dasWYNKpeLKlSuV+tw+ffoU+Y5VpD59+qBSqRg8eHCRxzSf3f/93/8ZfG5ZvqO3b99mwYIF+Pr6Ymtri42NDS1atGDUqFEcPHhQ79zw8HDGjRtH8+bNsbCwwMnJCT8/P6ZMmUJqamqZ+vPRRx/Rrl071Gp1mfrwf//3f0U+i9zcXP73v//xwAMP4ODggJWVFU2aNGH48OFs3ry5yPuj+WVqaoqjoyMPPPAA06ZN48yZMyW2dfr06ahUKoYNG2bw8fPnz2NmZkZoaGiZ+i5qLgl2RJ1y7do1/Pz82L17N9OnT2f79u188803jB49muPHj3P58uUyX2v16tUAJCUlsWXLlkpqMbzwwgscPny40q5vSFBQEHPmzKnVwc7QoUM5fPgw7u7u1d2UCrFz50727dtXrueU9h3Nz8/H39+fefPm8eSTT/Lzzz/zyy+/MG3aNG7dusVff/2lPffEiRN06dKFs2fP8sEHH7Bjxw5WrlzJ0KFD2blzJ0lJSaW25/r16yxcuJCPPvoII6N7v72MGzeOV199lb59+7Ju3Tp+//133nvvPUxMTNi5c2eR81999VUOHz7MwYMHWbt2LSNGjGDr1q34+vqyaNEig6+Rm5vLunXrANixYwcxMTFFzmndujVjx45l2rRp99wXUUMoQtRgH374oQIoCQkJBh9v37690rt3b+3PH3zwgQIoly9fNnh+fn5+mV43NjZWMTExUfr166dYWFgoAwcOLHfbS3P79u0Kv6YhkZGRCqB8++232mOLFi1SACUyMrJCXysjI0NRq9UVes3K8O2335a5/71799b7jlW03r17K61bt1aaN2+udOnSRe/903x2ixYtKvK8snxH9+3bpwDKN998Y/C1df8+PPfcc4q1tbWSmppq8NyyfK4zZsxQGjVqpHfdkvqgKEW/i5cvX1YA5YMPPii1zSVdOyMjQxk8eLACKNu3by/y+M8//6wAytChQxVAmTdvnsHXCw4OVgDln3/+KbbfouaTkR1RpyQmJmJkZISLi4vBx8v6v83vvvuOvLw8pk2bxhNPPMHevXuJiooqcl5KSgrPP/88Dg4ONGjQgKFDh3L58mVUKhWzZ8/WnqeZqgoNDeXJJ5/E3t6eFi1a6D1W2IYNG+jRowcNGjSgQYMGdOrUSfs/eYCmTZvy3//+t8jzSpt2mT17Nm+99RYAzZo1004BHDhwAKBI24t7Pc1U0K5du5g4cSLOzs5YWVlpp1N+/PFHevTogbW1NQ0aNGDQoEGcOHGi2HZBwbSEiYkJn3zySZHHDh06hEql4ueff9Z7/cJTUXv27KF///7Y2tpiZWVFr1692Lt3b4mvC6AoCgsXLqRJkyZYWFjg5+fHn3/+WeQ8tVrN3Llz8fb2xtLSEjs7Ozp27Mjnn39e6msUx9TUlHnz5hESEsKPP/5YpueU5TuamJgIUOzol+7fh8TERGxtbWnQoIHBc0ubas3JyWH16tWMGTPmvkZ1ytPmklhaWrJ69WpMTU0Nju6sXr0aMzMzvv32W7y8vPj2229RDOyL3aVLF9q2bcvKlSvL0QtR00iwI+qUHj16oFareeKJJ9i5c2eZ8wwK++abb3B3d2fIkCFMnDgRtVqtl/MCBTe9xx57jA0bNvD222+zefNmunfvbjD3QuOJJ56gZcuW/PzzzyX+4/nBBx8wduxYPDw8WLNmDZs3b2b8+PEGA67yeuGFF3j11VcB2LRpE4cPH+bw4cP4+fnd0/UmTpyIqakpa9eu5ZdffsHU1JT58+czevRo2rVrx08//cTatWtJS0vj4Ycf5uzZs8Veq2nTpjz++OOsXLmS/Px8vceWL1+Oh4cH//nPf4p9/rp16/D398fW1pbvvvuOn376CQcHBwYNGlRqwDNnzhzefvttBg4cyJYtWwgICODFF18kIiJC77yFCxcye/ZsRo8ezbZt2/jxxx95/vnn73tK8Omnn6ZLly6899575Obmlnp+Wb6jXbt2xdTUlNdff53169cTGxtb7PV69OhBbGwsY8eO5eDBg2RmZpar/UePHiUxMZG+ffuW63mFtW3bFjs7O+bMmcOqVavuKSdLw8PDgy5duhAUFEReXp72+LVr19i1axfDhw/H2dmZ8ePHc/HiRQ4dOmTwOn369OHPP/80GAyJWqKaR5aEKFF5p7HUarXy0ksvKUZGRgqgqFQqpW3btsq0adPKPGVz6NAhBVDeeecd7TWbNWumNGnSRG8of9u2bQqgBAYG6j3/k08+UQDlww8/LNIPQ0Pzmsc0Ll++rBgbGytjx44tsZ1NmjRRxo8fX+R44WmX8k5jFW57ca+nmQp67rnn9M67evWqYmJiorz66qt6x9PS0hQ3Nzdl1KhRJfZr//79CqBs3rxZeywmJkYxMTFR5syZU+T1NX24ffu24uDgoDz22GN618vPz1d8fX2Vbt26Ffvc5ORkxcLCQvnPf/6j99x//vlHAfTez2HDhimdOnUqsQ/l0bt3b6V9+/aKoijKnj17FEBZtmyZoijFT9OU9TuqKIqyevVqpUGDBgqgAIq7u7vy3HPPKYcOHdI7LysrSxkxYoT2PGNjY6Vz587KrFmzlPj4+FL7sWDBAgVQ4uLi9I6XdxpLUQr+bjk5OWnb4ujoqDz11FPK1q1by3VtRVGUp59+WgGUGzduaI999NFHCqDs2LFDUZSCv3MqlUoZN26cwWt89dVXCqCEh4eX+B6ImktGdkSdolKpWLlyJZcvX+bLL79kwoQJ5ObmsmTJEtq3b19k9YkhmqmiiRMnaq/53//+l6ioKL3RAc21Ro0apff80aNHF3vtkSNHlvr6u3fvJj8/n1deeaXUc2uCwn3auXMneXl5PPfcc+Tl5Wl/WVhY0Lt3b+10WXH69OmDr68vK1as0B5buXIlKpWKSZMmFfu8oKAgkpKSGD9+vN7rqtVqBg8ezPHjx7l9+7bB5x4+fJisrCzGjh2rd7xnz540adJE71i3bt04efIkkydPvq/RQ0P69++Pv78/H330EWlpacWeV9bvqOaca9eusWHDBl577TW8vLxYt24dvXv31pveMTc3Z/PmzZw9e5YlS5bwzDPPkJCQwLx582jbtm2REa7Crl+/jkqlwsnJ6V67r/Xoo49y9epVNm/ezJtvvkn79u3ZsmULjz/+OFOmTCnXtZRCozGKominrgYOHAgUTOf26dOHX3/91eDnqZkWN5TELGoHCXZEjWZiYgJQZEpDIy8vD1NT0yLHmzRpQkBAAKtXr+bChQv8+OOPZGVlaXNVipOWlsbPP/9Mt27dcHZ2JiUlhZSUFP7zn/+gUqn0cmYSExMxMTHBwcFB7xqurq7FXr8sK4cSEhIA8PT0LPXcmqBwn27cuAHAAw88gKmpqd6vH3/8sdgyArpee+019u7dS0REBLm5uXz11Vc8+eSTuLm5Ffsczes++eSTRV53wYIFKIpS7IoiTZ6IoesXPjZz5kz+7//+jyNHjjBkyBAcHR3p378/wcHBpfarLBYsWMDNmzeLXap
|
2023-01-12 21:45:54 +02:00
|
|
|
"text/plain": [
|
|
|
|
|
"<Figure size 640x480 with 1 Axes>"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"output_type": "display_data"
|
|
|
|
|
}
|
|
|
|
|
],
|
2023-01-12 20:25:12 +02:00
|
|
|
"source": [
|
2023-01-13 20:54:15 +02:00
|
|
|
"plt.plot(corn_year,corn_bu_per_acre,\"o\",markerfacecolor=\"none\",label=\"corn\")\n",
|
2023-01-12 20:25:12 +02:00
|
|
|
"plt.xlabel(\"time (calendar year)\")\n",
|
|
|
|
|
"plt.ylabel(\"Corn, bu/acre\")\n",
|
|
|
|
|
"plt.title(\"US Agriculture yields, NASS (USDA)\")\n",
|
|
|
|
|
"plt.legend()\n",
|
2023-01-13 20:54:15 +02:00
|
|
|
"plt.grid()\n",
|
2023-01-12 20:25:12 +02:00
|
|
|
"\n",
|
|
|
|
|
"#plt.show()\n",
|
|
|
|
|
"#plt.savefig(\"bear-quadratic.pdf\")"
|
|
|
|
|
]
|
2023-01-12 16:13:39 +02:00
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2023-01-12 21:45:54 +02:00
|
|
|
"execution_count": 5,
|
2023-01-12 20:25:12 +02:00
|
|
|
"id": "9c6ffb7f",
|
2023-01-12 16:13:39 +02:00
|
|
|
"metadata": {},
|
2023-01-12 21:45:54 +02:00
|
|
|
"outputs": [
|
|
|
|
|
{
|
|
|
|
|
"data": {
|
|
|
|
|
"text/plain": [
|
2023-01-13 20:54:15 +02:00
|
|
|
"<matplotlib.legend.Legend at 0x255cc993610>"
|
2023-01-12 21:45:54 +02:00
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
"execution_count": 5,
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"output_type": "execute_result"
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"data": {
|
2023-01-13 20:54:15 +02:00
|
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAjcAAAHFCAYAAAAOmtghAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAA9hAAAPYQGoP6dpAACUY0lEQVR4nOzdd1hU19YH4N/Q29C7NBtNFEXUYIm9a2LUaJRrS1GxJOpnitHEcjVGTdQYI8bYEiXGEtuNBqMINrDQYgMERUEcQHpRabO/P/AcZ5ihCkxhvc/Dc69nzpzZe5jMWey99toCxhgDIYQQQoia0FB0AwghhBBCGhMFN4QQQghRKxTcEEIIIUStUHBDCCGEELVCwQ0hhBBC1AoFN4QQQghRKxTcEEIIIUStUHBDCCGEELVCwQ0hhBBC1AoFN0RprFixAgKBAFlZWXIf9/LyQr9+/aSOpaamYs6cOXB1dYW+vj7Mzc3RsWNHfPTRR0hNTa3za588eRICgQAWFhYoKSl5nW7UiOtjU3n48CEEAgH27t3LHwsPD8eKFSuQl5fXZK/b1Pbu3QuBQICHDx826XP79esn8xlrTP369YNAIMCwYcNkHuN+d999953c59blM1pcXIx169bB29sbxsbGEAqFaNu2LSZMmIALFy5InRsXF4cpU6agTZs20NPTg6WlJXx8fDBv3jwUFBTUqT+rVq2Cp6cnxGJxnfrw3XffyfwuysrK8PPPP6Nbt24wNzeHgYEBnJ2d8fbbb+PYsWMy7w/3o62tDQsLC3Tr1g0LFy7EnTt3amzrokWLIBAIMGrUKLmP37t3Dzo6OoiOjq5T34lyo+CGqKzHjx/Dx8cHZ8+exaJFi3D69Gns3r0bkyZNwo0bN/DgwYM6X2vXrl0AgJycHBw/fryJWgx8+OGHiIiIaLLryxMeHo6VK1eqdHAzcuRIREREwM7OTtFNaRRnzpzB+fPn6/Wc2j6jFRUVGDJkCNasWYPx48fj8OHDOHLkCBYuXIj8/HxcunSJPzcmJgZdu3bF3bt38fXXXyM4OBjbt2/HyJEjcebMGeTk5NTanidPnmD9+vVYtWoVNDQafiuZMmUK5s+fj/79+2P//v343//+h2XLlkFLSwtnzpyROX/+/PmIiIjAhQsXsG/fPowZMwYnT56Et7c3NmzYIPc1ysrKsH//fgBAcHAw0tLSZM5xdXWFv78/Fi5c2OC+ECXCCFESy5cvZwDY06dP5T7eoUMH1rdvX/7fX3/9NQPAHjx4IPf8ioqKOr2uSCRiWlpabMCAAUxPT48NHjy43m2vTXFxcaNfU57k5GQGgO3Zs4c/tmHDBgaAJScnN+prPXv2jInF4ka9ZlPYs2dPnfvft29fqc9YY+vbty9zdXVlbdq0YV27dpV6/7jf3YYNG2SeV5fP6Pnz5xkAtnv3brmvLfnfw9SpU5mhoSErKCiQe25dfq+fffYZa9WqldR1a+oDY7KfxQcPHjAA7Ouvv661zTVd+9mzZ2zYsGEMADt9+rTM44cPH2YA2MiRIxkAtmbNGrmvFxkZyQCwK1euVNtvohpo5IaorOzsbGhoaMDa2lru43X9a/LXX39FeXk5Fi5ciLFjxyIkJASPHj2SOS8vLw8ffPABzM3NYWRkhJEjR+LBgwcQCARYsWIFfx439RQdHY3x48fDzMwMbdu2lXqsqt9//x1+fn4wMjKCkZEROnfuzP+lDgAuLi6YPn26zPNqm0ZZsWIFPv30UwBA69at+SH9sLAwAJBpe3Wvx03t/PPPP3j//fdhZWUFAwMDfnrk4MGD8PPzg6GhIYyMjDB06FDExMRU2y6gcppBS0sLa9eulXns4sWLEAgEOHz4sNTrV51aOnfuHAYOHAhjY2MYGBigV69eCAkJqfF1AYAxhvXr18PZ2Rl6enrw8fHB33//LXOeWCzG6tWr4ebmBn19fZiamqJTp0744Ycfan2N6mhra2PNmjWIiorCwYMH6/ScunxGs7OzAaDa0S3J/x6ys7NhbGwMIyMjuefWNnVaWlqKXbt2YfLkya81alOfNtdEX18fu3btgra2ttzRm127dkFHRwd79uyBo6Mj9uzZAyZnz+iuXbvCw8MD27dvr0cviDKi4IaoLD8/P4jFYowdOxZnzpypc55AVbt374adnR2GDx+O999/H2KxWCpnBai8yY0ePRq///47Pv/8cxw7dgw9evSQmzvBGTt2LNq1a4fDhw/X+GX59ddfw9/fH/b29ti7dy+OHTuGadOmyQ2w6uvDDz/E/PnzAQBHjx5FREQEIiIi4OPj06Drvf/++9DW1sa+fftw5MgRaGtr45tvvsGkSZPg6emJQ4cOYd++fSgsLESfPn1w9+7daq/l4uKCt956C9u3b0dFRYXUY1u3boW9vT3eeeedap+/f/9+DBkyBMbGxvj1119x6NAhmJubY+jQobUGOCtXrsTnn3+OwYMH4/jx4wgICMBHH32EhIQEqfPWr1+PFStWYNKkSTh16hQOHjyIDz744LWn+CZOnIiuXbti2bJlKCsrq/X8unxGfX19oa2tjU8++QRBQUEQiUTVXs/Pzw8ikQj+/v64cOECnj9/Xq/2X7t2DdnZ2ejfv3+9nleVh4cHTE1NsXLlSuzYsaNBOVUce3t7dO3aFeHh4SgvL+ePP378GP/88w/efvttWFlZYdq0aUhKSsLFixflXqdfv374+++/5QY/RIUoeOSIEF59p6XEYjGbNWsW09DQYACYQCBgHh4ebOHChXWegrl48SIDwL744gv+mq1bt2bOzs5SQ/OnTp1iAFhgYKDU89euXcsAsOXLl8v0Q95QO/cY58GDB0xTU5P5+/vX2E5nZ2c2bdo0meNVp1HqOy1Vte3VvR43tTN16lSp81JSUpiWlhabP3++1PHCwkJma2vLJkyYUGO/QkNDGQB27Ngx/lhaWhrT0tJiK1eulHl9rg/FxcXM3NycjR49Wup6FRUVzNvbm3Xv3r3a5+bm5jI9PT32zjvvSD33ypUrDIDU+zlq1CjWuXPnGvtQH3379mUdOnRgjDF27tw5BoD9+OOPjLHqp13q+hlljLFdu3YxIyMjBoABYHZ2dmzq1Kns4sWLUue9ePGCjRkzhj9PU1OTdenShS1dupRlZmbW2o9169YxACw9PV3qeH2npRir/G/L0tKSb4uFhQV799132cmTJ+t1bcYYmzhxIgPAMjIy+GOrVq1iAFhwcDBjrPK/OYFAwKZMmSL3Gr/88gsDwOLi4mp8D4hya9EjNxcvXsTo0aNhb28PgUBQ70RSboqh6o+hoWHTNJhIEQgE2L59Ox48eIBt27ZhxowZKCsrw6ZNm9ChQweZ1SHycFM/77//Pn/N6dOn49GjR1J//XPXmjBhgtTzJ02aVO21x40bV+vrnz17FhUVFZg7d26t5yqDqn06c+YMysvLMXXqVJSXl/M/enp66Nu3Lz/9VZ1+/frB29sbP/30E39s+/btEAgEmDlzZrXPCw8PR05ODqZNmyb1umKxGMOGDcONGzdQXFws97kRERF48eIF/P39pY737NkTzs7OUse6d++Of//9F3PmzHmt0UF5Bg4ciCFDhmDVqlUoLCys9ry6fka5cx4/fozff/8dH3/8MRwdHbF//3707dtXarpGV1cXx44dw927d7Fp0ya89957ePr0KdasWQMPDw+ZEayqnjx5AoFAAEtLy4Z2nzdixAikpKTg2LFjWLx4MTp06IDjx4/jrbfewrx58+p1LVZltIUxxk9FDR48GEDl9Gy/fv3w559/yv19ctPc8pKOiepo0cFNcXExvL29sXXr1gY9f/HixRCJRFI/np6eePfddxu5pS2DlpYWAMhMUXDKy8uhra0tc9zZ2RkBAQHYtWsXEhMTcfDgQbx48YLPNalOYWEhDh8+jO7du8PKygp5eXnIy8vDO++8A4FAIJXzkp2dDS0tLZibm0tdw8bGptrr12Vlz9OnTwEADg4OtZ6rDKr2KSMjAwDQrVs3aGtrS/0cPHiw2mX9kj7++GOEhIQgISEBZWVl+OWXXzB+/HjY2tpW+xzudcePHy/zuuvWrQNjrNo
|
2023-01-12 21:45:54 +02:00
|
|
|
"text/plain": [
|
|
|
|
|
"<Figure size 640x480 with 1 Axes>"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"output_type": "display_data"
|
|
|
|
|
}
|
|
|
|
|
],
|
2023-01-12 16:13:39 +02:00
|
|
|
"source": [
|
2023-01-13 20:54:15 +02:00
|
|
|
"plt.plot(corn_year,corn_kcal_per_acre,\"o\",markerfacecolor=\"none\",label=\"corn\")\n",
|
2023-01-12 16:13:39 +02:00
|
|
|
"plt.xlabel(\"time (calendar year)\")\n",
|
2023-01-12 20:25:12 +02:00
|
|
|
"plt.ylabel(\"kcal/acre\")\n",
|
2023-01-12 16:13:39 +02:00
|
|
|
"plt.title(\"US Agriculture yields, NASS (USDA)\")\n",
|
|
|
|
|
"plt.legend()\n",
|
|
|
|
|
"\n",
|
|
|
|
|
"#plt.show()\n",
|
|
|
|
|
"#plt.savefig(\"bear-quadratic.pdf\")"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "markdown",
|
2023-01-12 20:25:12 +02:00
|
|
|
"id": "fa891f0e",
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"source": [
|
|
|
|
|
"## Wheat"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2023-01-12 21:45:54 +02:00
|
|
|
"execution_count": 6,
|
2023-01-12 20:25:12 +02:00
|
|
|
"id": "263ea341",
|
|
|
|
|
"metadata": {},
|
2023-01-12 21:45:54 +02:00
|
|
|
"outputs": [
|
|
|
|
|
{
|
|
|
|
|
"data": {
|
|
|
|
|
"text/plain": [
|
|
|
|
|
"array([[ nan, 2022. , nan, nan, nan, nan, nan, nan,\n",
|
|
|
|
|
" nan, nan, nan, nan, nan, 0. , nan, nan,\n",
|
|
|
|
|
" nan, nan, 46.5, nan],\n",
|
|
|
|
|
" [ nan, 2021. , nan, nan, nan, nan, nan, nan,\n",
|
|
|
|
|
" nan, nan, nan, nan, nan, 0. , nan, nan,\n",
|
|
|
|
|
" nan, nan, 44.3, nan],\n",
|
|
|
|
|
" [ nan, 2020. , nan, nan, nan, nan, nan, nan,\n",
|
|
|
|
|
" nan, nan, nan, nan, nan, 0. , nan, nan,\n",
|
|
|
|
|
" nan, nan, 49.7, nan]])"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
"execution_count": 6,
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"output_type": "execute_result"
|
|
|
|
|
}
|
|
|
|
|
],
|
2023-01-12 20:25:12 +02:00
|
|
|
"source": [
|
|
|
|
|
"csv=np.genfromtxt ('./data/wheat.csv',\n",
|
|
|
|
|
" skip_header=1,\n",
|
|
|
|
|
" delimiter=\",\")\n",
|
|
|
|
|
"csv[0:3]\n",
|
|
|
|
|
"#csv[0:3].astype(float)"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2023-01-12 21:45:54 +02:00
|
|
|
"execution_count": 7,
|
2023-01-12 20:25:12 +02:00
|
|
|
"id": "5aa38243",
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [],
|
|
|
|
|
"source": [
|
|
|
|
|
"wheat_year=csv[:,1]\n",
|
|
|
|
|
"wheat_bu_per_acre=csv[:,-2]\n",
|
|
|
|
|
"\n",
|
|
|
|
|
"wheat_lbs_per_bu=60\n",
|
|
|
|
|
"#grams_per_lbs\n",
|
|
|
|
|
"wheat_kcal_per_gram=327.0/100.0\n",
|
|
|
|
|
"\n",
|
|
|
|
|
"wheat_kcal_per_acre=wheat_bu_per_acre*wheat_lbs_per_bu*grams_per_lbs*wheat_kcal_per_gram"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2023-01-12 21:45:54 +02:00
|
|
|
"execution_count": 8,
|
2023-01-12 20:25:12 +02:00
|
|
|
"id": "5dd77124",
|
|
|
|
|
"metadata": {},
|
2023-01-12 21:45:54 +02:00
|
|
|
"outputs": [
|
|
|
|
|
{
|
|
|
|
|
"data": {
|
2023-01-13 20:54:15 +02:00
|
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAjoAAAHFCAYAAAD7ZFORAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAA9hAAAPYQGoP6dpAAC7+UlEQVR4nOydeVwU9f/HX8uNcl+CyqUoqIiIN5qI921maoZnpb/UMu2wTC31K1lY2YWYpllGZodXad6gpigioikIigioIMghiAiyO78/thlmd2cv2GUXeD8fDx7AnJ+ZHfi85n2KGIZhQBAEQRAE0QQxMfQACIIgCIIg9AUJHYIgCIIgmiwkdAiCIAiCaLKQ0CEIgiAIoslCQocgCIIgiCYLCR2CIAiCIJosJHQIgiAIgmiykNAhCIIgCKLJQkKHIAiCIIgmCwkdwihZtWoVRCIRHjx4ILg+MDAQgwYNklmWm5uLBQsWoGPHjrC2toaTkxO6du2KuXPnIjc3V+Nz79+/HyKRCM7OzqiqqqrPZaiEvUZ9cfv2bYhEImzfvp1bdvbsWaxatQqlpaV6O6++2b59O0QiEW7fvq3XfQcNGqTwjOmSQYMGQSQSYeTIkQrr2M/u008/FdxXk2e0oqICn3zyCbp16wY7OzvY2tqiffv2mDJlCk6ePCmzbVpaGmbMmIF27drBysoKLi4uCAkJwWuvvYaysjKNrmfNmjXo3LkzJBKJRtfw6aefKnwWT58+xbfffotevXrByckJLVq0gLe3NyZMmIA9e/Yo3B/2y9zcHM7OzujVqxeWLFmCa9euqRzrm2++CZFIhLFjxwquz8jIgIWFBZKTkzW6dsK4IaFDNAnu3LmDkJAQHD16FG+++SYOHjyIbdu2Ydq0abhw4QJu3bql8bG2bt0KACguLsbevXv1NGLglVdeQUJCgt6OL8TZs2exevXqRi10xowZg4SEBHh4eBh6KDrh8OHDOHHihFb7qHtGxWIxhg8fjsjISDz//PP47bff8Pvvv2PJkiV4+PAhTp8+zW176dIl9OjRA6mpqfjggw9w6NAhbNq0CWPGjMHhw4dRXFysdjz37t1DVFQU1qxZAxOTuk8rM2bMwOuvv47w8HD89NNP+PPPP7FixQqYmZnh8OHDCtu//vrrSEhIwMmTJ7Fjxw48++yz2L9/P7p164b169cLnuPp06f46aefAACHDh3C3bt3Fbbp2LEjIiIisGTJkjpfC2FEMARhhHz44YcMAKawsFBwfZcuXZiwsDDu9w8++IABwNy6dUtwe7FYrNF58/LyGDMzM2bw4MGMlZUVM2zYMK3Hro6KigqdH1OIrKwsBgDz/fffc8vWr1/PAGCysrJ0eq7Hjx8zEolEp8fUB99//73G1x8WFibzjOmasLAwpmPHjky7du2YHj16yNw/9rNbv369wn6aPKMnTpxgADDbtm0TPDf/72HmzJlMy5YtmbKyMsFtNflcly5dyrRp00bmuKqugWEUn8Vbt24xAJgPPvhA7ZhVHfvx48fMyJEjGQDMwYMHFdb/9ttvDABmzJgxDAAmMjJS8HxJSUkMAObMmTNKr5toHJBFh2gSFBUVwcTEBG5uboLrNX3L/OGHH1BTU4MlS5bgueeew/Hjx5Gdna2wXWlpKV5++WU4OTnBxsYGY8aMwa1btyASibBq1SpuO9Y9lZycjOeffx6Ojo5o3769zDp5fv75Z/Tr1w82NjawsbFBcHAw9wYPAD4+Ppg9e7bCfupcLatWrcI777wDAPD19eXM/vHx8QCgMHZl52PdP0eOHMFLL70EV1dXtGjRgnOh7Nq1C/369UPLli1hY2ODESNG4NKlS0rHBUhdEWZmZli3bp3CulOnTkEkEuG3336TOb+8++nYsWMYMmQI7Ozs0KJFC/Tv3x/Hjx9XeV4AYBgGUVFR8Pb2hpWVFUJCQvD3338rbCeRSLB27Vr4+/vD2toaDg4OCAoKwpdffqn2HMowNzdHZGQkLl68iF27dmm0jybPaFFREQAotXrx/x6KiopgZ2cHGxsbwW3VuVerq6uxdetWvPjii/Wy5mgzZlVYW1tj69atMDc3F7TqbN26FRYWFvj+++/h6emJ77//HoxAb+sePXqgU6dO2LRpkxZXQRgjJHSIJkG/fv0gkUjw3HPP4fDhwxrHFcizbds2eHh4YNSoUXjppZcgkUhkYlwA6YQ3btw4/Pzzz3j33XexZ88e9OnTRzDWguW5556Dn58ffvvtN5X/OD/44ANERESgdevW2L59O/bs2YNZs2YJii1teeWVV/D6668DAHbv3o2EhAQkJCQgJCSkTsd76aWXYG5ujh07duD333+Hubk5PvroI0ybNg2dO3fGr7/+ih07dqC8vBzPPPMMUlNTlR7Lx8cH48ePx6ZNmyAWi2XWffPNN2jdujUmTpyodP+ffvoJw4cPh52dHX744Qf8+uuvcHJywogRI9SKndWrV+Pdd9/FsGHDsHfvXsyfPx9z585Fenq6zHZRUVFYtWoVpk2bhgMHDmDXrl14+eWX6+0GnDp1Knr06IEVK1bg6dOnarfX5Bnt2bMnzM3N8cYbbyA2NhZ5eXlKj9evXz/k5eUhIiICJ0+eRGVlpVbjP3/+PIqKihAeHq7VfvJ06tQJDg4OWL16NTZv3lynGCyW1q1bo0ePHjh79ixqamq45Xfu3MGRI0cwYcIEuLq6YtasWbh58yZOnToleJxBgwbh77//FhRCRCPCwBYlghBEW9eVRCJh/u///o8xMTFhADAikYjp1KkTs2TJEo3dNKdOnWIAMO+99x53TF9fX8bb21vGfH/gwAEGABMTEyOz/7p16xgAzIcffqhwHULmeHYdy61btxhTU1MmIiJC5Ti9vb2ZWbNmKSyXd7Vo67qSH7uy87Hun5kzZ8psl5OTw5iZmTGvv/66zPLy8nLG3d2dmTJlisrriouLYwAwe/bs4ZbdvXuXMTMzY1avXq1wfvYaKioqGCcnJ2bcuHEyxxOLxUy3bt2Y3r17K923pKSEsbKyYiZOnCiz75kzZxgAMvdz7NixTHBwsMpr0IawsDCmS5cuDMMwzLFjxxgAzNdff80wjHLXjKbPKMMwzNatWxkbGxsGAAOA8fDwYGbOnMmcOnVKZrsnT54wzz77LLedqakp0717d2b58uVMQUGB2uv45JNPGABMfn6+zHJtXVcMI/3bcnFx4cbi7OzMTJ48mdm/f79Wx2YYhpk6dSoDgLl//z63bM2aNQwA5tChQwzDSP/mRCIRM2PGDMFjbNmyhQHApKWlqbwHhHFDFh2iSSASibBp0ybcunULGzduxJw5c/D06VNs2LABXbp0UcgyEYJ1D7300kvcMWfPno3s7GwZqwB7rClTpsjsP23aNKXHnjRpktrzHz16FGKxGAsXLlS7rTEgf02HDx9GTU0NZs6ciZqaGu7LysoKYWFhnItMGYMGDUK3bt0QHR3NLdu0aRNEIhHmzZundL+zZ8+iuLgYs2bNkjmvRCLByJEjceHCBVRUVAjum5CQgCdPniAiIkJmeWhoKLy9vWWW9e7dG5cvX8aCBQvqZTUUYsiQIRg+fDjWrFmD8vJypdtp+oyy29y5cwc///wzFi1aBE9PT/z0008ICwuTcelYWlpiz549SE1NxYYNG/DCCy+gsLAQkZGR6NSpk4JlS5579+5BJBLBxcWlrpfPMXr0aOTk5GDPnj14++230aVLF+zduxfjx4/Ha6+9ptWxGDkrDMMwnLtq2LBhAKQu3EGDBuGPP/4Q/DxZV7hQwDLReCChQxglZmZmAKDgxmCpqamBubm5wnJvb2/Mnz8fW7duxY0bN7Br1y48efKEi01RRnl5OX777Tf07t0brq6uKC0tRWlpKSZOnAiRSCQTI1NUVAQzMzM4OTnJHKNVq1ZKj69JhlBhYSEAoG3btmq3NQbkr+n+/fsAgF69esHc3Fzma9euXUpLBfBZtGgRjh8/jvT0dDx9+hRbtmzB888/D3d3d6X7sOd9/vnnFc77ySefgGEYpZlDbFyI0PHlly1btgyffvopzp07h1GjRsHZ2RlDhgxBUlKS2uvShE8
|
2023-01-12 21:45:54 +02:00
|
|
|
"text/plain": [
|
|
|
|
|
"<Figure size 640x480 with 1 Axes>"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"output_type": "display_data"
|
|
|
|
|
}
|
|
|
|
|
],
|
2023-01-12 20:25:12 +02:00
|
|
|
"source": [
|
2023-01-13 20:54:15 +02:00
|
|
|
"plt.plot(corn_year,corn_bu_per_acre,\"o\",markerfacecolor=\"none\",label=\"Corn (BU/ACRE)\")\n",
|
2023-01-12 20:25:12 +02:00
|
|
|
"plt.plot(wheat_year,wheat_bu_per_acre,\"x\",label=\"Wheat (BU/ACRE)\")\n",
|
|
|
|
|
"plt.xlabel(\"time (calendar year)\")\n",
|
|
|
|
|
"plt.ylabel(\"Average per acre production\")\n",
|
|
|
|
|
"plt.title(\"US Agriculture yields, NASS (USDA)\")\n",
|
|
|
|
|
"plt.legend()\n",
|
|
|
|
|
"plt.yscale(\"log\")\n",
|
|
|
|
|
"#plt.show()\n",
|
|
|
|
|
"#plt.savefig(\"bear-quadratic.pdf\")"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "markdown",
|
|
|
|
|
"id": "89c8f489",
|
2023-01-12 16:13:39 +02:00
|
|
|
"metadata": {},
|
|
|
|
|
"source": [
|
|
|
|
|
"## Potatoes"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2023-01-12 21:45:54 +02:00
|
|
|
"execution_count": 9,
|
2023-01-12 20:25:12 +02:00
|
|
|
"id": "c31828d0",
|
2023-01-12 16:13:39 +02:00
|
|
|
"metadata": {},
|
2023-01-12 21:45:54 +02:00
|
|
|
"outputs": [
|
|
|
|
|
{
|
|
|
|
|
"data": {
|
|
|
|
|
"text/plain": [
|
|
|
|
|
"array([[ nan, 2022., nan, nan, nan, nan, nan, nan, nan,\n",
|
|
|
|
|
" nan, nan, nan, nan, 0., nan, nan, nan, nan,\n",
|
|
|
|
|
" 438., nan],\n",
|
|
|
|
|
" [ nan, 2021., nan, nan, nan, nan, nan, nan, nan,\n",
|
|
|
|
|
" nan, nan, nan, nan, 0., nan, nan, nan, nan,\n",
|
|
|
|
|
" 444., nan],\n",
|
|
|
|
|
" [ nan, 2020., nan, nan, nan, nan, nan, nan, nan,\n",
|
|
|
|
|
" nan, nan, nan, nan, 0., nan, nan, nan, nan,\n",
|
|
|
|
|
" 461., nan]])"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
"execution_count": 9,
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"output_type": "execute_result"
|
|
|
|
|
}
|
|
|
|
|
],
|
2023-01-12 16:13:39 +02:00
|
|
|
"source": [
|
|
|
|
|
"csv=np.genfromtxt ('./data/potatoes.csv',\n",
|
|
|
|
|
" skip_header=1,\n",
|
|
|
|
|
" delimiter=\",\")\n",
|
|
|
|
|
"csv[0:3]\n",
|
|
|
|
|
"#csv[0:3].astype(float)"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2023-01-12 21:45:54 +02:00
|
|
|
"execution_count": 10,
|
2023-01-12 20:25:12 +02:00
|
|
|
"id": "86091fed",
|
2023-01-12 16:13:39 +02:00
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [],
|
|
|
|
|
"source": [
|
|
|
|
|
"potatoes_year=csv[:,1]\n",
|
2023-01-12 20:25:12 +02:00
|
|
|
"potatoes_cwt_per_acre=csv[:,-2]\n",
|
|
|
|
|
"\n",
|
|
|
|
|
"lbs_per_cwt=100\n",
|
|
|
|
|
"#grams_per_lbs\n",
|
|
|
|
|
"potatoes_kcal_per_gram=77.0/100.0\n",
|
|
|
|
|
"potatoes_kcal_per_acre=potatoes_cwt_per_acre*lbs_per_cwt*grams_per_lbs*potatoes_kcal_per_gram"
|
2023-01-12 16:13:39 +02:00
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2023-01-12 21:45:54 +02:00
|
|
|
"execution_count": 11,
|
2023-01-12 20:25:12 +02:00
|
|
|
"id": "5ea8b2ac",
|
2023-01-12 16:13:39 +02:00
|
|
|
"metadata": {},
|
2023-01-12 21:45:54 +02:00
|
|
|
"outputs": [
|
|
|
|
|
{
|
|
|
|
|
"data": {
|
2023-01-13 20:54:15 +02:00
|
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAjoAAAHFCAYAAAD7ZFORAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAA9hAAAPYQGoP6dpAADitElEQVR4nOydeVhUdd+H72FHWQQEcUFEEUVx31Ir3DW1rNwyzWzzDW0zy6LVejRLM21RDHN5NDPSsix71FAxSxQV0RRcQUAF2RUBUZnz/jGd4wwzAzM47L/7uuZy5ne27xkG58N3VUmSJCEQCAQCgUBQB7GqbgMEAoFAIBAIKgshdAQCgUAgENRZhNARCAQCgUBQZxFCRyAQCAQCQZ1FCB2BQCAQCAR1FiF0BAKBQCAQ1FmE0BEIBAKBQFBnEUJHIBAIBAJBnUUIHYFAIBAIBHUWIXQENZK5c+eiUqnIysoyuD0oKIgBAwborKWmpjJjxgwCAgJwdHTE3d2dTp068dxzz5Gammrytbdu3YpKpcLDw4Pi4uK7uY0yke+xsrhw4QIqlYq1a9cqa/v372fu3Lnk5eVV2nUrm7Vr16JSqbhw4UKlHjtgwAC9z5glGTBgACqVihEjRuhtk392n376qcFjTfmMFhQU8Mknn9ClSxdcXFxwdnamTZs2TJgwgb179+rsm5CQwBNPPEHr1q1xcHCgcePGdO/enRdeeIFr166ZdD8ffvghHTp0QK1Wm3QPn376qd7P4tatW3z99df06tULd3d3GjRogK+vL2PGjGHLli1674/8sLW1xcPDg169ejFr1ixOnjxZpq2vvvoqKpWK0aNHG9x+5swZ7OzsiI2NNeneBTUbIXQEdYKLFy/SvXt3/vjjD1599VV+//13Vq9ezaRJkzh06BCJiYkmn2vVqlUA5OTk8PPPP1eSxfDss88SHR1daec3xP79+/nggw9qtdAZNWoU0dHRNG3atLpNsQg7duxg9+7dZh1T3me0pKSEYcOGMX/+fMaNG8emTZvYvHkzs2bN4urVq+zbt0/Z9+jRo/To0YP4+Hjee+89tm/fzooVKxg1ahQ7duwgJyenXHsuX77MwoUL+fDDD7GyqvjXyhNPPMGLL77IwIED+fbbb/n111955513sLGxYceOHXr7v/jii0RHR7N3717Wr1/Pww8/zNatW+nSpQuLFi0yeI1bt27x7bffArB9+3YuXbqkt09AQACTJ09m1qxZFb4XQQ1CEghqIO+//74ESJmZmQa3d+zYUQoODlZev/feexIgJSYmGty/pKTEpOumpaVJNjY20qBBgyQHBwdp6NChZtteHgUFBRY/pyGSkpIkQFqzZo2ytmjRIgmQkpKSLHqtwsJCSa1WW/SclcGaNWtMvv/g4GCdz5ilCQ4OlgICAqTWrVtLPXr00Hn/5J/dokWL9I4z5TO6e/duCZBWr15t8Nravw9Tp06VGjZsKF27ds3gvqb8XOfMmSM1b95c57xl3YMk6X8WExMTJUB67733yrW5rHMXFhZKI0aMkADp999/19u+adMmCZBGjRolAdL8+fMNXu/w4cMSIP39999G71tQOxAeHUGdIDs7GysrK7y8vAxuN/WvzP/+97/cvn2bWbNm8eijj7Jr1y6Sk5P19svLy+OZZ57B3d0dJycnRo0aRWJiIiqVirlz5yr7yeGp2NhYxo0bh5ubG23atNHZVprvvvuOvn374uTkhJOTE127dlX+ggdo1aoV06ZN0zuuvFDL3Llzef311wHw8/NT3P5RUVEAerYbu54c/tm5cydPP/00np6eNGjQQAmhRERE0LdvXxo2bIiTkxPDhw/n6NGjRu0CTSjCxsaGBQsW6G37888/UalUbNq0Sef6pcNPkZGRDB48GBcXFxo0aED//v3ZtWtXmdcFkCSJhQsX4uvri4ODA927d+d///uf3n5qtZp58+bRrl07HB0dadSoEZ07d+bzzz8v9xrGsLW1Zf78+Rw5coSIiAiTjjHlM5qdnQ1g1Oul/fuQnZ2Ni4sLTk5OBvctL7x68+ZNVq1axeOPP35X3hxzbC4LR0dHVq1aha2trUGvzqpVq7Czs2PNmjX4+PiwZs0aJAOzrXv06EFgYCArVqww4y4ENREhdAR1gr59+6JWq3n00UfZsWOHyXkFpVm9ejVNmzblgQce4Omnn0atVuvkuIDmC+/BBx/ku+++44033mDLli306dPHYK6FzKOPPoq/vz+bNm0q8z/O9957j8mTJ9OsWTPWrl3Lli1bePLJJw2KLXN59tlnefHFFwH46aefiI6OJjo6mu7du1fofE8//TS2trasX7+ezZs3Y2try0cffcSkSZPo0KEDP/zwA+vXryc/P5/77ruP+Ph4o+dq1aoVDz30ECtWrKCkpERn21dffUWzZs145JFHjB7/7bffMmzYMFxcXPjvf//LDz/8gLu7O8OHDy9X7HzwwQe88cYbDB06lJ9//pmQkBCee+45Tp8+rbPfwoULmTt3LpMmTWLbtm1ERETwzDPP3HUYcOLEifTo0YN33nmHW7dulbu/KZ/Rnj17Ymtry8svv8yGDRtIS0szer6+ffuSlpbG5MmT2bt3L0VFRWbZf/DgQbKzsxk4cKBZx5UmMDCQRo0a8cEHHxAeHl6hHCyZZs2a0aNHD/bv38/t27eV9YsXL7Jz507GjBmDp6cnTz75JOfOnePPP/80eJ4BAwbwv//9z6AQEtQiqtmjJBAYxNzQlVqtlv7v//5PsrKykgBJpVJJgYGB0qxZs0wO0/z5558SIL355pvKOf38/CRfX18d9/22bdskQAoLC9M5fsGCBRIgvf/++3r3YcgdL2+TSUxMlKytraXJkyeXaaevr6/05JNP6q2XDrWYG7oqbbux68nhn6lTp+rsl5KSItnY2Egvvviiznp+fr7k7e0tTZgwocz72rNnjwRIW7ZsUdYuXbok2djYSB988IHe9eV7KCgokNzd3aUHH3xQ53wlJSVSly5dpN69exs9Njc3V3JwcJAeeeQRnWP//vtvCdB5P0ePHi117dq1zHswh+DgYKljx46SJElSZGSkBEhffvmlJEnGQzOmfkYlSZJWrVolOTk5SYAESE2bNpWmTp0q/fnnnzr73bhxQ3r44YeV/aytraVu3bpJb7/9tpSRkVHufXzyyScSIKWnp+usmxu6kiTN71bjxo0VWzw8PKTx48dLW7duNevckiRJEydOlADpypUrytqHH34oAdL27dslSdL8zqlUKumJJ54weI6VK1dKgJSQkFDmeyCo2QiPjqBOoFKpWLFiBYmJiSxfvpynnnqKW7dusWTJEjp27KhXZWIIOTz09NNPK+ecNm0aycnJOl4B+VwTJkzQOX7SpElGzz127Nhyr//HH39QUlLCzJkzy923JlD6nnbs2MHt27eZOnUqt2/fVh4ODg4EBwcrITJjDBgwgC5durBs2TJlbcWKFahUKqZPn270uP3795OTk8OTTz6pc121Ws2IESM4dOgQBQUFBo+Njo7mxo0bTJ48WWe9X79++Pr66qz17t2bY8eOMWPGjLvyGhpi8ODBDBs2jA8//JD8/Hyj+5n6GZX3uXjxIt999x0vvfQSPj4+fPvttwQHB+uEdOzt7dmyZQvx8fEsWbKExx57jMzMTObPn09gYKCeZ6s0ly9fRqVS0bhx44revsLIkSNJSUlhy5YtvPbaa3Ts2JGff/6Zhx56iBdeeMGsc0mlvDCSJCnhqqFDhwKaEO6AAQP48ccfDf485VC4oYRlQe1BCB1BjcTGxgZAL4whc/v2bWxtbfXWfX19CQkJYdWqVZw9e5aIiAhu3Lih5KYYIz8/n02bNtG7d288PT3Jy8sjLy+PRx55BJVKpZMjk52djY2NDe7u7jrnaNKkidHzm1IhlJmZCUCLFi3K3bcmUPqerly5AkCvXr2wtbXVeURERBhtFaDNSy+9xK5duzh9+jS3bt1i5cqVjBs3Dm9vb6PHyNcdN26c3nU/+eQTJEkyWjkk54UYOn/ptdDQUD799FMOHDjAAw88gIeHB4MHD+bw4cPl3pcpfPLJJ2RlZRktxzb
|
2023-01-12 21:45:54 +02:00
|
|
|
"text/plain": [
|
|
|
|
|
"<Figure size 640x480 with 1 Axes>"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"output_type": "display_data"
|
|
|
|
|
}
|
|
|
|
|
],
|
2023-01-12 20:25:12 +02:00
|
|
|
"source": [
|
2023-01-13 20:54:15 +02:00
|
|
|
"plt.plot(corn_year,corn_bu_per_acre,\"o\",markerfacecolor=\"none\",label=\"Corn (BU/ACRE)\")\n",
|
2023-01-12 20:25:12 +02:00
|
|
|
"plt.plot(wheat_year,wheat_bu_per_acre,\"x\",label=\"Wheat (BU/ACRE)\")\n",
|
|
|
|
|
"plt.plot(potatoes_year,potatoes_cwt_per_acre,\"+\",label=\"Potatoes (CWT / ACRE)\")\n",
|
|
|
|
|
"plt.xlabel(\"time (calendar year)\")\n",
|
|
|
|
|
"plt.ylabel(\"Average per acre production\")\n",
|
|
|
|
|
"plt.title(\"US Agriculture yields, NASS (USDA)\")\n",
|
|
|
|
|
"plt.legend()\n",
|
|
|
|
|
"plt.yscale(\"log\")\n",
|
|
|
|
|
"#plt.show()\n",
|
|
|
|
|
"#plt.savefig(\"bear-quadratic.pdf\")"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "markdown",
|
|
|
|
|
"id": "0201b588",
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"source": [
|
|
|
|
|
"## Soybeans"
|
|
|
|
|
]
|
2023-01-12 16:13:39 +02:00
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2023-01-12 21:45:54 +02:00
|
|
|
"execution_count": 12,
|
2023-01-12 20:25:12 +02:00
|
|
|
"id": "ce223c5b",
|
2023-01-12 16:13:39 +02:00
|
|
|
"metadata": {},
|
2023-01-12 21:45:54 +02:00
|
|
|
"outputs": [
|
|
|
|
|
{
|
|
|
|
|
"data": {
|
|
|
|
|
"text/plain": [
|
|
|
|
|
"array([[ nan, 2022. , nan, nan, nan, nan, nan, nan,\n",
|
|
|
|
|
" nan, nan, nan, nan, nan, 0. , nan, nan,\n",
|
|
|
|
|
" nan, nan, 50.2, nan],\n",
|
|
|
|
|
" [ nan, 2021. , nan, nan, nan, nan, nan, nan,\n",
|
|
|
|
|
" nan, nan, nan, nan, nan, 0. , nan, nan,\n",
|
|
|
|
|
" nan, nan, 51.7, nan],\n",
|
|
|
|
|
" [ nan, 2020. , nan, nan, nan, nan, nan, nan,\n",
|
|
|
|
|
" nan, nan, nan, nan, nan, 0. , nan, nan,\n",
|
|
|
|
|
" nan, nan, 51. , nan]])"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
"execution_count": 12,
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"output_type": "execute_result"
|
|
|
|
|
}
|
|
|
|
|
],
|
2023-01-12 16:13:39 +02:00
|
|
|
"source": [
|
2023-01-12 20:25:12 +02:00
|
|
|
"csv=np.genfromtxt ('./data/soybeans.csv',\n",
|
|
|
|
|
" skip_header=1,\n",
|
|
|
|
|
" delimiter=\",\")\n",
|
|
|
|
|
"csv[0:3]\n",
|
|
|
|
|
"#csv[0:3].astype(float)"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2023-01-12 21:45:54 +02:00
|
|
|
"execution_count": 13,
|
2023-01-12 20:25:12 +02:00
|
|
|
"id": "12fe395c",
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [],
|
|
|
|
|
"source": [
|
|
|
|
|
"soybeans_year=csv[:,1]\n",
|
2023-01-12 21:45:54 +02:00
|
|
|
"soybeans_bu_per_acre=csv[:,-2]\n",
|
|
|
|
|
"\n",
|
|
|
|
|
"soybeans_lbs_per_bu=60\n",
|
|
|
|
|
"#grams_per_lbs\n",
|
|
|
|
|
"soybeans_kcal_per_gram=446.0/100.0\n",
|
|
|
|
|
"soybeans_kcal_per_acre=soybeans_bu_per_acre*soybeans_lbs_per_bu*grams_per_lbs*soybeans_kcal_per_gram"
|
2023-01-12 20:25:12 +02:00
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2023-01-12 21:45:54 +02:00
|
|
|
"execution_count": 14,
|
2023-01-12 20:25:12 +02:00
|
|
|
"id": "ae3f6ad4",
|
|
|
|
|
"metadata": {},
|
2023-01-12 21:45:54 +02:00
|
|
|
"outputs": [
|
|
|
|
|
{
|
|
|
|
|
"data": {
|
2023-01-13 20:54:15 +02:00
|
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAjoAAAHFCAYAAAD7ZFORAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOydeVhU5duA72FHNlnFBXFFQERF09QSF1xyrVzILLPNLzUrtSyyhcolNdMWxTTUn0tKWpplqaFiLihuuIHigogKgmwiArKc749pTjPMDDAIsvje13UuZ97znvc85zDjeeZZFZIkSQgEAoFAIBDUQYyqWwCBQCAQCASCqkIoOgKBQCAQCOosQtERCAQCgUBQZxGKjkAgEAgEgjqLUHQEAoFAIBDUWYSiIxAIBAKBoM4iFB2BQCAQCAR1FqHoCAQCgUAgqLMIRUcgEAgEAkGdRSg6ghpJcHAwCoWC27dv69zv4+NDr169NMYSExOZNGkSHh4eWFpa4uDgQLt27Xj99ddJTEws97m3bduGQqHA0dGR/Pz8B7mMUlFdY1Vx9epVFAoFq1evlscOHTpEcHAwmZmZVXbeqmb16tUoFAquXr1apcf26tVL6zNWmfTq1QuFQsHAgQO19qn+dl999ZXOY8vzGc3JyWHevHm0b98eW1tbbGxsaNmyJaNHj2bfvn0ac2NjY3nxxRdp0aIFFhYWODk54efnx5tvvsmdO3fKdT2ff/453t7eFBcXl+savvrqK62/RUFBAT/88AOPPfYYDg4O1KtXD3d3d4YPH86WLVu07o9qMzU1xdHRkccee4ypU6dy7ty5UmWdNm0aCoWCIUOG6NwfFxeHmZkZJ06cKNe1C2o2QtER1AmuX7+On58ff//9N9OmTePPP/9k5cqVjBkzhqNHj3LlypVyrxUaGgpAeno6W7durSKJ4bXXXiMyMrLK1tfFoUOH+Oyzz2q1ojN48GAiIyNp2LBhdYtSKezcuZM9e/YYdExZn9GioiL69+/P7NmzGTlyJJs2bWLz5s1MnTqVrKws9u/fL889efIknTp1IiYmhk8++YQdO3awbNkyBg8ezM6dO0lPTy9Tnps3bzJ//nw+//xzjIwq/lh58cUXmTJlCr1792bdunX8/vvvfPTRR5iYmLBz506t+VOmTCEyMpJ9+/axdu1ann76abZt20b79u1ZsGCBznMUFBSwbt06AHbs2MGNGze05nh4eDB27FimTp1a4WsR1CAkgaAG8umnn0qAlJqaqnN/27ZtJX9/f/n9J598IgHSlStXdM4vKioq13mTkpIkExMTqU+fPpKFhYXUr18/g2Uvi5ycnEpfUxfx8fESIK1atUoeW7BggQRI8fHxlXque/fuScXFxZW6ZlWwatWqcl+/v7+/xmessvH395c8PDykFi1aSJ06ddK4f6q/3YIFC7SOK89ndM+ePRIgrVy5Uue51b8P48aNk6ysrKQ7d+7onFuev+uMGTOkxo0ba6xb2jVIkvZn8cqVKxIgffLJJ2XKXNra9+7dkwYOHCgB0p9//qm1f9OmTRIgDR48WAKk2bNn6zzfsWPHJEA6ePCg3usW1A6ERUdQJ0hLS8PIyAgXFxed+8v7K/N///sfhYWFTJ06lWeffZbdu3eTkJCgNS8zM5NXX30VBwcHrK2tGTx4MFeuXEGhUBAcHCzPU7mnTpw4wciRI7G3t6dly5Ya+0ry008/0a1bN6ytrbG2tqZDhw7yL3iAZs2aMX78eK3jynK1BAcH89577wHQvHlz2ewfEREBoCW7vvOp3D+7du3ilVdewdnZmXr16skulLCwMLp164aVlRXW1tYMGDCAkydP6pULlK4IExMT5s6dq7Xvn3/+QaFQsGnTJo3zl3Q/hYeH07dvX2xtbalXrx49evRg9+7dpZ4XQJIk5s+fj7u7OxYWFvj5+fHXX39pzSsuLmbWrFm0adMGS0tL6tevj6+vL998802Z59CHqakps2fP5vjx44SFhZXrmPJ8RtPS0gD0Wr3Uvw9paWnY2tpibW2tc25Z7tX79+8TGhrK888//0DWHENkLg1LS0tCQ0MxNTXVadUJDQ3FzMyMVatW4ebmxqpVq5B09Lbu1KkTXl5eLFu2zICrENREhKIjqBN069aN4uJinn32WXbu3FnuuIKSrFy5koYNG/LUU0/xyiuvUFxcrBHjAsoH3tChQ/npp594//332bJlC127dtUZa6Hi2WefpVWrVmzatKnU/zg/+eQTxo4dS6NGjVi9ejVbtmzhpZde0qlsGcprr73GlClTAPj111+JjIwkMjISPz+/Cq33yiuvYGpqytq1a9m8eTOmpqbMmTOHMWPG4O3tzc8//8zatWvJzs7mySefJCYmRu9azZo1Y9iwYSxbtoyioiKNfd9//z2NGjXimWee0Xv8unXr6N+/P7a2tvzvf//j559/xsHBgQEDBpSp7Hz22We8//779OvXj61btzJx4kRef/11Lly4oDFv/vz5BAcHM2bMGLZv305YWBivvvrqA7sBAwMD6dSpEx999BEFBQVlzi/PZ7Rz586Ympry9ttvs379epKSkvSu161bN5KSkhg7diz79u0jNzfXIPmPHDlCWloavXv3Nui4knh5eVG/fn0+++wzli9fXqEYLBWNGjWiU6dOHDp0iMLCQnn8+vXr7Nq1i+HDh+Ps7MxLL73EpUuX+Oeff3Su06tXL/766y+dipCgFlHNFiWBQCeGuq6Ki4ul//u//5OMjIwkQFIoFJKXl5c0derUcrtp/vnnHwmQPvjgA3nN5s2bS+7u7hrm++3bt0uAFBISonH83LlzJUD69NNPta5DlzletU/FlStXJGNjY2ns2LGlyunu7i699NJLWuMlXS2Guq5Kyq7vfCr3z7hx4zTmXbt2TTIxMZGmTJmiMZ6dnS25urpKo0ePLvW69u7dKwHSli1b5LEbN25IJiYm0meffaZ1ftU15OTkSA4ODtLQoUM11isqKpLat28vdenSRe+xGRkZkoWFhfTMM89oHHvw4EEJ0LifQ4YMkTp06FDqNRiCv7+/1LZtW0mSJCk8PFwCpO+++06SJP2umfJ+RiVJkkJDQyVra2sJkACpYcOG0rhx46R//vlHY15eXp709NNPy/OMjY2ljh07SjNnzpRSUlLKvI558+ZJgJScnKwxbqjrSpKU3y0nJydZFkdHR2nUqFHStm3bDFpbkiQpMDBQAqRbt27JY59//rkESDt27JAkSfmdUygU0osvvqhzjRUrVkiAFBsbW+o9ENRshEVHUCdQKBQsW7aMK1eusHTpUl5++WUKCgpYtGgRbdu21coy0YXKPfTKK6/Ia44fP56EhAQNq4BqrdGjR2scP2bMGL1rjxgxoszz//333xQVFTF58uQy59YESl7Tzp07KSwsZNy4cRQWFsqbhYUF/v7+sotMH7169aJ9+/YsWbJEHlu2bBkKhYIJEyboPe7QoUOkp6fz0ksvaZy3uLiYgQMHcvToUXJycnQeGxkZSV5eHmPHjtUY7969O+7u7hpjXbp04dSpU0yaNOmBrIa66Nu3L/379+fzzz8nOztb77zyfkZVc65fv85PP/3EW2+9hZubG+vWrcPf31/DpWNubs6WLVuIiYlh0aJFPPfcc6SmpjJ79my8vLy0LFsluXnzJgqFAicnp4pevsygQYO4du0aW7Zs4d1336Vt27Zs3bqVYcOG8eabbxq0llTCCiNJkuyu6tevH6B04fbq1YtffvlF599T5QrXFbAsqD0IRUdQIzExMQHQcmOoKCwsxNTUVGvc3d2diRMnEhoaysWLFwkLCyMvL0+OTdFHdnY2mzZtokuXLjg7O5OZmUlmZibPPPMMCoVCI0YmLS0NExMTHBwcNNZo0KCB3vXLkyGUmpoKQJMmTcqcWxMoeU23bt0C4LHHHsPU1FRjCwsL01sqQJ233nqL3bt3c+HCBQoKClixYgUjR47E1dVV7zGq844cOVLrvPPmzUOSJL2ZQ6q4EF3rlxwLCgriq6++4vDhwzz11FM4OjrSt29fjh07VuZ1lYd58+Zx+/ZtvenYhnxGVdjZ2TFmzBi
|
2023-01-12 21:45:54 +02:00
|
|
|
"text/plain": [
|
|
|
|
|
"<Figure size 640x480 with 1 Axes>"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"output_type": "display_data"
|
|
|
|
|
}
|
|
|
|
|
],
|
2023-01-12 20:25:12 +02:00
|
|
|
"source": [
|
2023-01-13 20:54:15 +02:00
|
|
|
"plt.plot(corn_year,corn_bu_per_acre,\"o\",markerfacecolor=\"none\",label=\"Corn (BU/ACRE)\")\n",
|
2023-01-12 20:25:12 +02:00
|
|
|
"plt.plot(wheat_year,wheat_bu_per_acre,\"x\",label=\"Wheat (BU/ACRE)\")\n",
|
|
|
|
|
"plt.plot(potatoes_year,potatoes_cwt_per_acre,\"+\",label=\"Potatoes (CWT / ACRE)\")\n",
|
2023-01-13 20:54:15 +02:00
|
|
|
"plt.plot(soybeans_year,soybeans_bu_per_acre,\"*\",markerfacecolor=\"none\",label=\"Soybeans (BU/ACRE)\")\n",
|
2023-01-12 16:13:39 +02:00
|
|
|
"plt.xlabel(\"time (calendar year)\")\n",
|
2023-01-12 20:25:12 +02:00
|
|
|
"plt.ylabel(\"Average per acre production\")\n",
|
2023-01-12 16:13:39 +02:00
|
|
|
"plt.title(\"US Agriculture yields, NASS (USDA)\")\n",
|
|
|
|
|
"plt.legend()\n",
|
2023-01-12 20:25:12 +02:00
|
|
|
"plt.yscale(\"log\")\n",
|
|
|
|
|
"#plt.show()\n",
|
|
|
|
|
"#plt.savefig(\"bear-quadratic.pdf\")"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "markdown",
|
|
|
|
|
"id": "c68984fa",
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"source": [
|
|
|
|
|
"## Sugarbeets"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2023-01-12 21:45:54 +02:00
|
|
|
"execution_count": 15,
|
2023-01-12 20:25:12 +02:00
|
|
|
"id": "6e74b0b8",
|
|
|
|
|
"metadata": {},
|
2023-01-12 21:45:54 +02:00
|
|
|
"outputs": [
|
|
|
|
|
{
|
|
|
|
|
"data": {
|
|
|
|
|
"text/plain": [
|
|
|
|
|
"array([[ nan, 2022. , nan, nan, nan, nan, nan, nan,\n",
|
|
|
|
|
" nan, nan, nan, nan, nan, 0. , nan, nan,\n",
|
|
|
|
|
" nan, nan, 29.1, nan],\n",
|
|
|
|
|
" [ nan, 2021. , nan, nan, nan, nan, nan, nan,\n",
|
|
|
|
|
" nan, nan, nan, nan, nan, 0. , nan, nan,\n",
|
|
|
|
|
" nan, nan, 33.2, nan],\n",
|
|
|
|
|
" [ nan, 2020. , nan, nan, nan, nan, nan, nan,\n",
|
|
|
|
|
" nan, nan, nan, nan, nan, 0. , nan, nan,\n",
|
|
|
|
|
" nan, nan, 29.4, nan]])"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
"execution_count": 15,
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"output_type": "execute_result"
|
|
|
|
|
}
|
|
|
|
|
],
|
2023-01-12 20:25:12 +02:00
|
|
|
"source": [
|
|
|
|
|
"csv=np.genfromtxt ('./data/sugarbeets.csv',\n",
|
|
|
|
|
" skip_header=1,\n",
|
|
|
|
|
" delimiter=\",\")\n",
|
|
|
|
|
"csv[0:3]\n",
|
|
|
|
|
"#csv[0:3].astype(float)"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2023-01-12 21:45:54 +02:00
|
|
|
"execution_count": 16,
|
2023-01-12 20:25:12 +02:00
|
|
|
"id": "ab73b355",
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [],
|
|
|
|
|
"source": [
|
|
|
|
|
"sugarbeets_year=csv[:,1]\n",
|
|
|
|
|
"sugarbeets_tons_per_acre=csv[:,-2]"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2023-01-12 21:45:54 +02:00
|
|
|
"execution_count": 17,
|
2023-01-12 20:25:12 +02:00
|
|
|
"id": "cbaa454a",
|
|
|
|
|
"metadata": {},
|
2023-01-12 21:45:54 +02:00
|
|
|
"outputs": [
|
|
|
|
|
{
|
|
|
|
|
"data": {
|
2023-01-13 20:54:15 +02:00
|
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAjoAAAHFCAYAAAD7ZFORAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOydd1RURxuHf0vvHQQEAUWaCIjY0AgoirF3JNiSKIk1dkUTRT9L1KhJLBgMYuzETkyCCoKigIAIRkFAREAF6UWawM73x2ZvdtldYJUmznPOPezOzJ07s4X77ltZhBACCoVCoVAolE6IRHsvgEKhUCgUCqW1oIIOhUKhUCiUTgsVdCgUCoVCoXRaqKBDoVAoFAql00IFHQqFQqFQKJ0WKuhQKBQKhULptFBBh0KhUCgUSqeFCjoUCoVCoVA6LVTQoVAoFAqF0mmhgg6lQ+Lj4wMWi4WCggKh/dbW1nB2duZry87OxsKFC2FmZgZ5eXloaGigd+/emD9/PrKzs5t97aCgILBYLGhqaqKmpuZ9ttEo3D22Fs+fPweLxcKxY8eYtsjISPj4+KCkpKTVrtvaHDt2DCwWC8+fP2/Vc52dnQU+Yy2Js7MzWCwWRo0aJdDHfe9++OEHoec25zNaUVGBnTt3wtbWFioqKlBWVkaPHj0wffp03Lp1i29scnIyZs2ahe7du0NOTg5aWlqwt7fH4sWLUVZW1qz9bNmyBVZWVmCz2c3aww8//CDwXtTW1uKXX35Bv379oKGhAQUFBRgZGWHChAm4dOmSwOvDPaSlpaGpqYl+/fph+fLlePz4caNrXbFiBVgsFsaOHSu0PzU1FTIyMoiPj2/W3ikdGyroUDoFL168gL29PW7cuIEVK1bgr7/+wtGjR+Hh4YHY2Fg8e/as2XP5+/sDAIqKinD58uVWWjEwb948REVFtdr8woiMjMTmzZs/aEFnzJgxiIqKgp6eXnsvpUW4du0abt68KdY5TX1G6+vrMXLkSGzbtg1Tp07FuXPncP78eSxfvhylpaWIiIhgxj548AB9+/ZFUlISNm7ciODgYBw+fBhjxozBtWvXUFRU1OR6Xr16hV27dmHLli2QkHj328qsWbOwZMkSuLi44OTJk/jjjz/w7bffQkpKCteuXRMYv2TJEkRFReHWrVs4ceIEJk6ciKCgINja2mL37t1Cr1FbW4uTJ08CAIKDg/Hy5UuBMWZmZvD09MTy5cvfeS+UDgShUDogmzZtIgBIfn6+0P5evXoRJycn5vnGjRsJAPLs2TOh4+vr65t13ZycHCIlJUWGDRtG5OTkyIgRI8Ree1NUVFS0+JzCyMjIIABIQEAA07Z7924CgGRkZLTotSorKwmbzW7ROVuDgICAZu/fycmJ7zPW0jg5OREzMzPSvXt30rdvX77Xj/ve7d69W+C85nxGb968SQCQo0ePCr027/dh9uzZRFFRkZSVlQkd25z3dc2aNaRr16588za2B0IEP4vPnj0jAMjGjRubXHNjc1dWVpJRo0YRAOSvv/4S6D937hwBQMaMGUMAkG3btgm9XlxcHAFA7t69K3LflA8DqtGhdAoKCwshISEBHR0dof3N/ZX522+/oa6uDsuXL8fkyZMRGhqKzMxMgXElJSX48ssvoaGhASUlJYwZMwbPnj0Di8WCj48PM45rnoqPj8fUqVOhrq6OHj168PU15PTp0xg0aBCUlJSgpKQEOzs75hc8ABgbG2Pu3LkC5zVlavHx8cHq1asBACYmJozaPzw8HAAE1i7qelzzz/Xr1/HFF19AW1sbCgoKjAklMDAQgwYNgqKiIpSUlODm5oYHDx6IXBfAMUVISUlhx44dAn23b98Gi8XCuXPn+K7f0PwUEhKC4cOHQ0VFBQoKChg8eDBCQ0MbvS4AEEKwa9cuGBkZQU5ODvb29vj7778FxrHZbGzduhXm5uaQl5eHmpoabGxs8NNPPzV5DVFIS0tj27ZtuH//PgIDA5t1TnM+o4WFhQAgUuvF+30oLCyEiooKlJSUhI5tyrz69u1b+Pv747PPPnsvbY44a24MeXl5+Pv7Q1paWqhWx9/fHzIyMggICIChoSECAgJAhNS27tu3LywtLXH48GExdkHpiFBBh9IpGDRoENhsNiZPnoxr164126+gIUePHoWenh4+/fRTfPHFF2Cz2Xw+LgDnhjdu3DicPn0aa9euxaVLlzBgwAChvhZcJk+eDFNTU5w7d67Rf5wbN26Ep6cn9PX1cezYMVy6dAlz5swRKmyJy7x587BkyRIAwMWLFxEVFYWoqCjY29u/03xffPEFpKWlceLECZw/fx7S0tLYvn07PDw8YGVlhd9//x0nTpxAeXk5PvnkEyQlJYmcy9jYGOPHj8fhw4dRX1/P13fgwAHo6+tj0qRJIs8/efIkRo4cCRUVFfz222/4/fffoaGhATc3tyaFnc2bN2Pt2rUYMWIELl++jAULFmD+/PlISUnhG7dr1y74+PjAw8MDf/75JwIDA/Hll1++txnQ3d0dffv2xbfffova2tomxzfnM+rg4ABpaWl88803OHXqFHJyckTON2jQIOTk5MDT0xO3bt1CVVWVWOu/d+8eCgsL4eLiItZ5DbG0tISamho2b94MPz+/d/LB4qKvr4++ffsiMjISdXV1TPuLFy9w/fp1TJgwAdra2pgzZw6ePn2K27dvC53H2dkZf//9t1BBiPIB0c4aJQpFKOKarthsNvnqq6+IhIQEAUBYLBaxtLQky5cvb7aZ5vbt2wQAWbduHTOniYkJMTIy4lPf//nnnwQA8fX15Tt/x44dBADZtGmTwD6EqeO5fVyePXtGJCUliaenZ6PrNDIyInPmzBFob2hqEdd01XDtoq7HNf/Mnj2bb1xWVhaRkpIiS5Ys4WsvLy8nurq6ZPr06Y3uKywsjAAgly5dYtpevnxJpKSkyObNmwWuz91DRUUF0dDQIOPGjeObr76+ntja2pL+/fuLPLe4uJjIycmRSZMm8Z179+5dAoDv9Rw7diyxs7NrdA/i4OTkRHr16kUIISQkJIQAIPv37yeEiDbNNPczSggh/v7+RElJiQAgAIienh6ZPXs2uX37Nt+46upqMnHiRGacpKQk6dOnD9mwYQPJy8trch87d+4kAEhubi5fu7imK0I43y0tLS1mLZqammTatGkkKChIrLkJIcTd3Z0AIK9fv2batmzZQgCQ4OBgQgjnO8discisWbOEznHkyBECgCQnJzf6GlA6NlSjQ+kUsFgsHD58GM+ePcOhQ4fw+eefo7a2Fvv27UOvXr0EokyEwTUPffHFF8ycc+fORWZmJp9WgDvX9OnT+c738PAQOfeUKVOavP6NGzdQX1+PRYsWNTm2I9BwT9euXUNdXR1mz56Nuro65pCTk4OTkxNjIhOFs7MzbG1tcfDgQabt8OHDYLFY8PLyEnleZGQkioqKMGfOHL7rstlsjBo1CrGxsaioqBB6blRUFKqrq+Hp6cnX7ujoCCMjI762/v37IzExEQsXLnwvraEwhg8fjpEjR2LLli0oLy8XOa65n1HumBcvXuD06dNYunQpDA0NcfLkSTg5OfGZdGRlZXHp0iUkJSVh3759mDFjBvLz87Ft2zZYWloKaLYa8urVK7BYLGhpab3r9hlGjx6NrKwsXLp0CatWrUKvXr1w+fJljB8/HosXLxZrLtJAC0MIYcxVI0aMAMAx4To7O+PChQtC30+uKVyYwzLlw4EKOpQOiZSUFAAImDG41NXVQVpaWqDdyMgICxYsgL+/P9LS0hAYGIjq6mrGN0UU5eXlOHfuHPr37w9tbW2UlJSgpKQEkyZNAovF4vORKSwshJSUFDQ0NPjm6NKli8j5mxMhlJ+fDwAwMDBocmxHoOGeXr9+DQDo168fpKWl+Y7AwECRqQJ4Wbp0KUJDQ5GSkoLa2locOXIEU6dOha6urshzuNedOnWqwHV37twJQojIyCGuX4iw+Ru2eXt744cffkB0dDQ+/fRTaGpqYvjw4YiLi2tyX81h586dKCgoEBmOLc5nlIuqqio8PDzw008/4d69e3j48CG6dOmCDRs2CJjcLC0
|
2023-01-12 21:45:54 +02:00
|
|
|
"text/plain": [
|
|
|
|
|
"<Figure size 640x480 with 1 Axes>"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"output_type": "display_data"
|
|
|
|
|
}
|
|
|
|
|
],
|
2023-01-12 20:25:12 +02:00
|
|
|
"source": [
|
2023-01-13 20:54:15 +02:00
|
|
|
"plt.plot(corn_year,corn_bu_per_acre,\"o\",markerfacecolor=\"none\",label=\"Corn (BU/ACRE)\")\n",
|
2023-01-12 20:25:12 +02:00
|
|
|
"plt.plot(wheat_year,wheat_bu_per_acre,\"x\",label=\"Wheat (BU/ACRE)\")\n",
|
|
|
|
|
"plt.plot(potatoes_year,potatoes_cwt_per_acre,\"+\",label=\"Potatoes (CWT / ACRE)\")\n",
|
2023-01-13 20:54:15 +02:00
|
|
|
"plt.plot(soybeans_year,soybeans_bu_per_acre,\"*\",markerfacecolor=\"none\",label=\"Soybeans (BU/ACRE)\")\n",
|
2023-01-12 20:25:12 +02:00
|
|
|
"plt.plot(sugarbeets_year,sugarbeets_tons_per_acre,\".\",label=\"Sugarbeets (TONS/ACRE)\")\n",
|
|
|
|
|
"plt.xlabel(\"time (calendar year)\")\n",
|
|
|
|
|
"plt.ylabel(\"Average per acre production\")\n",
|
|
|
|
|
"plt.title(\"US Agriculture yields, NASS (USDA)\")\n",
|
|
|
|
|
"plt.legend()\n",
|
|
|
|
|
"plt.yscale(\"log\")\n",
|
|
|
|
|
"#plt.show()\n",
|
|
|
|
|
"#plt.savefig(\"bear-quadratic.pdf\")"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "markdown",
|
|
|
|
|
"id": "884a7b05",
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"source": [
|
|
|
|
|
"## Sunflowers"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2023-01-12 21:45:54 +02:00
|
|
|
"execution_count": 18,
|
2023-01-12 20:25:12 +02:00
|
|
|
"id": "846722a4",
|
|
|
|
|
"metadata": {},
|
2023-01-12 21:45:54 +02:00
|
|
|
"outputs": [
|
|
|
|
|
{
|
|
|
|
|
"data": {
|
|
|
|
|
"text/plain": [
|
|
|
|
|
"array([[ nan, 2022., nan, nan, nan, nan, nan, nan, nan,\n",
|
|
|
|
|
" nan, nan, nan, nan, 0., nan, nan, nan, nan,\n",
|
|
|
|
|
" 1782., nan],\n",
|
|
|
|
|
" [ nan, 2021., nan, nan, nan, nan, nan, nan, nan,\n",
|
|
|
|
|
" nan, nan, nan, nan, 0., nan, nan, nan, nan,\n",
|
|
|
|
|
" 1529., nan],\n",
|
|
|
|
|
" [ nan, 2020., nan, nan, nan, nan, nan, nan, nan,\n",
|
|
|
|
|
" nan, nan, nan, nan, 0., nan, nan, nan, nan,\n",
|
|
|
|
|
" 1790., nan]])"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
"execution_count": 18,
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"output_type": "execute_result"
|
|
|
|
|
}
|
|
|
|
|
],
|
2023-01-12 20:25:12 +02:00
|
|
|
"source": [
|
|
|
|
|
"csv=np.genfromtxt ('./data/sunflowers.csv',\n",
|
|
|
|
|
" skip_header=1,\n",
|
|
|
|
|
" delimiter=\",\")\n",
|
|
|
|
|
"csv[0:3]\n",
|
|
|
|
|
"#csv[0:3].astype(float)"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2023-01-12 21:45:54 +02:00
|
|
|
"execution_count": 19,
|
2023-01-12 20:25:12 +02:00
|
|
|
"id": "c79d0d11",
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [],
|
|
|
|
|
"source": [
|
|
|
|
|
"sunflowers_year=csv[:,1]\n",
|
2023-01-12 21:45:54 +02:00
|
|
|
"sunflowers_lbs_per_acre=csv[:,-2]\n",
|
|
|
|
|
"\n",
|
|
|
|
|
"#grams_per_lbs\n",
|
|
|
|
|
"sunflowers_kcal_per_gram=584.0/100.0\n",
|
|
|
|
|
"sunflowers_kcal_per_acre=sunflowers_lbs_per_acre*grams_per_lbs*sunflowers_kcal_per_gram"
|
2023-01-12 20:25:12 +02:00
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2023-01-13 20:54:15 +02:00
|
|
|
"execution_count": 34,
|
2023-01-12 20:25:12 +02:00
|
|
|
"id": "569111b9",
|
|
|
|
|
"metadata": {},
|
2023-01-12 21:45:54 +02:00
|
|
|
"outputs": [
|
|
|
|
|
{
|
|
|
|
|
"data": {
|
2023-01-13 20:54:15 +02:00
|
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAjsAAAHFCAYAAAAUpjivAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAA9hAAAPYQGoP6dpAADMGklEQVR4nOzdd3hUxdfA8e+mJ6SQXiBC6CShB5CioRdBQFAUUAFRFBAEFARBRaQIqKCo+FMpiiC+KiAq0qv0jpDQQ08IhDRISNv7/nHdzW6yS7KQnvN5nn12M/fuvTObhBxmzsxoFEVREEIIIYQoo6yKuwJCCCGEEIVJgh0hhBBClGkS7AghhBCiTJNgRwghhBBlmgQ7QgghhCjTJNgRQgghRJkmwY4QQgghyjQJdoQQQghRpkmwI4QQQogyTYIdUaJNmTIFjUbDrVu3TB4PDQ2lTZs2RmVXrlxh+PDh1KpVC0dHRzw8PKhXrx6vvPIKV65cyfe916xZg0ajwdPTk7S0tIdpxn3p2lhYLl68iEajYcmSJfqy3bt3M2XKFBISEgrtvoVtyZIlaDQaLl68WKjvbdOmTa6fsYLUpk0bNBoNXbp0yXVM9737+OOPTb43Pz+jd+/eZdasWTRo0ABXV1dcXFyoXr06ffv2Zfv27UbnRkZG8sILL1CtWjUcHBzw8vKicePGvP766yQlJeWrPVOnTiU4OBitVpuvNnz88ce5vhcZGRn873//o2nTpnh4eODk5ESVKlXo2bMnq1atyvX56B62trZ4enrStGlTxowZw8mTJ+9b17Fjx6LRaOjevbvJ42fOnMHOzo7Dhw/nq+2i5JJgR5QpV69epXHjxmzcuJGxY8eydu1aFi1aRL9+/Thw4AAXLlzI97UWLlwIwO3bt1m9enUh1Rhefvll9uzZU2jXN2X37t188MEHpTrY6datG3v27MHf37+4q1Ig1q9fz5YtWyx6T14/o1lZWXTq1Inp06fz9NNP88svv/Drr78yZswYEhMT2blzp/7cI0eO0KRJEyIiInjvvfdYt24dX3/9Nd26dWP9+vXcvn07z/pcv36d2bNnM3XqVKysHvzPywsvvMDIkSNp27YtP/74I3/88QeTJ0/GxsaG9evX5zp/5MiR7Nmzh+3bt7N06VJ69erFmjVraNCgAXPmzDF5j4yMDH788UcA1q1bx7Vr13KdU6tWLQYMGMCYMWMeuC2ihFCEKMHef/99BVBu3rxp8nhISIgSHh6u//q9995TAOXChQsmz8/KysrXfaOjoxUbGxulXbt2ioODg9KxY0eL656Xu3fvFvg1TYmKilIAZfHixfqyOXPmKIASFRVVoPdKSUlRtFptgV6zMCxevDjf7Q8PDzf6GSto4eHhSq1atZRq1aopTZo0Mfr8dN+7OXPm5Hpffn5Gt2zZogDKokWLTN7b8PfhxRdfVCpUqKAkJSWZPDc/39fx48crlSpVMrru/dqgKLl/Fi9cuKAAynvvvZdnne937ZSUFKVLly4KoKxduzbX8V9++UUBlG7duimAMn36dJP3O3jwoAIou3btMttuUfJJz44oU+Li4rCyssLHx8fk8fz+b/P7778nMzOTMWPG0Lt3bzZv3sylS5dynZeQkMCQIUPw8PDA2dmZbt26ceHCBTQaDVOmTNGfpxuqOnz4ME8//TTu7u5Ur17d6FhOy5cvp0WLFjg7O+Ps7EzDhg31/5MHqFq1KoMGDcr1vryGXaZMmcK4ceMACAoK0g8BbNu2DSBX3c3dTzcUtGHDBl566SW8vb1xcnLSD6f8/PPPtGjRggoVKuDs7Eznzp05cuSI2XqBOixhY2PDzJkzcx3bsWMHGo2GX375xej+OYeiNm3aRPv27XF1dcXJyYlWrVqxefPm+94XQFEUZs+eTZUqVXBwcKBx48b8/fffuc7TarVMmzaN2rVr4+joSMWKFalfvz6fffZZnvcwx9bWlunTp3Po0CF+/vnnfL0nPz+jcXFxAGZ7vwx/H+Li4nB1dcXZ2dnkuXkNtaanp7Nw4UL69+//UL06ltT5fhwdHVm4cCG2trYme3cWLlyInZ0dixcvJjAwkMWLF6OY2Be7SZMm1K1bl6+//tqCVoiSRoIdUaa0aNECrVZL7969Wb9+fb7zDHJatGgR/v7+dO3alZdeegmtVmuU8wLqH70nn3yS5cuX8/bbb7Nq1SqaN29uMvdCp3fv3tSoUYNffvnlvv94vvfeewwYMICAgACWLFnCqlWrGDhwoMmAy1Ivv/wyI0eOBGDlypXs2bOHPXv20Lhx4we63ksvvYStrS1Lly7l119/xdbWlhkzZtCvXz+Cg4P5v//7P5YuXUpycjKPPfYYERERZq9VtWpVevTowddff01WVpbRsS+++IKAgACeeuops+//8ccf6dSpE66urnz//ff83//9Hx4eHnTu3DnPgOeDDz7g7bffpmPHjqxevZphw4bxyiuvcPr0aaPzZs+ezZQpU+jXrx9//fUXP//8M0OGDHnoIcFnn32WJk2aMHnyZDIyMvI8Pz8/o2FhYdja2vLGG2+wbNkyoqOjzV6vRYsWREdHM2DAALZv305qaqpF9d+3bx9xcXG0bdvWovflVLduXSpWrMgHH3zAN99880A5WToBAQE0adKE3bt3k5mZqS+/evUqGzZsoGfPnnh7ezNw4EDOnTvHjh07TF6nTZs2/P333yaDIVFKFHPPkhD3ZekwllarVV599VXFyspKARSNRqPUrVtXGTNmTL6HbHbs2KEAyoQJE/TXDAoKUqpUqWLUlf/XX38pgLJgwQKj98+cOVMBlPfffz9XO0x1zeuO6Vy4cEGxtrZWBgwYcN96VqlSRRk4cGCu8pzDLpYOY+Wsu7n76YaCXnzxRaPzLl++rNjY2CgjR440Kk9OTlb8/PyUvn373rddW7duVQBl1apV+rJr164pNjY2ygcffJDr/ro23L17V/Hw8FCefPJJo+tlZWUpDRo0UJo1a2b2vfHx8YqDg4Py1FNPGb13165dCmD0eXbv3l1p2LDhfdtgifDwcCUkJERRFEXZtGmTAijz589XFMX8ME1+f0YVRVEWLlyoODs7K4ACKP7+/sqLL76o7Nixw+i8e/fuKb169dKfZ21trTRq1EiZNGmSEhsbm2c7Zs2apQBKTEyMUbmlw1iKov5ueXl56evi6empPPPMM8qaNWssuraiKMqzzz6rAMqNGzf0ZVOnTlUAZd26dYqiqL9zGo1GeeGFF0xe49tvv1UAJTIy8r6fgSi5pGdHlCkajYavv/6aCxcu8NVXXzF48GAyMjKYO3cuISEhuWafmKIbKnrppZf01xw0aBCXLl0y6h3QXatv375G7+/Xr5/Za/fp0yfP+2/cuJGsrCxGjBiR57klQc42rV+/nszMTF588UUyMzP1DwcHB8LDw/XDZea0adOGBg0a8OWXX+rLvv76azQaDUOHDjX7vt27d3P79m0GDhxodF+tVkuXLl04cOAAd+/eNfnePXv2cO/ePQYMGGBU3rJlS6pUqWJU1qxZM44dO8bw4cMfqvfQlPbt29OpUyemTp1KcnKy2fPy+zOqO+fq1assX76cUaNGERgYyI8//kh4eLjR8I69vT2rVq0iIiKCuXPn8txzz3Hz5k2mT59O3bp1c/Vw5XT9+nU0Gg1eXl4P2ny9J554gsuXL7Nq1SreeustQkJCWL16NT169OD111+36FpKjt4YRVH0Q1cdO3YE1OHcNm3a8Ntvv5n8fuqGxU0lMYvSQYIdUaLZ2NgA5BrS0MnMzMTW1jZXeZUqVRg2bBgLFy7k7Nmz/Pzzz9y7d0+fq2JOcnIyv/zyC82aNcPb25uEhAQSEhJ46qmn0Gg0RjkzcXFx2NjY4OHhYXQNX19fs9fPz8yhmzdvAlC5cuU8zy0Jcrbpxo0bADRt2hRbW1ujx88//2x2GQFDo0aNYvPmzZw+fZqMjAy+/fZbnn76afz8/My+R3ffp59+Otd9Z82ahaIoZmcU6fJETF0/Z9nEiRP5+OOP2bt3L127dsXT05P27dtz8ODBPNuVH7NmzeLWrVtmp2pb8jOq4+bmRr9
|
|
|
|
|
"text/plain": [
|
|
|
|
|
"<Figure size 640x480 with 1 Axes>"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"output_type": "display_data"
|
|
|
|
|
}
|
|
|
|
|
],
|
|
|
|
|
"source": [
|
|
|
|
|
"plt.plot(corn_year,corn_bu_per_acre,\"o\",markerfacecolor=\"none\",label=\"Corn (BU/ACRE)\")\n",
|
|
|
|
|
"#plt.plot(wheat_year,wheat_bu_per_acre,\"x\",label=\"Wheat (BU/ACRE)\")\n",
|
|
|
|
|
"plt.plot(potatoes_year,potatoes_cwt_per_acre,\"+\",label=\"Potatoes (CWT/ACRE)\")\n",
|
|
|
|
|
"#plt.plot(soybeans_year,soybeans_bu_per_acre,\"*\",markerfacecolor=\"none\",label=\"Soybeans (BU/ACRE)\")\n",
|
|
|
|
|
"#plt.plot(sunflowers_year,sunflowers_lbs_per_acre,\"d\",markerfacecolor=\"none\",label=\"Sunflowers (LB/ACRE)\")\n",
|
|
|
|
|
"#plt.plot(sugarbeets_year,sugarbeets_tons_per_acre,\".\",label=\"Sugarbeets (TONS/ACRE)\")\n",
|
|
|
|
|
"plt.xlabel(\"time (calendar year)\")\n",
|
|
|
|
|
"plt.ylabel(\"Average per acre production\")\n",
|
|
|
|
|
"plt.title(\"US Agriculture yields, NASS (USDA)\")\n",
|
|
|
|
|
"plt.legend()\n",
|
|
|
|
|
"#plt.yscale(\"log\")\n",
|
|
|
|
|
"plt.grid()\n",
|
|
|
|
|
"#plt.show()\n",
|
|
|
|
|
"plt.savefig(\"corn_potatoes_raw_production_per_acre.pdf\")"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
|
|
|
|
"execution_count": 31,
|
|
|
|
|
"id": "536bea26",
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [
|
|
|
|
|
{
|
|
|
|
|
"data": {
|
|
|
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAjoAAAHFCAYAAAD7ZFORAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOydd1hT1//H3zfMMAIYQECRIQICIkvcAlXBra0Td63ar6PO1lGtWuuoWrXLUSvugXWW2rplKSg4QAWUraggyJZNcn5/5JfbXJIgWBDF83qe+5ice+4ZN4n3w2cyhBACCoVCoVAolGYIr6kXQKFQKBQKhdJYUEGHQqFQKBRKs4UKOhQKhUKhUJotVNChUCgUCoXSbKGCDoVCoVAolGYLFXQoFAqFQqE0W6igQ6FQKBQKpdlCBR0KhUKhUCjNFiroUCgUCoVCabZQQYfyTrJq1SowDIOXL18qPO/k5ARvb29OW0ZGBmbOnAlbW1vw+Xy0aNECHTp0wLRp05CRkVHnuYOCgsAwDIRCISoqKv7LNmpFusfGIj09HQzDYN++fWxbREQEVq1ahYKCgkabt7HZt28fGIZBenp6o17r7e0t9x1rSLy9vcEwDPr16yd3TvrZ/fDDDwqvrct3tKSkBBs2bEDHjh0hEAigq6uLtm3bYtSoUQgNDeX0TUhIwIQJE2BtbQ1NTU0YGhrCzc0Ns2fPRlFRUZ32s3r1ajg4OEAsFtdpDz/88IPcZ1FVVYXffvsNnTp1QosWLaClpQULCwsMHToUp0+flrs/0kNNTQ1CoRCdOnXC/PnzERcXV+taFyxYAIZhMGjQIIXnExMToa6ujjt37tRp75R3GyroUJoFT58+hZubGy5duoQFCxbgn3/+wZ49e+Dv74/o6GikpqbWeayAgAAAQF5eHs6cOdNIKwamTp2KyMjIRhtfEREREfj222/fa0Fn4MCBiIyMhKmpaVMvpUG4cOECrl69Wq9rXvcdFYlE8PX1xdq1azFixAgcP34cJ06cwPz581FYWIjw8HC27927d+Hu7o74+HisWLEC58+fx86dOzFw4EBcuHABeXl5r13P8+fPsXHjRqxevRo83ps/ViZMmIAvvvgCPj4+OHToEP766y8sX74cqqqquHDhglz/L774ApGRkQgNDcXBgwcxbNgwBAUFoWPHjti0aZPCOaqqqnDo0CEAwPnz5/Hs2TO5Pra2thg3bhzmz5//xnuhvEMQCuUdZOXKlQQAycnJUXje0dGReHl5se9XrFhBAJDU1FSF/UUiUZ3mzczMJKqqquSjjz4impqapG/fvvVe++soKSlp8DEVkZaWRgCQvXv3sm2bNm0iAEhaWlqDzlVaWkrEYnGDjtkY7N27t8779/Ly4nzHGhovLy9ia2tLrK2tibu7O+f+ST+7TZs2yV1Xl+/o1atXCQCyZ88ehXPL/h4mTpxItLW1SVFRkcK+dflcFy1aRFq1asUZt7Y9ECL/XUxNTSUAyIoVK1675trGLi0tJf369SMAyD///CN3/vjx4wQAGThwIAFA1q5dq3C+W7duEQDk+vXrSvdNeT+gGh1KsyA3Nxc8Hg/GxsYKz9f1r8z9+/ejuroa8+fPxyeffIIrV67g8ePHcv0KCgrw2WefoUWLFtDR0cHAgQORmpoKhmGwatUqtp/UPHXnzh2MGDECBgYGaNu2LedcTY4cOYKuXbtCR0cHOjo6cHFxYf+CBwBLS0tMnjxZ7rrXmVpWrVqFr776CgBgZWXFqv1DQkIAQG7tyuaTmn8uXryIKVOmwMjICFpaWqwJ5dixY+jatSu0tbWho6MDPz8/3L17V+m6AIkpQlVVFevXr5c7FxYWBoZhcPz4cc78Nc1Ply9fRu/evSEQCKClpYXu3bvjypUrtc4LAIQQbNy4ERYWFtDU1ISbmxvOnTsn108sFmPNmjWws7MDn8+Hvr4+nJ2d8dNPP712DmWoqalh7dq1uH37No4dO1ana+ryHc3NzQUApVov2d9Dbm4uBAIBdHR0FPZ9nXm1srISAQEBGDt27H/S5tRnzbXB5/MREBAANTU1hVqdgIAAqKurY+/evTA3N8fevXtBFNS2dnd3R/v27bFz58567ILyLkIFHUqzoGvXrhCLxfjkk09w4cKFOvsV1GTPnj0wNTVF//79MWXKFIjFYo6PCyB54A0ePBhHjhzB4sWLcfr0aXTu3Fmhr4WUTz75BDY2Njh+/Hit/3GuWLEC48aNg5mZGfbt24fTp09j0qRJCoWt+jJ16lR88cUXAIBTp04hMjISkZGRcHNze6PxpkyZAjU1NRw8eBAnTpyAmpoa1q1bB39/fzg4OOCPP/7AwYMHUVxcjJ49eyI+Pl7pWJaWlhgyZAh27twJkUjEOffrr7/CzMwMH3/8sdLrDx06BF9fXwgEAuzfvx9//PEHWrRoAT8/v9cKO99++y0WL16Mvn374syZM5gxYwamTZuGR48ecfpt3LgRq1atgr+/P/7++28cO3YMn3322X82A44ePRru7u5Yvnw5qqqqXtu/Lt9RDw8PqKmpYe7cuTh8+DAyMzOVjte1a1dkZmZi3LhxCA0NRVlZWb3Wf/PmTeTm5sLHx6de19Wkffv20NfXx7fffotdu3a9kQ+WFDMzM7i7uyMiIgLV1dVs+9OnT3Hx4kUMHToURkZGmDRpEpKTkxEWFqZwHG9vb5w7d06hIER5j2hijRKFopD6mq7EYjH5/PPPCY/HIwAIwzCkffv2ZP78+XU204SFhREAZMmSJeyYVlZWxMLCgqO+//vvvwkAsmPHDs7169evJwDIypUr5fahSB0vPSclNTWVqKiokHHjxtW6TgsLCzJp0iS59pqmlvqarmquXdl8UvPPxIkTOf2ePHlCVFVVyRdffMFpLy4uJiYmJmTUqFG17is4OJgAIKdPn2bbnj17RlRVVcm3334rN790DyUlJaRFixZk8ODBnPFEIhHp2LEj8fT0VHptfn4+0dTUJB9//DHn2uvXrxMAnPs5aNAg4uLiUuse6oOXlxdxdHQkhBBy+fJlAoD88ssvhBDlppm6fkcJISQgIIDo6OgQAAQAMTU1JRMnTiRhYWGcfuXl5WTYsGFsPxUVFeLq6kqWLVtGsrOzX7uPDRs2EAAkKyuL015f0xUhkt+WoaEhuxahUEhGjhxJgoKC6jU2IYSMHj2aACAvXrxg21avXk0AkPPnzxNCJL85hmHIhAkTFI7x+++/EwAkISGh1ntAebehGh1Ks4BhGOzcuROpqanYvn07Pv30U1RVVWHr1q1wdHSUizJRhNQ8NGXKFHbMyZMn4/HjxxytgHSsUaNGca739/dXOvbw4cNfO/+lS5cgEokwa9as1/Z9F6i5pwsXLqC6uhoTJ05EdXU1e2hqasLLy4s1kSnD29sbHTt2xLZt29i2nTt3gmEYTJ8+Xel1ERERyMvLw6RJkzjzisVi9OvXD9HR0SgpKVF4bWRkJMrLyzFu3DhOe7du3WBhYcFp8/T0RGxsLGbOnPmftIaK6N27N3x9fbF69WoUFxcr7VfX76i0z9OnT3HkyBHMmTMH5ubmOHToELy8vDgmHQ0NDZw+fRrx8fHYunUrxowZg5ycHKxduxbt27eX02zV5Pnz52AYBoaGhm+6fZYBAwbgyZMnOH36NL788ks4OjrizJkzGDJkCGbPnl2vsUgNLQwhhDVX9e3bF4DEhOvt7Y2TJ08q/DylpnBFDsuU9wcq6FDeSVRVVQFAzowhpbq6GmpqanLtFhYWmDFjBgICApCUlIRjx46hvLyc9U1RRnFxMY4fPw5PT08YGRmhoKAABQUF+Pjjj8EwDMdHJjc3F6qqqmjRogVnjJYtWyodvy4RQjk5OQCA1q1bv7bvu0DNPb148QIA0KlTJ6ipqXGOY8eOKU0VIMucOXNw5coVPHr0CFVVVfj9998xYsQImJiYKL1GOu+IESPk5t2wYQMIIUojh6R+IYrGr9m2dOlS/PDDD7hx4wb69+8PoVCI3r1749atW6/dV13YsGEDXr58qTQcuz7fUSl6enrw9/fHTz/9hJs3b+LevXto2bIlli1bJmdya9+
|
2023-01-12 21:45:54 +02:00
|
|
|
"text/plain": [
|
|
|
|
|
"<Figure size 640x480 with 1 Axes>"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"output_type": "display_data"
|
|
|
|
|
}
|
|
|
|
|
],
|
2023-01-12 20:25:12 +02:00
|
|
|
"source": [
|
2023-01-13 20:54:15 +02:00
|
|
|
"plt.plot(corn_year,corn_bu_per_acre,\"o\",markerfacecolor=\"none\",label=\"Corn (BU/ACRE)\")\n",
|
2023-01-12 20:25:12 +02:00
|
|
|
"plt.plot(wheat_year,wheat_bu_per_acre,\"x\",label=\"Wheat (BU/ACRE)\")\n",
|
2023-01-13 20:54:15 +02:00
|
|
|
"plt.plot(potatoes_year,potatoes_cwt_per_acre,\"+\",label=\"Potatoes (CWT/ACRE)\")\n",
|
|
|
|
|
"plt.plot(soybeans_year,soybeans_bu_per_acre,\"*\",markerfacecolor=\"none\",label=\"Soybeans (BU/ACRE)\")\n",
|
|
|
|
|
"plt.plot(sunflowers_year,sunflowers_lbs_per_acre,\"d\",markerfacecolor=\"none\",label=\"Sunflowers (LB/ACRE)\")\n",
|
2023-01-12 20:25:12 +02:00
|
|
|
"plt.plot(sugarbeets_year,sugarbeets_tons_per_acre,\".\",label=\"Sugarbeets (TONS/ACRE)\")\n",
|
|
|
|
|
"plt.xlabel(\"time (calendar year)\")\n",
|
|
|
|
|
"plt.ylabel(\"Average per acre production\")\n",
|
|
|
|
|
"plt.title(\"US Agriculture yields, NASS (USDA)\")\n",
|
|
|
|
|
"plt.legend()\n",
|
|
|
|
|
"plt.yscale(\"log\")\n",
|
2023-01-13 20:54:15 +02:00
|
|
|
"plt.grid()\n",
|
2023-01-12 16:13:39 +02:00
|
|
|
"#plt.show()\n",
|
2023-01-13 20:54:15 +02:00
|
|
|
"plt.savefig(\"raw_production_per_acre.pdf\")"
|
2023-01-12 16:13:39 +02:00
|
|
|
]
|
2023-01-12 20:25:12 +02:00
|
|
|
},
|
2023-01-12 21:45:54 +02:00
|
|
|
{
|
|
|
|
|
"cell_type": "markdown",
|
|
|
|
|
"id": "5828e06b",
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"source": [
|
|
|
|
|
"### All together, kcal/acre"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
|
|
|
|
"execution_count": null,
|
|
|
|
|
"id": "9d58d630",
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [],
|
|
|
|
|
"source": []
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
|
|
|
|
"execution_count": 21,
|
|
|
|
|
"id": "7903705c",
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [
|
|
|
|
|
{
|
|
|
|
|
"data": {
|
2023-01-13 20:54:15 +02:00
|
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAjoAAAHFCAYAAAD7ZFORAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOzdd1zV1f/A8ddl770VEBFQxIXmTHFrucqZmWlDSy1Lzcys1NIcDcu+SmmOfmZqmqmluQHNhTgTcONmb9nj/P643o9clqAgjvN8PHwIn3nu5XJ533Pe57xVQgiBJEmSJEnSE0inphsgSZIkSZJUXWSgI0mSJEnSE0sGOpIkSZIkPbFkoCNJkiRJ0hNLBjqSJEmSJD2xZKAjSZIkSdITSwY6kiRJkiQ9sWSgI0mSJEnSE0sGOpIkSZIkPbFkoPOIWLhwISqVCj8/v5puyiOnTp069O7du9R9YWFhqFQqVq5cqbV9x44ddO/eHRcXFwwNDXFxcaFjx47MnTu3Uvfu378/KpWKd955536bXyF16tRh5MiR1Xb9GTNmoFKptLYtXry4xPP2uOnYsSMdO3as1nOvXLlS6mvsfowcORIzM7MHvs6Trk6dOqhUKt5+++0S+4KDg1GpVGzYsKHUcyvyXpqYmMjUqVPx9fXF1NQUS0tL6tevz/Dhwzl9+rTWsUeOHOHFF1/Ezc0NQ0NDHB0dadOmDZMmTarw43n99dfp2bNnhR/DO++8U+L3NSMjg3nz5tGkSRMsLCwwNzfH09OTwYMHExISUuLamn8GBgbY29vTrl07pk2bxtWrV8tt673e8/bs2YOZmRk3b96s6MOvcTLQeUQsX74cgPDwcI4cOVLDrXm8/fjjj/Ts2RMLCwv+97//sWPHDubNm0eDBg3KfGMpTVxcHH///TcAq1evJjs7u7qazJ9//smnn35abdcvzZMQ6CxevJjFixfXdDOkarJs2TLOnTtXqXPu9V56+/ZtWrduzcqVK3nzzTfZsmULq1evZvTo0URFRXHy5Enl2K1bt9K2bVvS0tKYP38+O3fu5Pvvv6ddu3asW7euQu05ceIEv/zyC7NmzarU4yiqoKCA7t27M3v2bAYOHMj69evZsGEDEyZMIDU1lf3795c458svv+TQoUMEBQWxbNkyOnbsyPLly2nQoAGrV68u9T4Vec/r0qULLVu25OOPP77vx/PQCanGHT16VACiV69eAhCjRo166G0oLCwUmZmZD/2+FeHu7i569epV6j7Nc7dixQplm5ubm+jQoUOpxxcUFFT4vl999ZXWz2X16tWVandFPKznfPr06aL4r3vDhg1FQEBAld8rIyOjyq9ZHQICAir0+KOiokq8xu7XiBEjhKmp6QNf53FyP68Hd3d30aZNG2FpaSn69++vtS8oKEgAYv369SXOq8h76fLlywUg9u7dW+q9i75HdOjQQXh6eoq8vLxyjyvP4MGDRevWrSv8GIQQYty4cVq/r3v37hWAWL58+T3bXN61ExMTRbNmzYSenp44ffp0if0Vfc/bsGGD0NXVFdeuXSt1/6NG9ug8ApYtWwbA3Llzadu2LWvXriUzMxOAvLw8HBwcGD58eInzUlJSMDY2ZuLEicq2tLQ0PvjgAzw8PDAwMKBWrVq8//77ZGRkaJ2r6Zr88ccfadCgAYaGhvzyyy8AzJw5k1atWmFjY4OFhQX+/v4sW7YMUaz+a05ODpMmTcLJyQkTExM6dOjAsWPHSh2GiYmJ4a233qJ27doYGBjg4eHBzJkzyc/Pf+Dnr7jExEScnZ1L3aejU/GX/PLly3F0dOSXX37B2NhY+aRY3L///kubNm0wMjKiVq1afPrpp/z888+oVCquXLmiHKcZgtu4cSPNmjXDyMiImTNnKvuKP2cpKSlMmjSJunXrYmhoiIODA88//zxnz54F7nZRBwcHa51XkaGWOnXqEB4eTkhIiNLFXadOHQBWrlxZou1l3a9jx474+fmxb98+2rZti4mJCa+//jpQ8ddicV988QV6enpcv369xL7XX38dW1tb5ZNmacNPubm5zJo1i/r162NoaIi9vT2vvfYa8fHx5d4X4NatWwwePBhzc3MsLS0ZMmQIMTExJY67fPkyL730kjI06ujoSJcuXbR6A8oTHh5Oly5dMDU1xd7ennfeeUf5nQf1p+b69euX+J0TQlCvXj169epV5rXfeOMNbGxstK6n0blzZxo2bKh1vcWLF9O0aVOMjY2xtrZm4MCBXL58Weu8Xbt20a9fP2rXro2RkRH16tXjrbfeIiEhQes4zRDp8ePHGThwINbW1nh6elboOSnOxsaGjz76iI0bN3L48OEKnVPee6lGYmIiQIXeIxITE7Gzs0NPT6/c48oSGxvLn3/+Wer7d2VUps3lsbGx4aeffiI/P58FCxaU2F/R97w+ffpgZmbG0qVLK/gIaljNxllSZmamsLS0FM8884wQQoiff/5ZAGLlypXKMRMmTBDGxsYiNTVV69zFixcLQInMMzIyRNOmTYWdnZ349ttvxe7du8X3338vLC0tRefOnUVhYaFyLiBq1aolGjduLH777Texd+9ecebMGSGEECNHjhTLli0Tu3btErt27RJffPGFMDY2FjNnztS6/9ChQ4WOjo746KOPxM6dO8V3330nXF1dhaWlpRgxYoRyXHR0tHB1dRXu7u7ip59+Ert37xZffPGFMDQ0FCNHjrznc1TZHp2uXbsKPT09MX36dHHy5EmRn59/z3sUd+DAAQGIyZMnCyGEeOWVV4RKpRKXL1/WOu7UqVPCyMhING7cWKxdu1Zs2bJFPP/886JOnToCEFFRUVqPw9nZWdStW1csX75cBAUFidDQUGVf0ecsLS1NNGzYUJiamorPP/9c7NixQ/zxxx/ivffeUz6Jaj65BQUFabWptB6I4j06x48fF3Xr1hXNmjUThw4dEocOHRLHjx8XQgixYsWKEm0v634BAQHCxsZGuLq6ih9++EEEBQWJkJCQSr0Wi4uNjRWGhoZi2rRpWtsTExOFsbGx8jPR3L9or0xBQYHo2bOnMDU1FTNnzhS7du0SP//8s6hVq5bw9fXV6kErfm5mZqZo0KCBsLS0FD/88IPYsWOHGD9+vHBzcyvxfPr4+Ih69eqJVatWiZCQEPHHH3+ISZMmlfhZFDdixAhhYGAg3NzcxOzZs8XOnTvFjBkzhJ6enujdu7dy3ObNmwUgdu3apXX+1q1bBSC2bt1a5j1OnTolALF06VKt7eHh4QIQixYtUraNGjVK6Ovri0mTJont27eL3377TdSvX184OjqKmJgY5bjAwEAxZ84csWXLFhESEiJ++eUX0aRJE+Hj4yNyc3OV4zSvM3d3dzFlyhSxa9cusWnTpnKfk9JofuczMzNFrVq1RPv27ZV9ZfVYVOS9VAgh/v33XwGIZ555Rvz5558iISGhzHa8+eabAhDvvvuuOHz4sNZjrYj/+7//E4CIiIjQ2l7ZHp2oqCihr68vvL29xa+//ipu3bpV5j3vdW0hhHB2dhaenp5a2yr6nqfx3HPPCX9//zLv8SiRgU4N0/wi/Pjjj0IIIdLT04WZmZnWL/bp06cFIJYsWaJ1bsuWLUXz5s2V7+fMmSN0dHTE0aNHtY7bsGGDAMS2bduUbYCwtLQUSUlJ5bavoKBA5OXlic8//1zY2toqf6A0b5pTpkzROn7NmjUC0Pqj/dZbbwkzMzNx9epVrWO//vprAYjw8PBy21DZQOfixYvCz89PAAIQxsbGokuXLuJ///tfhd+oXn/9dQGIyMhIIcTdN49PP/1U67hBgwYJU1NTER8fr2wrKCgQvr6+pQY6urq64ty5c6U+xqLP2eeff17qH7qiHiTQEaLsoavKBjqA2LNnj9axlXktlmbEiBHCwcFB5OTkKNvmzZsndHR0tNpVPFjRvP7++OMPretpXieLFy8u89zAwEABiM2bN2udO2rUKK3nMyEhQQDiu+++K/cxlPW4APH9999rbZ89e7YAxL///iuEUL+G6tatK/r166d13HPPPSc8PT3LDRQ1j61p06Za28aMGSMsLCxEenq6EEKIQ4c
|
2023-01-12 21:45:54 +02:00
|
|
|
"text/plain": [
|
|
|
|
|
"<Figure size 640x480 with 1 Axes>"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"output_type": "display_data"
|
|
|
|
|
}
|
|
|
|
|
],
|
|
|
|
|
"source": [
|
2023-01-13 20:54:15 +02:00
|
|
|
"plt.plot(corn_year,corn_kcal_per_acre,\"o\",markerfacecolor=\"none\",label=\"Corn\")\n",
|
2023-01-12 21:45:54 +02:00
|
|
|
"plt.plot(wheat_year,wheat_kcal_per_acre,\"x\",label=\"Wheat\")\n",
|
|
|
|
|
"plt.plot(potatoes_year,potatoes_kcal_per_acre,\"+\",label=\"Potatoes\")\n",
|
2023-01-13 20:54:15 +02:00
|
|
|
"plt.plot(soybeans_year,soybeans_kcal_per_acre,\"*\",markerfacecolor=\"none\",label=\"Soybeans\")\n",
|
2023-01-12 21:45:54 +02:00
|
|
|
"#plt.plot(sugarbeets_year,sugarbeets_kcal_per_acre,\".\",label=\"Sugarbeets\")\n",
|
2023-01-13 20:54:15 +02:00
|
|
|
"plt.plot(sunflowers_year,sunflowers_kcal_per_acre,\"d\",markerfacecolor=\"none\",label=\"Sunflowers\")\n",
|
2023-01-12 21:45:54 +02:00
|
|
|
"plt.xlabel(\"time (calendar year)\")\n",
|
|
|
|
|
"plt.ylabel(\"kcal/acre production\")\n",
|
|
|
|
|
"plt.title(\"Average US Agriculture yields by year, NASS (USDA)\")\n",
|
|
|
|
|
"plt.legend()\n",
|
|
|
|
|
"plt.yscale(\"log\")\n",
|
|
|
|
|
"#plt.show()\n",
|
|
|
|
|
"#plt.savefig(\"bear-quadratic.pdf\")"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
2023-01-13 20:54:15 +02:00
|
|
|
"execution_count": 29,
|
2023-01-12 21:45:54 +02:00
|
|
|
"id": "db78da0c",
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [
|
|
|
|
|
{
|
|
|
|
|
"data": {
|
2023-01-13 20:54:15 +02:00
|
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAjMAAAHFCAYAAAAHcXhbAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOzdd3zN1//A8dfN3pE9JESQZQa1R+xdm6JqtF9taLWotkorVAdtVWml1db4KaUU1dIYEatoEKEkRoidJYnsfc/vj/ReuRkkJJFxno+Hh9zzWefce3PvO59zzvsohBACSZIkSZKkakrrWVdAkiRJkiTpachgRpIkSZKkak0GM5IkSZIkVWsymJEkSZIkqVqTwYwkSZIkSdWaDGYkSZIkSarWZDAjSZIkSVK1JoMZSZIkSZKqNRnMSJIkSZJUrclgphKtWLEChUJB06ZNn3VVqhwXFxcGDRpU7LbTp0+jUChYt26dRvnevXvp06cPjo6O6Ovr4+joiI+PD5999lmZrj18+HAUCgWvv/76k1a/VFxcXJg0aVKFnd/Pzw+FQqFRtmrVqiLPW3Xj4+ODj49PhR5748aNYt9jT2LSpEmYmJg89XlqOhcXFxQKBa+99lqRbYcOHUKhULBt27Zijy3NZ2l8fDxz587Fy8sLY2NjzM3N8fDwYMKECZw/f15j33/++Ydhw4ZRr1499PX1sbOzo0OHDsyePbvU7ZkyZQr9+vUrdRtef/31Ir+vaWlpLFmyhBYtWmBmZoapqSkNGzZk9OjRHD58uMi5Vf/09PSwsbGhU6dOzJs3j5s3bz6yro/7zAsMDMTExIS7d++WtvnPnAxmKtGaNWsAuHjxIv/8888zrk319t1339GvXz/MzMz45ptv2Lt3L0uWLMHT07PED4/ixMbG8ueffwKwceNGMjMzK6rK7Nixgw8++KDCzl+cmhDMrFq1ilWrVj3rakgV5KeffuLy5ctlOuZxn6Wpqam0b9+edevW8corr7Br1y42btzI1KlTiYyMJDQ0VL3v7t276dixI8nJySxdupR9+/bx9ddf06lTJ7Zs2VKq+pw9e5b169ezePHiMrWjoLy8PPr06cPHH3/MyJEj2bp1K9u2bWPmzJkkJSVx9OjRIsd88sknnDhxgqCgIH766Sd8fHxYs2YNnp6ebNy4sdjrlOYzr2fPnrRt25b333//idtT6YRUKU6dOiUAMXDgQAGI//3vf5VeB6VSKdLT0yv9uqVRv359MXDgwGK3qZ67tWvXqsvq1asnunbtWuz+eXl5pb7u559/rvG6bNy4sUz1Lo3Kes4XLFggCv9KN2nSRHTr1q3cr5WWllbu56wI3bp1K1X7IyMji7zHntTEiROFsbHxU5+nOnmS90P9+vVFhw4dhLm5uRg+fLjGtqCgIAGIrVu3FjmuNJ+la9asEYA4ePBgsdcu+BnRtWtX0bBhQ5GTk/PI/R5l9OjRon379qVugxBCTJ8+XeP39eDBgwIQa9aseWydH3Xu+Ph44e3tLXR0dMT58+eLbC/tZ962bduEtra2uHXrVrHbqxp5Z6aS/PTTTwB89tlndOzYkc2bN5Oeng5ATk4Otra2TJgwochxDx48wNDQkFmzZqnLkpOTefvtt2nQoAF6enrUrVuXt956i7S0NI1jVbcRv/vuOzw9PdHX12f9+vUALFy4kHbt2mFpaYmZmRmtWrXip59+QhRadzQrK4vZs2djb2+PkZERXbt25cyZM8V2mURHR/Pqq6/i5OSEnp4eDRo0YOHCheTm5j7181dYfHw8Dg4OxW7T0ir923rNmjXY2dmxfv16DA0N1X/xFXbs2DE6dOiAgYEBdevW5YMPPuDHH39EoVBw48YN9X6q7rLt27fj7e2NgYEBCxcuVG8r/Jw9ePCA2bNn4+rqir6+Pra2tgwYMIBLly4BD28nHzp0SOO40nSLuLi4cPHiRQ4fPqy+He3i4gLAunXritS9pOv5+PjQtGlTjhw5QseOHTEyMmLKlClA6d+LhX300Ufo6Ohw+/btItumTJmClZWV+i/G4rqKsrOzWbx4MR4eHujr62NjY8PkyZOJi4t75HUB7t27x+jRozE1NcXc3JwxY8YQHR1dZL/r16/zwgsvqLsx7ezs6Nmzp8Zf9Y9y8eJFevbsibGxMTY2Nrz++uvq33nI/+vXw8OjyO+cEIJGjRoxcODAEs/98ssvY2lpqXE+lR49etCkSRON861atYqWLVtiaGiIhYUFI0eO5Pr16xrH7d+/nyFDhuDk5ISBgQGNGjXi1Vdf5f79+xr7qbozQ0JCGDlyJBYWFjRs2LBUz0lhlpaWvPfee2zfvp2TJ0+W6phHfZaqxMfHA5TqMyI+Ph5ra2t0dHQeuV9JYmJi2LFjR7Gf32VRljo/iqWlJd9//z25ubl89dVXRbaX9jNv8ODBmJiY8MMPP5SyBc/Ys42laof09HRhbm4unnvuOSGEED/++KMAxLp169T7zJw5UxgaGoqkpCSNY1etWiUAdYSdlpYmWrZsKaytrcWyZcvEgQMHxNdffy3Mzc1Fjx49hFKpVB8LiLp164rmzZuLTZs2iYMHD4oLFy4IIYSYNGmS+Omnn8T+/fvF/v37xUcffSQMDQ3FwoULNa4/duxYoaWlJd577z2xb98+sXz5cuHs7CzMzc3FxIkT1ftFRUUJZ2dnUb9+ffH999+LAwcOiI8++kjo6+uLSZMmPfY5KuudmV69egkdHR2xYMECERoaKnJzcx97jcL+/vtvAYg5c+YIIYR48cUXhUKhENevX9fY79y5c8LAwEA0b95cbN68WezatUsMGDBAuLi4CEBERkZqtMPBwUG4urqKNWvWiKCgIBEcHKzeVvA5S05OFk2aNBHGxsZi0aJFYu/eveK3334Tb775pvovStVfYEFBQRp1Ku5OQuE7MyEhIcLV1VV4e3uLEydOiBMnToiQkBAhhBBr164tUveSrtetWzdhaWkpnJ2dxcqVK0VQUJA4fPhwmd6LhcXExAh9fX0xb948jfL4+HhhaGiofk1U1y94dyUvL0/069dPGBsbi4ULF4r9+/eLH3/8UdStW1d4eXlp3AkrfGx6errw9PQU5ubmYuXKlWLv3r1ixowZol69ekWeT3d3d9GoUSOxYcMGcfjwYfHbb7+J2bNnF3ktCps4caLQ09MT9erVEx9//LHYt2+f8PPzEzo6OmLQoEHq/X7//XcBiP3792scv3v3bgGI3bt3l3iNc+fOCUD88MMPGuUXL14UgPj222/VZf/73/+Erq6umD17tggICBCbNm0SHh4ews7OTkRHR6v38/f3F59++qnYtWuXOHz4sFi/fr1o0aKFcHd3F9nZ2er9VO+z+vXri3fffVfs379f7Ny585HPSXFUv/Pp6emibt26okuXLuptJd15KM1nqRBCHDt2TADiueeeEzt27BD3798vsR6vvPKKAMQbb7whTp48qdHW0vi///s/AYiwsDCN8rLemYmMjBS6urrCzc1N/Pzzz+LevXslXvNx5xZCCAcHB9GwYUONstJ+5qn0799ftGrVqsRrVCUymKkEqjf7d999J4QQIiUlRZiYmGj88p4/f14AYvXq1RrHtm3bVrRu3Vr9+NNPPxVaWlri1KlTGvtt27ZNAGLPnj3qMkCYm5uLhISER9YvLy9P5OTkiEWLFgkrKyv1l5Dqg/Hdd9/V2P+XX34RgMYX86uvvipMTEzEzZs3Nfb94osvBCAuXrz4yDqUNZiJiIgQTZs2FYAAhKGhoejZs6f45ptvSv1hNGXKFAGI8PBwIcTDD4gPPvhAY79Ro0YJY2NjERcXpy7Ly8sTXl5exQYz2tra4vLly8W2seBztmjRomK/zAp6mmBGiJK7mcoazAAiMDBQY9+yvBeLM3HiRGFrayuysrLUZUuWLBFaWloa9SockKjef7/99pvG+VTvk1WrVpV4rL+/vwDE77//rnHs//73P43n8/79+wIQy5cvf2QbSmoXIL7++muN8o8//lgA4tixY0KI/PeQq6urGDJkiMZ+/fv3Fw0bNnxkMKhqW8uWLTXKfH19hZmZmUhJSRFCCHH
|
2023-01-12 21:45:54 +02:00
|
|
|
"text/plain": [
|
|
|
|
|
"<Figure size 640x480 with 1 Axes>"
|
|
|
|
|
]
|
|
|
|
|
},
|
|
|
|
|
"metadata": {},
|
|
|
|
|
"output_type": "display_data"
|
|
|
|
|
}
|
|
|
|
|
],
|
|
|
|
|
"source": [
|
2023-01-13 20:54:15 +02:00
|
|
|
"plt.plot(corn_year,corn_kcal_per_acre/1.0e6,\"o\",markerfacecolor=\"none\",label=\"Corn\")\n",
|
|
|
|
|
"plt.plot(wheat_year,wheat_kcal_per_acre/1.0e6,\"x\",label=\"Wheat\")\n",
|
|
|
|
|
"plt.plot(potatoes_year,potatoes_kcal_per_acre/1.0e6,\"+\",label=\"Potatoes\")\n",
|
|
|
|
|
"plt.plot(soybeans_year,soybeans_kcal_per_acre/1.0e6,\"*\",markerfacecolor=\"none\",label=\"Soybeans\")\n",
|
2023-01-12 21:45:54 +02:00
|
|
|
"#plt.plot(sugarbeets_year,sugarbeets_kcal_per_acre,\".\",label=\"Sugarbeets\")\n",
|
2023-01-13 20:54:15 +02:00
|
|
|
"plt.plot(sunflowers_year,sunflowers_kcal_per_acre/1.0e6,\"d\",markerfacecolor=\"none\",label=\"Sunflowers\")\n",
|
2023-01-12 21:45:54 +02:00
|
|
|
"plt.xlabel(\"time (calendar year)\")\n",
|
2023-01-13 20:54:15 +02:00
|
|
|
"plt.ylabel(\"average yield (million kcal / acre)\")\n",
|
2023-01-12 21:45:54 +02:00
|
|
|
"plt.title(\"Average US Agriculture yields by year, NASS (USDA)\")\n",
|
|
|
|
|
"plt.legend()\n",
|
|
|
|
|
"plt.grid()\n",
|
|
|
|
|
"#plt.yscale(\"log\")\n",
|
|
|
|
|
"#plt.show()\n",
|
2023-01-13 20:54:15 +02:00
|
|
|
"plt.savefig(\"kcal_per_acre_yields.pdf\")"
|
2023-01-12 21:45:54 +02:00
|
|
|
]
|
|
|
|
|
},
|
2023-01-12 20:25:12 +02:00
|
|
|
{
|
|
|
|
|
"cell_type": "code",
|
|
|
|
|
"execution_count": null,
|
2023-01-12 21:45:54 +02:00
|
|
|
"id": "26e9a397",
|
2023-01-12 20:25:12 +02:00
|
|
|
"metadata": {},
|
|
|
|
|
"outputs": [],
|
|
|
|
|
"source": []
|
2023-01-12 16:13:39 +02:00
|
|
|
}
|
|
|
|
|
],
|
|
|
|
|
"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.9.13"
|
|
|
|
|
}
|
|
|
|
|
},
|
|
|
|
|
"nbformat": 4,
|
|
|
|
|
"nbformat_minor": 5
|
|
|
|
|
}
|
