[BioC] Differentially expressed gene list

carotenuto@igb.cnr.it carotenuto at igb.cnr.it
Fri Oct 7 16:11:28 CEST 2005



Hi,when I perform SAM on  my array data(siggenes)I have some problems in 
retrieving the separate lists of up regulated and down regulated genes.
When I write:
fold<-function(x){				
gruppi<-split(x,controllo)
geni1<-abs(mean(gruppi[[2]])-mean(gruppi[[1]]))
return(geni1)
}
fold<-esApply(expr.contr.tratt.4,1,fold)
filtro.geni.4_prova<-expr.contr.tratt.4[which(fold>=1),]
filtro.geni.4_prova

everything is working well and I can obtain the correct HTML output,but if I 
write the same using fold change<1 to obtain downregulated genes,I have the 
following  warning message:

###Warning messages:
####1: the condition has length > 1 and only the first element will be used in: 
if (!nchar(text <- getText(object))) return("") 
####2: the condition has length > 1 and only the first element will be used in: 
if (!nchar(text <- getText(object))) return("") 
####3: the condition has length > 1 and only the first element will be used in: 
if (!nchar(text <- getText(object))) return("") 
####4: the condition has length > 1 and only the first element will be used in: 
if (!nchar(text <- getText(object))) return("") 
####5: the condition has length > 1 and only the first element will be used in: 
if (!nchar(text <- getText(object))) return("")

Where is the mistake?
Thank you for your help
Sara



More information about the Bioconductor mailing list