| bai.in | R Documentation |
Convert multiple ring-width series to basal area increment (i.e., ring area) going from the pith to the bark.
bai.in(rwl, d2pith = NULL)
rwl |
a |
d2pith |
an optional |
This converts ring-width series (mm) to ring-area series (mm squared)
(aka basal area increments) based on the distance between the
innermost measured ring and the pith of the tree. It is related to
bai.out, which calculates each ring area starting from
the outside of the tree and working inward. Both methods assume a
circular cross section (Biondi 1999). See the references below for
further details.
An object of class c("bai", "data.frame") containing the ring
areas for each series with the column names, row names and dimensions
of rwl. It carries the provenance record of
rwl, if it has one, and records how it was made in
attr(x, "dplR.bai"): the function (fun) and whether
d2pith was given (d2pith). See as.bai
for which functions take it.
DendroLab website: https://dendrolaborg.wordpress.com/
Code by Andy Bunn based on work from DendroLab, University of Nevada Reno, USA. Patched and improved by Mikko Korpela.
Biondi, F. (1999) Comparing tree-ring chronologies and repeated timber inventories as forest monitoring tools. Ecological Applications, 9(1), 216–227.
Biondi, F. and Qeadan, F. (2008) A theory-driven approach to tree-ring standardization: Defining the biological trend from expected basal area increment. Tree-Ring Research, 64(2), 81–96.
bai.out, as.bai
library(graphics)
library(stats)
library(utils)
## Toy
n <- 100
## Make three fake tree-ring series to show that these funcs work on rwl objects
base.series <- 0.75 + exp(-0.2 * 1:n)
rwl <- data.frame(x1 = base.series + abs(rnorm(n, 0, 0.05)),
x2 = base.series + abs(rnorm(n, 0, 0.05)),
x3 = base.series + abs(rnorm(n, 0, 0.05)))
## The inside out method
foo <- bai.in(rwl = rwl)
## The outside in method
bar <- bai.out(rwl = rwl)
## Identical
head(bar)
head(foo)
## Use gp data
data(gp.rwl)
data(gp.d2pith)
foo <- bai.in(rwl = gp.rwl, d2pith = gp.d2pith)
foo.crn <- chron(foo)
yrs <- time(foo.crn)
plot(yrs, foo.crn[, 1], type = "n",
xlab = "Year", ylab = expression(mm^2))
lines(yrs, foo.crn[, 1], col = "grey", lty = "dashed")
lines(yrs, caps(foo.crn[, 1], nyrs = 32), col = "red", lwd = 2)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.