[BioC] scan date information
Rob Dunne
rob.dunne at csiro.au
Tue Sep 22 00:42:22 CEST 2009
Thanks for that.
I am not using the development version yet but I will look out for the
new slot.
Saroj, your method doesn't work for me, perhaps your cel file is ascii?
strings SB_20D.CEL | grep DatHeader
However, I have found
$ grep --text d.a.t.e. SB_20D.CEL
text/plainaffymetrix-scan-date(2008-04-03T04:45:53Z
the "--text" option makes grep read a binary file as thought it was
text. I am not sure why I need the dots in date.
Bye
Rob
Patrick Aboyoun wrote:
> Robert,
> The answer depends on which version of R and BioC are you using. If you
> are using R <= 2.9, BioC <= 2.4, you will need to devise your own
> method; one of which was given by Saroj. If you are using R-devel and
> BioC 2.5 (devel), the eSet abstract class and its derived classes such
> as ExpressionSet contain a new slot called protocolData that contains an
> AnnotatedDataFrame object. This slot is to be populated by metadata
> contained in microarray data files. In BioC 2.5 (devel) the
> read.affybatch from affy and read.celfile from affyio add a ScanData
> column to the protocolData slot with the metadata you are looking for.
>
>
> Cheers,
> Patrick
>
>
> Saroj K Mohapatra wrote:
>> Hi Rob:
>>
>> I have a file called _16.CEL. I want to find out the date information
>> in its header. The following gives me:
>>
>> $ strings _16.CEL | grep DatHeader
>> DatHeader=[2..65534] _16:CLS=7365 RWS=7365 XIN=1 YIN=1 VE=30
>> 2.0 10/27/06 10:57:45 50207590 M10 I find a date 10/27/06. Is this
>> what you are looking for?
>>
>> Best,
>>
>> Saroj
>>
>>
>> Robert Dunne wrote:
>>> Hi List,
>>>
>>> I apologise for what may be a very simple question. How can I retrieve
>>> the scan date information from cel files?
>>>
>>> I can find the information using some editors, kate under linux shows
>>> "a f f y m e t r i x - s c a n - d a t e ( 2 0 0 8 - 0 4 - 0 3"
>>> but I can't find it all all using vi or emacs. I suppose this is
>>> something to do with encoding.
>>> Also "string file.cel | grep "d a t e"" does not work.
>>>
>>> I have tried the affxparser library but
>>> readCelHeader("file.cel")
>>> does not pick up the date.
>>>
>>> Unfortunately in many experiments the scan date turns out to be the
>>> major effect.
>>>
>>> Bye
>>> Rob
>>>
>>> _______________________________________________
>>> 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
>>>
>>
>> ------------------------------------------------------------------------
>>
>> _______________________________________________
>> 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
>
> _______________________________________________
> 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