Newer
Older
% Please edit documentation in R/CST_BiasCorrection.R
\name{CST_BiasCorrection}
\alias{CST_BiasCorrection}
\title{Bias Correction of a CSTools object based on the mean and standard deviation adjustment following Torralba et al. (2017) }
\item{data}{a CSTools object (an s2dverification object as output by the \code{Load} function from the s2dverification package).}
\code{$mod} {a CSTools object (s2dverification object) with the bias corrected forecasts (provided in $mod) with the same dimensions as data$mod.}
}
\description{
This function applies the simple bias adjustment technique described in Torralba et al. (2017). The adjusted forecasts have an equivalent standard deviation and mean to that of the reference dataset.
}
\examples{
# Creation of sample s2dverification objects. These are not complete
# s2dverification objects though. The Load function returns complete objects.
# Example
mod1 <- 1 : (1 * 3 * 4 * 5 * 6 * 7)
dim(mod1) <- c(dataset = 1, member = 3, sdate = 4, ftime = 5, lat = 6, lon = 7)
obs1 <- 1 : (1 * 1 * 4 * 5 * 6 * 7)
dim(obs1) <- c(dataset = 1, member = 1, sdate = 4, ftime = 5, lat = 6, lon = 7)
lon <- seq(0, 30, 5)
lat <- seq(0, 25, 5)
data1 <- list(mod = mod1, obs = obs1, lat = lat, lon = lon)
str(a)
dim(mod1) <-
c(
dataset = 1,
member = 3,
sdate = 4,
ftime = 5,
lat = 6,
lon = 7
)
dim(obs1) <-
c(
dataset = 1,
member = 1,
sdate = 4,
ftime = 5,
lat = 6,
lon = 7
)
lon <- seq(0, 30, 5)
lat <- seq(0, 25, 5)
data1 <- list(
mod = mod1,
obs = obs1,
lat = lat,
lon = lon
)
mod2 <- mod1
mod2[1, 2, 1, 1, 1, 1] <- NA
data2 <- list(
mod = mod2,
obs = obs1,
lat = lat,
lon = lon
)
obs2 <- obs1
obs2[1, 1, 2, 1, 1, 1] <- NA
data3 <- list(
mod = mod1,
obs = obs2,
lat = lat,
lon = lon
)
data4 <- list(
mod = mod2,
obs = obs2,
lat = lat,
lon = lon
)
lat2 <- lat
lat2[3] <- NA
data5 <- list(
mod = mod1,
obs = obs1,
lat = lat2,
lon = lon
)
lon2 <- lon
lon2[5] <- NA
data6 <- list(
mod = mod1,
obs = obs1,
lat = lat,
lon = lon2
)
data7 <- list(
mod = mod1,
obs = obs1,
lat = lat2,
lon = lon2
)
data8 <- list(
mod = mod1,
obs = obs2,
lat = lat2,
lon = lon2
)
data9 <- list(
mod = mod2,
obs = obs1,
lat = lat2,
lon = lon2
)
data10 <- list(
mod = mod2,
obs = obs2,
lat = lat2,
lon = lon2
)
}
\references{
Torralba, V., Doblas-Reyes, F. J., MacLeod, D., Christel, I., & Davis, M. (2017). Seasonal climate prediction: A new source of information for the management of wind energy resources. Journal of Applied Meteorology and Climatology, 56(5), 1231-1247.
}
\author{
Verónica Torralba, \email{veronica.torralba@bsc.es}
}