WACSvalid {WACS}

Performs validations of WACS simulations

Description

The validation is based on different types of statistics computed on WACS data, WACS parameters and WACS simulations.

Usage

WACSvalid(
  what = "Sim",
  wacsdata = NULL,
  wacspar = NULL,
  wacssimul = NULL,
  varname = NULL,
  varname2 = NULL,
  base = 0,
  above = TRUE,
  months = 1:12
)

Arguments

what

Type of validation. Possible choices are:

="Sim"	Compares a simulation run to data
="Rain"	qq-plots of rainfall, per season
="MeanSd"	Compares monthly mean and standard deviations
="BiVar"	Compares monthly variate correlations
="CorTemp"	Compares monthly temporal correlations
="SumBase"	Compares sums above a threshold
="Persistence"	Compares persistence of a variable above (or below) a threshold

wacsdata

WACS data obtained when calling WACSdata

wacspar

WACS parameters estimated when calling WACSestim

wacssimul

WACS simulation obtained when calling WACSsimul

varname

Variable on which the validation is performed

varname2

Second variable on which validation is performed (only needed if what=BiVar)

base

Threshold used for "SumBase" and "Persistence"

above

Boolean value used for "Persistence": TRUE if data is considered above threshold; FALSE otherwise

months

Months to which the analysis should be restricted (only for "SumBase" and "Persistence")

Value

A list containing all information needed for plots; contains also the type of validation, as a class

Note

If what=sim, data and simulations are displayed as a function of the day of the year, from 1 to 365. Smoothed versions of daily average and daily envelopes (defined by average +/- 2. standard deviations) are also displayed.

If what=rain, qq-plots and superimposition of histograms and models of rain are produced for each season.

If what=MeanSd, boxplots of monthly means and monthly standard deviations are compared. The median value of the monthly mean, resp. monthly standard deviation, of the data are displayed on top of the boxplots computed on the simulations.

#' If what=BiVar, boxplots of monthly correlations coefficients between varname and varname2 are compared. The median value of the correlation coefficient computed on the data is displayed on top of the boxplots computed on the simulations.

If what=CorTemp, boxplots of monthly 1-day auto correlation are compared. The median value of the auto-correlation coefficient computed on the data is displayed on top of the boxplots computed on the simulations.

If what=SumBase, boxplots of the sum of the variable varname above a given threshold, base, is computed during the months provided in the variable months.

If what=Persistence, histograms of consecutive days of the variable varname above (or below) a given threshold, base, are compared. If above=TRUE, consecutive days above the threshold are computed, whereas days below the threshold are computed if above=FALSE. Months can be selected with the variable months.

Examples

## Not run: 
  ## Simple example
  data(ClimateSeries)
  ThisData = WACSdata(ClimateSeries)
  ThisPar  = WACSestim(ThisData)
  ThisSim  = WACSsimul(ThisPar, from="1995-01-01", to="2012-12-31")
  Val1 = WACSvalid(what="Sim",wacsdata = ThisData, 
                   wacspar = ThisPar, wacssimul = ThisSim, varname="tmin")
  Val2 = WACSvalid(what="MeanSd",wacsdata = ThisData, 
                   wacssimul = ThisSim, varname="RG")
  Val3 = WACSvalid(what="SumBase", wacsdata = ThisData,  
                   wacssimul = ThisSim, varname="tmoy", base=5, month=2:5)
  Val4 = WACSvalid(what="Persistence",wacsdata = ThisData,  
                   wacssimul = ThisSim, varname="tmin", base=0, above=FALSE)

## End(Not run)