#' @name conditional
#' @aliases conditional
#' @aliases conditionalRecal
#'
#' @title
#' Bias Conditional on Forecast
#'
#' @description
#' Compute calibration for biases that are conditional on ensemble mean forecast
#'
#' @param ... arguments passed to \code{\link{linmod}}
#'
#' @seealso linmod
#'
#' @examples
#' ## initialise forcast observation pairs
#' seasonal <- sin(seq(0,4,length=215))
#' signal <- outer(seasonal, rnorm(30), '+')
#' fcst <- array(rnorm(215*30*51), c(215, 30, 15)) +
#' 2*c(signal)
#' obs <- array(rnorm(215*30, mean=2), c(215, 30)) +
#' signal
#' fc.time <- outer(1:215, 1981:2010, function(x,y) as.Date(paste0(y, '-11-01')) - 1 + x)
#' fcst.debias <- biascorrection:::conditional(fcst[,1:20,],
#' obs[,1:20], fcst.out=fcst, fc.time=fc.time[,1:20], fcout.time=fc.time, span=0.5)
#'
#' @rdname conditional
#' @keywords util
conditional <- function(...){
return(linmod(..., formula=obs ~ fcst + fcst:poly(lead,3) + fcst:exp(-lead/5), recal=FALSE))
}
#' @rdname conditional
conditionalRecal <- function(...){
return(linmod(..., formula=obs ~ fcst + fcst:poly(lead,3) + fcst:exp(-lead/5), recal=TRUE))
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.