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


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

--

Legend:

Unmodified
Added
Removed
Modified
  • udg/ecoms/RPackage/biascorrection

    v11 v12  
    11= 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 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]//
     5--------------
    26
    37The 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).
     
    610
    711Note 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
    813
    914||= R code =||= Output =||
     
    9196method <- "eqm" # Empirical quantile mapping
    9297pr.threshold <- 1 # The minimum value that is considered as a non-zero 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.
     98window <- 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.
    9499}}}
    95100
     
    101106variables <- c("ps","wss","huss","tas","tasmax","tasmin","tp","rsds","rlds")
    102107for (v in 1:length(variables)){
    103    obs <- NULL
    104    prd <- NULL
    105108   obs <- loadECOMS(dataset = "WFDEI",
    106109                    var = variables[v],