pointer.rgc: Calculate pointer years using the relative growth change...

Description Usage Arguments Details Value Author(s) References Examples

View source: R/pointer.rgc.R

Description

The function calculates event and pointer years on a data.frame with tree-ring series using the relative growth change method, described as abrupt growth change method in Schweingruber et al. (1990). This method relates tree growth in year i to the average growth of n preceding years. Thresholds for event- and pointer-year calculations can be adjusted.

Usage

1
2
pointer.rgc(data, nb.yrs = 4, rgc.thresh.pos = 60, rgc.thresh.neg = 40, 
            series.thresh = 75)

Arguments

data

a data.frame with raw tree-ring series as columns and years as rows (e.g., output of read.rwl of package dplR).

nb.yrs

an integer specifying the number of preceding years to be used in calculating relative growth changes. Defaults to 4.

rgc.thresh.pos

a numeric specifying the threshold above which a relative growth change (in percentage) for a specific tree and year is considered a positive event year. Defaults to 60.

rgc.thresh.neg

a numeric specifying the threshold below which a relative growth change (in percentage) for a specific tree and year is considered a negative event year. Defaults to 40.

series.thresh

a numeric specifying the minimum percentage of trees that should display a positive (or negative) event year for that year to be considered as positive (or negative) pointer year. Defaults to 75.

Details

The function calculates the ratio of tree growth in year i and the average growth of n preceding years for individual trees. Resulting relative growth changes are used to identify event years for trees, and these event years to define pointer years for the site.

Following Schweingruber et al. (1990), nb.yrs, rgc.thresh.pos, rgc.thresh.neg and series.thresh are set to 4, 60, 40 and 75 respectively, meaning that a positive or negative pointer year will be defined when at least 75% of the tree-ring series display an event year with a growth increase or decrease of at least 60% or 40%, respectively, relative to the average growth in the 4 preceding years.

Note that the resulting time series are truncated by nb.yrs at the beginning inherent to the calculation methods.

Value

The function returns a list containing the following components:

rgc

a matrix with relative growth changes for individual tree-ring series

EYvalues

a matrix indicating positive (1), negative (-1) and non-event years (0) for individual tree-ring series

out

a data.frame containing the following columns:

year - time stamp

nb.series - number of series considered

perc.pos - percentage of trees showing a positive event year

perc.neg - percentage of trees showing a negative event year

nature - number indicating whether the year is a positive (1), negative (-1) or no pointer year (0)

dev_mean - mean growth deviation in percentage over the available series

dev_sd - standard deviation of the growth deviation

spec.param

a data.frame specifying the arguments used in the calculation

Author(s)

Marieke van der Maaten-Theunissen and Ernst van der Maaten.

References

Schweingruber, F.H., Eckstein, D., Serre-Bachet, F. and Br<c3><a4>ker, O.U. (1990) Identification, presentation and interpretation of event years and pointer years in dendrochronology. Dendrochronologia 8: 9-38.

In writing the function, the code of the dplR function pointer (Pierre M<c3><a9>rian) was used as a reference.

Examples

1
2
3
4
5
## Calculate pointer years on tree-ring series
data(s033)
py <- pointer.rgc(s033, nb.yrs = 4, rgc.thresh.pos = 60, rgc.thresh.neg = 40, 
                  series.thresh = 75)
py$out

pointRes documentation built on May 1, 2019, 8:45 p.m.