smooth_scale | R Documentation |
Computes mean de-biasing with loess smoothing and adjusts variance. smooth_scale
scales the ensemble deviations from the climatology, whereas smooth_scalespread
only scales the ensemble anomalies (change to the ensemble spread).
smooth_scale(fcst, obs, fcst.out = fcst, span = min(1, 31/nrow(fcst)), ...)
smooth_scalespread(
fcst,
obs,
fcst.out = fcst,
span = min(1, 31/nrow(fcst)),
...
)
fcst |
n x m x k array of n lead times, m forecasts, of k ensemble members |
obs |
n x m matrix of veryfing observations |
fcst.out |
array of forecast values to which bias correction
should be applied (defaults to |
span |
the parameter which controls the degree of smoothing (see |
... |
additional arguments for compatibility with other bias correction methods |
The standard deviation (1/n definition) of the observations and simulations is computed against the smoothed climatology in order to get consistent results should the forecast ensemble collapse on the observations. The disadvantage of this approach is the dependence of the scaling on the climatological fit.
## initialise forcast observation pairs
fcst <- array(rnorm(215*30*51, mean=3, sd=0.2), c(215, 30, 51)) +
0.5*sin(seq(0,4,length=215))
obs <- array(rnorm(215*30, mean=2), c(215, 30)) +
sin(seq(0,4, length=215))
fcst.debias <- biascorrection:::smooth_scale(fcst[,1:20,], obs[,1:20], fcst.out=fcst, span=0.5)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.