ext.risk: Create an ext.risk Object

In BruceKendall/PVA: Population viability analysis in R

Description

Create an object of class ext.risk.

Usage

ext.risk(extrisk, year, nboot, ext.boot, alpha, ext.ci, model.type, model.params)

Arguments

extrisk	The cumulative probability of quasi-extinction projected at various times in the future
year	The times associated with the risks in extrisk
nboot	Number of bootstrap replicates used to calculate the confidence intervals
ext.boot	Extinction risk curves for each of the bootstrap replicates
alpha	One minus the confidence level
ext.ci	Confidence intervals of cumulative extinction risk
model.type	Name of the model used to perform the PVA
model.params	Estimated parameters of the model used to perform the PVA

Details

If extrisk and year are vectors of length n, then ext.boot has dimension (n, nboot) and ext.ci has dimension (n, 2). Each column of ext.boot represents a single bootstrap replicate.

If there are no bootstrap replicates, set nboot, ext.boot, alpha, and ext.ci all to NULL.

Value

An object of class ext.risk

Note

This function was developed in order to document the ext.risk class. In the normal course of events, there should be no need for the user to call this function.

Author(s)

Bruce E. Kendall ([e-mail])

See Also

print.ext.risk and plot.ext.risk for printing and plotting ext.risk objects; these objects are generated, for example, by count.DI.PVA

Examples

## Load the grizzly bear data and run a DI PVA
data(grizzly)
grizzly.ext <- count.DI.PVA(grizzly$N, Nx=20, Nc=50, nboot=1000)

## Print the confidence intervals for all 100 years
cbind(grizzly.ext$year, grizzly.ext$ext.ci)

## Plot a selection of bootstrap replicates
matplot(grizzly.ext$year, grizzly.ext$ext.boot[,1:100], type='l', col='black', 
    lty=1, xlab="Year", ylab="Cumulative extinction probability")

BruceKendall/PVA documentation built on Jan. 23, 2021, 2:56 a.m.