{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": false }, "outputs": [], "source": [ "from cameo.strain_design.heuristic import ReactionKnockoutOptimization\n", "from cameo import load_model\n", "from cameo.strain_design.heuristic.objective_functions import biomass_product_coupled_yield, product_yield, number_of_knockouts\n", "from cameo.flux_analysis.simulation import fba\n", "import cameo\n", "import inspyred" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": false }, "outputs": [], "source": [ "model = load_model(\"../tests/data/iJO1366.xml\")" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "collapsed": false }, "outputs": [], "source": [ "of1 = product_yield(model.reactions.EX_ac_lp_e_rp_.id, model.reactions.EX_glc_lp_e_rp_.id)\n", "of2 = number_of_knockouts()\n", "of3 = biomass_product_coupled_yield(model.reactions.Ec_biomass_iJO1366_core_53p95M.id,\n", " model.reactions.EX_ac_lp_e_rp_.id,\n", " model.reactions.EX_glc_lp_e_rp_.id)" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": false }, "outputs": [], "source": [ "ko = ReactionKnockoutOptimization(model=model, objective_function=[of1, of2],\n", " simulation_method=fba, heuristic_method=inspyred.ec.emo.NSGA2)" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Starting optimization at Mon, 18 May 2015 00:08:14\n" ] }, { "data": { "text/html": [ " \n", " \n", " \n", " " ] }, "metadata": {}, "output_type": "display_data" }, { "name": "stdout", "output_type": "stream", "text": [ "Using saved session configuration for http://localhost:5006/\n", "To override, pass 'load_from_config=False' to Session\n" ] }, { "data": { "text/html": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/html": [ " 0%" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "application/javascript": [ "jQuery('#progress-bar-80219fc9-dd8d-445c-aa13-a2aa9700fb44').val('100')" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "application/javascript": [ "jQuery('#perc-progress-bar-80219fc9-dd8d-445c-aa13-a2aa9700fb44').html('0%')" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "application/javascript": [ "jQuery('#progress-bar-80219fc9-dd8d-445c-aa13-a2aa9700fb44').val('200')" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "application/javascript": [ "jQuery('#perc-progress-bar-80219fc9-dd8d-445c-aa13-a2aa9700fb44').html('0%')" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "application/javascript": [ "jQuery('#progress-bar-80219fc9-dd8d-445c-aa13-a2aa9700fb44').val('300')" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "application/javascript": [ "jQuery('#perc-progress-bar-80219fc9-dd8d-445c-aa13-a2aa9700fb44').html('0%')" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "application/javascript": [ "jQuery('#progress-bar-80219fc9-dd8d-445c-aa13-a2aa9700fb44').val('400')" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "application/javascript": [ "jQuery('#perc-progress-bar-80219fc9-dd8d-445c-aa13-a2aa9700fb44').html('0%')" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "application/javascript": [ "jQuery('#progress-bar-80219fc9-dd8d-445c-aa13-a2aa9700fb44').val('500')" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "application/javascript": [ "jQuery('#perc-progress-bar-80219fc9-dd8d-445c-aa13-a2aa9700fb44').html('1%')" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "application/javascript": [ "jQuery('#progress-bar-80219fc9-dd8d-445c-aa13-a2aa9700fb44').val('600')" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "application/javascript": [ "jQuery('#perc-progress-bar-80219fc9-dd8d-445c-aa13-a2aa9700fb44').html('1%')" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "application/javascript": [ "jQuery('#progress-bar-80219fc9-dd8d-445c-aa13-a2aa9700fb44').val('700')" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "application/javascript": [ "jQuery('#perc-progress-bar-80219fc9-dd8d-445c-aa13-a2aa9700fb44').html('1%')" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "application/javascript": [ "jQuery('#progress-bar-80219fc9-dd8d-445c-aa13-a2aa9700fb44').val('800')" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "application/javascript": [ "jQuery('#perc-progress-bar-80219fc9-dd8d-445c-aa13-a2aa9700fb44').html('1%')" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "application/javascript": [ "jQuery('#progress-bar-80219fc9-dd8d-445c-aa13-a2aa9700fb44').val('900')" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "application/javascript": [ "jQuery('#perc-progress-bar-80219fc9-dd8d-445c-aa13-a2aa9700fb44').html('1%')" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "application/javascript": [ "jQuery('#progress-bar-80219fc9-dd8d-445c-aa13-a2aa9700fb44').val('1000')" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "application/javascript": [ "jQuery('#perc-progress-bar-80219fc9-dd8d-445c-aa13-a2aa9700fb44').html('2%')" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "application/javascript": [ "jQuery('#progress-bar-80219fc9-dd8d-445c-aa13-a2aa9700fb44').val('1100')" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "application/javascript": [ "jQuery('#perc-progress-bar-80219fc9-dd8d-445c-aa13-a2aa9700fb44').html('2%')" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "application/javascript": [ "jQuery('#progress-bar-80219fc9-dd8d-445c-aa13-a2aa9700fb44').val('1200')" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "application/javascript": [ "jQuery('#perc-progress-bar-80219fc9-dd8d-445c-aa13-a2aa9700fb44').html('2%')" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "ko.run(max_evaluations=50000, n=10, mutation_rate=0.3, populations_size=100, crossover_rate=0.2)" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true }, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 2", "language": "python", "name": "python2" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 2 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython2", "version": "2.7.9" } }, "nbformat": 4, "nbformat_minor": 0 }