# path to local PedonPC back-end DB
dsn <- "S:/Service_Center/NRCS/pedon/pedon.accdb"
# get routinely used soil data SoilProfileCollection object
f <- fetchPedonPC(dsn)
# plot only those profiles with densic contact
plot(f[which(f$densic.contact), ], name='hzname')
Run the code above in your browser using DataLab