Splits a projection matrix into transition and fertility matrices where
A = T + F.
a projection matrix
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
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 (and T is just A-F)
Caswell, H. 2001. Matrix population models: construction, analysis, and interpretation, Second edition. Sinauer, Sunderland, Massachusetts, USA.
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[] 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")
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.