I just find out that the simulation results from my BacArena system are the same with different values of tstep, i.e. the time (in h) per iteration.
What I did is
1) initialize the system:
arena <- BacArena::Arena(n=100, m=100, stir=TRUE, Lx=0.025, Ly=0.025, tstep=stephours)
2) set modles:
bac <- BacArena::Bac(model=Bsubtilis, deathrate=0.05, minweight=0.05, growtype="exponential", limit_growth=TRUE, lyse=TRUE)
arena <- BacArena::addOrg(arena, bac, amount=500)
3) simulate the system for 5 iterations:
sim <- BacArena::simEnv(arena, time=5)
I tried to set stephours=0.5, 1, 2, 4. However, the simulation results are almost the same.
Hi @jotech ,
I just find out that the simulation results from my BacArena system are the same with different values of tstep, i.e. the time (in h) per iteration. What I did is 1) initialize the system:
arena <- BacArena::Arena(n=100, m=100, stir=TRUE, Lx=0.025, Ly=0.025, tstep=stephours) 2) set modles: bac <- BacArena::Bac(model=Bsubtilis, deathrate=0.05, minweight=0.05, growtype="exponential", limit_growth=TRUE, lyse=TRUE) arena <- BacArena::addOrg(arena, bac, amount=500) 3) simulate the system for 5 iterations: sim <- BacArena::simEnv(arena, time=5)
I tried to set stephours=0.5, 1, 2, 4. However, the simulation results are almost the same.
Did you know the reason? Thanks!
Best, Xuejian