mwStatsGrid: 'Dynamic window' moving average, median and variance of...

In astrochron: A Computational Tool for Astrochronology

'Dynamic window' moving average, median and variance of stratigraphic series, using evenly spaced spatial/temporal grid

Description

'Dynamic window' moving average, median and variance of stratigraphic series. This routine adjusts the number of data points in the window so it has a constant duration in time or space, for use with unevenly sampled data. The results are output on an evenly spaced spatial/temporal grid (this contrasts with mwStats).

Usage

mwStatsGrid(dat,cols=NULL,win=NULL,step=NULL,start=NULL,end=NULL,output=T,norm=F,
            palette=6,ncolors=100,genplot=1,verbose=T)

Arguments

dat	Your data frame containing stratigraphic data; any number of columns (variables) are permitted, but the first column should be a location identifier (e.g., depth, height, time).
cols	A vector that identifies the variable column to be extracted (first column automatically extracted).
win	Moving window size, in units of space or time.
step	Moving window step size, in units of space or time.
start	Starting point for analysis, in units of space or time.
end	Ending point for analysis, in units of space or time.
norm	Normalize density estimates to maximum value? (T or F). If false, density estimates are normalized to unit area.
output	Output results? (T or F)
palette	What color palette would you like to use? (1) rainbow, (2) grayscale, (3) blue, (4) red, (5) blue-white-red, (6) viridis
ncolors	Number of colors to use in plot.
genplot	Generate summary plots? (0=none, 1=all time series, 2=kernel density estimates for each window, 3=kernel density estimates with median, 4=kernel density estimates with mean)
verbose	Verbose output? (T or F)

Value

A data frame with four columns: Center of window, Average, Median, Variance

Examples

# generate example series from ar1 noise, 5 kyr sampling interval
ex = ar1(npts=1001,dt=5)

# jitter sampling times
ex[1]=ex[1]+rnorm(1001,sd=1)
# sort
ex = ex[order(ex[,1],na.last=NA,decreasing=FALSE),]

# run mwStats
mwStatsGrid(ex,win=100)

astrochron documentation built on Sept. 30, 2024, 9:14 a.m.