Entering edit mode
Monther Alhamdoosh
▴
40
@monther-alhamdoosh-10001
Last seen 5.5 years ago
Australia/Melbourne/CSL Limited
Hi Gabriel,
I was trying to use your package to read a Genbank file but it failed with the following error:
> gb = readGenBank(gbFile, verbose = TRUE) Error in validObject(.Object) : invalid class “GRanges” object: 'mcols(x)' is not parallel to 'x' > traceback() 10: stop(msg, ": ", errors, domain = NA) 9: validObject(.Object) 8: initialize(value, ...) 7: initialize(value, ...) 6: new(Class, seqnames = seqnames, ranges = ranges, strand = strand, elementMetadata = mcols, seqinfo = seqinfo) 5: new_GRanges("GRanges", seqnames = seqnames, ranges = ranges, strand = strand, mcols = mcols, seqlengths = seqlengths, seqinfo = seqinfo) 4: GRanges(seqnames = stk$seqnames, ranges = IRanges(start = stk$start, end = stk$end), strand = stk$strand) 3: fill_stack_df(resthang$FEATURES[typs == "source"]) 2: parseGenBank(text = text, partial = partial, verbose = verbose, ret.seq = TRUE) 1: readGenBank(gbFile, verbose = TRUE) > sessionInfo() R version 3.4.2 (2017-09-28) Platform: x86_64-apple-darwin15.6.0 (64-bit) Running under: macOS Sierra 10.12.6 Matrix products: default BLAS: /System/Library/Frameworks/Accelerate.framework/Versions/A/Frameworks/vecLib.framework/Versions/A/libBLAS.dylib LAPACK: /Library/Frameworks/R.framework/Versions/3.4/Resources/lib/libRlapack.dylib locale: [1] en_AU.UTF-8/en_AU.UTF-8/en_AU.UTF-8/C/en_AU.UTF-8/en_AU.UTF-8 attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] genbankr_1.6.0 BiocInstaller_1.28.0 loaded via a namespace (and not attached): [1] Rcpp_0.12.14 compiler_3.4.2 GenomeInfoDb_1.14.0 [4] XVector_0.18.0 GenomicFeatures_1.30.0 prettyunits_1.0.2 [7] bitops_1.0-6 tools_3.4.2 zlibbioc_1.24.0 [10] progress_1.1.2 biomaRt_2.34.0 digest_0.6.12 [13] bit_1.1-12 BSgenome_1.46.0 lattice_0.20-35 [16] RSQLite_2.0 memoise_1.1.0 tibble_1.3.4 [19] rlang_0.1.4 Matrix_1.2-12 DelayedArray_0.4.1 [22] DBI_0.7 parallel_3.4.2 GenomeInfoDbData_0.99.1 [25] rtracklayer_1.38.2 stringr_1.2.0 knitr_1.17 [28] Biostrings_2.46.0 S4Vectors_0.16.0 IRanges_2.12.0 [31] grid_3.4.2 stats4_3.4.2 bit64_0.9-7 [34] Biobase_2.38.0 R6_2.2.2 AnnotationDbi_1.40.0 [37] XML_3.98-1.9 RMySQL_0.10.13 BiocParallel_1.12.0 [40] blob_1.1.0 magrittr_1.5 matrixStats_0.52.2 [43] GenomicAlignments_1.14.1 Rsamtools_1.30.0 BiocGenerics_0.24.0 [46] GenomicRanges_1.30.0 SummarizedExperiment_1.8.0 assertthat_0.2.0 [49] stringi_1.1.6 RCurl_1.95-4.8 VariantAnnotation_1.24.2 >
Do you know what would be the reason for this error?
Thanks,
Monther