popbio (version 2.7)

stage.vector.plot: Plot stage vector projections

Description

Plots short-term dynamics and convergence to stage stage distribution using stage vector projections

Usage

stage.vector.plot(stage.vectors, proportions = TRUE,
  legend.coords = "topright", ylim = NULL, xlab = "Years",
  ylab = NULL, col = rainbow(8), ...)

Value

A plot of stage or age class projections

Arguments

stage.vectors

a matrix listing stage class vectors in columns

proportions

plot proportional changes or total numbers, defaults to proportions

legend.coords

a legend keyword or vector of x,y coordinates, defaults to top-right corner

ylim

the y limits of the plot, defaults to min and max values in stage.vectors

xlab

a label for the x axis

ylab

a label for the y axis

col

vector of line colors, defaults to rainbow(8)

...

additional options are passed to plot function

Author

Chris Stubben

Details

see section 2.2 in Caswell 2001

See Also

see pop.projection

Examples

## matrix from Example 2.1 in Caswell
A <- matrix2(c(
0, 0.3,   0,
1,   0, 0.5,
5,   0,   0
), 1:3)
n <- c(1,0,0)
p <- pop.projection(A,n,60)
## Plots in Figure 2.3
stage.vector.plot(p$stage.vector[,1:15], col='black', las=1, prop=FALSE)
stage.vector.plot(p$stage.vector[,1:40], col=2:4, las=1)
## log-scale with custom y-axis
stage.vector.plot(p$stage.vector, col=2:4, prop=FALSE,
ylim=c(.01, 10), log='y', legend="bottomright", yaxt='n')
pwrs <- -2:1
# major ticks
axis(2, at = 10^pwrs, labels=parse(text=paste("10^", pwrs, sep = "")),
las=1, tcl= -.6)
# minor ticks
axis(2, at = 1:9 * rep(10^pwrs[-1] / 10, each = 9),
    tcl = -0.3, labels = FALSE)