mortality.dbh {CTFS} | R Documentation |
Computes annual mortality rate for DBH classes.It is provided as a
shortcut for using using mortality
with a split1
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 DBH class.
mortality.dbh(census1, census2, dbhclass = c(10, 100, 300), alivecode = c("A", "AB", "AS"))
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 |
vector of 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". |
See CTFS.mortality
and mortality
for details on the computation of
mortality rates and associated functions.
for
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.dbh
can be organized into
dataframes with the use of
assemble.demography
.
mortality.dbh
returns a list of arrays with 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 |
Rick Condit
## Not run: 1. Default use of mortality mort.dbh.out <- mortality.dbh(tst.bci90.full) 2. User defined DBH classes dbhclass=c(10, 20, 30, 40, 50, 60, 70, 80, 90, 100, 200, 300) mort.dbh.out <- mortality.dbh(tst.bci90.full, classbreak = dbhclass) ## End(Not run)