Saving/loading TxDb objects
1
1
Entering edit mode
@james-w-macdonald-5106
Last seen 18 hours ago
United States

What is the recommended method for saving/loading a TxDb object?

> example(makeTxDb)

mkTxDb> transcripts <- data.frame(
mkTxDb+                    tx_id=1:3,
mkTxDb+                    tx_chrom="chr1",
mkTxDb+                    tx_strand=c("-", "+", "+"),
mkTxDb+                    tx_start=c(1, 2001, 2001),
mkTxDb+                    tx_end=c(999, 2199, 2199))

mkTxDb> splicings <-  data.frame(
mkTxDb+                    tx_id=c(1L, 2L, 2L, 2L, 3L, 3L),
mkTxDb+                    exon_rank=c(1, 1, 2, 3, 1, 2),
mkTxDb+                    exon_start=c(1, 2001, 2101, 2131, 2001, 2131),
mkTxDb+                    exon_end=c(999, 2085, 2144, 2199, 2085, 2199),
mkTxDb+                    cds_start=c(1, 2022, 2101, 2131, NA, NA),
mkTxDb+                    cds_end=c(999, 2085, 2144, 2193, NA, NA))

mkTxDb> txdb <- makeTxDb(transcripts, splicings)
Warning message:
In .normarg_makeTxDb_chrominfo(chrominfo, transcripts$tx_chrom,  :
  chromosome lengths and circularity flags are not available for this TxDb object

> txdb
TxDb object:
# Db type: TxDb
# Supporting package: GenomicFeatures
# transcript_nrow: 3
# exon_nrow: 4
# cds_nrow: 4
# Db created by: GenomicFeatures package from Bioconductor
# Creation time: 2015-06-09 07:13:11 -0700 (Tue, 09 Jun 2015)
# GenomicFeatures version at creation time: 1.20.1
# RSQLite version at creation time: 1.0.0
# DBSCHEMAVERSION: 1.1


> save(txdb, file = "tmp.Rdata")
> rm(txdb)
> load("tmp.Rdata")
> txdb
TxDb object:
Error in sqliteSendQuery(con, statement, bind.data) :
  expired SQLiteConnection

 

genomicfeatures • 3.2k views
ADD COMMENT
4
Entering edit mode
@martin-morgan-1513
Last seen 17 days ago
United States

saveDb() / loadDb() in AnnotationDbi.

ADD COMMENT
0
Entering edit mode

Nice! Thanks Martin.

ADD REPLY

Login before adding your answer.

Traffic: 624 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6