Function for chemotaxis of bacteria to their prefered substrate
The generic function
chemotaxis implements a bacterial movement in the Moore neighbourhood to the highest substrate concentration.
1 2 3 4
An object of class Bac.
An object of class Arena.
The number of the iteration of interest.
Bacteria move to a position in the Moore neighbourhood which has the highest concentration of the prefered substrate, which is not occupied by other individuals. The prefered substance is given by slot
chem in the
Bac object. If there is no free space the individuals stays in the same position. If the concentration in the Moore neighbourhood has the same concentration in every position, then random movement is implemented.
1 2 3 4 5 6 7
data(Ec_core, envir = environment()) #get Escherichia coli core metabolic model bac <- Bac(Ec_core,deathrate=0.05, chem = "EX_o2(e)", minweight=0.05,growtype="exponential") #initialize a bacterium arena <- Arena(n=20,m=20) #initialize the environment arena <- addOrg(arena,bac,amount=10) #add 10 organisms arena <- addSubs(arena,40) #add all possible substances chemotaxis(bac,arena,1)