Fixing the weighted crps.

atmoswing · Jun 8, 2017 · df490e8 · df490e8
1 parent 317c2b2
commit df490e8
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/R/performanceScore.R b/R/performanceScore.R
@@ -23,12 +23,12 @@ crps <- function(x, x0, a=0.44, b=0.12, w=NA) {
 
   n <- length(x)
 
-  if (is.na(w)) {
+  if (length(w) == 1) {
     x.s <- sort(x)
     r <- 1:n
     Fx <- (r - a) / (n + b)
   } else {
-    assert_that(length(Fx) == n)
+    assert_that(length(w) == n)
     sorted <- sort.int(x, index.return=TRUE)
     x.s <- sorted$x
     w <- w[sorted$ix]