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

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

View source: R/splitA.R

Description

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

Usage

1
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 (and T is just A-F)

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 internally to split the matrix

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
data(whale)
splitA(whale)
# teasel -fertilitiles in last column
data(teasel)
splitA(teasel, r=1:6, c=6)
# hudsonia - fertilities in first two columns
data(hudsonia)
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")

popbio documentation built on May 4, 2018, 1:04 a.m.