Commit 8e2c2628 authored by Manuela's avatar Manuela

edits to resolve notes

parent 4b36dac5
......@@ -9,14 +9,14 @@ Authors@R: c(person("Manuela", "Brunner", role = c("aut", "cre"),
person("Reinhard", "Furrer", role = c("aut"),
email = "reinhard.furrer@math.uzh.ch",
comment = c(ORCID = "0000-0002-6319-2332")))
Author: Manuela Brunner [aut, cre] (<https://orcid.org/0000-0001-8824-877X>),
Author: Manuela Brunner [aut, cre] (<https://orcid.org/0000-0001-8824-877X>),
Reinhard Furrer [aut] (<https://orcid.org/0000-0002-6319-2332>)
Maintainer: Manuela Brunner <manuela.i.brunner@gmail.com>
Description: Provides a simulation framework to simulate streamflow time series with similar
main characteristics as observed data. These characteristics include the distribution of daily
streamflow values and their temporal correlation as expressed by short- and long-range
dependence. The approach is based on the randomization of the phases of the Fourier
transform or the phases of the wavelet transform. The function prsim() is applicable to singel site simulation and uses the Fourier transform.
transform or the phases of the wavelet transform. The function prsim() is applicable to single site simulation and uses the Fourier transform.
The function prsim.wave() extends the approach to multiple sites and is based on the complex wavelet transform. We further use the flexible four-parameter Kappa distribution, which allows
for the extrapolation to yet unobserved low and high flows. Alternatively, the empirical or any other distribution can be used. A detailed description of
the simulation approach for single sites and an application example can be found
......
File added
......@@ -190,8 +190,8 @@ prsim <- function(data, station_id="Qobs", number_sim=1, win_h_length=15,
data_random <- rCDF(n=length(data_window), theta)
# density_gengam[[d]] <- density(data_gengam)
hist(data_window)
hist(data_random,add=T,col="red")
# hist(data_window)
# hist(data_random,add=T,col="red")
if (tolower(GoFtest)=="ks"){
p_vals[d] <- ks.test(data_window,data_random)$p.value
}
......
......@@ -57,7 +57,7 @@ Manuela Brunner
}
\examples{
data(runoff_multi_sites)
prsim.wave(runoff_multi_sites, "Qobs", 1, suppWarn=TRUE)
prsim.wave(runoff_multi_sites[c(1:2)], "Qobs", 1, suppWarn=TRUE)
# warnings() # as a follow-up to `suppWarn=TRUE`
## Specifying particular CDFs:
......
......@@ -27,7 +27,8 @@ The data is included to illustrate the validation and visualization routines in
\source{
The data has been generated with
\code{prsim.wave(data=runoff_multi_sites, number_sim=10, marginal="kappa", GoFtest = "KS",out_dir="~/PRSim-devel/data/simulations_multi_sites")}
\code{prsim.wave(data=runoff_multi_sites, number_sim=10, marginal="kappa",
GoFtest = "KS",out_dir="~/PRSim-devel/data/simulations_multi_sites")}
(default values for all other arguments).
}
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment