|
@@ -173,9 +173,10 @@ CALDER(contact_file_hic=hic_file,
|
|
|
|
|
|
### Example four: run CALDER on other genomes
|
|
|
```
|
|
|
+## prepare feature_track
|
|
|
library(rtracklayer)
|
|
|
-feature_track = import('ENCFF934YOE.bigWig') ## from ENCODE
|
|
|
-feature_track = data.table::as.data.table(feature_track)[, c(1:3, 6)]
|
|
|
+feature_track_raw = import('ENCFF934YOE.bigWig') ## from ENCODE https://www.encodeproject.org/files/ENCFF934YOE/@@download/ENCFF934YOE.bigWig
|
|
|
+feature_track = data.table::as.data.table(feature_track_raw)[, c(1:3, 6)]
|
|
|
|
|
|
chrs = c(21:22)
|
|
|
contact_file_dump = as.list(system.file("extdata", sprintf("mat_chr%s_10kb_ob.txt.gz", chrs),
|