wux (version 2.1-1)

plotAnnualCycle: Plots the annual cycle

Description

plotAnnualCycle plots the monthly or seasonal annual cycle of indicated models and the box-whisker plots of the underlying distribution.

This plotting routine extracts all the information from the input data frame which has to be 'WUX-style' (see models2wux).

Usage

plotAnnualCycle(datain.df, var.name = NULL, subreg.subset = NULL,
season.subset = NULL, plot.quantiles = NULL, quantile.method = 7,
mark.df = NULL, plot.legend = FALSE, cex.names = 1.2, cex.lab = 1.2,
ylab = NULL, main = NULL, out.file.directory = NULL, out.file.name =
NULL, copyright = FALSE, ...)

Arguments

datain.df

WUX data frame obtained from models2wux.

var.name

Character string of parameter in WUX dataset.

subreg.subset

Vector of subregions to be plotted (e.g. c("EU.ENS", "GAR")).

season.subset

Vector of seasons to be plotted (e.g. c("MAM", "DJF")).

plot.quantiles

5 element vector indictaing the quantiles to be plotted (e.g. c(0.02,0.25,0.5,0.75,0.98)).

quantile.method

An integer between 1 and 9 selecting one of the nine quantile types in quantiles with default 7.

mark.df

Subset of WUX data frame indicating the models to be marked.

plot.legend

Boolean. Indicating if a plot legend indicating the models of mark.df and sample size should be plotted. Default is FALSE.

cex.names

Expansion factor for numeric axis labels in bxp. Default is 1.2.

cex.lab

Expansion factor for axis names (bar labels) in bxp. Default is 1.2.

ylab

Label for y-axis.

main

Main title.

out.file.directory

String of the directory where the plots are exported (e.g. "/tmp/plots/").

out.file.name

Prefix of the file names of the plots. Files will be stored as out.file.name_subreg_season.eps, where subreg is one realization of the subreg.subset argument and season is one realization of

copyright

Boolean. If a copyright message should be plotted. Default is FALSE.

...

Further optional arguments passed to bxp.

Examples

## load WUX and read WUX test data
require(wux)
data(ensembles)

wuxtest.df <- subset(ensembles, subreg == "GAR")

## set data frame for model marks
mark.df <- subset(wuxtest.df, acronym %in% c("ICTP-REGCM3", "MPI-M-REMO"))
mark.df2 <- gdata::drop.levels(mark.df)

plotAnnualCycle(wuxtest.df, "perc.delta.precipitation_amount", mark.df =
mark.df, plot.legend = TRUE, boxfill = "light yellow", notch =
FALSE,  boxwex = 0.5, ylim = c(-60,60), plot.quantiles =
c(0.02,0.25,0.5,0.75,0.98), boxcol = "red", ylab = "Precipitation
Amount [%]", main = "Annual cycle ", las = 1, copyright = TRUE)