ercrema / cultureallee Goto Github PK

We need to find an analytical solution for the equilibrium condition of the competition model. This is not Km,Kn, can be obtained using the simulation and it should be possible to solve this with the equations (possibly using Wolfram Alpha?).

I adedd some of the equations to the LaTeX manuscript here: fdd1ffa
We have to confirm that these are the final equations, and add the rest.

The current equation compares the intrinsic growth rate between the two populations. The intrinsic growth rate is given by these lines of code:

https://github.com/ercrema/CultureAllee/blob/56a663f39068f26528ff64ead65ba56653330c75/src/allee.R#L17-18

The problem is that before Rm and Rn were internally multiplied with m and n:
d497e27#L17-18

Now with the old version this guaranteed that Rm=0 when m=0, and Rn=0 when n=0. But with the current implementation this doesn't happen. So you might have a negative infinite growth rate. Try the following settings:

r.m= 0.4946667
r.n= 0.05
z = 0.03
An = 0.56
Am =0.1
Kn= 1.06
Km = 0.6
ini.m= 0.8615385
ini.n= 0
b=2

with all the competition coefficient set to 0.

Possible bug in the basin() function, results do not coincide with independent individual runs with the same parameter settings. Possible error generated by the parallelisation?

-[ ] Complete Analysis for Scenario 1
-[ ] Complete Analysis for Scenario 2
-[ ] Complete Analysis for Scenario 3
-[ ] Add discussion in the manuscript