genes.up <- c(1:10)
genes.down <- c(21:30)
set.seed(123)
scores <- matrix(rnorm(200), ncol=2)
## the first gene set receives increased /
## decreased scores in the first experiment
scores[genes.up,1] <- scores[genes.up ,1] + 1
scores[genes.down,1] <- scores[genes.down,1] - 1
## significantly greater
signedRankSumTest( statistics = scores[,1],
index.up = genes.up,
index.down = genes.down)
## not signficant
signedRankSumTest( statistics = scores[,2],
index.up = genes.up,
index.down = genes.down)
Run the code above in your browser using DataLab