bug in xvErrorMeasures for simulated output corrected

DESCRIPTION

 Package: gmGeostats
-Version: 0.10-7.9004
+Version: 0.10-7.9005
 Date: 2020-10-05
 Title: Geostatistics for Compositional Analysis
 Authors@R: c(person(given = "Raimon", 

NEWS.md

+# gmGeostats 0.10.7.9005
+
+* (2021-06-30) bugs in xvErrorMeasures() for simulated data corrected
+* (2021-06-30) bugs in accuracy() and precision() corrected
+
 # gmGeostats 0.10.7.9004
 
 * (2021-06-29) bugs in accuracy() and precision() corrected

R/accuracy.R

@@ -387,7 +387,17 @@ xvErrorMeasures.DataFrameStack = function(x, observed, output="ME",
   
   outputs = c("ME","MSE")
   output = outputs[pmatch(output, outputs, duplicates.ok = TRUE)]
-  mn = gmApply(x, FUN=colMeans)
-  class(mn) = "data.frame"
-  xvErrorMeasures(mn, observed, output, univariate)
+  dims = noStackDim(x)
+  dims = c(ifelse(is.numeric(dims), 1, "loc"), dims)
+  mn = gmApply(x, MARGIN=dims, FUN=mean)
+  # class(mn) = "data.frame"
+  # xvErrorMeasures(mn, observed, output, univariate=TRUE)
+  
+  resids = mn-observed
+  myfun = function(output){
+    if(output=="ME") return(colMeans(resids, na.rm=TRUE))
+    if(output=="MSE") return(mean(rowSums(resids^2, na.rm=T), na.rm=TRUE))
+  }
+  
+  return(sapply(outputs, myfun))
 }