Error in H5Fopen with kallisto .h5 files
0
0
Entering edit mode
uso2015 • 0
@uso2015-20314
Last seen 5.1 years ago

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

tximport kallisto rhdf5 • 2.3k views
ADD COMMENT
0
Entering edit mode

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 REPLY
0
Entering edit mode

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 REPLY
0
Entering edit mode

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 REPLY
0
Entering edit mode

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

ADD REPLY

Login before adding your answer.

Traffic: 794 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6