Merge pull request #62 from alyosama/package

Fixing F-distribution degree of freedom bug
JEFworks-Lab · Nov 4, 2024 · 1d3bcca · 1d3bcca
2 parents 7b0e0a4 + b5f69b2
commit 1d3bcca
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/R/import.R b/R/import.R
@@ -167,8 +167,8 @@ getOverdispersedGenes <- function(counts,
   }
   df$res <- -Inf
   df$res[vi] <- resid(m,type='response')
-  n.cells <- ncol(mat)
-  n.obs <- nrow(mat)
+  n.cells <- nrow(mat)
+  n.obs <- ncol(mat)
   df$lp <- as.numeric(pf(exp(df$res),n.obs,n.obs,lower.tail=FALSE,log.p=TRUE))
   df$lpa <- bh.adjust(df$lp,log=TRUE)
   df$qv <- as.numeric(qchisq(df$lp, n.cells-1, lower.tail = FALSE,log.p=TRUE)/n.cells)