Error with autoplot of transcriptDB object

Mark Dunning <[email protected]> Wed, 3 Sep 2014 11:16:06 +0100
Newsgroups gmane.science.biology.informatics.conductor
Message-ID <CAGpynuk-bEWHxVMkG_Z0H8bkGjohaHnGa72bS5M9H-8jRsg5ZQ@mail.gmail.com>
Hi all,

I've having problems running some code that was previously Ok for me. Can
anyone see what is going wrong? Did ggbio change how it plots transcripts?

Regards,

Mark

>library(TxDb.Hsapiens.UCSC.hg19.knownGene)
>txdb <- TxDb.Hsapiens.UCSC.hg19.knownGene

>exons <- exonsBy(txdb, "gene")
>ibrary(ggbio)

>autoplot(txdb,which=exons[["49"]])

Parsing transcripts...
Parsing exons...
Parsing cds...
Parsing utrs...
------exons...
------cdss...
------introns...
------utr...
aggregating...
Done
"gap" is not matching to following arbitrary model terms"cds CDS Cds exon
EXON Exon utr UTR Utr"
Constructing graphics...
Error in `[[<-`(`*tmp*`, name, value = 1L) :
  1 elements in value to replace 0 elements




> sessionInfo()
R version 3.1.1 (2014-07-10)
Platform: x86_64-pc-linux-gnu (64-bit)

locale:
 [1] LC_CTYPE=en_GB.UTF-8       LC_NUMERIC=C
LC_TIME=en_GB.UTF-8        LC_COLLATE=en_GB.UTF-8
 [5] LC_MONETARY=en_GB.UTF-8    LC_MESSAGES=en_GB.UTF-8
LC_PAPER=en_GB.UTF-8       LC_NAME=C
 [9] LC_ADDRESS=C               LC_TELEPHONE=C
LC_MEASUREMENT=en_GB.UTF-8 LC_IDENTIFICATION=C

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

other attached packages:
 [1] XVector_0.4.0
ggbio_1.12.10
 [3] ggplot2_1.0.0
TxDb.Hsapiens.UCSC.hg19.knownGene_2.14.0
 [5] GenomicFeatures_1.16.2
AnnotationDbi_1.26.0
 [7] Biobase_2.24.0
GenomicRanges_1.16.4
 [9] GenomeInfoDb_1.0.2
IRanges_1.22.10
[11] BiocGenerics_0.10.0

loaded via a namespace (and not attached):
 [1] BatchJobs_1.3            BBmisc_1.7
BiocParallel_0.6.1       biomaRt_2.20.0
 [5] Biostrings_2.32.1        biovizBase_1.12.3
bitops_1.0-6             brew_1.0-6
 [9] BSgenome_1.32.0          checkmate_1.3
cluster_1.15.2           codetools_0.2-9
[13] colorspace_1.2-4         DBI_0.2-7
dichromat_2.0-0          digest_0.6.4
[17] fail_1.2                 foreach_1.4.2
Formula_1.1-2            GenomicAlignments_1.0.6
[21] grid_3.1.1               gridExtra_0.9.1
gtable_0.1.2             Hmisc_3.14-4
[25] iterators_1.0.7          lattice_0.20-29
latticeExtra_0.6-26      MASS_7.3-34
[29] munsell_0.4.2            plyr_1.8.1
proto_0.3-10             RColorBrewer_1.0-5
[33] Rcpp_0.11.2              RCurl_1.95-4.1
reshape2_1.4             Rsamtools_1.16.1
[37] RSQLite_0.11.4           rtracklayer_1.24.2
scales_0.2.4             sendmailR_1.1-2
[41] splines_3.1.1            stats4_3.1.1
stringr_0.6.2            survival_2.37-7
[45] tools_3.1.1              VariantAnnotation_1.10.5
XML_3.98-1.1             zlibbioc_1.10.0

	[[alternative HTML version deleted]]

_______________________________________________
Bioconductor mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/bioconductor
Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor