| eigen.analysis {popbio} | R Documentation |
Calculate population growth rate and other demographic parameters from a projection matrix model using matrix algebra
eigen.analysis(A, zero=TRUE)
A |
a projection matrix |
zero |
Set sensitivities for unobserved transitions to zero |
A list with 6 items
lambda |
dominant eigenvalue |
stable.stage |
proportional stable stage distribution |
sensitivities |
matrix of eigenvalue sensitivities |
elasticities |
matrix of eigenvalue elasticities |
repro.value |
reproductive value scaled so v[1]=1 |
damping.ratio |
dominant eigenvalue divided by second largest eigenvalue |
Original code by James Holland Jones, Stanford University, Department of Anthropological Sciences, 12 August 2005 at http://popstudies.stanford.edu/summer_course/
Caswell, H. 2001. Matrix population models: construction, analysis, and interpretation, Second edition. Sinauer, Sunderland, Massachusetts, USA.
projection.matrix to construct matrix model
data(whale)
whaleA<-whale$T+whale$F
a<-eigen.analysis(whaleA)
a
ymax<-max(a$repro.value)*1.25
barplot(a$repro.value, col="blue", ylim=c(0, ymax ), xpd=FALSE,
ylab="Reproductive value", xlab="Killer whale stage class")
box()
## display sensitititivies of "impossible" transitions
## (postreprod -> juvenile, yearling -> mature, etc)
eigen.analysis(whaleA, zero=FALSE)$sensitivities
data(teasel)
teaselA<-teasel$T + teasel$F
a<-eigen.analysis(teaselA)
a
barplot(a$stable.stage, col="green", ylim=c(0,1),
ylab="Stable stage proportion", xlab="Teasel stage class")
box()