[BioC] XPS with Arabidopsis Gene 1.1 ST Array Strip

cstrato cstrato at aon.at
Wed Jul 9 21:41:28 CEST 2014


Dear Benedikt,

I have just downloaded from Affymetrix strip array data for the 
HuGene-1_1_st array and can confirm your error messages. This means that 
xps does currently only support plate array data but does not support 
strip array data for the same type of expression array.

Best regards,
Christian


On 7/9/14 8:19 PM, cstrato wrote:
> Dear Benedikt,
>
> In principle, xps should support this array.
>
> Can you show me how you have created the scheme file 'scheme.Ara'? The
> code should be something like:
>
> scheme.aragene11st <- import.exon.scheme("aragene11st", filedir =
> file.path(scmdir, "na34"),
> file.path(libdir, "AraGene-1_1-st.clf"),
> file.path(libdir, "AraGene-1_1-st.pgf"),
> file.path(anndir, "Version14May",
> "AraGene-1_1-st-v1.na34.tair10.probeset.csv"),
> file.path(anndir, "Version14May",
> "AraGene-1_1-st-v1.na34.tair10.transcript.csv"))
>
> The error messages indicate that the CEL-files do not belong to the
> scheme or that they are in an unknown format.
>
> The GEO profile number for the Affymetrix AraGene-1_1-st array seems to
> be GPL17133. Sadly, there does not exist any data samples so that I
> cannot test this array. Thus could you please send me some CEL-files for
> this array so that I can test it.
>
> Best regards,
> Christian
> _._._._._._._._._._._._._._._._._._
> C.h.r.i.s.t.i.a.n   S.t.r.a.t.o.w.a
> V.i.e.n.n.a           A.u.s.t.r.i.a
> e.m.a.i.l:        cstrato at aon.at
> _._._._._._._._._._._._._._._._._._
>
>
>
>
> On 7/9/14 4:30 PM, Benedikt [guest] wrote:
>> Hi,
>>
>> I would like to analyse data from the Arabidopsis Gene 1.1 ST Array
>> Strip with the XPS package. I created a scheme and tried to load some
>> example data with with the function import.data(). This error message
>> will appear:
>>
>> -------------------------
>> dataAll <- import.data(scheme.Ara, "AraGene1_1", celdir=celdir,
>> celfiles=celfiles)
>> Opening file
>> </home/bathmer/R_projekte/monika/xps_test/libraryfiles/schemes/Scheme_Ara11ST_NA34.root>
>> in <READ> mode...
>> Creating new file
>> </home/bathmer/R_projekte/monika/xps_test/AraGene1_1_cel.root>...
>> Importing
>> </home/bathmer/R_projekte/monika/xps_test/cel/AraGene-1_1-st_hybridization-rep_strip-1_A01.cel>
>> as <AraGene_1_1_st_hybridization_rep_strip_1_A01.cel>...
>> Error: Number of data groups is not 1!
>> Error: CEL-file with version/magic number
>> </home/bathmer/R_projekte/monika/xps_test/cel/AraGene-1_1-st_hybridization-rep_strip-1_A01.cel>
>> is not supported.
>> Fehler in import.data(scheme.Ara, "AraGene1_1", celdir = celdir,
>> celfiles = celfiles) :
>>    error in function ‘ImportData’
>> Zusätzlich: Warnmeldung:
>> In import.data(scheme.Ara, "AraGene1_1", celdir = celdir, celfiles =
>> celfiles) :
>>    characters [](){}.:# etc in 'celnames' will be replaced  with '_'
>> ----------------------------
>>
>> Does XPS support this kind of cel files?
>>
>> Regards,
>>
>> Benedikt
>>
>>   -- output of sessionInfo():
>>
>> sessionInfo()
>> R version 3.0.2 (2013-09-25)
>> Platform: x86_64-pc-linux-gnu (64-bit)
>>
>> locale:
>>   [1] LC_CTYPE=de_DE.UTF-8       LC_NUMERIC=C
>>   [3] LC_TIME=en_US.UTF-8        LC_COLLATE=de_DE.UTF-8
>>   [5] LC_MONETARY=en_US.UTF-8    LC_MESSAGES=de_DE.UTF-8
>>   [7] LC_PAPER=en_US.UTF-8       LC_NAME=C
>>   [9] LC_ADDRESS=C               LC_TELEPHONE=C
>> [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
>>
>> attached base packages:
>> [1] stats     graphics  grDevices utils     datasets  methods   base
>>
>> other attached packages:
>> [1] xps_1.24.0
>>
>> loaded via a namespace (and not attached):
>> [1] tools_3.0.2
>>> traceback()
>> Kein traceback vefügbar
>>
>> --
>> Sent via the guest posting facility at bioconductor.org.
>>
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at r-project.org
> 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