Question: tximport error reading salmon quantification files...
gravatar for nschurch
19 months ago by
nschurch10 wrote:

I'm trying to use tximport with salmon quantifications. Initially I had some problems installing tximport because of the version of R (3.3.1) and bioconductor (3.2) installed on our Rstudio server so (while I wait for the guys who maintain that to fix it) I figured I'd install it directly from the bioconductor github site:


This went just fine, but when I then try to load some data I get the following error:

tximp_files <- file.path(exp[[3]], "quant.sf")
txi <- tximport(files, type="salmon", tx2gene=t2g)

reading in files with read_tsv

Error in guess_header_(datasource, tokenizer, locale) : basic_string::resize
stop(structure(list(message = "basic_string::resize", call = guess_header_(datasource, tokenizer, locale), cppstack = NULL), .Names = c("message", "call", "cppstack"), class = c("std::length_error", "C++Error", "error", "condition")))
guess_header_(datasource, tokenizer, locale)
guess_header(ds_header, tokenizer, locale)
col_spec_standardise(data, skip = skip, comment = comment, guess_max = guess_max, col_names = col_names, col_types = col_types, tokenizer = tokenizer, locale = locale)
read_delimited(file, tokenizer, col_names = col_names, col_types = col_types, locale = locale, skip = skip, comment = comment, n_max = n_max, guess_max = guess_max, progress = progress)
readr::read_tsv(x, progress = FALSE, col_types = readr::cols())
tximport(files, type = "salmon", tx2gene = t2g)

The quant.sf data were generated with salmon 0.5.1, the file paths are correct and the files look like:

# salmon (mapping-based) v0.5.1
# [ program ] => salmon 
# [ command ] => quant 
# [ index ] => { /Projects/Arabidopsis_annotations/TAIR10/ensembl/release_26/Sequence/SalmonIndex/transcripts_index_salmon0.5.1_quasi }
# [ libType ] => { ISR }
# [ mates1 ] => { /Projects/salmon/0.5.1_bootstraps/Col-1/R1.fastq }
# [ mates2 ] => { /Projects/salmon/0.5.1_bootstraps/Col-1/R2.fastq }
# [ output ] => { /Projects/salmon/0.5.1_bootstraps/Col-1/TAIR10 }
# [ threads ] => { 8 }
# [ numBootstraps ] => { 100 }
# [ mapping rate ] => { 66.1342% }
# Name    Length    TPM    NumReads
ATMG00010.1    462    1.99872    65
ATMG00030.1    324    4334.17    73792.3
ATMG00040.1    948    11.5158    1018.12
ATMG00050.1    396    0.551205    13.7941
ATMG00060.1    542    564.503    23531.1
ATMG00070.1    573    318.244    14399.1


ADD COMMENTlink modified 19 months ago by James W. MacDonald48k • written 19 months ago by nschurch10
gravatar for James W. MacDonald
19 months ago by
United States
James W. MacDonald48k wrote:

If you already have mis-matched R and Bioconductor versions (R-3.3.1 should be running BioC 3.4, not Bioc 3.2), the fix is probably not to install devel versions of other packages.

Anyway, the error you are getting is from Hadley's readr package, which is not a Bioc package at all. You can always use read.table instead, or alternatively you could install the release version of readr, which is certainly working for me.

ADD COMMENTlink written 19 months ago by James W. MacDonald48k


Can you post the error when installing with biocLite()? This is the installation route that we commit to provide support for, going off that trail you're a bit on your own.

ADD REPLYlink written 19 months ago by Michael Love19k

sure, here it is:

> source("")
Bioconductor version 3.2 (BiocInstaller 1.20.3), ?biocLite for help
> biocLite("tximport")
Using Bioconductor 3.2 (BiocInstaller 1.20.3), R 3.3.1 (2016-06-21).
Installing package(s) ‘tximport’
Old packages: 'commonmark', 'mclust', 'R.filesets', 'abind', 'acepack',
  'backports', 'BBmisc', 'BH', 'bookdown', 'broom', 'car', 'checkmate',
  'chron', 'coda', 'colorspace', 'curl', 'data.table', 'DBI', 'DEoptimR',
  'deSolve', 'digest', 'e1071', 'effects', 'energy', 'expm', 'fastcluster',
  'fields', 'futile.logger', 'future', 'GGally', 'ggthemes', 'ggvis', 'git2r',
  'globals', 'gstat', 'httr', 'irlba', 'jsonlite', 'ks', 'Lahman', 'lambda.r',
  'lawstat', 'logcondens', 'loo', 'lubridate', 'MALDIquant', 'maps',
  'matrixStats', 'mclust', 'mime', 'mnormt', 'msm', 'mstate', 'multcomp',
  'multicool', 'mvtnorm', 'NeuralNetTools', 'nlme', 'numDeriv', 'nycflights13',
  'openssl', 'packrat', 'permute', 'plotly', 'plotrix', 'pracma', 'proto',
  'PSCBS', 'psych', 'quantreg', 'R.devices', 'R.filesets', 'R.oo', 'R.rsp',
  'R.utils', 'R6', 'rbokeh', 'Rcpp', 'RcppArmadillo', 'RcppEigen',
  'RefManageR', 'reshape', 'reshape2', 'rgl', 'robustbase', 'roxygen2',
  'rsconnect', 'RSNNS', 'RSQLite', 'rstan', 'rversions', 'scatterplot3d',
  'shiny', 'snow', 'sp', 'spacetime', 'spam', 'SparseM', 'StanHeaders',
  'stringi', 'stringr', '', 'tidyr', 'UpSetR', 'vegan', 'VGAM',
  'VGAMdata', 'XML', 'xml2', 'yaml', 'zoo', 'BiocInstaller', 'cluster',
  'codetools', 'foreign', 'lattice', 'Matrix', 'mgcv', 'nlme', 'survival'
Update all/some/none? [a/s/n]: 

Update all/some/none? [a/s/n]:
Warning message:
package ‘tximport’ is not available (for R version 3.3.1)
> library(tximport)
Error in library(tximport) : there is no package called ‘tximport’

It's just a warning because of the version mismatch and a failure to install. Looks like I'll just have to suck it up and wait for the sys-admins to upgrade the bioconductor version tot he right one...

ADD REPLYlink written 19 months ago by nschurch10

Yeah, Bioc 3.2 is from October 2015, so is 1.5 years out of date. tximport was introduced in Bioc 3.3 in April 2016. This information is on the package landing page:

"In Bioconductor since BioC 3.3"
ADD REPLYlink written 19 months ago by Michael Love19k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 180 users visited in the last hour