Feasible environment for given 2 strains to coexist.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "2cf0cb26",
   "metadata": {},
   "source": [
    "### Given parameters g1, K1 and g2, K2 of the two strains, calculate the feasible dilution factor & initial resource concentration where the 2 strains coexist. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "7a35a019",
   "metadata": {},
   "outputs": [],
   "source": [
    "import tkinter\n",
    "import math\n",
    "\n",
    "#Here you can change the in-program parameters. stepD is how much greater next dilution coefficient checked is than the previous, measured in fractions of Dmax-1\n",
    "#(with dilution coefficient greater than Dmax coexistence is surely impossible)\n",
    "stepD = 0.01\n",
    "#stepy is how much greater next dimensionless resource concentration value is than the previous, measured in fractions of y0 \n",
    "#(with value of y0 of less coexistence is surely impossible).\n",
    "stepy = 0.1\n",
    "\n",
    "#In line 38 you can appoint a maximum value for dimensionless resource concentration value y.\n",
    "\n",
    "def destruction():\n",
    "    window.destroy()\n",
    "    return()\n",
    "def checkData():\n",
    "    global stepD\n",
    "    global stepy\n",
    "    g1 = float(entryg1.get())\n",
    "    g2 = float(entryg2.get())\n",
    "    k1 = float(entryk1.get())\n",
    "    k2 = float(entryk2.get())\n",
    "    if(g2>g1):g1, g2, k1, k2 = g2, g1, k2, k1\n",
    "    if (g1/g2)>=(k1/k2):\n",
    "        errorWindow = tkinter.Tk()\n",
    "        errorLabel = tkinter.Label(errorWindow, text = 'Error: ratio of growth rates g\\u2081/g\\u2082 is equal to or greater than ratio of Monod constants K\\u2081/K\\u2082.\\nIt is impossible for species with these characteristics to coexist on the same food source.')\n",
    "        errorLabel.grid(row = 0, column = 0)\n",
    "        errorButton = tkinter.Button(errorWindow, text = 'OK', command = errorDestruction)\n",
    "        errorButton.grid(row = 1, column = 0)\n",
    "    else:\n",
    "        possible = []\n",
    "        possibleD = []\n",
    "        possibler0 = []\n",
    "        gamma = g2/g1 \n",
    "        k = k1/k2 \n",
    "        y0 = (gamma*k-1)/((1-gamma)*k)\n",
    "        y = (1+stepy)*y0\n",
    "        ymax = 10*k\n",
    "        while y<ymax:\n",
    "            m = gamma*(k-1)/((1-gamma)*(y*k+1)-gamma*(k-1))\n",
    "            Dmax = (y*k+1)**m\n",
    "            D = 1+stepD*(Dmax-1)\n",
    "            while D<=Dmax and D<11000:\n",
    "                logD = math.log(D)\n",
    "                less = (y*k*D/(D-1)+1)*logD/((y*k*D/(D-1)+k)*logD+(k-1)*math.log(y*k+1))\n",
    "                more = ((y*k*D/(D-1)+1)*logD-(k-1)*math.log(y+1))/((y*k*D/(D-1)+k)*logD)\n",
    "                if less<gamma and more>gamma:\n",
    "                    r0 = y*k1\n",
    "                    possible+=[[D, r0]]\n",
    "                    possibleD+=[D]\n",
    "                    possibler0+=[r0]\n",
    "                elif len(possible)>0 and possible[len(possible)-1][1]==y*k1:\n",
    "                    break\n",
    "                D+=stepD*(Dmax-1)\n",
    "            y+=stepy*y0\n",
    "        for value in window.children.values():\n",
    "            value.grid_remove()\n",
    "        canvas = tkinter.Canvas(window, width = 480, height = 480, bg = 'white')\n",
    "        canvas.grid(row = 0, column = 0)\n",
    "        canvas.create_text(240, 470, text = 'Starting concentration')\n",
    "        canvas.create_text(10, 240, angle = 90, text = 'Dilution coefficient (logarithmic scale)')\n",
    "        canvas.create_line(40, 40, 40, 440)\n",
    "        canvas.create_line(40, 440, 440, 440)\n",
    "        f = open('data/bacteria_data.txt', 'w')\n",
    "        canvas.create_text(30, 440, text = '1', angle = 90)\n",
    "        startR1 = str(y0*k1//1)\n",
    "        startR2 = str(y0*k1)[0:6]\n",
    "        if len(startR1)>len(startR2):\n",
    "            startR = startR1\n",
    "        else:\n",
    "            startR = startR2\n",
    "        maxD1 = str(max(possibleD)//1)\n",
    "        maxD2 = str(max(possibleD))[0:6]\n",
    "        if len(maxD1)>len(maxD2):\n",
    "            maxD = maxD1\n",
    "        else:\n",
    "            maxD = maxD2\n",
    "        maxR1 = str(max(possibler0)//1)\n",
    "        maxR2 = str(max(possibler0))[0:6]\n",
    "        if len(maxR1)>len(maxR2):\n",
    "            maxR = maxR1\n",
    "        else:\n",
    "            maxR = maxR2\n",
    "        canvas.create_text(40, 450, text = startR)\n",
    "        canvas.create_text(30, 40, text = maxD, angle = 90)\n",
    "        canvas.create_text(440, 450, text = maxR)\n",
    "        for dataset in possible:\n",
    "            Daxis = 400*(math.log(dataset[0]))/(math.log(max(possibleD)))\n",
    "            Raxis = 400*(dataset[1]-y0*k1)/(max(possibler0)-y0*k1)\n",
    "            canvas.create_oval(39+Raxis, 439-Daxis, Raxis+41, 441-Daxis, outline = 'black', fill = 'black')\n",
    "            f.write(str(dataset[0])+', '+str(dataset[1])+'\\n')\n",
    "        f.close()\n",
    "        labelText = 'This plot is an approximation.\\n'\n",
    "        labelText+= 'Possible combinations written to file bacteria_data.txt.\\n'\n",
    "        labelText+= 'Every line of that file consists of two numbers\\n'\n",
    "        labelText+= 'separated by a comma and a space, where\\n'\n",
    "        labelText+= 'the first number is dilution coefficient,\\n'\n",
    "        labelText+= 'the second number is the corresponding\\n'\n",
    "        labelText+= 'starting concentration.\\n'\n",
    "        labelText+= 'The plot is clickable; if you click on a point in\\n'\n",
    "        labelText+= 'the plot area, corresponding points will be marked\\n'\n",
    "        labelText+= 'and labeled on the axes.'\n",
    "        label = tkinter.Label(window, text = labelText)\n",
    "        label.grid(row = 1, column = 0)\n",
    "        button = tkinter.Button(window, text = 'Exit', command = destruction)\n",
    "        button.grid(row = 2, column = 0)\n",
    "        def canvasClick(event):\n",
    "            if 40<=event.x and 440>=event.x and 40<=event.y and 440>=event.y:\n",
    "                canvas.delete('text')\n",
    "                canvas.create_line (event.x, 445, event.x, 435, tag = 'text')\n",
    "                canvas.create_line (35, event.y, 45, event.y, tag = 'text')\n",
    "                DaxisClick = 440-event.y \n",
    "                RaxisClick = event.x-40\n",
    "                Raxis = RaxisClick*(max(possibler0)-y0*k1)/400+y0*k1\n",
    "                Daxis = DaxisClick*math.log(max(possibleD))/400\n",
    "                D = math.e**Daxis\n",
    "                D1 = str(D//1)\n",
    "                D2 = str(D)[0:6]\n",
    "                if len(D1)>len(D2):\n",
    "                    textD = D1\n",
    "                else:\n",
    "                    textD = D2\n",
    "                canvas.create_text(30, event.y, text = textD, angle = 90, tag = 'text')\n",
    "                R1 = str(Raxis//1)\n",
    "                R2 = str(Raxis)[0:6]\n",
    "                if len(R1)>len(R2):\n",
    "                    textR = R1\n",
    "                else:\n",
    "                    textR = R2\n",
    "                canvas.create_text(event.x, 450, text = textR, tag = 'text')\n",
    "            return()\n",
    "        canvas.bind('<Button-1>', canvasClick)\n",
    "    return()\n",
    "\n",
    "window = tkinter.Tk()\n",
    "label = tkinter.Label(window, text = 'Enter the characteristics of your bacteria species.')\n",
    "label.grid(row = 0, column = 0, columnspan = 2)\n",
    "label1 = tkinter.Label(window, text = 'Species 1')\n",
    "label1.grid(row = 1, column = 0, columnspan = 2)\n",
    "labelg1 = tkinter.Label(window, text = 'Growth rate (g\\u2081): ')\n",
    "labelg1.grid(row = 2, column = 0)\n",
    "entryg1 = tkinter.Entry(window)\n",
    "entryg1.grid(row = 2, column = 1)\n",
    "labelk1 = tkinter.Label(window, text = 'Monod constant (K\\u2081): ')\n",
    "labelk1.grid(row = 3, column = 0)\n",
    "entryk1 = tkinter.Entry(window)\n",
    "entryk1.grid(row = 3, column = 1)\n",
    "label2 = tkinter.Label(window, text = 'Species 2')\n",
    "label2.grid(row = 4, column = 0, columnspan = 2)\n",
    "labelg2 = tkinter.Label(window, text = 'Growth rate (g\\u2082): ')\n",
    "labelg2.grid(row = 5, column = 0)\n",
    "entryg2 = tkinter.Entry(window)\n",
    "entryg2.grid(row = 5, column = 1)\n",
    "labelk2 = tkinter.Label(window, text = 'Monod constant (K\\u2082): ')\n",
    "labelk2.grid(row = 6, column = 0)\n",
    "entryk2 = tkinter.Entry(window)\n",
    "entryk2.grid(row = 6, column = 1)\n",
    "enterButton = tkinter.Button(window, text = 'OK', command = checkData)\n",
    "enterButton.grid(row = 7, column = 0)\n",
    "exitButton = tkinter.Button(window, text = 'Exit', command = destruction)\n",
    "exitButton.grid(row = 7, column = 1)\n",
    "window.mainloop()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "4e0b8a69",
   "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.8.8"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}