Question: PrepareAnnotationRefseq2 function error!
gravatar for yeh234
3.2 years ago by
yeh2340 wrote:
Hello everyone,

I follow the "customProDB"  tutorial for prepare annotation files and use "PGA package"  PrepareAnnotationRefseq2 function. But I got an error message and fail to generate depentent RData.

Am I current run R input below,

> options(width=70)
> transcript_ids <- c("NM_001126112", "NM_033360", "NR_073499", "NM_004448",
+ "NM_000179", "NR_029605", "NM_004333", "NM_001127511")
> pepfasta <- system.file("extdata", "refseq_pro_seq.fasta",
+ package="customProDB")
> CDSfasta <- system.file("extdata", "refseq_coding_seq.fasta",
+ package="customProDB")
> annotation_path <- tempdir()
> PrepareAnnotationRefseq2(genome='hg19', CDSfasta, pepfasta, annotation_path,
+ dbsnp = NULL, transcript_ids=transcript_ids,
+ splice_matrix=FALSE, COSMIC=FALSE)
Build TranscriptDB object (txdb.sqlite) ... 
Download the refGene table ... OK
Download the hgFixed.refLink table ... OK
Extract the 'transcripts' data frame ... OK
Extract the 'splicings' data frame ... OK
Download and preprocess the 'chrominfo' data frame ... OK
Prepare the 'metadata' data frame ... OK
Make the TxDb object ... OK
Prepare gene/transcript/protein id mapping information (ids.RData) ... Error in normArgTable(value, x) : unknown table name 'refLink'

> sessionInfo()
R version 3.3.1 (2016-06-21)
Platform: x86_64-apple-darwin13.4.0 (64-bit)
Running under: OS X 10.11.6 (El Capitan)

PGA_1.2.2; customProDB_1.12.0    

> traceback()
9: stop("unknown table name '", name, "'")
8: normArgTable(value, x)
7: .local(x, ..., value)
6: `tableName<-`(`*tmp*`, value = "refLink")
5: `tableName<-`(`*tmp*`, value = "refLink")
4: .local(x, ...)
3: ucscTableQuery(session, "refGene", table = "refLink", names = refGene[, 
2: ucscTableQuery(session, "refGene", table = "refLink", names = refGene[, 
1: PrepareAnnotationRefseq2(genome = "hg19", CDSfasta, pepfasta, 
       annotation_path, dbsnp = NULL, transcript_ids = transcript_ids, 
       splice_matrix = FALSE, COSMIC = FALSE)

Does anyone can help me to solve this error.


Thank you,


customprodb pga • 554 views
ADD COMMENTlink written 3.2 years ago by yeh2340
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 176 users visited in the last hour