nignatiadis

ihw_deseq <- function(deseq_res, filter, alpha, adjustment_type){

    if (missing(filter)) {
    filter <- deseq_res$baseMean
  }
  stopifnot(length(filter) == nrow(deseq_res))
  
  ihw_res <- ihw(deseq_res$pvalue, filter, alpha=alpha, adjustment_type=adjustment_type)
  
  deseq_res$padj <- adj_pvalues(ihw_res)