teasel {popbio}R Documentation

T and F matrix for teasel

Description

Transition T and Fertility F matrices for the plant teasel

Usage

data(teasel)

Format

A list with T and F matrices

Source

Example 5.2

References

Caswell, H. 2001. Matrix population models: construction, analysis, and interpretation, Second edition. Sinauer, Sunderland, Massachusetts, USA.

Examples

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=""))


[Package popbio version 1.1.2 Index]