Using NB distribution to calculate p-value for each gene as well as adjust p-value
the method for adjusting p-value, default is 'BH'. For details, see
switch for DE detection through fold-change, which will use a normal distribution (N(0,sd)) to test the significance of log2 fold-change. The sd is estimated through a quantile function of gamma distribution at
This function firstly calls p-value used
pnbinom to call pvalue based on sum of counts difference between two
groups or used
pnorm to call pvalue via log2 fold-change, then adjusts the pvalues via
p.adjust method. In addition, it also shrink the log2 fold-change towards a common dispersion
after pvalue calling.
ABSDataSet object with additional elements: shrinked log2 fold-change, pvalue and adjusted p-value,
denoted by foldChange pvalue and adj-pvalue, respectively. Use the
results method to get access it.
this function should run after
1 2 3 4 5 6