compcolM {primer} | R Documentation |
Multi-species competition colonization model, with habitat
destruction, after Nee and May (1992). For use with ode
in the deSolve
package.
compcolM(t, y, params)
t |
Argument for each time point |
y |
A vector for the populations |
params |
Vector or list of parameters |
Returns a list of length one, for use with ode
in the deSolve
package.
Component 1 |
vector of the state variables, y. |
Hank Stevens <HStevens@muohio.edu>
Nee, S. and May, R.M. (1992) Dynamics of metapopulations: habitat destruction and competitive coexistence. Journal of Animal Ecology, 61, 37–40.
Stevens. M.H.H. (2009) A Primer of Ecology with R. Use R! Series. Springer.
library(deSolve) S <- 10 ci <- 2^seq(-5, 5, length=S) m <- rep(.1, S) params <- list(ci=ci, m=m, S=S, D=0) init.N <- rep(0.01, S); t=seq(1, 200, .1) cc.out <- ode(init.N, t, compcolM, params) matplot(t, cc.out[, -1], type="l", ylab="Proportion of Habitat", xlab="Years")