# \donttest{
if(require(soilDB) &
require(aqp) &
require(latticeExtra)) {
# sample data, an SPC
data(loafercreek, package='soilDB')
# convert to logical matrix
hp <- multinominal2logical(loafercreek, 'hillslopeprof')
# join-in to site data
site(loafercreek) <- hp
# variable names
v <- c('lithic.contact', 'paralithic.contact',
'argillic.horizon', 'toeslope', 'footslope',
'backslope', 'shoulder', 'summit')
# visualize with some other diagnostic features
x <- diagnosticPropertyPlot(loafercreek, v, k = 5,
grid.label = 'bedrckkind', dend.label = 'pedon_id')
}
# }
Run the code above in your browser using DataLab