Question on Affxparser (createCel)
1
0
Entering edit mode
Guido Hooiveld ★ 3.0k
@guido-hooiveld-2020
Last seen 5 days ago
Wageningen University, Wageningen, the …
As you may recall, I am trying to create CEL files from an AffyBatch object (after having my original CEL files be corrected by Harshlight). http://article.gmane.org/gmane.science.biology.informatics.conductor/1 87 57 Since my R-skills are very basic, for this I mainly rely on understanding examples available in vignettes or BioC mailing list... As suggested, I would like to use the 'createCel' and 'updateCel' commands from Affxparser. However, there is one thing in the 'createCel' help page i don't understand, which has to do with the type of the CEL file (binary,ASCII): why should the version of the original input CEL files (that are used as template to create the new, updated files) explicity be of version 3 (type = ASCII?)? After reading the help page I understand that the OUTPUT cannot be anything else than v4 (type = binary), but what about the input? See below. Thanks, Guido <copy paste="" from="" help="" page="" 'createcel'=""> # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - # Read the CEL header # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - hdr <- readCelHeader(file) stopifnot(hdr$version == 3) ------------------------------------------------ Guido Hooiveld, PhD Nutrition, Metabolism & Genomics Group Division of Human Nutrition Wageningen University Biotechnion, Bomenweg 2 NL-6703 HD Wageningen the Netherlands tel: (+)31 317 485788 fax: (+)31 317 483342 internet: http://nutrigene.4t.com <http: nutrigene.4t.com=""/> email: guido.hooiveld@wur.nl [[alternative HTML version deleted]]
affxparser affxparser • 526 views
ADD COMMENT
0
Entering edit mode
@henrik-bengtsson-4333
Last seen 4 months ago
United States
Hi. On Tue, Aug 5, 2008 at 6:46 AM, Hooiveld, Guido <guido.hooiveld at="" wur.nl=""> wrote: > > As you may recall, I am trying to create CEL files from an AffyBatch > object (after having my original CEL files be corrected by Harshlight). > http://article.gmane.org/gmane.science.biology.informatics.conductor /187 > 57 > > Since my R-skills are very basic, for this I mainly rely on > understanding examples available in vignettes or BioC mailing list... > As suggested, I would like to use the 'createCel' and 'updateCel' > commands from Affxparser. > > However, there is one thing in the 'createCel' help page i don't > understand, which has to do with the type of the CEL file > (binary,ASCII): > why should the version of the original input CEL files (that are used as > template to create the new, updated files) explicity be of version 3 > (type = ASCII?)? > > After reading the help page I understand that the OUTPUT cannot be > anything else than v4 (type = binary), but what about the input? See > below. The template CEL header (argument 'header') can be of any version, i.e. v3 (ASCII), v4 (binary;XDA), or v1 (binary;Calvin). I have clarified this in the help(createCel) in the next release/the developers version. The easiest is to "reuse" one you get from readCelHeader() of an existing CEL file The stopifnot(hdr$version == 3) expression of example(createCel) is just to assert that the generic code above found an ASCII file. It is just to catch stupid mistakes when R CMD check runs. In the next release/the developers version, I have added a comment explaining this. Hope this helps Cheers Henrik > > Thanks, > Guido > > <copy paste="" from="" help="" page="" 'createcel'=""> > # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - > # Read the CEL header > # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - > hdr <- readCelHeader(file) > stopifnot(hdr$version == 3) > > > ------------------------------------------------ > Guido Hooiveld, PhD > Nutrition, Metabolism & Genomics Group > Division of Human Nutrition > Wageningen University > Biotechnion, Bomenweg 2 > NL-6703 HD Wageningen > the Netherlands > tel: (+)31 317 485788 > fax: (+)31 317 483342 > internet: http://nutrigene.4t.com <http: nutrigene.4t.com=""/> > email: guido.hooiveld at wur.nl > > > [[alternative HTML version deleted]] > > _______________________________________________ > 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 >
ADD COMMENT

Login before adding your answer.

Traffic: 413 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6