# NOT RUN {
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