Description Usage Arguments Details Value Note Author(s) References See Also Examples
Convert multiple ring-width series to basal area increment (i.e., ring area) going from the pith to the bark.
1 |
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.
A data.frame
containing the ring areas for each series with
column names, row names and dimensions of rwl
.
DendroLab website: http://dendrolab.org/
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.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 | 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, ffcsaps(foo.crn[, 1], nyrs = 32), col = "red", lwd = 2)
|
x1 x2 x3
1 7.923204 7.889906 8.380657
2 21.614002 21.199569 21.995239
3 30.485703 30.733817 31.762276
4 40.661058 37.593654 39.292309
5 45.222284 43.714669 45.807184
6 48.594832 48.562144 50.973255
x1 x2 x3
1 7.923204 7.889906 8.380657
2 21.614002 21.199569 21.995239
3 30.485703 30.733817 31.762276
4 40.661058 37.593654 39.292309
5 45.222284 43.714669 45.807184
6 48.594832 48.562144 50.973255
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.