bio8: bio8: Calculates mean temperature of Wettest quarter

In ilyamaclean/climvars: Functions for calculating biologically meaningful climate variables

Description

bio8 is used to calculate the mean temperature in the wettest quarter of the year

Usage

bio8(temps, prec, tme1, tme2, method = "")

Arguments

temps	a vector of temperatures, normally for one year (see details).
prec	a vector of precipitation values, normally for one year (see details).
tme1	a POSIXlt object representing the date and time of each temps value.
tme2	a POSIXlt object representing the date and time of each prec value.
method	An optional character string describing the methods used to calculate mean temperature in the wettest quarter. Options are "anuclim" or unspecified (see details).

Details

If method is "anuclim", mean monthly temperature and total monthly precipitation is calculated and then spline interpolated to a weekly time period. Precipitation is calculated for each 13-week period and the mean temperature for the wettest period returned. If data spans more than one year, data are aggregated by unique month irrespective of year and one value returned. If method is unspecified, the mean temperature of the wettest quarter is calculated using all temps values and precipitation per quarter is calculated using the time interval for measurements. If data span more than one year, calculations are performed on all data and a single value returned.

Value

a single numeric value of mean temperature of the wettest quarter of the year.

See Also

the tmecreate() function can be used to create a POSIXlt object.

Examples

temps <- 10 * sin(c(0:1459) / (pi * 150)) + rnorm(1460)
tme1 <- tmecreate(2010, 6)
plot(temps~as.POSIXct(tme1), type = "l", xlab = "Month", ylab = "Temperature")
prec <- (10 * sin(c(0:364) * (pi / -360)) + rnorm(365) + 12)
tme2 <- tmecreate(2010, 24)
plot(prec~as.POSIXct(tme2), type = "l", xlab = "Month", ylab = "Precipitation")
bio8(temps, prec, tme1, tme2)
bio8(temps, prec, tme1, tme2, method = "anuclim")

ilyamaclean/climvars documentation built on June 19, 2019, 2:22 p.m.