fixed post processing issue

scripts/post_process_imputations.R

@@ -12,7 +12,7 @@ df_gose              <- readRDS(gosefile)
 
 posteriors <- list.files(modelimputations_dir, pattern = "*.rds")
 
-# load and combine all posteriros over different (baseline) imputations
+# load and combine all posteriors over different (baseline) imputations
 df_imputations <- c()
 for (i in 1:length(posteriors)) {
   df_imputations <- rbind(
@@ -40,11 +40,12 @@ df_imputations <- df_imputations %>%
       {
         if (length(.) == 0) { # fix all probs. NA case
           return(NA)
-        } else {
-          return(.)
         }
+        if (length(.) > 1) { # need to resolve ambiguous gose predictions !
+          return(sample(., 1))
+        }
+        return(.)
       } %>%
-      sample(1) %>% # need to resolve ambiguous gose predictions !
       rep(n()),
     GOSE = factor(
       GOSE,