Fixed issue for cases where there are no features retained after filt…

…ering
biobakery · Apr 28, 2023 · c07170d · c07170d
1 parent ba54e9f
commit c07170d
Show file tree

Hide file tree

Showing 2 changed files with 28 additions and 0 deletions.
diff --git a/NEWS b/NEWS
@@ -1,6 +1,9 @@
 
 # MTXmodel History #
 
+## Changes in v1.2.3 (2023-04-28) ##
+* Fixed issue if no features left after filtering
+
 ## Changes in v1.2.2 (2022-11-14) ##
 * Fixed warning on if conditions
 * Fixed issue for sample matching

diff --git a/R/MTXmodel.R b/R/MTXmodel.R
@@ -1006,6 +1006,31 @@ MTXmodel <-
             }
         }
 
+		if (ncol(filtered_data_norm) < 1) { 
+            all_outfile <- file.path(output, "all_results.tsv")
+            mycolnames = c( 
+                "feature",
+                "metadata",
+                "value",
+                "coef",
+                "stderr",
+                "N", 
+                "N.not.0",
+                "pval",
+                "qval")
+            write.table(
+                t(mycolnames),
+                file = all_outfile,
+                sep = "\t",
+                quote = FALSE,
+                col.names = FALSE,
+                row.names = FALSE)
+
+            residuals_file = file.path(output, "residuals.rds")
+            file.create(residuals_file)
+            return(NULL)
+        }		
+
         # apply the method to the data with the correction
         logging::loginfo(
             "Running selected analysis method: %s", analysis_method)