Skip to content

Instantly share code, notes, and snippets.

@jotech
Last active August 1, 2023 03:39
Show Gist options
  • Star 1 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save jotech/cd6c073ee99856b2ad0742aad9b8a31c to your computer and use it in GitHub Desktop.
Save jotech/cd6c073ee99856b2ad0742aad9b8a31c to your computer and use it in GitHub Desktop.
Adding oxygen after some time
library(BacArena)
data("Ec_core")
bac <- Bac(Ec_core)
arena <- Arena(n=20, m=20)
arena <- addOrg(arena,bac,amount=30)
arena <- addSubs(arena, smax=0.05, mediac="EX_glc(e)", unit="mM")
arena <- addSubs(arena, smax=1, mediac=c("EX_pi(e)", "EX_h2o(e)", "EX_nh4(e)"), unit="mM")
# Aenaerobic simulation for 6h
sim <- simEnv(arena, time=6, with_shadow = T)
# Add oxygen and continue simulation for another 6 hours
sim <- addSubs(sim, smax=0.05, mediac="EX_o2(e)", unit="mM")
sim <- simEnv(sim, time=6)
plotGrowthCurve(sim)
# anaerobic growth with production of fermentation products (acetate,formate),
# which stops after oxygen is added but continous with declining oxygen levels again.
plotSubCurve(sim, c("EX_glc(e)", "EX_o2(e)", "EX_for(e)", "EX_ac(e)"))[[2]]
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment