Changes between Version 11 and Version 12 of udg/ecoms/RPackage/biascorrection
 Timestamp:
 May 13, 2016 3:07:39 PM (6 years ago)
Legend:
 Unmodified
 Added
 Removed
 Modified

udg/ecoms/RPackage/biascorrection
v11 v12 1 1 = Bias correction of seasonal forecasting data = 2 3  4 //**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 postprocessing. 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/ApplicationtoSeasonalForecasts downscaleR's wiki]// 5  2 6 3 7 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/BiasCorrectionandModelOutputStatistics(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). … … 6 10 7 11 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. 12 8 13 9 14 = R code == Output = … … 91 96 method < "eqm" # Empirical quantile mapping 92 97 pr.threshold < 1 # The minimum value that is considered as a nonzero precipitation. 93 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.98 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. 94 99 }}} 95 100 … … 101 106 variables < c("ps","wss","huss","tas","tasmax","tasmin","tp","rsds","rlds") 102 107 for (v in 1:length(variables)){ 103 obs < NULL104 prd < NULL105 108 obs < loadECOMS(dataset = "WFDEI", 106 109 var = variables[v],