damping.ratio {popbio} | R Documentation |
Calculate the damping ratio of a projection matrix
damping.ratio(A)
A |
A projection matrix |
see section 4.7 in Caswell (2001)
Damping ratio
The damping ratio is calculated by dividing the dominant eigenvalue by the eigenvalue with the second largest magnitude.
Chris Stubben
Caswell, H. 2001. Matrix population models: construction, analysis, and interpretation, Second edition. Sinauer, Sunderland, Massachusetts, USA.
## whale converges slowly to stable stage distribution
matplot2(pop.projection(whale, c(1,1,1,1), 60)$stage.vectors,
prop=TRUE, legend=NA,
main=paste("whale damping ratio = ", round(damping.ratio(whale),3) ) )
# Calathea - compare to Table 12 in Horvitz and Schemske (1995)
x <- sapply(calathea[-17], damping.ratio)
x <- matrix(x, nrow=4, byrow=TRUE, dimnames= list(paste("plot", 1:4), 1982:1985))
x
matplot2(x, type='b', ylab="Damping ratio", main="Calathea")