Description Usage Arguments Details Value Author(s) References See Also Examples
Compute genewise exact tests for differences in the means between two groups of negative-binomially distributed counts.
1 2 3 4 5 6 | exactTest(object, pair=1:2, dispersion="auto", rejection.region="doubletail",
big.count=900, prior.count=0.125)
exactTestDoubleTail(y1, y2, dispersion=0, big.count=900)
exactTestBySmallP(y1, y2, dispersion=0)
exactTestByDeviance(y1, y2, dispersion=0)
exactTestBetaApprox(y1, y2, dispersion=0)
|
object |
an object of class |
pair |
vector of length two, either numeric or character, providing the pair of groups to be compared; if a character vector, then should be the names of two groups (e.g. two levels of |
dispersion |
either a numeric vector of dispersions or a character string indicating that dispersions should be taken from the data object.
If a numeric vector, then can be either of length one or of length equal to the number of genes.
Allowable character values are |
rejection.region |
type of rejection region for two-sided exact test. Possible values are |
big.count |
count size above which asymptotic beta approximation will be used. |
prior.count |
average prior count used to shrink log-fold-changes. Larger values produce more shrinkage. |
y1 |
numeric matrix of counts for the first 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 |
y2 |
numeric 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 |
The functions test for differential expression between two groups of count libraries.
They implement the exact test proposed by Robinson and Smyth (2008) for a difference in mean between two groups of negative binomial random variables.
The functions accept two groups of count libraries, and a test is performed for each row of data.
For each row, the test is conditional on the sum of counts for that row.
The test can be viewed as a generalization of the well-known exact binomial test (implemented in binomTest
) but generalized to overdispersed counts.
exactTest
is the main user-level function, and produces an object containing all the necessary components for downstream analysis.
exactTest
calls one of the low level functions exactTestDoubleTail
, exactTestBetaApprox
, exactTestBySmallP
or exactTestByDeviance
to do the p-value computation.
The low level functions all assume that the libraries have been normalized to have the same size, i.e., to have the same expected column sum under the null hypothesis.
exactTest
equalizes the library sizes using equalizeLibSizes
before calling the low level functions.
The functions exactTestDoubleTail
, exactTestBySmallP
and exactTestByDeviance
correspond to different ways to define the two-sided rejection region when the two groups have different numbers of samples.
exactTestBySmallP
implements the method of small probabilities as proposed by Robinson and Smyth (2008).
This method corresponds exactly to binomTest
as the dispersion approaches zero, but gives poor results when the dispersion is very large.
exactTestDoubleTail
computes two-sided p-values by doubling the smaller tail probability.
exactTestByDeviance
uses the deviance goodness of fit statistics to define the rejection region, and is therefore equivalent to a conditional likelihood ratio test.
Note that rejection.region="smallp"
is no longer recommended.
It is preserved as an option only for backward compatiblity with early versions of edgeR.
rejection.region="deviance"
has good theoretical statistical properties but is relatively slow to compute.
rejection.region="doubletail"
is just slightly more conservative than rejection.region="deviance"
, but is recommended because of its much greater speed.
For general remarks on different types of rejection regions for exact tests see Gibbons and Pratt (1975).
exactTestBetaApprox
implements an asymptotic beta distribution approximation to the conditional count distribution.
It is called by the other functions for rows with both group counts greater than big.count
.
exactTest
produces an object of class DGEExact
containing the following components:
table |
data frame containing columns for the log2-fold-change ( |
comparison |
character vector giving the names of the two groups being compared. |
genes |
data frame containing annotation for each gene; taken from |
The low-level functions, exactTestDoubleTail
etc, produce a numeric vector of genewise p-values, one for each row of y1
and y2
.
Mark Robinson, Davis McCarthy, Gordon Smyth
Robinson MD and Smyth GK (2008). Small-sample estimation of negative binomial dispersion, with applications to SAGE data. Biostatistics, 9, 321-332.
Gibbons, JD and Pratt, JW (1975). P-values: interpretation and methodology. The American Statistician 29, 20-25.
1 2 3 4 5 6 7 8 9 10 | # generate raw counts from NB, create list object
y <- matrix(rnbinom(80,size=1/0.2,mu=10),nrow=20,ncol=4)
d <- DGEList(counts=y, group=c(1,1,2,2), lib.size=rep(1000,4))
de <- exactTest(d, dispersion=0.2)
topTags(de)
# same p-values using low-level function directly
p.value <- exactTestDoubleTail(y[,1:2], y[,3:4], dispersion=0.2)
sort(p.value)[1:10]
|
Loading required package: limma
Comparison of groups: 2-1
logFC logCPM PValue FDR
4 1.6476983 13.71491 0.05281761 0.7394864
2 1.3510744 13.63245 0.11455220 0.7394864
20 -1.3006595 13.38629 0.14738543 0.7394864
7 1.0675633 13.96003 0.18838665 0.7394864
6 -1.0922585 13.57473 0.20686359 0.7394864
19 -1.0374747 13.63245 0.22600371 0.7394864
1 -1.0518389 13.24543 0.25882024 0.7394864
8 0.7274741 13.54499 0.41330869 0.7792405
12 -0.8129144 13.00442 0.42612512 0.7792405
16 0.6358437 13.57473 0.47721648 0.7792405
[1] 0.05281761 0.11455220 0.14738543 0.18838665 0.20686359 0.22600371
[7] 0.25882024 0.41330869 0.42612512 0.47721648
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.