diagnostic_hz-set: Add data to the diagnostic slot

In aqp: Algorithms for Quantitative Pedology

Add data to the diagnostic slot

Description

Diagnostic data in an object inheriting from data.frame can easily be added via merge (LEFT JOIN). There must be one or more same-named columns containing profile ID on the left and right hand side to facilitate the join: diagnostic_hz(spc) <- newdata

Usage

diagnostic_hz(object) <- value

Arguments

object	A SoilProfileCollection
value	An object inheriting data.frame

Examples

# 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)

aqp documentation built on April 12, 2022, 1:05 a.m.