[BioC] rtracklayer - export.wig issue
Lavinia Gordon
lavinia.gordon at mcri.edu.au
Fri Apr 1 05:34:47 CEST 2011
Dear Bioconductors,
I am trying to export a RangedData object as a WIG track, but I am
getting this error:
> wig <- export.wig(B18mtrack, mycon, dataFormat = "auto")
Error in export.ucsc(object, con, subformat, append, ...) :
Track not compatible with WIG/bedGraph: Overlapping features must be
on separate strands and every feature width must be positive
> B18mtrack
RangedData with 33608 rows and 2 value columns across 1 space
space ranges | chr score
<character> <IRanges> | <character> <numeric>
1 1 [ 15865, 15866] | chr1 0.003188034
2 1 [ 18827, 18828] | chr1 0.013269155
3 1 [ 29407, 29408] | chr1 0.006626525
4 1 [ 29425, 29426] | chr1 0.004264605
5 1 [ 29435, 29436] | chr1 -0.004953431
6 1 [ 69591, 69592] | chr1 0.013326089
7 1 [135252, 135253] | chr1 0.035002515
8 1 [533950, 533951] | chr1 0.002628425
9 1 [534242, 534243] | chr1 0.020180509
... ... ... ... ... ...
33600 1 [46361270, 46361271] | chr10 0.0051111176
33601 1 [46616291, 46616292] | chr10 -0.0518829821
33602 1 [46762634, 46762635] | chr10 0.0082029002
33603 1 [46762703, 46762704] | chr10 0.0017677983
33604 1 [46763197, 46763198] | chr10 -0.0079985873
33605 1 [46763211, 46763212] | chr10 -0.0001168199
33606 1 [46776799, 46776800] | chr10 0.0929834798
33607 1 [46915981, 46915982] | chr10 -0.0077882567
33608 1 [46951477, 46951478] | chr10 -0.0025202364
However, if I take away the last line of this object, it works:
> B18mtrack
RangedData with 33607 rows and 2 value columns across 1 space
space ranges | chr score
<character> <IRanges> | <character> <numeric>
1 1 [ 15865, 15866] | chr1 0.003188034
2 1 [ 18827, 18828] | chr1 0.013269155
3 1 [ 29407, 29408] | chr1 0.006626525
4 1 [ 29425, 29426] | chr1 0.004264605
5 1 [ 29435, 29436] | chr1 -0.004953431
6 1 [ 69591, 69592] | chr1 0.013326089
7 1 [135252, 135253] | chr1 0.035002515
8 1 [533950, 533951] | chr1 0.002628425
9 1 [534242, 534243] | chr1 0.020180509
... ... ... ... ... ...
33599 1 [46344329, 46344330] | chr10 -0.0091757753
33600 1 [46361270, 46361271] | chr10 0.0051111176
33601 1 [46616291, 46616292] | chr10 -0.0518829821
33602 1 [46762634, 46762635] | chr10 0.0082029002
33603 1 [46762703, 46762704] | chr10 0.0017677983
33604 1 [46763197, 46763198] | chr10 -0.0079985873
33605 1 [46763211, 46763212] | chr10 -0.0001168199
33606 1 [46776799, 46776800] | chr10 0.0929834798
33607 1 [46915981, 46915982] | chr10 -0.0077882567
(The .wig file loads perfectly into the UCSC genome browser)
Can anyone tell me what is wrong?
I thought it might be related to this post:
https://stat.ethz.ch/pipermail/bioconductor/2009-July/028470.html
but I can't see anything like this in my data.
With thanks for your time,
Lavinia Gordon.
R version 2.12.0 (2010-10-15)
Platform: x86_64-unknown-linux-gnu (64-bit)
locale:
[1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
[3] LC_TIME=en_US.UTF-8 LC_COLLATE=en_US.UTF-8
[5] LC_MONETARY=C LC_MESSAGES=en_US.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] IRanges_1.8.9 rtracklayer_1.10.6 RCurl_1.5-0 bitops_1.0-4.1
loaded via a namespace (and not attached):
[1] Biobase_2.10.0 Biostrings_2.18.4 BSgenome_1.18.3
[4] GenomicRanges_1.2.3 XML_3.2-0
--
Senior Bioinformatics Officer
Murdoch Childrens Research Institute
Royal Children's Hospital
Flemington Road
Parkville
Victoria 3052
Australia
www.mcri.edu.au
______________________________________________________________________
This email has been scanned by the MessageLabs Email Security System.
For more information please visit http://www.messagelabs.com/email
More information about the Bioconductor
mailing list