rankSeries: Create lithofacies rank series from bed thickness data

In astrochron: A Computational Tool for Astrochronology

Create lithofacies rank series from bed thickness data

Description

Create lithofacies rank series from bed thickness data.

Usage

rankSeries(dat,dt,start=0,genplot=T,verbose=T)

Arguments

dat	First column should be bed thickness, and second column should bed lithofacies rank.
dt	Sampling interval for piecewise linear interpolation. By default a grid spacing that is 5 times smaller than the thinnest bed is used. If dt is set to zero, interpolation is skipped.
start	Start at what time/depth/height value?
genplot	Generate summary plots? (T or F)
verbose	Verbose output? (T or F)

Examples

# generate example series with random bed thicknesses
exThick=rnorm(n=20,mean=10,sd=2)
# assign alternating rank of 1 and 2
rank=double(20)
rank[seq(from=1,to=19,by=2)] <- 1
rank[seq(from=2,to=20,by=2)] <- 2

# combine into a dataframe
ex=cb(exThick,rank)

# generate lithofacies rank series
rankSeries(ex)

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