Nothing
qtimedelay <- function(timedelay, deltim, size_frac_future) {
# Computes the fraction of runoff in future time steps
# Author: Claudia Vitolo
#
# Args:
# timedelay: mparam$timedelay
# deltim: time step in days
# size_frac_future: maximum number of future time steps
#
# Returns:
# Fraction of runoff in future time steps.
frac_future <- rep(0, size_frac_future)
alpha <- 2.5 # shape parameter
alamb <- alpha/timedelay # quantile (in the original fortran code this is defined scale parameter, why?)
psave <- 0.0 # cumulative probability at jtim-1
# loop through time steps and compute the fraction of runoff in future time steps
for (jtim in seq(along = frac_future)) {
tfuture <- jtim * deltim # future time (units of days)
#cumprob <- pgamma(alamb*tfuture,alpha) # cumulative probability at jtim (original formula in fortran)
cumprob <- pgamma(alamb*tfuture,shape=alpha,scale=1) # cumulative probability at jtim
frac_future[jtim] <- max(0, cumprob-psave) # probability between jtim-1 and jtim
psave <- cumprob # cumulative probability at jtim-1
}
# ensure that the fractions sum to 1.0 (account for rounding errors, and not enough bins)
frac_future <- frac_future / sum(frac_future)
return(frac_future)
}
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.