stage.vector.plot: Plot stage vector projections

In cstubben/popbio: Construction and Analysis of Matrix Population Models

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,
  ...
)

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

Details

see section 2.2 in Caswell 2001

Value

A plot of stage or age class projections

Author(s)

Chris Stubben

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)

cstubben/popbio documentation built on April 2, 2024, 4:18 a.m.