added degrees of freedom to scLANE output -- closes #267

jr-leary7 · Oct 30, 2024 · 712c7f8 · 712c7f8
1 parent 507a552
commit 712c7f8
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 1 deletion.
diff --git a/R/getResultsDE.R b/R/getResultsDE.R
@@ -40,7 +40,7 @@ getResultsDE <- function(test.dyn.res = NULL,
                                      function(x) {
                                        purrr::map_dfr(x,
                                                       function(y) {
-                                                        as.data.frame(rbind(y[seq(14)])) %>%
+                                                        as.data.frame(rbind(y[seq(15)])) %>%
                                                           dplyr::mutate(dplyr::across(tidyselect::everything(), \(z) unname(unlist(z))))
                                                       })
                                      }) %>%

diff --git a/R/testDynamic.R b/R/testDynamic.R
@@ -324,6 +324,7 @@ testDynamic <- function(expr.mat = NULL,
                                                        ifelse(gee.test == "wald", "Wald", "Score"), 
                                                        "LRT"),
                                Test_Stat_Note = NA_character_,
+                               Degrees_Freedom = NA_real_, 
                                P_Val = NA_real_,
                                LogLik_MARGE = marge_sumy$ll_marge,
                                LogLik_Null = null_sumy$null_ll,
@@ -366,6 +367,7 @@ testDynamic <- function(expr.mat = NULL,
                                            ifelse(gee.test == "wald", test_res$Wald_Stat, test_res$Score_Stat), 
                                            test_res$LRT_Stat)
      lineage_list[[j]]$Test_Stat_Note <- test_res$Notes
+     lineage_list[[j]]$Degrees_Freedom <- test_res$DF 
      lineage_list[[j]]$P_Val <- test_res$P_Val
     }
     names(lineage_list) <- paste0("Lineage_", LETTERS[seq(n_lineages)])
@@ -393,6 +395,7 @@ testDynamic <- function(expr.mat = NULL,
                                      ifelse(gee.test == "wald", "Wald", "Score"), 
                                      "LRT"),
              Test_Stat_Note = NA_character_,
+             Degrees_Freedom = NA_real_, 
              P_Val = NA_real_,
              LogLik_MARGE = NA_real_,
              LogLik_Null = NA_real_,