Nothing
#' Cumulative uncertainty(DEPRECATED)
#'
#' This function performs uncertainty decomposition based on the cumulative uncertainty.
#' @param data a data frame containing scenarios(factor or character) for each stages and the variable of interest(numeric).
#' data should contain all combinations of scenarios.
#' @param var_name the name of the variable of interest
#' @param stages names of the stages in the modeling chain. should be ordered by the order of the modeling chain
#' @param U a function that returns uncertainty such as range and variance of a given numeric vector.
#' This package have built-in uncertainty functions var0() and drange(). Default is var0().
#' @return summary of uncertainties
#' @import stats
#' @export
#' @examples
#' stage1 <- LETTERS[1:3]
#' stage2 <- LETTERS[1:2]
#' stage3 <- LETTERS[1:4]
#' y <- rnorm(3*2*4)
#' data <- expand.grid(stage1=stage1,
#' stage2=stage2,
#' stage3=stage3)
#' data <- cbind(data, y)
#' # cum_uncertainty() is deprecated. Use UD_cum_stage()
#' # cum_uncertainty(data,"y", names(data)[-4])
#' # cum_uncertainty(data,"y", names(data)[-4],drange)
cum_uncertainty<-function(data, var_name,
stages=setdiff(names(data),var_name),
U=var0){
.Deprecated(old = "cum_uncertainty", new = "UD_cum_stage")
nstage <- length(stages)
data[stages] <- lapply(data[stages], as.factor)
cu <- rep(0, nstage)
for(k in 1:(nstage-1)) {
given_stages <- paste(stages[(k+1):nstage], collapse="+")
formula0 <- as.formula(paste(var_name, "~" , given_stages, sep=""))
agg <- aggregate(formula0, data=data, U)
cu[k] <- mean(agg[[ var_name ]])
}
cu[nstage] <- U(data[[ var_name ]])
u_cum = c(cu[1], diff(cu))
prop = u_cum/cu[nstage]
result <- data.frame(stage=stages,
cum_uncer=cu,
stage_uncer=u_cum,
Prop=prop)
return(result)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.