There's an assumption in the last model (Metapopulations and propagule rain - Tutorial in EcoVirtual) that we can get rid off: the assumption that the propagule rain is constant. Instead, we can make the colonization rate a function of the percentage of occupied patches. In a simple formulation of this model, the propagule source is just from within the system (closed system) and the probability of colonization increases linearly with the proportion of occupied patches.

This way, our model does not have a constant colonization rate $p_i$, but a probability of colonization that depends on the number of occupied patches:

$$p_i = if$$

here $i$ is a constant that indicates how much the colonization probability increases with every new occupied patch. This way, the more occupied patches we have, the greater is the probability that new patches will be occupied. Substituting this new formula for $p_i$ in our old equation, we have:

$$\frac{df}{dt}=if(1-f)- p_e f $$

The fraction of the occupied patches in equilibrium(F) changes to:

$$ F=1-\frac{p_e}{i} $$

Let's see this first-hand by simulating this situation. As in the last exercise, we have created an R function that runs the simulation for us. As before, this function simply draws random colonization and extinction events for every patch at every time step, following the model's rules. It then shows the graph and an occupancy matrix for every time interval.

The parameters are:

Now you can simulate the model with the values you chose. For example,

tmax = 100
ncol = 10 
nrow = 10  
f0 = 0.1  
ci = 1
pe = 0.5

tmax = 100 
cl = 10
rw = 10  
f0 = 0.1 
ci = 0.5
pe = 0.5 

• Gotelli, N. 2007. Ecologia. Londrina, Ed. Planta. Capítulo 4.
• Stevens, M. H. 2009. A primer of ecology with R. New York. Springer.Capítulo 4.
• Gotelli, N. 1991. Metapopulation models: the rescue effect, the propagule rain, and the core-satellite hypothesis. The American Naturalist, 138: 768-776. pdf no site do autor