Skip to content

Issue with bins in imported single resolution .cool data when trying to use getCompartments() #4

Description

@meisterpeter

I like OHCA (introduced to me by Martial).
I imported single resolution (50k) Hi-C matrix produced by HiC-Pro, converted into .cool format. Import goes well, bins are correctly imported.

`bins(hic)
GRanges object with 2009 ranges and 2 metadata columns:
                      seqnames            ranges strand |    bin_id     weight
                         <Rle>         <IRanges>  <Rle> | <numeric>  <numeric>
            I_1_50000        I           1-50000      * |         0 0.00979230
       I_50001_100000        I      50001-100000      * |         1 0.00891737
      I_100001_150000        I     100001-150000      * |         2 0.00703530
      I_150001_200000        I     150001-200000      * |         3 0.00997862
      I_200001_250000        I     200001-250000      * |         4 0.00947075
                  ...      ...               ...    ... .       ...        ...
  X_17500001_17550000        X 17500001-17550000      * |      2004 0.00902023
  X_17550001_17600000        X 17550001-17600000      * |      2005 0.01092203
  X_17600001_17650000        X 17600001-17650000      * |      2006 0.00740209
  X_17650001_17700000        X 17650001-17700000      * |      2007 0.00957678
  X_17700001_17718942        X 17700001-17718942      * |      2008        NaN`

However, when trying to call compartments with getCompartments(hic, genome=genome), even when defining the genome, the function returns
Erreur : BiocParallel errors 1 remote errors, element index: 1 6 unevaluated and other errors first remote error: Error in rhdf5::h5read(file, name = path, index = list(idx), ..., bit64conversion = "double"): Object '/resolutions/50000/bins' does not exist in this HDF5 file.

Is there something I can debug myself or are the bins not passed correctly?

Activity

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Projects

    No projects

      Milestone

      No milestone

      Relationships

      None yet

      Development

      No branches or pull requests

      Issue actions