timeseries: Time series

In brasmus/preciparea: Precipitation area ananlysis.

Description

Plots the time series

Usage

timeseries(x, is = 1, denominator='total area', ylab=NULL,...)

Arguments

x	a time series object
is	index - which series to plot
denominator	Denominator
ylab	see link{plot}
...	arguments passed on to plot

Author(s)

[e-mail]

Examples

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (x, is = 2, ...) 
{
    par(bty = "n")
    y <- x[, is]
    plot(y, main = attr(x, "longname")[is], plot.type = "single", 
        col = c("red", "blue", "grey"), sub = paste("TMRR ", 
            paste(abs(range(lat(x))), collapse = "S-"), "N: ", 
            min(year(y)), " - ", max(year(y)), sep = ""), ylab = attr(x, 
            "unit")[is], xlab = "", ylim = c(0, max(y)), ...)
    for (i in 1:length(is)) {
        lines(trend(y[, i]))
        t1 <- trend(y[, i])[1]
        t2 <- trend(y[, i])[length(index(y))]
        lines(c(as.Date("1990-01-01"), index(y)[1]), rep(t1, 
            2), col = rgb(0.5, 0.5, 0.5, 0.5), lwd = 2)
        lines(c(as.Date("1990-01-01"), index(y)[length(index(y))]), 
            rep(t2, 2), col = rgb(0.5, 0.5, 0.5, 0.5), lwd = 2)
        text(as.Date("1997-08-01"), t1, round(t1), cex = 0.7)
        text(as.Date("1997-08-01"), t2, round(t2), cex = 0.7)
    }
    grid()
  }

brasmus/preciparea documentation built on May 7, 2019, 11:11 a.m.