# need these
library(soilDB)
library(sharpshootR)
# soils of interest
s.list <- c('hornitos', 'perkins', 'argonaut', 'inks', 'mokelumne', 'dunstone', 'auburn', 'pentz', 'pardee', 'peters', 'amador', 'laniger')
# fetch and convert data into an SPC
s <- fetchOSD(s.list)
# estimate soil depth based on horizon designation
s$soil.depth <- profileApply(s, estimateSoilDepth, name='hzname', top='top', bottom='bottom')
# plot profiles, order by soil depth
par(mar=c(0,0,0,0))
plot(s, name='hzname', plot.order=order(s$soil.depth), cex.names=0.85, axis.line.offset=-5)
# plot dendrogram + profiles
SoilTaxonomyDendrogram(s, cex.taxon.labels=0.8)
This document is based on aqp
version 1.7-4 and soilDB
version 1.3-4.