Question: warning trying to create a TxDB for S.pombe
1
gravatar for Lila M
2.5 years ago by
Lila M 790
UK
Lila M 790 wrote:

Hi! I want to annotate the genes with a peak for a S.pombe ChiPseq experiment. I always do that using ChIPseeker. As there is no TxDB genome for S.pombe, I've created it as follow

txdb_pombe <- makeTxDbFromGFF("Schizosaccharomyces_pombe.ASM294v2.34.gff3.gz", format="gff3", organism = "Schizosaccharomyces pombe", taxonomyId = "4896")

But I have the follow warning:

Import genomic features from the file as a GRanges object ... OK
Prepare the 'metadata' data frame ... OK
Make the TxDb object ... OK
Warning messages:
1: Named parameters not used in query: internal_chrom_id, chrom, length, is_circular 
2: Named parameters not used in query: internal_id, name, type, chrom, strand, start, end 
3: Named parameters not used in query: internal_id, name, chrom, strand, start, end 
4: Named parameters not used in query: internal_id, name, chrom, strand, start, end 
5: Named parameters not used in query: internal_tx_id, exon_rank, internal_exon_id, internal_cds_id 
6: Named parameters not used in query: gene_id, internal_tx_id

And when I've tried to run the program

for (i in fileList){
   peak <- readPeakFile(i, header=F)
  peakAnno <- annotatePeak(peak, tssRegion=c(-500, 2000), TxDb= "txdb_pombe")
  print(peakAnno)
    write.table(peakAnno, paste("annot_", i, sep="") , sep="\t", col.names=T, row.names = F, quote=F)
}

The follow error is reported

Error in (function (classes, fdef, mtable)  : 
  unable to find an inherited method for function ‘metadata’ for signature ‘"character"’

Any suggestion?

Many thanks!

ADD COMMENTlink modified 2.5 years ago by ddiez1.8k • written 2.5 years ago by Lila M 790
2
gravatar for ddiez
2.5 years ago by
ddiez1.8k
Japan
ddiez1.8k wrote:

The help for annotatePeak says this about the TxDb argument:

TxDb TxDb object

But you are passing a character string "txdb_pombe". Try this instead:

annotatePeak(peak, tssRegion = c(-500, 2000), TxDb = txdb_pombe)
ADD COMMENTlink written 2.5 years ago by ddiez1.8k

Totally right, problem solved!!!

ADD REPLYlink written 2.5 years ago by Lila M 790
Please log in to add an answer.

Help
Access

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