Learn R Programming

IPMpack (version 1.6)

makeListFmatrix: Builds list of F (sexual reproduction) matrices from lists of growth, survival and fecundity objects.

Description

Uses posteriors (or other) lists of growth, fecundity and survival objects to create a list of F (sexual reproduction) matrices.

Usage

makeListFmatrix(fecObjList, nBigMatrix, minSize, 
	maxSize, cov = FALSE, envMat = NULL, integrateType = "midpoint",
		correction = "none")

Arguments

fecObjList
list of fecundity objects.
nBigMatrix
number of meshpoints.
minSize
minimum size.
maxSize
maximum size.
cov
level of covariate.
envMat
environmental matrix, defaults to NULL.
integrateType
integration type, defaults to "midpoint" (which uses probability density function); other option is "cumul" (which uses the cumulative density function).
correction
correction type, defaults to "none"; option is "constant" which will multiply every column of the IPM by a constant sufficient to adjust values to those predicted for survival at that size.

Value

  • list of Fmatrices.

References

Yang, Jongejans, Yang & Bishop. 2011. The effect of consumers and mutualists of Vaccinum membranaceum at Mount St. Helens: dependence on successional context. PLoS One 10, p1-11.

See Also

makeListPmatrix, makeListIPMs

Examples

Run this code
dff <- generateData()
fvlist <- makePostFecObjs(dataf = dff,burnin=100, nitt = 500)

# Only first 10 for speed
Flist <- makeListFmatrix(fecObjList = fvlist[1:10], 
	nBigMatrix = 20, minSize = -5, maxSize = 20, cov = FALSE, 
		envMat = NULL)

# Plots four F matrices at random:
par(mfrow=c(2,2))
image(Flist[[1]]@meshpoints, Flist[[1]]@meshpoints, 
	t(Flist[[1]]), xlab = "Continuous (e.g. size) stage at time t", 
		ylab = "Continuous (e.g. size) stage at time t+1")

image(Flist[[2]]@meshpoints, Flist[[2]]@meshpoints, 
	t(Flist[[2]]), xlab = "Continuous (e.g. size) stage at time t", 
	ylab = "Continuous (e.g. size) stage at time t+1")

image(Flist[[3]]@meshpoints, Flist[[3]]@meshpoints, 
	t(Flist[[3]]), xlab = "Continuous (e.g. size) stage at t", 
	ylab = "Continuous (e.g. size) stage at t+1")

image(Flist[[4]]@meshpoints, Flist[[4]]@meshpoints, 
	t(Flist[[4]]), xlab = "Continuous (e.g. size) stage at time t", 
	ylab = "Continuous (e.g. size) stage at time t+1")

Run the code above in your browser using DataLab