succniche: A Four-state model of Successional Dynamics
In primer: Functions and Data for the Book, a Primer of Ecology with R

Description Usage Arguments Value Author(s) References See Also Examples

This is the five-state, two-species model of the succession-niche model, after Pacala and Rees (1998). For use with ode in the deSolve package.

1	succniche(t, y, params)

`t`	Argument for the time point at integration
`y`	A vector of length four, for states, E, M, S, and R.
`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 variable (a scalar for the proportion of sites occupied).

Hank Stevens <HStevens@muohio.edu>

Pacala, S.W. and Rees, M. (1998) Models suggesting field experiments to test two hypotheses explaining successional diversity. The American Naturalist, 152, 729–737.

Stevens, M.H.H. (2009) A Primer of Ecology with R. Use R! Series. Springer.

levins, compcol, compcolM

params.suc <- c(a=7, c=0.2, g=.1, m=0.04, D=0)
t=seq(0,50,.1)
init.suc <- c(S=0, E=0.5, M=0.5, R=0.00) # Free space is implicit, F=1-(S+E+M+R).
ccg.out <- data.frame(ode(init.suc, t, succniche, params.suc))
matplot(t, ccg.out[,-1], type="l", ylab="Relative Frequency",
         xlab="Time", ylim=c(0,1) )
legend("right", colnames(ccg.out)[5:2], lty=4:1,  bty="n")