[BioC] how to save the expression values for two color data after normalization

Mark Cowley m.cowley at garvan.org.au
Tue Sep 28 01:23:43 CEST 2010


dat2$genes is probably a data.frame or matrix, not a vector

On 28/09/2010, at 12:03 AM, Steve Lianoglou wrote:

> Hi,
> 
> On Mon, Sep 27, 2010 at 7:09 AM, Budhayash Gautam
> <budhayashgautam at gmail.com> wrote:
>> Dear all
>> 
>> i am working on some smd data which is basically two color data. i
>> have normalize the data by the following command
>> 
>> 
>>> dat2<-normalizeWithinArrays(dat, method="loess","normexp", offset=50)
>> 
>> but as now i am trying to save its expression values by the following
>> command its is not giving any genes names with respect to rows without
>> which further analysis (annotation etc ) could not be done .
>> 
>>> dat.m<-dat2$M
>>> rownames(dat.m)<-dat2$genes$ProbeName
>>> write.table(dat.m, "tgondiiexpressionvalues.txt", sep="\t", row.names=T, col.names=T, quote=F)
> 
> This gives you probe id's, right? You can "reverse engineer" which
> probes map to which genes using several methods. Depending on the
> array you are using, one way would be to query biomart (using the
> biomaRt bioc package) for a list of probe id's for your chip and have
> it return the probe id & the gene symbol, or entrez id, or whatever
> you like.
> 
> Another way would be to use one of the many annotation packages
> available from bioconductor.
> 
>>> rownames(dat.m) <-dat2$genes
>> 
>> following error is coming .....
>> 
>> Error in dimnames(x) <- dn :
>>  length of 'dimnames' [1] not equal to array extent
> 
> That's because the number of rows in your `dat.m` matrix does not
> equal the length of your dat2$genes vector, so R can't make much sense
> of the assignment you are trying to make.
> 
> Look at the result of:
> 
> R> nrow(dat.m) == length(dat2$genes)
> 
> Where did you get the dat2$genes vector? Perhaps you can show us what
> you get from:
> 
> R> head(dat2$genes)
> 
> is it possible that the probe id's are the names() of this vector?
> That could make your life a bit easier.
> 
> -- 
> Steve Lianoglou
> Graduate Student: Computational Systems Biology
>  | Memorial Sloan-Kettering Cancer Center
>  | Weill Medical College of Cornell University
> Contact Info: http://cbio.mskcc.org/~lianos/contact
> 
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor



More information about the Bioconductor mailing list