wq (version 0.3-2)

utilities: Miscellaneous utility functions

Description

A variety of small utilities used in other functions. decyear2date: Converts decimal year to object of class "Date".

layerMean: Acts on a matrix or data frame with depth in the first column and observations for different variables (or different sites, or different times) in each of the remaining columns. The trapezoidal mean over the given depths is calculated for each of the variables. Replicate depths are averaged, and missing values or data with only one unique depth are handled. Data are not extrapolated to cover missing values at the top or bottom of the layer. The result can differ markedly from the simple mean even for equal spacing of depths, because the top and bottom values are weighted by 0.5 in a trapezoidal mean. leapYear: TRUE if x is a leap year, FALSE otherwise. meanSub: Mean of a subset of a vector. monthNum: Converts dates to the corresponding numeric month. tsSub: Drops seasons from a matrix or vector time series.

years: Converts dates to the corresponding numeric years.

Usage

decyear2date(x)
layerMean(d)
leapYear(x)
meanSub(x, sub, na.rm = FALSE)
monthNum(y)
tsSub(x1, seas = 1:frequency(x1))
years(y)

Rdversion

1.1

Examples

decyear2date(c(1996.001, 1999.999, 2004.162, 2005.162))

z = c(1,2,3,5,10)  # 5 depths
x = matrix(rnorm(30), nrow = 5)  # 6 variables at 5 depths
layerMean(cbind(z, x))

leapYear(seq(1500, 2000, 100))
leapYear(c(1996.9, 1997))

## Aggregate monthly time series over Feb-Apr only.
aggregate(sfbayChla, 1, meanSub, sub=2:4)

monthNum(as.Date(c('2007-03-17', '2003-06-01')))

## Ignore certain seasons in a Seasonal Kendall test.
c27 <- sfbayChla[, 's27']
seaKen(tsSub(c27))  # Aug and Dec missing the most key data
seaKen(tsSub(c27, seas = c(1:7, 9:11)))

y = Sys.time()
years(y)