Question: Error in H5Fopen with kallisto .h5 files
0
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
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.

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.

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?