Fix the logic of the check

LieberInstitute · Jan 9, 2024 · cba82cf · cba82cf
1 parent b6bb2ed
commit cba82cf
Showing 1 changed file with 5 additions and 0 deletions.
diff --git a/R/DEqual.R b/R/DEqual.R
@@ -54,6 +54,11 @@ DEqual <- function(DE) {
         stop("Error: Row names of DE are null.")
       }
 
+    # Check if any rowname in degradation_tstats is in rownames DE
+    if (!any(rownames(DE) %in% rownames(qsvaR::degradation_tstats) | rownames(DE) %in% (gsub('\\..*', '', rownames(qsvaR::degradation_tstats))))) {
+      stop("The rownames of the DE output and the rownames of the  degradation t-statistic do not match")
+    }
+
     ## Locate common transcripts
     common <- intersect(rownames(qsvaR::degradation_tstats), rownames(DE))
     #stopifnot(length(common) > 0)