Question: Error in H5Fopen with kallisto .h5 files
0
gravatar for uso2015
7 weeks ago by
uso20150
uso20150 wrote:

I used the kallisto for RNA-seq, then import the result "*.h5" by the tximport tools, but I get the following error: Error in H5Fopen(file, "H5FACCRDONLY", native = native) : HDF5. File accessibilty. Unable to open file.

I am trying to import count data from Kallisto to DESeq2 using the tximport package following the instructions here. there is the same error.

sessionInfo()

R version 3.5.1 (2018-07-02) Platform: x8664-condacos6-linux-gnu (64-bit) Running under: CentOS Linux 7 (Core)

Matrix products: default BLAS/LAPACK: /software/anaconda2/lib/R/lib/libRblas.so

locale: [1] LCCTYPE=enUS.UTF-8 LCNUMERIC=C [3] LCTIME=enUS.UTF-8 LCCOLLATE=enUS.UTF-8 [5] LCMONETARY=enUS.UTF-8 LCMESSAGES=enUS.UTF-8 [7] LCPAPER=enUS.UTF-8 LCNAME=C [9] LCADDRESS=C LCTELEPHONE=C [11] LCMEASUREMENT=enUS.UTF-8 LC_IDENTIFICATION=C

attached base packages: [1] stats graphics grDevices utils datasets methods base

other attached packages: [1] rhdf5_2.26.2

loaded via a namespace (and not attached): [1] compiler3.5.1 tools3.5.1 Rhdf5lib_1.4.2

rhdf5 tximport kallisto • 99 views
ADD COMMENTlink modified 7 weeks ago • written 7 weeks ago by uso20150

This sounds like the issue at reported at https://github.com/pachterlab/sleuth/issues/120 and https://github.com/grimbough/Rhdf5lib/issues/15.

Are you using R / Bioconductor installed via bioconda? It would be helpful if you could update your question to include the output of the command sessionInfo() so we can see the versions of package that you're using. It's also useful to tag the question with rhdf5 since that's where the error is ultimately coming from, and then the maintainer of that package will be notified when post something.

ADD REPLYlink written 7 weeks ago by Mike Smith3.5k

yes, I installed R via bioconda. I saw the issues, the reason seems to the Rhdf5lib, but I still don't know how to solve it.

ADD REPLYlink modified 7 weeks ago • written 7 weeks ago by uso20150

Does updating to Rhdf5lib 1.4.3 resolve this? That version was created specifically to try and fix problems with the bioconda build system.

Can you also try running the command Sys.setenv(HDF5_USE_FILE_LOCKING = "FALSE") and then re-running whatever command is failing.

If that doesn't work, is it possible for you to share the .h5 file so I can try and dig a bit deeper into the HDF5 error message?

ADD REPLYlink modified 7 weeks ago • written 7 weeks ago by Mike Smith3.5k

I have solved the problem by reinstalling Rhdf5lib and rhdf5, thanks.

ADD REPLYlink written 6 weeks ago by uso20150
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: 135 users visited in the last hour