aqp (version 1.4)

get.ml.hz: Determine ML Horizon Boundaries

Description

This function accepts input from slab() along with a vector of horizon names, and returns a data.frame of the most likely horizonat boundaries.

Usage

get.ml.hz(x, hz.names)

Arguments

x

output from slab.SPC

hz.names

vector of horizon names

Value

• A dataframe with the following columns:
• hz
  horizon names
• top
  top boundary
• bottom
  bottom boundary
• confidence
  integrated probability / ML horizon thickness, rounded to the nearest integer

Details

This function is expecting that x is a data.frame containing columns identified in hz.names.

See Also

slab.SPC

Examples

data(sp1)

# normalize horizon names: result is a factor
sp1$name <- generalize.hz(sp1$name, 
  new=c('O','A','B','C'), 
  pat=c('O', '^A','^B','C'))

# compute slice-wise probability so that it sums to contributing fraction, from 0-150
a <- slab(sp1, fm= ~ name, class_prob_mode=1, seg_vect=0:150)

# generate table of ML horizonation
get.ml.hz(a, c('O','A','B','C'))