Plot Two Scores With Confidence Intervals In A Common Plot
Plots two input variables having the same dimensions in a common plot.
One plot for all experiments.
Input variables should have dimensions (nexp/nmod, nltime).
Plot2VarsVsLTime(
var1,
var2,
toptitle = "",
ytitle = "",
monini = 1,
freq = 12,
nticks = NULL,
limits = NULL,
listexp = c("exp1", "exp2", "exp3"),
listvars = c("var1", "var2"),
biglab = FALSE,
hlines = NULL,
leg = TRUE,
siglev = FALSE,
sizetit = 1,
show_conf = TRUE,
fileout = "output_plot2varsvsltime.eps",
width = 8,
height = 5,
size_units = "in",
res = 100,
...
)var1 |
Matrix of dimensions (nexp/nmod, nltime). |
var2 |
Matrix of dimensions (nexp/nmod, nltime). |
toptitle |
Main title, optional. |
ytitle |
Title of Y-axis, optional. |
monini |
Starting month between 1 and 12. Default = 1. |
freq |
1 = yearly, 12 = monthly, 4 = seasonal, ... Default = 12. |
nticks |
Number of ticks and labels on the x-axis, optional. |
limits |
c(lower limit, upper limit): limits of the Y-axis, optional. |
listexp |
List of experiment names, up to three, optional. |
listvars |
List of names of input variables, optional. |
biglab |
TRUE/FALSE for presentation/paper plot. Default = FALSE. |
hlines |
c(a, b, ...) Add horizontal black lines at Y-positions a, b,
... |
leg |
TRUE/FALSE if legend should be added or not to the plot. Default = TRUE. |
siglev |
TRUE/FALSE if significance level should replace confidence
interval. |
sizetit |
Multiplicative factor to change title size, optional. |
show_conf |
TRUE/FALSE to show/not confidence intervals for input variables. |
fileout |
Name of output file. Extensions allowed: eps/ps, jpeg, png,
pdf, bmp and tiff. |
width |
File width, in the units specified in the parameter size_units (inches by default). Takes 8 by default. |
height |
File height, in the units specified in the parameter size_units (inches by default). Takes 5 by default. |
size_units |
Units of the size of the device (file or window) to plot in. Inches ('in') by default. See ?Devices and the creator function of the corresponding device. |
res |
Resolution of the device (file or window) to plot in. See ?Devices and the creator function of the corresponding device. |
... |
Arguments to be passed to the method. Only accepts the following
graphical parameters: |
Examples of input:
——————
RMSE error for a number of experiments and along lead-time: (nexp, nltime)
History:
1.0 - 2013-03 (I. Andreu-Burillo, isabel.andreu-burillo@ic3.cat)
- Original code
# Load sample data as in Load() example:
example(Load)
clim <- Clim(sampleData$mod, sampleData$obs)
ano_exp <- Ano(sampleData$mod, clim$clim_exp)
ano_obs <- Ano(sampleData$obs, clim$clim_obs)
runmean_months <- 12
dim_to_smooth <- 4 # Smooth along lead-times
smooth_ano_exp <- Smoothing(ano_exp, runmean_months, dim_to_smooth)
smooth_ano_obs <- Smoothing(ano_obs, runmean_months, dim_to_smooth)
dim_to_mean <- 2 # Mean along members
required_complete_row <- 3 # Discard start dates that contain NA along lead-times
leadtimes_per_startdate <- 60
rms <- RMS(Mean1Dim(smooth_ano_exp, dim_to_mean),
Mean1Dim(smooth_ano_obs, dim_to_mean),
compROW = required_complete_row,
limits = c(ceiling((runmean_months + 1) / 2),
leadtimes_per_startdate - floor(runmean_months / 2)))
smooth_ano_exp_m_sub <- smooth_ano_exp - InsertDim(Mean1Dim(smooth_ano_exp, 2,
narm = TRUE), 2, dim(smooth_ano_exp)[2])
spread <- Spread(smooth_ano_exp_m_sub, c(2, 3))
Plot2VarsVsLTime(InsertDim(rms[, , , ], 1, 1), spread$sd,
toptitle = 'RMSE and spread', monini = 11, freq = 12,
listexp = c('CMIP5 IC3'), listvar = c('RMSE', 'spread'),
fileout = 'plot2vars.eps')Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.