https://github.com/cran/aqp
Raw File
Tip revision: 490762247abee87bbc3158f73ed5620d5377b55a authored by Dylan Beaudette on 28 November 2023, 05:20:10 UTC
version 2.0.2
Tip revision: 4907622
denormalize.Rd
% Generated by roxygen2: do not edit by hand
% Please edit documentation in R/denormalize.R
\name{denormalize}
\alias{denormalize}
\title{Create a (redundant) horizon-level attribute from a site-level attribute}
\usage{
denormalize(object, attr)
}
\arguments{
\item{object}{A SoilProfileCollection}

\item{attr}{Site-level attribute name (character string) to denormalize to horizon.}
}
\value{
A vector of values of equal length to the number of rows in the horizon table of the input SPC.
}
\description{
Create a (redundant) horizon-level attribute from a site-level attribute. Specify a SoilProfileCollection and a site-level attribute from that SPC (by name) to receive a vector of length equal to the number of horizons containing the site-level values. This vector is directly usable with the SoilProfileCollection horizon setter.

\code{denormalize} is the inverse operation for the formula interface that "normalizes" a horizon level variable to site level:

\code{site(object) <- ~ horizonvar}
}
\details{
"Denormalization" is the process of trying to improve the read performance of a database, at the expense of losing some write performance, by adding redundant copies of data or by grouping data. Sometimes it is beneficial to have site-level attributes denormalized for grouping of horizon-level data in analyses. \code{denormalize} achieves this result for SoilProfileCollections.
}
\examples{

data(sp1)

# create a SoilProfileCollection from horizon data
depths(sp1) <- id ~ top + bottom

# create random site-level attribute `sitevar` with a binary (0/1) outcome
sp1$sitevar <- round(runif(length(sp1)))

# use denormalize() to create a mirror of sitevar in the horizon table
# name the attribute something different (e.g. `hz.sitevar`) to 
# prevent collision with the site attribute
# the attributes can have the same name but you will then need 
# site() or horizons() to access explicitly
sp1$hz.sitevar <- denormalize(sp1, 'sitevar')

# compare number of profiles to number of sitevar assignments
length(sp1)
table(sp1$sitevar)

# compare number of horizons to number of horizon-level copies of sitevar `hz.'sitevar`
nrow(sp1)
table(sp1$hz.sitevar)
}
\author{
Andrew G. Brown, Dylan Beaudette
}
back to top