R: Monthly Maximum 1-day Precipitation

climdex.rx1day {climdex.pcic} | R Documentation

## Monthly Maximum 1-day Precipitation

### Description

This function computes the climdex index Rx1day.

### Usage

    climdex.rx1day(ci, freq = c("monthly", "annual"))

### Arguments

`ci` | Object of type climdexInput.
`freq` | Time frequency to aggregate to.

### Details

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

### Value

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

### Note

The default behaviour of climdex.rx5day differs somewhat from fclimdex, as fclimdex and climdex.pcic differ on the definition of Rx5day. The running sum series computed by fclimdex is off by 2 days, and the first day a running sum can be computed for is left out entirely. The behaviour of fclimdex can be replicated by setting center.mean.on.last.day to TRUE.

### References

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

### See Also

`climdexInput.raw`, `climdexInput.csv`.

### Examples

    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 1-day precipitation.
    rx1day <- climdex.rx1day(ci)

[Package climdex.pcic version 1.1-11 Index]