Fix missing value handling in moderators

jasp-stats · Nov 8, 2023 · 3aff05e · 3aff05e
1 parent 23e9834
commit 3aff05e
Show file tree

Hide file tree

Showing 4 changed files with 251 additions and 1 deletion.
diff --git a/R/classicProcess.R b/R/classicProcess.R
@@ -674,8 +674,14 @@ ClassicProcess <- function(jaspResults, dataset = NULL, options) {
       # Return matrix with dummy coding for each factor
       return(conMat[, colIdx])
     }
+    # If listwise missing value deletion, only use complete cases for quantiles
+    if (options$naAction == "listwise") {
+      isComplete <- complete.cases(dataset)
+    } else {
+      isComplete <- !logical(nrow(dataset))
+    }
     # If not factor return quantiles of continuous moderator
-    return(quantile(dataset[[nms]], probs = probeVals))
+    return(quantile(dataset[[nms]][isComplete], probs = probeVals, na.rm = TRUE))
   })
 
   names(modProbes) <- modVars

diff --git a/...hat/_snaps/classic-process-integration-general/conceptual-path-plot-missing.svg b/...hat/_snaps/classic-process-integration-general/conceptual-path-plot-missing.svg