Entering edit mode
21 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:
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