Description Usage Arguments Details Value Examples
View source: R/nest_confects.R
Find sets of discoveries for a range of effect sizes, controlling the False Discovery Rate (FDR) for each set.
1  nest_confects(n, pfunc, fdr = 0.05, step = 0.001, full = FALSE)

n 
Number of items being tested. 
pfunc 
A 
fdr 
False Discovery Rate to control for. 
step 
Granularity of effect sizes to test. 
full 
If TRUE, also include FDRadjusted pvalue that effect size is nonzero. Note that this is against the spirit of the topconfects approach. 
This is a general purpose function, which can be applied to any method of calculting pvalues (supplied as a function argument) for the null hypothesis that the effect size is smaller than a given amount.
A "Topconfects" object, containing a table of results and various associated information.
The most important part of this object is the $table element, a data frame with the following columns:
rank
 Ranking by confect
and for equal
confect
by pvalue at that effect size.
index
 Number of
the test, between 1 and n.
confect
 CONfident efFECT size.
The usage is as follows: To find a set of tests which have effect size
greater than x with the specified FDR, take the rows with abs(confect)
>= x
. Once the set is selected, the confect values provide confidence bounds
on the effect size with False Coveragestatement Rate (FCR) at the same level
as the FDR.
One may essentially take the top however many rows of the data frame and
these will be the best set of results of that size to dependably have an
effect size that is as large as possible. However if some genes have the same
abs(confect)
, all or none should be selected.
Some rows in the output may be given the same confect
, even if
step
is made small. This is an expected behaviour of the algorithm.
(This is similar to FDR adjustment of pvalues sometimes resulting in a run
of the same adjusted pvalue, even if all the input pvalues are distinct.)
Some wrappers around this function may add a sign to the confect
column, if it makes sense to do so. They will also generally add an
effect
column, containing an estimate of the effect size that aims to
be unbiassed rather than a conservative lower bound.
1 2 3 4 5 6 7 
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.