Apply Tukey's HSD over 1-unit depth slices.
slicedHSD(object, fm, conf = 0.95)
SoilProfileCollection
object
a formula describing depth sequence, horizon attribute, and site (grouping) attribute. For example 0:100 ~ estimated_oc | taxonname
confidence applied in TukeyHSD