[BioC] probeAnno for tilingArray package.
Wolfgang Huber
whuber at embl.de
Sat Sep 11 14:14:54 CEST 2010
Dear Sohail
in R, the character string "NA" is different from the value NA. Try with
supplying the second one for the argument 'strands'. Also, please make
sure that the 'chr' argument is set to the right value, the default
(1:17) is likely not to work with S.pombe. You can see what the
chromosome names are in your organism by inspecting the probeAnno object.
Best wishes
Wolfgang
On Sep/8/10 10:33 PM, sohkhan at notes.cc.sunysb.edu wrote:
> Hi all,
>
> I am trying to use the "tilingArray" package to analyze Affymetrix S.pombe
> tiling arrays. I downloaded the .bpmap file from Affymetrix's website,
> and created a probeAnno object by bpmapToProbeAnno().
> However, I am getting the following error when I try to segment the data
> by segChrom().
>
>> segdf<-segChrom(cels,probeAnno,strands="NA")
> Running 'segment' on chromosome 1.NAError in probeAnno[w] :
> No mapping '1.NA.start' in this 'probeAnno' object.
>
>> sessionInfo()
> R version 2.11.1 (2010-05-31)
> x86_64-apple-darwin9.8.0
>
> locale:
> [1] en_US.UTF-8/en_US.UTF-8/C/C/en_US.UTF-8/en_US.UTF-8
>
> attached base packages:
> [1] grid stats graphics grDevices utils datasets methods
> base
>
> other attached packages:
> [1] davidTiling_1.2.12 GO.db_2.4.1 RSQLite_0.9-2 DBI_0.2-5
> AnnotationDbi_1.10.2 tilingArray_1.26.0 pixmap_0.4-10 Starr_1.4.4
>
> [9] affy_1.26.1 Ringo_1.12.0 Matrix_0.999375-39
> lattice_0.18-8 limma_3.4.4 RColorBrewer_1.0-2 Biobase_2.8.0
> affxparser_1.20.0
>
> loaded via a namespace (and not attached):
> [1] affyio_1.16.0 annotate_1.26.1 genefilter_1.30.0
> MASS_7.3-6 preprocessCore_1.10.0 pspline_1.0-14 splines_2.11.1
> strucchange_1.4-1
> [9] survival_2.35-8 tools_2.11.1 vsn_3.16.0 xtable_1.5-6
>
>
> I am wondering if I used the correct .bpmap file??
> Thanks for any suggestion/advice.
>
> -Sohail
> SBU, NY
>>
> [[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
--
Wolfgang Huber
EMBL
http://www.embl.de/research/units/genome_biology/huber
More information about the Bioconductor
mailing list