how to get genomic range for hi-c datasets in hdf5 file
0
0
Entering edit mode
13 months ago
rheab1230 ▴ 140

Hello everyone,

I have hi-c contact map in hdf5 file. I stored the data for chr9 in a variable. But it doesn't have genomic range. Does anyone know how to add colnames and rownames with respect to genomic range for the matrix.

BiocManager::install("rhd5f")
library(rhd5f)
h5f = H5Fopen("GSE77565_FBD_IC-heatmap-chr-40k.hdf5")
h5f
HDF5 FILE 
        name /
    filename 

        name       otype dclass         dim
0  (0, 0)    H5I_DATASET FLOAT  6232 x 6232
1  (1, 1)    H5I_DATASET FLOAT  6080 x 6080
2  (10, 10)  H5I_DATASET FLOAT  3376 x 3376
3  (9, 9)  H5I_DATASET FLOAT  3347 x 3347
4 _self_key H5I_DATASET STRING ( 0 )      

chr9 <- as.data.frame(h5read("GSE77565_FBD_IC-heatmap-chr-40k.hdf5", "(9, 9)"))

My data look like this:

enter image description here

In this case chr9 has interacting frequency but there is no rownames or colnames mentioning the genomic range which is interacting. Does anyone know how to incorporate that information

Thank you

R hi-c contact-map genomic-ranges • 307 views
ADD COMMENT

Login before adding your answer.

Traffic: 2629 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