mzytnicki/HiCDOC

Chromosomes attribute missing

cmdoret opened this issue · 1 comments

Hello,

I tried running HiCDOC on my samples, as explained in the readme. HiCDOCDataSetFromCool runs successfully, however I obtain the following error when running filterSmallChromosomes:

Error in eval_tidy(xs[[j]], mask) : 
  no slot of name "chromosomes" for this object of class "HiCDOCDataSet"
Calls: filterSmallChromosomes ... %>% -> eval -> eval -> tibble -> tibble_quos -> eval_tidy
Execution halted

If I skip filterSmallChromosomes and run filterWeakPositions directly, I get the same error:

Error in eval_tidy(xs[[j]], mask) : 
  no slot of name "chromosomes" for this object of class "HiCDOCDataSet"
Calls: filterWeakPositions ... %>% -> eval -> eval -> tibble -> tibble_quos -> eval_tidy

Inspecting the structure of the HiCDOCDataSet object reveals a chromosome attribute, instead of chromosomes. Could that be the issue (the missing s) ?

> print(str(object))
Formal class 'HiCDOCDataSet' [package "HiCDOC"] with 5 slots
  ..@ inputPath   : chr [1:10] "data/tmp/cool_compres/PM125.cool" "data/tmp/cool_compres/PM55.cool" "data/tmp/cool_compres/PM54.cool" "data/tmp/cool_compres/PM124.cool" ...
  ..@ interactions: tibble [38,525,930 × 6] (S3: tbl_df/tbl/data.frame)
  .. ..$ chromosome: Factor w/ 57 levels "NC_003197.2",..: 3 3 3 3 3 3 3 3 3 3 ...
  .. ..$ position.1: int [1:38525930(1d)] 2880000 2880000 2880000 2880000 2880000 2880000 2880000 2880000 2880000 2880000 ...
  .. ..$ position.2: int [1:38525930(1d)] 2880000 3040000 3200000 3360000 3520000 3680000 3840000 4000000 4160000 4320000 ...
  .. ..$ value     : int [1:38525930(1d)] 40 34 4 9 9 10 8 7 8 7 ...
  .. ..$ replicate : Factor w/ 10 levels "PM122_2","PM123_2",..: 4 4 4 4 4 4 4 4 4 4 ...
  .. ..$ condition : Factor w/ 2 levels "1","2": 1 1 1 1 1 1 1 1 1 1 ...
  ..@ replicates  : chr [1:10] "PM125_1" "PM55_1" "PM54_1" "PM124_1" ...
  ..@ conditions  : chr [1:10] "1" "1" "1" "1" ...
  ..@ binSize     : NULL

I figured my issue; I made a mistake and used a HiCDOCDataSet object instead of HiCDOCExp !