data(sp1)
# promote to SPC
depths(sp1) <- id ~ top + bottom
# set important metadata columns
hzdesgnname(sp1) <- "name"
hztexclname(sp1) <- "texture"
# set custom horizon property (clay content) column
hzmetaname(sp1, "clay") <- "prop"
# inspect metadata list
metadata(sp1)
# get horizon clay content column
hzmetaname(sp1, "clay")
# uses hzdesgname(), hztexclname(), hzmetaname(attr="clay") in function definition
estimatePSCS(sp1)
Run the code above in your browser using DataLab