sensitivity: Sensitivity analysis of a projection matrix

In popbio: Construction and Analysis of Matrix Population Models

Sensitivity analysis of a projection matrix

Description

Calculate the sensitivities of eigenvalues to changes in the projection matrix elements

Usage

sensitivity(A, zero = FALSE)

Arguments

A	A projection matrix
zero	Set sensitivities for unobserved transitions to zero, default is FALSE

Details

see section 9.1 in Caswell (2001)

Value

A sensitivity matrix

Author(s)

Chris Stubben

References

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

See Also

elasticity

Examples

sens <- sensitivity(teasel)
## IMAGE plot with smaller boxes
image2(sens, mar=c(1,3.5,5,1), box.offset=.1)
 title("Sensitivity matrix using image2", line=2.5)
## MATPLOT
matplot2(sens, log='y', type='b', yaxt='n', ltitle="Fate",
 ylab=expression(paste("Sensitivity of ",lambda)),
 main="Sensitivity matrix using matplot2")
pwrs <- -4:1
axis(2, 10^pwrs, parse(text=paste("10^", pwrs, sep = "")), las=1)

popbio documentation built on May 29, 2024, 4:35 a.m.