Question: Report error in the stage of "Prepare the 'metadata' data frame" when running exomepeak() or metdiff with
gravatar for 939930094
14 months ago by
9399300940 wrote:

Hi everyone,

I'm using exomePeak to call m6A peaks with toy data and my MeRIP-Seq data, They both report the same error:

Import genomic features from the file as a GRanges object ... OK
Prepare the 'metadata' data frame ... Error in get(name, envir = asNamespace(pkg), inherits = FALSE) :

These are my commands, as described in "An Introduction to exomePeak":

> library("exomePeak")
> gtf <- system.file("extdata", "example.gtf", package="exomePeak")
> f1 <- system.file("extdata", "IP1.bam", package="exomePeak")
> f2 <- system.file("extdata", "IP2.bam", package="exomePeak")
> f3 <- system.file("extdata", "IP3.bam", package="exomePeak")
> f4 <- system.file("extdata", "IP4.bam", package="exomePeak")
> f5 <- system.file("extdata", "Input1.bam", package="exomePeak") 
> f6 <- system.file("extdata", "Input2.bam", package="exomePeak")
> f7 <- system.file("extdata", "Input3.bam", package="exomePeak")
> result <- exomepeak(GENE_ANNO_GTF=gtf,IP_BAM=c(f1,f2,f3,f4),INPUT_BAM=c(f5,f6,f7))

The sessionInfo() output:

R version 3.4.2 (2017-09-28)
Platform: x86_64-redhat-linux-gnu (64-bit)
Running under: CentOS Linux 7 (Core)

Matrix products: default
BLAS/LAPACK: /usr/lib64/R/lib/

 [1] LC_CTYPE=zh_CN.utf8       LC_NUMERIC=C             
 [3] LC_TIME=zh_CN.utf8        LC_COLLATE=zh_CN.utf8    
 [5] LC_MONETARY=zh_CN.utf8    LC_MESSAGES=zh_CN.utf8   
 [7] LC_PAPER=zh_CN.utf8       LC_NAME=C                
 [9] LC_ADDRESS=C              LC_TELEPHONE=C           

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

other attached packages:
 [1] exomePeak_2.13.2           GenomicAlignments_1.14.2  
 [3] SummarizedExperiment_1.8.1 DelayedArray_0.4.1        
 [5] matrixStats_0.54.0         rtracklayer_1.38.3        
 [7] GenomicFeatures_1.30.3     AnnotationDbi_1.40.0      
 [9] Biobase_2.38.0             Rsamtools_1.30.0          
[11] Biostrings_2.46.0          XVector_0.18.0            
[13] GenomicRanges_1.30.3       GenomeInfoDb_1.16.0       
[15] IRanges_2.14.10            S4Vectors_0.19.19 
[17] BiocGenerics_0.26.0       

loaded via a namespace (and not attached):
 [1] Rcpp_0.12.18           compiler_3.4.2         prettyunits_1.0.2     
 [4] bitops_1.0-6           tools_3.4.2            zlibbioc_1.24.0       
 [7] progress_1.2.0         biomaRt_2.34.2         digest_0.6.17         
[10] bit_1.1-14             lattice_0.20-35        RSQLite_2.1.1         
[13] memoise_1.1.0          pkgconfig_2.0.2        rlang_0.2.2           
[16] Matrix_1.2-14          DBI_1.0.0              GenomeInfoDbData_1.0.0
[19] stringr_1.3.1          httr_1.3.1             hms_0.4.2             
[22] grid_3.4.2             bit64_0.9-7            R6_2.2.2              
[25] XML_3.98-1.16          RMySQL_0.10.15         BiocParallel_1.12.0   
[28] blob_1.1.1             magrittr_1.5           assertthat_0.2.0      
[31] stringi_1.2.4          RCurl_1.95-4.11        crayon_1.3.4

Here is the traceback:

6: get(name, envir = asNamespace(pkg), inherits = FALSE)
5: GenomeInfoDb:::.taxonomyId
4: .prepareGFFMetadata(file, dataSource, organism, taxonomyId, miRBaseBuild, 
3: makeTxDbFromGFF(PARAMETERS$GENE_ANNO_GTF, format = "gtf")
2: .read.gtf(PARAMETERS)
1: exomepeak(GENE_ANNO_GTF = gtf, IP_BAM = c(f1, f2, f3, f4), INPUT_BAM = c(f5, 
       f6, f7))

The same error also occured when using metdiff() with example data and commands. What's the problem and how can I fix it?  Thanks a lot!



exomepeak • 282 views
ADD COMMENTlink written 14 months ago by 9399300940
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: 456 users visited in the last hour