怎样使用 Juicer tools 的 dump 命令将.hic文件转换为交互矩阵matrix计数文件 (Windows)

创作日志： 万恶的生信...一个scHiC数据集没有提供处理好的计数文件，需要从.hic转换。Github一个个好长的文档看了好久才定位到 juicer tools 的dump命令，使用起来比想象中简单。

---

一、下载Juicer tools

注意：使用Juicer tools的前提是你的电脑里装了Java哦！

网址：Juicer tools jar包下载官网

我不知道有什么具体区别，下载的是第一个。下载之后也不需要有什么操作，直接就可以用。

二、dump命令讲解

1. juicer_tools dump 用法1------提取观测值或期望值

• 参数：

  observed/oe \] \[ NONE/VC/VC_SQRT/KR \] \[ hicFile(s) \] \[ chr1 \]\[:x1:x2\] \[ chr2 \]\[:y1:y2\] \[ BP/FRAG \] \[ binsize \] \[ outfile (可选)

• 参数解释：
  • observed/oe： 选择提取观测值 (observed) 或 观测/期望值 (oe)。
  • NONE/VC/VC_SQRT/KR： 选择归一化选项：

  NONE：无归一化

  VC：Vector Correction 归一化

  VC_SQRT：Square Root Vector Correction 归一化

  KR：Knight-Ruiz 归一化
  • hicFile(s)： 输入的 .hic 文件路径。
  • chr1:x1:x2： 第一个染色体及其范围（例如 chr1:0:100000）。
  • chr2:y1:y2： 第二个染色体及其范围（例如 chr2:0:100000），也可以是相同的染色体。
  • BP/FRAG： 选择单位：基对 (BP) 或 酶切片段 (FRAG)。
  • binsize： 分辨率，例如 10000 表示 10kb。
  • outfile： 输出文件路径（可选）。

2. juicer_tools dump 用法2------提取归一化或期望值

• 参数：

  norm/expected \] \[ NONE/VC/VC_SQRT/KR \] \[ hicFile(s) \] \[ chr \] \[ BP/FRAG \] \[ binsize \] \[ outfile (可选)

• 参数解释：
  • norm/expected：选择提取归一化值 (norm) 或 期望值 (expected)。
  • NONE/VC/VC_SQRT/KR：选择归一化选项（同上）。
  • hicFile(s)：输入的 .hic 文件路径。
  • chr：染色体。
  • BP/FRAG：选择单位：基对 (BP) 或 酶切片段 (FRAG)。
  • binsize：分辨率。
  • outfile：输出文件路径（可选）。

3. juicer_tools dump 用法3------提取染色质环loop或域domain信息

• 参数：

  loops/domains \] \[ hicFile URL\] \[ outfile (可选)

• 参数解释：
  • loops/domains：选择提取环 (loops) 或 域 (domains) 信息。
  • ：输入的 .hic 文件 URL。
  • outfile：输出文件路径（可选）。

三、使用实例

1. 打开Windows cmd
2. 输入java -jar 你的juicer_tools jar包安装路径 dump 以上所需参数 即可完成转换

拿我的举个例子：

最后在我的指定路径中生成了extract_matrix.txt文件：

打开内容是这样的：