[BioC] rtracklayer export method requires unique names

Samuel Younkin syounkin at jhsph.edu
Wed Jan 9 18:04:29 CET 2013


BioC,

It seems that the export method in rtracklayer for the GRanges class 
requires that the GRanges object have unique names.  The error message 
given when names are not unique is the following:

Error in export(object, FileForFormat(con), ...) :
   cannot export object of class 'GRanges'

As one can easily make a case that non-unique names may be troublesome, 
I am curious to know whether you intentionally disallowed non-unique names.

For future users, it would be great if the error message was not so 
misleading.  I imagine that others may label ranges by subject ID, as I 
did, and therefore are likely to create a GRanges object with non-unique 
names.

See code below for an example.

Thanks,

Sam

~~

 > foo <- GRangesForUCSCGenome( genome = "hg18", chrom = rep("chr1",10) 
, ranges = IRanges(start = 1:10, width = 5), strand = "*" )
 > names(foo) <- 1:10
 > export(foo, "foo.bed" )
 > foobar <- foo
 > names(foobar) <- c(rep(1,2),3:10)
 > export(foobar, "foobar.bed" )
Error in export(object, FileForFormat(con), ...) :
   cannot export object of class 'GRanges'
 > sessionInfo()
R Under development (unstable) (2013-01-08 r61589)
Platform: x86_64-unknown-linux-gnu (64-bit)

locale:
  [1] LC_CTYPE=en_US.iso885915       LC_NUMERIC=C
  [3] LC_TIME=en_US.iso885915        LC_COLLATE=en_US.iso885915
  [5] LC_MONETARY=en_US.iso885915    LC_MESSAGES=en_US.iso885915
  [7] LC_PAPER=C                     LC_NAME=C
  [9] LC_ADDRESS=C                   LC_TELEPHONE=C
[11] LC_MEASUREMENT=en_US.iso885915 LC_IDENTIFICATION=C

attached base packages:
[1] parallel  stats     graphics  grDevices datasets  utils     methods
[8] base

other attached packages:
  [1] rtracklayer_1.19.6    CleftCNVAssoc_0.1.0   MinimumDistance_1.3.3
  [4] GWASTools_1.5.5       gdsfmt_0.9.10         sandwich_2.2-9
  [7] zoo_1.7-9             GWASExactHW_1.01      ncdf_1.6.6
[10] Biobase_2.19.2        GenomicRanges_1.11.21 IRanges_1.17.24
[13] BiocGenerics_0.5.6

loaded via a namespace (and not attached):
  [1] affyio_1.27.1         annotate_1.37.3       AnnotationDbi_1.21.9
  [4] BiocInstaller_1.9.6   Biostrings_2.27.8     bit_1.1-9
  [7] bitops_1.0-5          BSgenome_1.27.1       codetools_0.2-8
[10] crlmm_1.17.12         DBI_0.2-5             DNAcopy_1.33.1
[13] ellipse_0.3-7         ff_2.2-10             foreach_1.4.0
[16] genefilter_1.41.1     grid_3.0.0            iterators_1.0.6
[19] lattice_0.20-13       lmtest_0.9-30         Matrix_1.0-10
[22] msm_1.1.4             mvtnorm_0.9-9994      oligoClasses_1.21.6
[25] preprocessCore_1.21.1 quantsmooth_1.25.0    RcppEigen_0.3.1.2
[28] RCurl_1.95-3          Rsamtools_1.11.14     RSQLite_0.11.2
[31] SNPchip_2.5.1         splines_3.0.0         stats4_3.0.0
[34] survival_2.37-2       tools_3.0.0           VanillaICE_1.21.10
[37] XML_3.95-0.1          xtable_1.7-0          zlibbioc_1.5.0
 >



More information about the Bioconductor mailing list