# Changes between Version 11 and Version 12 of udg/ecoms/RPackage/biascorrection

Ignore:
Timestamp:
May 13, 2016 3:07:39 PM (6 years ago)
Comment:

--

### Legend:

Unmodified
 v11 = Bias correction of seasonal forecasting data = ------------ //**NOTE**: The following examples may take more than 30 min of computation. They have been conceived to provide an overview of the loadeR.ECOMS capabilities in combination with the companion tools for climate data analysis and post-processing. For a quicker overview, consider reducing the spatial/temporal domain of the initial request. Faster bias correction examples applied to seasonal forecast predictions are available in the [https://github.com/SantanderMetGroup/downscaleR/wiki/Application-to-Seasonal-Forecasts downscaleR's wiki]// -------------- The seasonal forecasting data obtained from the loadeR.ECOMS package can be easily bias corrected (and downscaled) using the  [https://github.com/SantanderMetGroup/downscaleR/wiki downscaleR] package (see a description of the [https://github.com/SantanderMetGroup/downscaleR/wiki/Bias-Correction-and-Model-Output-Statistics-(MOS) bias correction functions]). This package has been developed in the framework of the SPECS and EUPORIAS projects  for bias correction and downscaling of daily climate model outputs (with special focus in seasonal forecasting). Note that, in order to facilitate the use of the resulting bias corrected data in different impact applications, the resulting bias corrected data can be easily exported to NetCDF format using the [https://github.com/SantanderMetGroup/loadeR.2nc loadeR.2nc] extension. ||= R code =||= Output =|| method <- "eqm" # Empirical quantile mapping pr.threshold <- 1 # The minimum value that is considered as a non-zero precipitation. window <- c(30,7) # Integer vector specifying the time window width and the time step used to calibrate. The window is centered on the target time frame. Default to \code{NULL}, which considers the whole period available. window <- c(30,10) # Integer vector specifying the time window width and the time step used to calibrate. The window is centered on the target time frame. Default to \code{NULL}, which considers the whole period available. }}} variables <- c("ps","wss","huss","tas","tasmax","tasmin","tp","rsds","rlds") for (v in 1:length(variables)){ obs <- NULL prd <- NULL obs <- loadECOMS(dataset = "WFDEI", var = variables[v],