inmoose.edgepy.exactTestDoubleTail

inmoose.edgepy.exactTestDoubleTail(y1, y2, dispersion=0, big_count=900)

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

This function computes two-sided p-values by doubling the smaller tail probability.

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.

  • big_count (int) – count size above which asymptotic beta approximation will be used.

Returns:

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

Return type:

ndarray