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 bark to the pith.

1 |

`rwl` |
a |

`diam` |
an optional |

This converts ring-width series (mm) to ring-area series (mm squared)
(aka basal area increments) based on the diameter of the tree and the
width of each ring moving towards the pith of the tree. It is related
to `bai.in`

, which calculates each ring area starting from
the inside of the tree and working outward. 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 31 32 33 34 35 | ```
library(graphics)
library(utils)
## Not run:
library(stats)
## 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)
## End(Not run)
## Use gp data
data(gp.rwl)
data(gp.dbh)
## dbh (minus the bark) from cm to mm
gp.dbh2 <- gp.dbh[, 1:2]
gp.dbh2[, 2] <- (gp.dbh[, 2] - gp.dbh[, 3]) * 10
bar <- bai.out(rwl = gp.rwl, diam = gp.dbh2)
bar.crn <- chron(bar)
yr <- as.numeric(rownames(bar.crn))
plot(yr, bar.crn[, 1], type = "n",
xlab = "Year", ylab = expression(mm^2))
lines(yr, bar.crn[, 1], col = "grey", lty = "dashed")
lines(yr, ffcsaps(bar.crn[, 1], nyrs = 32), col = "red", lwd = 2)
``` |

```
x1 x2 x3
1 7.770198 7.824806 8.029568
2 21.114120 21.185551 20.938930
3 32.502043 31.525539 31.693033
4 39.119675 39.752321 40.045688
5 44.495713 44.026240 45.698924
6 48.650308 51.379409 52.295201
x1 x2 x3
1 7.770198 7.824806 8.029568
2 21.114120 21.185551 20.938930
3 32.502043 31.525539 31.693033
4 39.119675 39.752321 40.045688
5 44.495713 44.026240 45.698924
6 48.650308 51.379409 52.295201
```

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.