怎样使用 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文件:

打开内容是这样的:

相关推荐
Red Red2 个月前
GEO数据库提取疾病样本和正常样本|GEO数据库区分疾病和正常样本|直接用|生物信息|生信
开发语言·数据库·笔记·学习·r语言·c#·生物信息
tRNA做科研2 个月前
Bio-Linux-shell详解-2-基本Shell命令快速掌握
linux·运维·服务器·生物信息·计算生物学
Red Red2 个月前
GEO数据的下载和处理|GEO数据转换为Gene symbol|GEO注释文件提取symbol|查看样本标签|查看GEO数据疾病或正常|生物信息基础
数据库·笔记·学习·r语言·生物信息·geo数据库
图灵生信3 个月前
python脚本:输入基因名,通过爬虫的方式获取染色体上的location。
生物信息
qq_214782614 个月前
新工具:轻松可视化基因组,部分功能超IGV~
机器学习·信息可视化·数据分析·生物信息
徐洲更hoptop5 个月前
使用ShinyCell展示你的单细胞数据
生物信息
qq_214782615 个月前
Python散点图矩阵代码模版
数据库·python·mysql·信息可视化·数据挖掘·数据分析·生物信息
陆沙6 个月前
gnuplot使用实例
开发语言·生物信息·gnuplot
丸丸丸子w7 个月前
生物信息学 GO、KEGG
生物信息·kegg