inmoose.edgepy.exactTestBySmallP

inmoose.edgepy.exactTestBySmallP(y1, y2, dispersion=0)

Compute genewise p-values for differences in the means between two groups of negative-binomially distributed counts.

This function implements the method of small probabilities as proposed by [Robinson2008]. This method corresponds exactly to binomTest() as the dispersion approaches zero, but gives poor results when the dispersion is very large.

See also

exactTest

Parameters:
  • y1 (matrix) – matrix of counts for the first of the two experimental groups to be tested for differences. Rows correspond to genes and columns to libraries. Libraries are assumed to be equal in size – e.g. adjusted pseudocounts from the output of equalizeLibSizes().

  • y2 (matrix) – matrix of counts for the second of the two experimental groups to be tested for differences. Rows correspond to genes and columns to libraries. Libraries are assumed to be equal in size – e.g. adjusted pseudocounts from the output of equalizeLibSizes().

  • dispersion (array_like of floats) – an array of dispersions, either of length one or of length equal to the number of genes.

Returns:

array of genewise p-values, one for each row of y1 and y2

Return type:

ndarray