calc_bins: Assign continuous values in a vector to discrete bins

In USGS-R/streamMetabolizer: Models for Estimating Aquatic Photosynthesis and Respiration

Assign continuous values in a vector to discrete bins

Description

Assigns each value in vec a new, discrete value corresponding to a bin. This function provides one interface to the functions 'base::cut', 'ggplot2::cut_interval', and 'ggplot2::cut_number'.

Usage

calc_bins(
  vec,
  method = c("bounds", "interval", "number", "width"),
  ...,
  bounds
)

Arguments

vec	the numeric vector whose values should be binned. log(discharge.daily) is a good candidate when using this function for pooling of K600 values.
method	a single character string indicating the automated bin selection method to use
...	other arguments (e.g. n, width) passed to the ggplot function corresponding to the value of cuts, if cuts is a character (otherwise ignored)
bounds	if method=='bounds', a numeric vector of bin boundaries

Examples

ln.disch <- log(rlnorm(100))

# for use in setting specs
brks <- calc_bins(ln.disch, 'width', width=0.8)$bounds
specs('b_Kb_oipi_tr_plrckm.stan', K600_lnQ_nodes_centers=brks)

# variations

# by 'number' method
bins_num <- calc_bins(ln.disch, 'number', n=5)
df_num <- data.frame(t=1:length(ln.disch), vec=ln.disch, bin=bins_num$names[bins_num$vec])
table(bins_num$vec)

# by 'interval' method
bins_int <- calc_bins(ln.disch, 'interval', n=5)
df_int <- data.frame(t=1:length(ln.disch), vec=ln.disch, bin=bins_int$names[bins_int$vec])
table(bins_int$vec)

# by 'width' method
bins_wid <- calc_bins(ln.disch, 'width', width=0.2, boundary=0)
df_wid <- data.frame(t=1:length(ln.disch), vec=ln.disch, bin=bins_wid$names[bins_wid$vec])
table(bins_wid$vec)

# choose your own arbitrary breaks
bins_arb <- calc_bins(ln.disch, bounds=seq(-4,4,by=1))
df_arb <- data.frame(t=1:length(ln.disch), vec=ln.disch, bin=bins_arb$names[bins_arb$vec])
table(bins_arb$vec)
## Not run: 
library(ggplot2)
ggplot(df_num, aes(x=t, y=vec, color=bin)) + geom_point() + 
  geom_hline(data=as.data.frame(bins_num['bounds']), aes(yintercept=bounds))
ggplot(df_int, aes(x=t, y=vec, color=bin)) + geom_point() +
  geom_hline(data=as.data.frame(bins_int['bounds']), aes(yintercept=bounds))
ggplot(df_wid, aes(x=t, y=vec, color=bin)) + geom_point() +
  geom_hline(data=as.data.frame(bins_wid['bounds']), aes(yintercept=bounds))
ggplot(df_arb, aes(x=t, y=vec, color=bin)) + geom_point() +
  geom_hline(data=as.data.frame(bins_arb['bounds']), aes(yintercept=bounds))

## End(Not run)

USGS-R/streamMetabolizer documentation built on Aug. 15, 2023, 7:50 a.m.