# NOT RUN {
# Example1: basic ideal case
obs <- 1:10
sim <- 1:10
KGE(sim, obs)
obs <- 1:10
sim <- 2:11
KGE(sim, obs)
##################
# Example2: Looking at the difference between 'method=2009' and 'method=2012'
# Loading daily streamflows of the Ega River (Spain), from 1961 to 1970
data(EgaEnEstellaQts)
obs <- EgaEnEstellaQts
# Simulated daily time series, initially equal to twice the observed values
sim <- 2*obs
# KGE 2009
KGE(sim=sim, obs=obs, method="2009", out.type="full")
# KGE 2012
KGE(sim=sim, obs=obs, method="2012", out.type="full")
##################
# Example3: KGE for simulated values equal to observations plus random noise
# on the first half of the observed values
# Randomly changing the first 1826 elements of 'sim', by using a normal distribution
# with mean 10 and standard deviation equal to 1 (default of 'rnorm').
sim <- obs
sim[1:1826] <- obs[1:1826] + rnorm(1826, mean=10)
# Computing the new 'KGE'
KGE(sim=sim, obs=obs)
# Randomly changing the first 2000 elements of 'sim', by using a normal distribution
# with mean 10 and standard deviation equal to 1 (default of 'rnorm').
sim[1:2000] <- obs[1:2000] + rnorm(2000, mean=10)
# Computing the new 'KGE'
KGE(sim=sim, obs=obs)
# }
Run the code above in your browser using DataLab