rcs {dplR}R Documentation

Regional Curve Standardization

Description

Detrend multiple ring-width series simultaneously using a regional curve.

Usage

  rcs(rwl,po,nyrs=NULL,f=0.5,biweight=TRUE,rc.out=FALSE,
  make.plot=TRUE,...)

Arguments

rwl a data.frame with series as columns and years as rows such as that produced by read.rwl
po a data.frame containing two variables. Variable one (series in the example below) gives the series ID as either characters or factors. These must exactly match colnames(rwl). Variable two (pith.offset in the example below) must be integers and give the years from the beginning of the core to the pith (or center) of the tree.
nyrs a number giving the rigidity of the smoothing spline, defaults to 0.1 of length of the maximum cambial age (i.e., the length of the regional curve) if nyrs is NULL.
f a number between 0 and 1 giving the frequency response or wavelength cutoff. Defaults to 0.5.
biweight logical flag. If TRUE then a robust mean is calculated using tbrm.
make.plot logical flag. Makes plots of the raw data and regional curve if TRUE.
rc.out logical flag. Returns the regional curve along with the ring-width indices if TRUE.
... other arguments passed to plot.

Details

This method detrends and standardizes tree-ring series by calculating an age-related growth curve specific to the rwl. The detrending is the estimation and removal of the tree's natural biological growth trend. The standardization is done by dividing each series by the growth trend to produce units in the dimensionless ring-width index (RWI).

The “spline” approach uses an n-year spline where the frequency response is 0.50 at a wavelength of 10 percent of the maximum cambial age unless specified differently using nyrs and f in the function ffcsaps.

This attempts to remove the low frequency variability that is due to biological or stand effects. See the references below for further details on detrending in general, and Biondi and Qeadan (2008) for an explanation of RCS.

Value

A data.frame containing the dimensionless and detrended ring widths indices with column names, row names and dimensions of rwl. If rc.out is TRUE then a list will be returned with a data.frame containing the detrended ring widths as above and a vector containing the regional curve.

Note

DendroLab website: http://dendrolab.org/

Author(s)

Code provided by DendroLab based on work by F. Qeadan, University of Nevada Reno, USA and adapted for dplR by Andy Bunn

References

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.

Cook, E.R. and Kairiukstis, L.A. (1990) Methods of Dendrochronology: Applications in the Environmental Sciences. Springer. ISBN-13: 978-0792305866.

Fritts, H.C. (2001) Tree Rings and Climate. Blackburn. ISBN-13: 978-1930665392.

See Also

detrend, chron, cms, ffcsaps

Examples

  data(gp.rwl)
  data(gp.po)
  gp.rwi <- rcs(rwl=gp.rwl,po=gp.po,biweight=TRUE,rc.out=TRUE,make.plot=FALSE)
  str(gp.rwi)
  gp.rwi <- rcs(rwl=gp.rwl,po=gp.po,biweight=TRUE,make.plot=TRUE,main='Regional Curve')

[Package dplR version 1.2.5 Index]