check length(xdat) == length(ds) in gevdfit

dc65bf7c · Laura Mack · 5a1b0e3a · dc65bf7c
Commit dc65bf7c authored Aug 04, 2020 by Laura Mack
Show whitespace changes
Inline Side-by-side

Showing with 3 additions and 1 deletion

R/gevdfit.R R/gevdfit.R +3 -1

No files found.
--- a/R/gevdfit.R
+++ b/R/gevdfit.R
@@ -77,7 +77,9 @@ gev.d.fit<-
           init.vals = NULL, theta_zero = FALSE,
           show = TRUE, method = "Nelder-Mead", maxit = 10000, ...)
  {
-    
+    if (length(xdat) != length(ds)) {
+      stop(paste0('The length of xdat is ',length(xdat),', but the length of ds is ',length(ds),'.'))
+    } 
    z <- list()
    # number of parameters (betas) to estimate for each parameter: 
    npmu <- length(mul) + 1