Last data update: 2014.03.03

R: Barplots of the ANOVA results
plot.wux.aovR Documentation

Barplots of the ANOVA results

Description

Barplots of the aovWux results displaying the relative or absolute contribution of the individual factors to the overall variance.

Usage

## S3 method for class 'wux.aov'
plot(x, ss.relative = TRUE, subreg.subset =
NULL, cex.names = 1.2, cex.lab = 1.2, legend.text = NULL, sd.text =
TRUE, sd.unit = "", ylim = NULL, ylab = NULL, main = NULL,
out.file.directory = NULL, out.file.name = NULL, copyright = FALSE, ...)

Arguments

x

Object of class wux.aov obtained from the ANOVA aovWux.

ss.relative

Boolean. Indicating if the relative contribution of the factors to the overall variance should be calculated. Default is TRUE.

subreg.subset

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

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.

legend.text

String vector of the factors (e.g. c("GCM", "RCM", "RES")).

sd.text

Boolean. Indicating if the overall standard deviation should be displayed. Default is TRUE.

sd.unit

Character string of the standard deviation unit with default "" (e.g. "K").

ylim

Range vector for the y-axis.

ylab

Label for y-axis.

main

Main title.

out.file.directory

String of the directory where the plots are exported (e.g. "/tmp/plots/"). If neither out.file.name nor out.file.directory are passed, the plot will be displayed on screen.

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 season.subset. For example: out.file.name = "Barplot" will store to the files to Barplot_EUROPE_DJF.eps and Barplot_EUROPE_JJA.eps. If neither out.file.name nor out.file.directory are passed, the plot will be displyed on screen.

copyright

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

...

Further optional arguments passed to barplot.

Author(s)

Georg Heinrich g.heinrich@uni-graz.at

Examples

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

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

## unique model acronyms are required for reconstruction
wuxtest.df$acronym <- factor(paste(wuxtest.df$institute, "_",
wuxtest.df$rcm, sep="")) 

## reconstruction of the  missing data
reconstructLES.df <- reconstruct(wuxtest.df, factor1.name =
"acronym", factor2.name = "gcm", data.name =
"perc.delta.precipitation_amount") 

## calculate ANOVA
anova.list <- aovWux(perc.delta.precipitation_amount ~ acronym +
gcm, reconstructLES.df)

## barplot of ANOVA results
## Not run: plot(anova.list, ss.relative = TRUE, las = 1,
sd.unit = "%", legend.text = c("RCM", "GCM", "RES"), mgp = c(2.5,1,0),
main = "ANOVA Barplot", ylim = c(0,110))

## End(Not run)

Results