From abc2bb4e4e0f20adb6dcf416999c9268299cc538 Mon Sep 17 00:00:00 2001 From: nguyenta Date: Thu, 8 Feb 2024 13:25:47 +0100 Subject: [PATCH] update vignettes --- vignettes/run_SWATplus_with_given_parametersets.R | 2 +- vignettes/run_SWATplus_with_given_parametersets.Rmd | 2 +- vignettes/run_SWATplus_with_given_parametersets.html | 2 +- 3 files changed, 3 insertions(+), 3 deletions(-) diff --git a/vignettes/run_SWATplus_with_given_parametersets.R b/vignettes/run_SWATplus_with_given_parametersets.R index 2ae4f64..deffa7c 100644 --- a/vignettes/run_SWATplus_with_given_parametersets.R +++ b/vignettes/run_SWATplus_with_given_parametersets.R @@ -116,7 +116,7 @@ knitr::opts_chunk$set( # # Now calculate the 95PPU (SUFI2 approach) for each output variables, example for variable 1 # # TODO: change variable and percentiles as your wish # variable <- 1 -# percentiles <- c(0.0025, 0.5, 0.975) +# percentiles <- c(0.025, 0.5, 0.975) # # for (itime in 1:nrow(output[[variable]])){ # ppu <- as.numeric(quantile(output[[variable]][itime,],percentiles)) diff --git a/vignettes/run_SWATplus_with_given_parametersets.Rmd b/vignettes/run_SWATplus_with_given_parametersets.Rmd index 5eec623..f101ffd 100644 --- a/vignettes/run_SWATplus_with_given_parametersets.Rmd +++ b/vignettes/run_SWATplus_with_given_parametersets.Rmd @@ -150,7 +150,7 @@ output <- readAllOutVar(workingFolder, numTimesteps) # Now calculate the 95PPU (SUFI2 approach) for each output variables, example for variable 1 # TODO: change variable and percentiles as your wish variable <- 1 -percentiles <- c(0.0025, 0.5, 0.975) +percentiles <- c(0.025, 0.5, 0.975) for (itime in 1:nrow(output[[variable]])){ ppu <- as.numeric(quantile(output[[variable]][itime,],percentiles)) diff --git a/vignettes/run_SWATplus_with_given_parametersets.html b/vignettes/run_SWATplus_with_given_parametersets.html index 0e6cf08..77a5f43 100644 --- a/vignettes/run_SWATplus_with_given_parametersets.html +++ b/vignettes/run_SWATplus_with_given_parametersets.html @@ -474,7 +474,7 @@

4. Processing model output

# Now calculate the 95PPU (SUFI2 approach) for each output variables, example for variable 1 # TODO: change variable and percentiles as your wish variable <- 1 -percentiles <- c(0.0025, 0.5, 0.975) +percentiles <- c(0.025, 0.5, 0.975) for (itime in 1:nrow(output[[variable]])){ ppu <- as.numeric(quantile(output[[variable]][itime,],percentiles))