Question: error with calculateRPKM in customProDB
0
gravatar for javamng
20 months ago by
javamng0
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/libRblas.so

locale:
 [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C
 [3] LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8
 [5] LC_MONETARY=en_US.UTF-8    LC_MESSAGES=en_US.UTF-8
 [7] LC_PAPER=en_US.UTF-8       LC_NAME=C
 [9] LC_ADDRESS=C               LC_TELEPHONE=C
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=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 • 680 views
ADD COMMENTlink written 20 months ago by javamng0
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 16.09
Traffic: 125 users visited in the last hour