[BioC] export Views of Rle in rtracklayer
Jason Lu
jasonlu68 at gmail.com
Tue Nov 22 18:39:24 CET 2011
Hi list,
I have a question and would like to get your help.
What I would like to get is to show the read coverage in the exon
regions (only) in the ucsc genome browser. My question is how to
export the Views so I can load into ucsc browser.
Here is what I have:
#
>cvg = coverage(ranges(ss))
>cvg.view = Views(cvg,ranges(exon.gr)) # exon.gr is a GRanges
>cvg.view
Views on a 22109556-length Rle subject
views:
start end width
[1] 22109317 22109688 372 [5 5 7 5 6 6 6 5 5 5 5 5 4 4 5 8 8 8 9 9 9 9 ...]
[2] 22084156 22084276 121 [4 4 4 3 3 3 3 3 3 3 3 3 3 4 4 4 4 4 4 4 4 4 ...]
[3] 22083039 22083195 157 [2 1 1 0 0 0 1 1 1 1 2 2 2 2 2 2 2 2 2 3 3 4 ...]
[4] 22079485 22079612 128 [0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 2 2 2 2 2 2 2 ...]
[5] 22079020 22079144 125 [6 6 6 6 7 7 7 7 7 7 7 7 7 7 6 6 6 6 6 6 5 4 ...]
I don't know how to modify this view such as adding chromosome info
and export this 'cvg.view' in rtracklayer.
Thanks in advance.
Jason
> sessionInfo()
R version 2.14.0 (2011-10-31)
Platform: x86_64-unknown-linux-gnu (64-bit)
locale:
[1] C
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] rtracklayer_1.14.3 RCurl_1.7-0 bitops_1.0-4.1
[4] Rsamtools_1.6.2 Biostrings_2.22.0 GenomicFeatures_1.6.4
[7] AnnotationDbi_1.16.4 Biobase_2.14.0 GenomicRanges_1.6.3
[10] IRanges_1.12.2 BiocInstaller_1.2.1
loaded via a namespace (and not attached):
[1] BSgenome_1.22.0 DBI_0.2-5 RSQLite_0.10.0 XML_3.4-3
[5] biomaRt_2.10.0 tools_2.14.0 zlibbioc_1.0.0
>
More information about the Bioconductor
mailing list