popbio (version 2.8)

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 = 1: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

...

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)