get p-values of t-test values for volcano

Description

get p-values of t-test values for volcano

Usage

1
getTValuesForVolcano(x, y, alternative = "two.sided")

Arguments

x

- one data matrix

y

- second data matrix

alternative

two.sided, less, greater

Value

list with three fields fchange (fold change) , pval and pvaladj

Examples

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 <- getTValuesForVolcano(a,b)
volcanoplot(res$fchange , res$pval)

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.