splitA: Split a projection matrix into separate T and F matrices
In popbio: Construction and Analysis of Matrix Population Models

Description Usage Arguments Details Value Note Author(s) References See Also Examples

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

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

`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

see section 5.1 in Caswell (2001)

A list with T and F matrices

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.

Chris Stubben

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

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

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")