teasel {popbio} | R Documentation |
Transition T and Fertility F matrices for the plant teasel
data(teasel)
A list with T and F matrices
Example 5.2
Caswell, H. 2001. Matrix population models: construction, analysis, and interpretation, Second edition. Sinauer, Sunderland, Massachusetts, USA.
data(teasel) A<-teasel$T + teasel$F A tea<-eigen.analysis(A, zero=FALSE) tea$lambda ## Summed elasticities? Growth in bottom-left triangle, stasis on diagonal, and # fertility in a prebreeding census for a monocarpic plant is in last column. elas<-tea$elasticities elas el<-c(fertility=sum(elas[,6]), stasis=sum(diag(elas)), growth=sum(elas[row(elas)>col(elas)])) el tea$sensitivities # Matrix plots like figure 9.4 in Caswell. No 3d barplot in R yet. Note log10 scale. op<-par(oma=c(1,2,4,0)) heatmap(log10(tea$sensitivities[6:1, ]), Rowv = NA, Colv = NA, scale = "none", margin = c(8,8), ylab = "Fate", xlab = "Stage", # col = gray(0:24 / 24)) col = heat.colors(24)) mtext("Sensitivity matrix for teasel", 3, outer=TRUE, line=1, cex=1.5) par(op) plot(log10(c(tea$sensitivities)), type="s", ylim=c(-5,2), xlab="Stage", ylab=expression(paste(Log[10], " sensitivity of ", lambda)), xaxt="n") axis(1, seq(1,36,6), 1:6) text(log10(c(tea$sensitivities)), cex=.7, pos=3, labels=paste(" ", 1:6, rep(1:6,each=6), sep=""))