plot\_NparOF.R: added checking for 'beh.thr'

NEWS

@@ -109,6 +109,8 @@ NEWS/ChangeLog for hydroPSO
                               -) argument 'sim' may also be 'integer' (before it had to be 'numeric', which does not allow 'integer' objects !)  
 
         o 'plot_NparOF'     : -) new argument 'beh.thr' to allow selecting only behavioral parameter sets
+
+        o 'params2ecdf'     : -) new arguments 'gofs', 'MinMax' and 'beh.thr' 
                             
         o new function 'pest2hydroPSO' for importing PEST input files into hydroPSO
         o new function 'hydroPSO2pest' for exporting hydroPSO input files to PEST

R/params2ecdf.R

@@ -36,6 +36,11 @@ params2ecdf <- function(params, ...) UseMethod("params2ecdf")
 ################################################################################
 params2ecdf.default <- function(params, 
                                 param.names=NULL,
+
+                                gofs,
+                                MinMax=NULL, 
+                                beh.thr=NA, 
+
                                 weights=NULL,                                                  
                                 byrow=FALSE, 
                                 plot=TRUE,

R/plot_NparOF.R

@@ -59,6 +59,12 @@ plot_NparOF <- function(params,
     # Setting 'MinMax' 
     MinMax <- match.arg(MinMax)
 
+    # Checking 'beh.thr'
+    if ( !is.na(beh.thr) ) {       
+      if ( is.null(gofs) )
+        stop("Missing argument: 'gofs' has to be provided before using 'beh.thr' !!")
+    } # IF end
+
     # Number of parameters that will be analysed
     npar <- length(param.names)