dplR (version 1.5.3)

series.rwl.plot: Plot Series and a Master

Description

Plots a tree-ring series with a master chronology and displays their fit, segments, and detrending options in support of the cross-dating functions.

Usage

series.rwl.plot(rwl, series, series.yrs = as.numeric(names(series)),
                seg.length = 100, bin.floor = 100, n = NULL,
                prewhiten = TRUE, biweight = TRUE, floor.plus1 = FALSE)

Arguments

rwl

a data.frame with series as columns and years as rows such as that produced by read.rwl.

series

a numeric vector. Usually a tree-ring series.

series.yrs

a numeric vector giving the years of series. Defaults to as.numeric(names(series)).

seg.length

an even integral value giving length of segments in years (e.g., 20, 50, 100 years).

bin.floor

a non-negative integral value giving the base for locating the first segment (e.g., 1600, 1700, 1800 AD). Typically 0, 10, 50, 100, etc.

n

NULL or an integral value giving the filter length for the hanning filter used for removal of low frequency variation.

prewhiten

logical flag. If TRUE each series is whitened using ar.

biweight

logical flag. If TRUE then a robust mean is calculated using tbrm.

floor.plus1

logical flag. If TRUE, one year is added to the base location of the first segment (e.g., 1601, 1701, 1801 AD).

Value

• A list containing the filtered vectors series and master.

Details

The function is typically invoked to produce four plots showing the effect of the detrending options n and prewhiten and the binning options seg.length and bin.floor. [object Object],[object Object],[object Object],[object Object] The series and master are returned as well. See help pages for corr.rwl.seg, corr.series.seg, and ccf.series.rwl for more information on these arguments.

See Also

corr.rwl.seg, corr.series.seg, ccf.series.rwl

Examples

data(co021)
dat <- co021
flagged <- dat$"646244"
names(flagged) <- rownames(dat)
dat$"646107" <- NULL
foo <- series.rwl.plot(rwl = dat, series = flagged, seg.length = 100,
                       n = 5)
## note effect of n on first year in the series
foo <- series.rwl.plot(rwl = dat, series = flagged, seg.length = 100,
                       n = 13, prewhiten = FALSE)
bar <- series.rwl.plot(rwl = dat, series = flagged, seg.length = 100,
                       n = 7, prewhiten = FALSE)
head(foo$series)
head(bar$series)