Description Usage Arguments Value Examples
View source: R/getValuesForVolcano.R
get p-values of wilcoxon rank sum test for volcano
1 | getWRValuesForVolcano(x, y, paired = FALSE, adjust = TRUE)
|
x |
- one data matrix |
y |
- second data matrix |
paired |
a logical indicating whether you want a paired t-test. |
adjust |
pvalues using Benjamin Hochberg |
list with two fields fchange (fold change) and pval
1 2 3 4 5 6 7 | a <- t(replicate(200,rnorm(20,runif(1,-3,3),1)))
b <- a[1:100,]
a <- a[101:200,]
boxplot(t(a[1:20,]))
boxplot(t(b[1:20,]))
res <- getWRValuesForVolcano(a,b)
volcanoplot(res$fchange , res$pval)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.