de Kroon, Plaisier, van Groenendael & Caswell. 1986. Elasticity: the relative contribution of demographic parameters to population growth rate. Ecology 67, p1427-1431.
de Kroon, van Groenendael & Ehrlen. 2000. Elasticities: a review of methods and model limitations. Ecology 81, p607-618.
R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> library(IPMpack)
Loading required package: Matrix
Loading required package: MASS
Loading required package: nlme
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/IPMpack/elas.Rd_%03d_medium.png", width=480, height=480)
> ### Name: elas
> ### Title: Estimates matrix element sensitivities and elasticities.
> ### Aliases: elas sens
>
> ### ** Examples
>
> dff <- generateData()
> Pmatrix <- makeIPMPmatrix(minSize = min(dff$size, na.rm=TRUE),
+ maxSize = max(dff$size, na.rm = TRUE), growObj = makeGrowthObj(dff),
+ survObj = makeSurvObj(dff))
>
> Fmatrix <- makeIPMFmatrix(minSize = min(dff$size, na.rm = TRUE),
+ maxSize = max(dff$size, na.rm = TRUE), fecObj = makeFecObj(dff))
[1] "Warning: fertility values < 0 exist in matrix, consider transforms. Negative values set to zero"
>
> IPM <- Pmatrix + Fmatrix
>
> par(mfrow = c(1, 2))
>
> senst <- sens(IPM)
> image(Pmatrix@meshpoints, Pmatrix@meshpoints,t(senst),
+ main = "Sensitivity", xlab = "Continuous (e.g. size) stage in t",
+ ylab = "Continuous (e.g. size) stage in t+1")
>
> elast <- elas(IPM)
> image(Pmatrix@meshpoints, Pmatrix@meshpoints, t(elast), main = "Elasticity",
+ xlab = "Continuous (e.g. size) stage in t",
+ ylab = "Continuous (e.g. size) stage in t+1")
>
>
>
>
>
> dev.off()
null device
1
>