[BioC] affylmGUI gene list from Venn diagram, was: (no subject)

Simon Kidd kidd at mail.rockefeller.edu
Sat Jun 19 19:55:42 CEST 2004


At 11:51am +1000 19/6/04, James Wettenhall wrote:
>On Fri, 18 Jun 2004, Simon Kidd wrote:
>>  >class(a) <- "matrix"
>>  >data.frame(geneIDs=geneIDs,a)
>
>>  This worked right upto the end where it gave me the following error message:
>>
>>  >  data.frame(geneIDs=geneIDs,a)
>>  Error in data.frame(geneIDs = geneIDs, a) :
>>	arguments imply differing number of rows: 14010, 8
>
>Sorry, my mistake.  What you need to use in the data.frame is
>the 'clas' object from classifyTestsF, not the 'a' object
>from vennCounts.
>
>clas <- classifyTestsF(tstat,df=10,p.value=0.05)
>class(clas) <- "matrix"
>data.frame(geneIDs=geneIDs,clas)

With this I get a slightly different error message, here is what I 
typed and the error message:
--------------
load("rma.lma")
tstat <- ContrastParameterizationList[[1]]$eb$t
tstat <- matrix(rt(300,df=10),100,3)
tstat[1:33,] <- tstat[1:33,]+2
clas <- classifyTestsF(tstat,df=10,p.value=0.05)
a <- vennCounts(clas)
geneIDs <- ls(getCdfInfo(RawAffyData))
class(clas) <- "matrix"
data.frame(geneIDs=geneIDs,clas)

Error in data.frame(geneIDs = geneIDs, clas) :	arguments imply 
differing number of rows: 14010, 100

-------------

I dont know if it matters but I have eight targets of which only four 
are used for contrasts, will that make a difference?

thanks for your help.

Simon



More information about the Bioconductor mailing list