# NOT RUN {
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)
# }
Run the code above in your browser using DataLab