climdex.pcic (version 1.1-11)

climdex.rx5day: Monthly Maximum 5-day Consecutive Precipitation

Description

This function computes the climdex index Rx5day.

Usage

climdex.rx5day(ci, freq = c("monthly", "annual"),
  center.mean.on.last.day = FALSE)

Arguments

ci

Object of type climdexInput.

freq

Time frequency to aggregate to.

center.mean.on.last.day

Whether to center the 5-day running mean on the last day of the window, instead of the center day.

Value

A vector containing the value of the index for each month of each year.

Details

This function takes a climdexInput object as input and computes the climdex index Rx5day: monthly or annual maximum 5-day consecutive precipitation.

References

http://etccdi.pacificclimate.org/list_27_indices.shtml

See Also

climdexInput.raw, climdexInput.csv.

Examples

# NOT RUN {
library(PCICt)

## Create a climdexInput object from some data already loaded in and
## ready to go.

## Parse the dates into PCICt.
tmax.dates <- as.PCICt(do.call(paste, ec.1018935.tmax[,c("year",
"jday")]), format="%Y %j", cal="gregorian")
tmin.dates <- as.PCICt(do.call(paste, ec.1018935.tmin[,c("year",
"jday")]), format="%Y %j", cal="gregorian")
prec.dates <- as.PCICt(do.call(paste, ec.1018935.prec[,c("year",
"jday")]), format="%Y %j", cal="gregorian")

## Load the data in.
ci <- climdexInput.raw(ec.1018935.tmax$MAX_TEMP,
ec.1018935.tmin$MIN_TEMP, ec.1018935.prec$ONE_DAY_PRECIPITATION,
tmax.dates, tmin.dates, prec.dates, base.range=c(1971, 2000))

## Create a timeseries of monthly maximum 5-day consecutive precipitation.
rx5day <- climdex.rx5day(ci)
# }