# load test data
data(sp2)
# promote to SPC
depths(sp2) <- id ~ top + bottom
# assign two profiles a zone related to the mollic epipedon
newdata <- data.frame(id = c("hon-1","hon-17"),
                      featkind = "fixed-depth surface sample",
                      featdept = 0,
                      featdepb = 18)
# do left join
diagnostic_hz(sp2) <- newdata
# inspect site table: newvalue TRUE only for horizons
#  with top depth equal to zero
diagnostic_hz(sp2)
Run the code above in your browser using DataLab