Question: error with calculateRPKM in customProDB
gravatar for javamng
2.3 years ago by
javamng0 wrote:

I am trying to create customProDB using the Human Ensembl 38.87 version.

Recently, when I run the claculateRPKM of the customProDB package after installing the related package, the same error is always output regardless of the data set.

By the way, I have seen the same error in the data that I have successfully executed before.

The error message is:

> RPKM <- calculateRPKM(bamFile, exon, proteincodingonly=TRUE, ids)
Error in keepSeqlevels(anno, seqlevels(galn)) :
  invalid seqlevels: KI270729.1, KI270743.1, GL000008.2, KI270747.1,
  KI270722.1, KI270742.1, KI270736.1, GL000224.1, KI270719.1, KI270712.1,
  KI270706.1, KI270725.1, KI270715.1, KI270749.1, KI270741.1, GL000221.1,
  KI270716.1, KI270751.1, KI270519.1, GL000214.1, KI270708.1, KI270730.1,
  KI270438.1, KI270737.1, KI270738.1, KI270748.1, KI270435.1, GL000208.1,
  KI270538.1, KI270756.1, KI270739.1, KI270757.1, KI270709.1, KI270746.1,
  KI270753.1, KI270589.1, KI270735.1, KI270745.1, KI270714.1, KI270732.1,
  KI270754.1, KI270710.1, KI270717.1, KI270724.1, KI270720.1, KI270723.1,
  KI270718.1, KI270317.1, KI270740.1, KI270755.1, KI270707.1, KI270579.1,
  KI270752.1, KI270512.1, KI270322.1, GL000226.1, KI270311.1, KI270366.1,
  KI270511.1, KI270448.1, KI270521.1, KI270581.1, KI270582.1, KI270515.1,
  KI270588.1, KI270591.1, KI270522.1, KI270507.1, KI270590.1, KI270584.1,
  KI270320.1, KI270382.1, KI270468.1, KI270467.1, KI270362.1, KI270517.1,
  KI270593.1, KI270528.1, KI270587.


>  sessionInfo()

R version 3.4.1 (2017-06-30)
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=en_US.UTF-8       LC_NUMERIC=C
 [3] LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8
 [7] LC_PAPER=en_US.UTF-8       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] customProDB_1.16.0   biomaRt_2.32.1       AnnotationDbi_1.38.2
[4] Biobase_2.36.2       IRanges_2.10.2       S4Vectors_0.14.3
[7] BiocGenerics_0.22.0

loaded via a namespace (and not attached):
 [1] Rcpp_0.12.12               plyr_1.8.4
 [3] compiler_3.4.1             GenomeInfoDb_1.12.2
 [5] XVector_0.16.0             GenomicFeatures_1.28.4
 [7] bitops_1.0-6               tools_3.4.1
 [9] zlibbioc_1.22.0            digest_0.6.12
[11] bit_1.1-12                 AhoCorasickTrie_0.1.0
[13] RSQLite_2.0                memoise_1.1.0
[15] tibble_1.3.3               lattice_0.20-35
[17] BSgenome_1.44.0            rlang_0.1.1
[19] Matrix_1.2-10              DelayedArray_0.2.7
[21] DBI_0.7                    GenomeInfoDbData_0.99.0
[23] stringr_1.2.0              rtracklayer_1.36.4
[25] Biostrings_2.44.2          bit64_0.9-7
[27] grid_3.4.1                 XML_3.98-1.9
[29] BiocParallel_1.10.1        magrittr_1.5
[31] blob_1.1.0                 Rsamtools_1.28.0
[33] matrixStats_0.52.2         GenomicRanges_1.28.4
[35] GenomicAlignments_1.12.1   SummarizedExperiment_1.6.3
[37] stringi_1.1.5              RCurl_1.95-4.8
[39] VariantAnnotation_1.22.3

software error customprodb • 753 views
ADD COMMENTlink written 2.3 years ago by javamng0
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: 169 users visited in the last hour