growth: Annual Growth Rates by Categories.

Description Usage Arguments Details Value See Also Examples

Description

The principle growth function, constructed like recruitment and mortality. It requires two complete datasets, one per census, with dbh, pom, and date for every individual of all species in at least 2 censuses (see Data Format).

Usage

1
2
3
growth(census1, census2, rounddown = FALSE, method = "I",
  stdev = FALSE, dbhunit = "mm", mindbh = 10, growthcol = "dbh",
  err.limit = 4, maxgrow = 75, split1 = NULL, split2 = NULL)

Arguments

census1

The R Analytical Table for a single census, either tree or stem.

census2

The matching R Analytical Table for a later census.

rounddown

If TRUE, all dbh < 55 are rounded down to the nearest multiple of 5.

method

Use 'I' to calculate annual dbh increment: (dbh2 - dbh1)/time, or 'E' to calculate the relative growth rate (log(dbh2) - log(dbh1))/time.

stdev

Logical. Default (FALSE) returns confidence limits, otherwise returns the SD in growth rate per group.

dbhunit

'cm' or 'mm', only used for basal area

mindbh

The minimum diameter above which the counts are done. Trees smaller than mindbh are excluded. If NULL, all living trees are included.

growthcol

defines how growth is measured, either 'dbh'or 'agb'(agb=biomass)

err.limit

A number. Numbers such as 10000 are high and will return all measures.

maxgrow

A number. Numbers such as 10000 are high and will return all measures.

split1

a vector of categories, one per individual

split2

another vector of categories, one per individual

Details

It calculates the mean growth rate in one or more categories defined by the split variables, split1 and split2. The column date is required for annualizing rates.

The columns status and stemID are both required, in order to determine which stems should have dbh change calculated.

The function trim.growth() handles all checks for trees to include; excluded are (a) cases where the stemID changes, (b) extreme values based on err.limit and maxgrow, (c) and trees below a minimum dbh in the first census.

Value

A list with components:

Pass the list to assemble.demography() (in utilities.r) with type = "g" to convert the list to a data.frame.

See Also

trim.growth().

Growth requires fill.dimension() in utilities.r.

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
## Not run: 
CTFSplot("bci", 56)
growth.data = growth(bciex::bci12t5mini, bciex::bci12t6mini)
growth.data$rate
growth.data = growth(
  bciex::bci12t5mini, bciex::bci12t6mini, 
  split1 = bciex::bci12t5mini$sp
)
growth.data$rate
assemble.demography(grow.data, type = 'g')

## End(Not run)

forestgeo/ctfs documentation built on May 3, 2019, 6:44 p.m.