wq (version 0.4.6)

seasonTrend: Determine seasonal trends

Description

Finds the trend for each season and each variable in a time series.

Usage

seasonTrend(x, plot = FALSE, type = c("slope", "relative"), pval = .05, ...)

Value

A data frame with the following fields:seriesseries namesseasonseason numbersen.slopeSen slope in original units per yearsen.slope.relSen slope divided by median for that specific season and seriespp-value for the trend according to the Mann-Kendall test.missingProportion of slopes joining first and last fifths of the data that are missing

Details

The Mann-Kendall test is applied for each season and series (in the case of a matrix). The actual and relative Sen slope (actual divided by median for that specific season and series); the p-value for the trend; and the fraction of missing slopes involving the first and last fifths of the data are calculated (see mannKen).

If plot = TRUE, each season for each series is represented by a bar showing the trend. The fill colour indicates whether $p < 0.05$ or not. If the fraction of missing slopes is 0.5 or more, the corresponding trends are omitted.

Parameters can be passed to the plotting function, in particular, to facet_wrap in ggplot2. The most useful parameters here are ncol (or nrow), which determines the number of columns (or rows) of plots, and scales, which can be set to "free_y" to allow the y-axis to change for each time series. Like all ggplot2 objects, the plot output can also be customized extensively by modifying and adding layers.

See Also

mannKen, plotSeason, facet_wrap

Examples

x <- sfbayChla
seasonTrend(x)
seasonTrend(x, plot = TRUE, ncol = 4)