mortality.eachspp {CTFS}R Documentation

Annual Mortality Rate for Species by DBH Class

Description

Computes annual mortality rate for each species by DBH classes. It is provided as a shortcut for using using mortality with a split1 vector of species names and split2 vector of DBH classes. Any two census datasets can be provided in order of date of census. The annual mortality rate is computed for each species for each DBH class.

Usage

mortality.eachspp(census1, census2, dbhclass = c(10, 100, 300), 
        alivecode = c("A", "AB", "AS"))

Arguments

census1 name of census datafile for first census, must be a dataframe, must be of same length as census2
census2 name of census datafile for second census, must be a dataframe, must be of same length as census1
dbhclass lower class boundaries for DBH classes, classes are nonoverlapping
alivecode character, codes of the variable status that indicate the tree is alive. The most general valid categories are: "A" and "AB" and "AS".

Details

See CTFS.mortality andmortality for details on the computation of mortality rates and associated functions.

Any two censuses on a datafile must be used. They do not have to be sequential, only that census1 has to be before census2.

The results of mortality.eachspp can be organized into dataframes with the use of assemble.demography.

Value

mortality returns a list of arrays with the species as the first dimension and DBH classes as the second dimension of the array. The array contains the following named components:

$N0 the number of living trees at the first census
$S the number of surviving trees from the first to second census
$D the number of trees that died between the first and second census
$rate the mortality in %/year
$lower the lower 95% confidence interval
$upper the upper 95% confidence interval
$dbhmean the mean dbh in mm for trees used to compute mortality
$meanyrs mean number of years between census for trees used in mortality rate computation.
$date0 mean date of first census
$date1 mean date of second census

Author(s)

Rick Condit

See Also

CTFS.mortality

Examples

## Not run: 

1.  Default use of mortality.eachspp()
mort.eachspp.out1<-mortality.eachspp(tst.bci90.full, tst.bci95.full)

2.  User defined DBH classes

dbhclass = c(10, 20, 30, 40, 50, 60, 70, 80, 90, 100, 200, 300)
mort.eachspp.out2 <- mortality.eachspp(tst.bci90.full, classbreak = dbhclass)

3.  This is equivalent to:
mort.out <- 
mortality(tst.bci90.full,tst.bci95.full,split1=tst.bci90.full$sp,split2=dbhclass)
## End(Not run)

[Package CTFS version 1.00 Index]