mortality.eachspp {CTFS} | R Documentation |
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.
mortality.eachspp(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 |
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
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
.
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 |
Rick Condit
## 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)