splitA: Split a projection matrix into separate T and F matrices

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

Split a projection matrix into separate T and F matrices

Description

Splits a projection matrix into transition and fertility matrices where A = T + F

Usage

splitA(A, r = 1, c = -1)

Arguments

A	a projection matrix
r	rows containing fertilities (default is first row) OR a logical matrix where TRUE is the location of a fertility value OR a complete fertility matrix
c	columns containing fertilities, default is all columns except first

Details

see section 5.1 in Caswell (2001)

Value

A list with T and F matrices

Note

By default, the fertility matrix will include elements in the first row (except first element). In some cases, it is not possible to split a projection matrix using only row and column indexes. Therefore, a logical matrix (where TRUE is the location of a fertility value) or the complete fertility matrix is also accepted.

Author(s)

Chris Stubben

References

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

See Also

functions like generation.time and net.reproductive.rate use splitA to split the matrix

Examples

splitA(whale)
# teasel -fertilitiles in last column
splitA(teasel, r=1:6, c=6)
# hudsonia - fertilities in first two columns
A <- hudsonia[[1]]
splitA(A, r=1:2)
## example using a logical matrix (if fertilities were in the upper diagonal)
splitA(A, row(A)<col(A))
# survival curves
x <- sapply(hudsonia, function(x) colSums(splitA(x, r=1:2)$T))
matplot2(t(x), legend="bottomright", ylab="Survival",
 main="Hudsonia survival curves")

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