Learn R Programming

wux (version 2.2-1)

plot.wux.aov: 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

"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.

Examples

Run this code
## 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)

Run the code above in your browser using DataLab