[BioC] Extracting expressions after a multtest procedure (part 2)

[JUAN MIGUEL MARIN DIAZARAQUE]

Hello,

In previous message I wonder how to extract just the expressions of 102 
genes under two conditions after a multtest procedure.

At last, I wrote this unelegant code that seems to work:

quedan <- resT[res$adjp[,"BH"]<0.05,]
busco <- dimnames(quedan)[[1]]

cosa <- NULL
for (i in 1:length(dimnames(exprs(esetSub))[[1]]))
	for (j in 1:length(busco))
		if (dimnames(exprs(esetSub[i,]))[[1]] == busco[j]) {
		cosa <- rbind(cosa, as.data.frame(exprs(esetSub))[i,]) }
	end
end

But it runs very slowly...

Best

jm~

_______________________________

     Juan Miguel Marin

http://www.est.uc3m.es/jmmarin

    Dep. of Statistics
University Carlos III of Madrid
        Spain (E.U.)