I'm trying to import salmon quant tables with tximeta.
I did the transcript quantification with salmon as below:
salmon quant -i ../gencode.v32_salmon0.14.1 -l A \
-1 ${samp}_1.fastq.gz \
-2 ${samp}_2.fastq.gz \
--numGibbsSamples 20 \
-p 8 --gcBias --validateMappings -o ../salmon_quant/${samp}_quant
When I'm trying to import it as a summarized experiment in R I get the following error.
> se <- tximeta(coldata)
importing quantifications
reading in files with read_tsv
1 2 3 4 5 6
couldn't find matching transcriptome, returning un-ranged SummarizedExperiment
Error in stopifnot(all(sapply(infReps, ncol) == nreps)) :
(list) object cannot be coerced to type 'integer'
I don't know what's going and appreciate the help. Thanks
> sessionInfo()
R version 3.6.1 (2019-07-05)
Platform: x86_64-apple-darwin15.6.0 (64-bit)
Running under: macOS Mojave 10.14.5
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.6/Resources/lib/libRlapack.dylib
locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
attached base packages:
[1] parallel stats4 stats graphics grDevices utils datasets
[8] methods base
other attached packages:
[1] tximeta_1.2.2 TSAtools_0.0.0.9000 GOfuncR_1.4.0
[4] vioplot_0.3.4 zoo_1.8-7 sm_2.2-5.6
[7] rtracklayer_1.44.4 GenomicRanges_1.36.1 GenomeInfoDb_1.20.0
[10] IRanges_2.18.3 S4Vectors_0.22.1 BiocGenerics_0.30.0
loaded via a namespace (and not attached):
[1] Biobase_2.44.0 httr_1.4.1
[3] jsonlite_1.6 bit64_0.9-7
[5] gtools_3.8.1 assertthat_0.2.1
[7] BiocFileCache_1.8.0 blob_1.2.0
[9] GenomeInfoDbData_1.2.1 Rsamtools_2.0.3
[11] progress_1.2.2 pillar_1.4.3
[13] RSQLite_2.2.0 backports_1.1.5
[15] lattice_0.20-38 glue_1.3.1
[17] digest_0.6.23 XVector_0.24.0
[19] Matrix_1.2-18 XML_3.98-1.20
[21] pkgconfig_2.0.3 biomaRt_2.40.5
[23] zlibbioc_1.30.0 purrr_0.3.3
[25] BiocParallel_1.18.1 tibble_2.1.3
[27] AnnotationFilter_1.8.0 SummarizedExperiment_1.14.1
[29] GenomicFeatures_1.36.4 lazyeval_0.2.2
[31] magrittr_1.5 crayon_1.3.4
[33] memoise_1.1.0 tools_3.6.1
[35] prettyunits_1.1.0 hms_0.5.3
[37] matrixStats_0.55.0 stringr_1.4.0
[39] DelayedArray_0.10.0 AnnotationDbi_1.46.1
[41] ensembldb_2.8.1 Biostrings_2.52.0
[43] compiler_3.6.1 mapplots_1.5.1
[45] rlang_0.4.2 grid_3.6.1
[47] RCurl_1.95-4.12 tximport_1.12.3
[49] rstudioapi_0.10 rappdirs_0.3.1
[51] bitops_1.0-6 tcltk_3.6.1
[53] DBI_1.1.0 curl_4.3
[55] R6_2.4.1 GenomicAlignments_1.20.1
[57] knitr_1.26 dplyr_0.8.3
[59] bit_1.1-15 zeallot_0.1.0
[61] ProtGenerics_1.16.0 readr_1.3.1
[63] stringi_1.4.5 Rcpp_1.0.3
[65] vctrs_0.2.1 dbplyr_1.4.2
[67] tidyselect_0.2.5 xfun_0.12
Hi Mike,
I hope you're doing good. You were partially correct. Two of my samples didn't have the inferential replicates. and that was triggering the:
I removed them from the coldata data.frame, but I still get the non-matchon Txom error.
Here is the content of one of the meta_info.json files
Thanks again, Omid.
Great, so there's no more error, the message above just means that your transcriptome is not a match for the records in tximeta 1.2.2 Which transcriptome are you using?
I'll work on some better messaging for when there is a mix of data with and without inf reps.
Thanks, Upgraded the tximeta to 1.4.2 and it identified the Gencode v32 correctly.