plot.wux.df: X - Y Scatterplot of climate change signals
In wux: Wegener Center Climate Uncertainty Explorer

Description Usage Arguments Author(s) Examples

plot.wux.df plots one or more scatterplots containing climate change signals of selected meteorological parameters.

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

## S3 method for class 'wux.df'
plot(x,
               var1.name = "delta.air_temperature",
               var2.name = "perc.delta.precipitation_amount",
               subreg.subset = NULL,
               season.subset = NULL, boxplots = TRUE,
               label.only.these.models = NULL, highlight.models = NULL,
                no.text = FALSE,
               vert.box.col = "cyan", horiz.box.col = "coral",
               zero.line.col = "gray80", median.line.col = "black",
               draw.legend = TRUE, draw.seperate.legend = FALSE,
               draw.median.lines = TRUE, use.rainbow.colors = TRUE,
               xlim = NULL, ylim = NULL, xlab = NULL, ylab = NULL,
               main = NULL, out.file.directory = NULL,
               out.file.name = NULL, copyright = FALSE, ...)

`x`	WUX data.frame (class `wux.df`) obtained from `models2wux`'.
`var1.name`	Character string of 1st parameter in WUX dataset. Default is temperature change.
`var2.name`	Character string of 2nd parameter in WUX dataset. Default is precipitation change.
`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")`).
`boxplots`	Boolean. Indicating if marginal boxplots for the two input parameters should be plotted. Default is `TRUE`.
`label.only.these.models`	Character vector of modelnames (acronyms) to be labeled in the scatterplot.
`highlight.models`	Character vector of modelnames (acronyms) to be highlighted in the scatterplot.
`no.text`	Boolean. Indicating if no models should be labeled. Default is `FALSE`.
`vert.box.col`	Color character for vertical boxplot. Default is `coral`.
`horiz.box.col`	Color character for horizontal boxplot. Default is `cyan`.
`zero.line.col`	Color character for the zero lines. Default is `gray80`.
`median.line.col`	Color character for the median lines. Default is `black`.
`use.rainbow.colors`	Boolean. Use `rainbow()` color palette if `TRUE`, otherwise a custom color palette with 17 colors is used. Default is `TRUE`.
`xlim`	Range vector for 1st parameter (x-axis).
`ylim`	Range vector for 2nd parameter (y-axis).
`xlab`	Label of 1st parameter (x-axis).
`ylab`	Label of 2nd parameter (y-axis).
`draw.legend`	Boolean. Indicating if legend with GCMs should be plotted. Default is `TRUE`.
`draw.seperate.legend`	Boolean. Should legend with GCMs be plotted on a seperate screen? Default is `FALSE`. Draws legend even if `draw.legend` is set `FALSE`.
`draw.median.lines`	Draw median lines for both parameters. Default is `TRUE`.
`main`	Main title.
`out.file.directory`	Directory where the plots shall be exported (e.g. "/tmp/plots/"). If neither `out.file.name` nor `out.file.directory` are passed, the plot will be displyed 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 = "scatterplot"` will store to the files `scatterplot_EUROPE_DJF.eps` and `scatterplot_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 to be passed to `plot`, such as graphical parameters (see `par`).

Thomas Mendlik thomas.mendlik@uni-graz.at and Georg Heinrich g.heinrich@uni-graz.at

require(wux)

### ENSEMBLES RCM analysis
data(ensembles)

## Not run: plot(ensembles, "perc.delta.precipitation_amount",
  "delta.air_temperature", boxplots = TRUE, xlim = c(-40,40),
  ylim = c(0, 4), label.only.these.models = c("ICTP-REGCM3", "MPI-M-REMO"), 
  xlab = "Precipitation Amount [%]", ylab = "2-m Air Temperature [K]",
  main = "Scatterplot", subreg.subset = c("GAR"))

## End(Not run)

### now see where ENSMEBLES GCMs lie within CMIP3 ensemble
data(ensembles_gcms) # GCMs for forcing of ENSEMBLES RCMs
data(cmip3_2050)     # GCMs of CMIP3 ensemble

ensembles.gcm.names <- levels(ensembles_gcms$acronym) #8 GCM names

cmip3_2050.sub <- subset(cmip3_2050, subreg %in% c("World", "EU.ENS")
                       & em.scn == "A1B") 
cmip3_2050.sub <- droplevels(cmip3_2050.sub)
ensembles_gcms.sub <- subset(ensembles_gcms, !acronym %in%
                                c("mpi_echam5-r3", "bccr_bcm2_0-r1",
                                  "ipsl_cm4-r2"))
ensembles_gcms.sub <- droplevels(ensembles_gcms.sub)
## combine cmip3 and ENSEMBLES GCMs in one data.frame
gcms.combined <- rbind(ensembles_gcms.sub, cmip3_2050.sub)

## Scatterplot
prec.range <- range(gcms.combined$perc.delta.precipitation_amount) + c(-1, 1)
tas.range <- range(gcms.combined$delta.air_temperature)
## Not run: plot(gcms.combined,
               "perc.delta.precipitation_amount", "delta.air_temperature",
               subreg.subset = "EU.ENS", draw.median.lines = FALSE,
               label.only.these.models = ensembles.gcm.names,
               xlim = prec.range,
               ylim = tas.range,
               main = "GCMs from ENSEMBLES project within CMIP3 SRESA1B ensemble",
               draw.seperate.legend = TRUE)
## End(Not run)