We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
1 parent d0f5372 commit 62e30b1Copy full SHA for 62e30b1
R/deg.R
@@ -565,7 +565,8 @@ deg.deseq2 <- function(se,
565
pval <- df$pvalue
566
pval[pval < 1e-300] <- 1e-300
567
pval[is.na(df$padj)] <- NA
568
- S4Vectors::metadata(se)$deg[[paste0(prefix, "_harmonic_mean_pvalue")]] <- 1./mean(1./pval, na.rm }
+ S4Vectors::metadata(se)$deg[[paste0(prefix, "_harmonic_mean_pvalue")]] <- 1./mean(1./pval, na.rm=TRUE)
569
+ }
570
return(se)
571
}
572
0 commit comments