set.seed(15)
yFactor <- as.factor( c(rep("treated", 5), rep("control", 5)) )
xData <- matrix(data=rnorm(length(letters)*10) ,nrow=length(letters), ncol=10)
rownames(xData) <- letters
geneSetABC15 <- GeneSet(geneIds=letters[1:15], setName="setABC15")
res <- npGSEA(x = xData, y = yFactor, set = geneSetABC15)
head( xSet(res) )
Run the code above in your browser using DataLab