astrochron (version 0.3.1)

rankSeries: Create lithofacies rank series from bed thickness data

Description

Create lithofacies rank series from bed thickness data.

Usage

rankSeries(dat,dt,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.

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)