timetostratratebin: Transform Patterns from Time into Stratigraphic Height (Binned)

Description

This function transforms a given binned temporal pattern into the corresponding stratigraphic pattern. The transformation is based on the age model derived from the given deposition rate. Essentially a wrapper of patterntransform(...,direction='time to height',depositionmodel='binned deposition rate',patternmode='binned') kept for backwards compatibility of the package. It is recommended to use patterntransform instead of this function.

Usage

timetostratratebin(binborder, depoval, signalval, pos = NULL)

Arguments

binborder

Vector of strictly increasing numerical values. Defines the borders of the bins in time

depoval

Strictly positive vector of length length(binborder)-1. Defines the deposition rate in the bins, i.e. depoval[i] is the deposition rate in the time interval between binborder[i] and binborder[i+1]

signalval

Positive vector of length length(binborder)-1. Defines the temporal pattern in the bins, i.e. signalval[i] is the value of the temporal pattern between the times binborder[i] and binborder[i+1]

pos

OPTIONAL. Vector of points in time that will be transformed into stratigraphic height, and where the stratigraphic pattern will evaluated. Default setting is a vector spanning the interval of interest with evenly spaced points (no. is chosen by the code)

Value

A list containing the following entries:

height

vector containing stratigraphic heights

val

vector containing the values of the stratigraphic pattern at the height given by height

NA is returned if the values coincide with a hiatus or when the deposition rate/temporal pattern are undefined.

References

Hohmann, Niklas. 2018. Quantifying the Effects of Changing Deposition Rates and Hiatii on the Stratigraphic Distribution of Fossils. <doi:10.13140/RG.2.2.23372.51841>

Examples

Run this code

# NOT RUN {
## Define deposition rate
binborder=1:6 #temporal bins for the deposition rate
depoval=c(5,4,3,1,2) #deposition rate in the bins
#plot deposition rate
depositionrate=approxfun(binborder,c(depoval,tail(depoval,1)),method="constant",yleft=NA,yright=NA)
plot(depositionrate(seq(from=min(binborder),to=max(binborder),length.out=100)),xlab='time',
  ylab='deposition rate', main='Deposition rate')
## Define temporal pattern
signalval=c(1,0,5,2,1)
#plot temporal pattern
temporalpat=approxfun(binborder,c(signalval,tail(signalval,1)),method="constant",
  yleft=NA,yright=NA)
plot(temporalpat(seq(from=min(binborder),to=max(binborder),length.out=100)),xlab='time',
  ylab='Intensity', main='Temporal Pattern')

## Transform temporal pattern into stratigraphic pattern
reslist=timetostratratebin(binborder,depoval,signalval)
#plot resulting stratigraphic pattern
plot(reslist$height,reslist$val,xlab='Stratigraphic Height',ylab='Intensity',
  main='Stratigraphic Pattern')
## With removal of sediment (hiatus)
depoval=c(5,4,-3,1,2) #erosion rate is 3 in the middle time bin
reslist=timetostratratebin(binborder,depoval,signalval)
#plot resulting stratigraphic pattern
plot(reslist$height,reslist$val,xlab='Stratigraphic Height',ylab='Intensity',
  main='Stratigraphic Pattern')
# }