[BioC] Cannot receive data tracks back rtracklayer session

tholzman at fhcrc.org tholzman at fhcrc.org
Fri Feb 6 02:02:05 CET 2009


Hi,

This is Ted Holzman.  I am new bioconductor user, and I have a Problem.

I am using the rtracklayer package within Bioconductor.  I am trying to
receive
a track back from a session.  The view generated by the session looks fine
in the browser.  But I get no data back in the track.

> conservation <- trackSet(session,"Conservation")
> conservation
ucscTrackSet (storageMode: lockedEnvironment)
assayData: 377 features, 1 samples
  element names: dataVals
phenoData
  sampleNames: score
  varLabels and varMetadata description: none
featureData
  featureNames: ., ..1, ..., ..376  (377 total)
  fvarLabels and fvarMetadata description:
    chrom: Chromosome ID
    start: Start position on chromosome
    ...: ...
    strand: DNA strand, sense (+) or antisense (-)
    (5 total)
experimentData: use 'experimentData(object)'
Annotation:  
genome: hg18
trackLine: track name=multiz30way description=multiz30way visibility=1

All this looks fine except the genome is incorrect.  It should be "mm9".

> cons <- trackData(conservation)
> head(cons)

    chrom    start      end name strand  featMid score
.   chr12 85105663 85105741    .     NA 85105702     0
..1 chr12 85105742 85105813    .     NA 85105778     0
..2 chr12 85105814 85106306    .     NA 85106060     0
..3 chr12 85106307 85106426    .     NA 85106366     0
..4 chr12 85106427 85106780    .     NA 85106604     0
..5 chr12 85106781 85106808    .     NA 85106794     0


All the scores are 0!

Can anybody help?

Thanks,
                -Ted

****

> sessionInfo()
R version 2.8.1 (2008-12-22)
x86_64-unknown-linux-gnu

locale:
LC_CTYPE=en_US.UTF-8;LC_NUMERIC=C;LC_TIME=en_US.UTF-8;LC_COLLATE=en_US.UTF-8;LC_MONETARY=C;LC_MESSAGES=en_US.UTF-8;LC_PAPER=en_US.UTF-8;LC_NAME=C;LC_ADDRESS=C;LC_TELEPHONE=C;LC_MEASUREMENT=en_US.UTF-8;LC_IDENTIFICATION=C

attached base packages:
[1] grid      tools     stats     graphics  grDevices utils    
datasets  methods   base     

other attached packages:
[1] GenomeGraphs_1.3.4 biomaRt_1.99.2     rtracklayer_1.2.2 
RCurl_0.94-0       Biobase_2.2.1     

loaded via a namespace (and not attached):
[1] Biostrings_2.10.15 IRanges_1.0.9      lattice_0.17-20   
Matrix_0.999375-17 rJava_0.6-2        XML_1.99-0



More information about the Bioconductor mailing list