[BioC] Error with autoplot of transcriptDB object
Mark Dunning
mark.dunning at gmail.com
Wed Sep 3 12:16:06 CEST 2014
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]]
More information about the Bioconductor
mailing list