Skip to content
GitLab
Select Git revision
  • develop
  • master default
Find file BlameHistoryPermalink
CreateInputsCrit.R 3.40 KiB
1 
CreateInputsCrit <-
2 
  function(FUN_CRIT,
3 
           InputsModel,
4 
           RunOptions,
5 
           Qobs,
6 
           BoolCrit = NULL,
7 
           transfo = "",
8 
           Ind_zeroes = NULL,
9 
           epsilon = NULL) {
10 
    ObjectClass <- NULL
11
12
13 
    ##check_FUN_CRIT
14 
    BOOL <- FALSE
15
16 
    if (identical(FUN_CRIT, ErrorCrit_NSE) | identical(FUN_CRIT, ErrorCrit_KGE) |
17 
        identical(FUN_CRIT, ErrorCrit_KGE2) | identical(FUN_CRIT, ErrorCrit_RMSE)) {
18 
      BOOL <- TRUE
19 
    }
20 
    if (!BOOL) {
21 
      stop("incorrect FUN_CRIT for use in CreateInputsCrit \n")
22 
      return(NULL)
23 
    }
24
25 
    ##check_arguments
26 
    if (inherits(InputsModel, "InputsModel") == FALSE) {
27 
      stop("InputsModel must be of class 'InputsModel' \n")
28 
      return(NULL)
29 
    }
30 
    if (inherits(RunOptions , "RunOptions") == FALSE) {
31 
      stop("RunOptions must be of class 'RunOptions' \n")
32 
      return(NULL)
33 
    }
34
35 
    LLL <- length(InputsModel$DatesR[RunOptions$IndPeriod_Run])
36
37 
    if (is.null(Qobs)) {
38 
      stop("Qobs is missing \n")
39 
      return(NULL)
40 
    }
41 
    if (!is.vector(Qobs)) {
42 
      stop(paste("Qobs must be a vector of numeric values \n", sep = ""))
43 
      return(NULL)
44 
    }
45 
    if (!is.numeric(Qobs)) {
46 
      stop(paste("Qobs must be a vector of numeric values \n", sep = ""))
47 
      return(NULL)
48 
    }
49 
    if (length(Qobs) != LLL) {
50 
      stop("Qobs and InputsModel series must have the same length \n")
51 
      return(NULL)
52 
    }
53 
    if (is.null(BoolCrit)) {
54 
      BoolCrit <- rep(TRUE, length(Qobs))
55 
    }
56 
    if (!is.logical(BoolCrit)) {
57 
      stop("BoolCrit must be a vector of boolean \n")
58 
      return(NULL)
59 
    }
60 
    if (length(BoolCrit) != LLL) {
61 
      stop("BoolCrit and InputsModel series must have the same length \n")
62 
      return(NULL)
63 
    }
64 
    if (is.null(transfo)) {
65 
      stop("transfo must be a chosen among the following: '', 'sqrt', 'log' or 'inv' or 'sort' \n")
66 
      return(NULL)
67 
    }
68 
    if (!is.vector(transfo)) {
69 
      stop("transfo must be a chosen among the following: '', 'sqrt', 'log' or 'inv' or 'sort' \n")
70 
      return(NULL)
7172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119
} if (length(transfo) != 1) { stop("transfo must be a chosen among the following: '', 'sqrt', 'log' or 'inv' or 'sort' \n") return(NULL) } if (!is.character(transfo)) { stop("transfo must be a chosen among the following: '', 'sqrt', 'log' or 'inv' or 'sort' \n") return(NULL) } if (transfo %in% c("", "sqrt", "log", "inv", "sort") == FALSE) { stop("transfo must be a chosen among the following: '', 'sqrt', 'log' or 'inv' or 'sort' \n") return(NULL) } if (!is.null(Ind_zeroes)) { if (!is.vector(Ind_zeroes)) { stop("Ind_zeroes must be a vector of integers \n") return(NULL) } if (!is.integer(Ind_zeroes)) { stop("Ind_zeroes must be a vector of integers \n") return(NULL) } } if (!is.null(epsilon)) { if (!is.vector(epsilon) | length(epsilon) != 1 | !is.numeric(epsilon)) { stop("epsilon must be single numeric value \n") return(NULL) } epsilon = as.double(epsilon) } ##Create_InputsCrit InputsCrit <- list(BoolCrit = BoolCrit, Qobs = Qobs, transfo = transfo, Ind_zeroes = Ind_zeroes, epsilon = epsilon) class(InputsCrit) <- c("InputsCrit", ObjectClass) return(InputsCrit) }