ErrorCrit_KGE2.R

ErrorCrit_KGE2 <- function(InputsCrit, OutputsModel, warnings = TRUE, verbose = TRUE) {
  
  
    ##Arguments_check________________________________
    if (inherits(InputsCrit, "InputsCrit") == FALSE) {
      stop("InputsCrit must be of class 'InputsCrit' \n")
      return(NULL)
    }
    if (inherits(OutputsModel, "OutputsModel") == FALSE) {
      stop("OutputsModel must be of class 'OutputsModel' \n")
      return(NULL)
    }
    
    
    ##Initialisation_________________________________
    CritName <- NA
    if (InputsCrit$transfo == "") {
      CritName <- "KGE'[Q]"
    }
    if (InputsCrit$transfo == "sqrt") {
      CritName <- "KGE'[sqrt(Q)]"
    }
    if (InputsCrit$transfo == "log") {
      CritName <- "KGE'[log(Q)]"
    }
    if (InputsCrit$transfo == "inv") {
      CritName <- "KGE'[1/Q]"
    }
    if (InputsCrit$transfo == "sort") {
      CritName <- "KGE'[sort(Q)]"
    }
    CritValue       <- NA
    CritBestValue   <- +1
    Multiplier      <- -1
    ### must be equal to -1 or +1 only
    
    
    ##Data_preparation_______________________________
    VarObs <- InputsCrit$Qobs
    VarObs[!InputsCrit$BoolCrit] <- NA
    VarSim <- OutputsModel$Qsim
    VarSim[!InputsCrit$BoolCrit] <- NA
    
    ##Data_transformation
    if ("Ind_zeroes" %in% names(InputsCrit) & "epsilon" %in% names(InputsCrit)) {
      if (length(InputsCrit$Ind_zeroes) > 0) {
        VarObs <- VarObs + InputsCrit$epsilon
        VarSim <- VarSim + InputsCrit$epsilon
      }
    }
    if (InputsCrit$transfo == "sqrt") {
      VarObs <- sqrt(VarObs)
      VarSim <- sqrt(VarSim)
    }
    if (InputsCrit$transfo == "log") {
      VarObs <- log(VarObs)
      VarSim <- log(VarSim)
      VarSim[VarSim      < -1e100] <- NA
    }
    if (InputsCrit$transfo == "inv") {
      VarObs <- 1 / VarObs
      VarSim <- 1 / VarSim
      VarSim[abs(VarSim) > 1e+100] <- NA
    }
    if (InputsCrit$transfo == "sort") {
      VarSim[is.na(VarObs)] <- NA
      VarSim <- sort(VarSim, na.last = TRUE)
      VarObs <- sort(VarObs, na.last = TRUE)
      InputsCrit$BoolCrit <- sort(InputsCrit$BoolCrit, decreasing = TRUE)
    }