View source: R/power_estimation.R
default_proper | R Documentation |
Recent reviewers of Najib's grants have taken an increased interest in knowing the statistical power of the various experiments. He queried Dr. Corrada-Bravo who suggested PROPER. I spent some time looking through it and, with some revervations, modified its workflow to (at least in theory) be able to examine any dataset. The workflow in question is particularly odd and warrants further discussion/analysis. This function invokes PROPER exactly as it was performed in their paper.
default_proper(
de_tables,
p = 0.05,
experiment = "cheung",
nsims = 20,
reps = c(3, 5, 7, 10),
de_method = "edger",
alpha_type = "fdr",
alpha = 0.1,
stratify = "expr",
target = "lfc",
add_coverage = TRUE,
filter = "none",
delta = 0.5
)
de_tables |
A set of differential expression results, presumably from EdgeR or DESeq2. |
p |
Cutoff |
experiment |
The default data set in PROPER is entitled 'cheung'. |
nsims |
Number of simulations to perform. |
reps |
Simulate these number of experimental replicates. |
de_method |
There are a couple choices here for tools which are pretty old, my version of this only accepts deseq or edger. |
alpha_type |
I assume p-adjust type. |
alpha |
Accepted fdr rate. |
stratify |
There are a few options here, I don't fully understand them. |
target |
Cutoff. |
add_coverage |
Add a line showing the actual coverage observed? |
filter |
Apply a filter? |
delta |
Not epsilon! (E.g. I forget what this does). |
List containing the various results and plots from proper.
[PROPER] DOI:10.1093/bioinformatics/btu640
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.