popbio (version 2.4.3)

splitA: 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

Value

A list with T and F matrices

Details

see section 5.1 in Caswell (2001)

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

# NOT RUN {
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")

# }