该问题已被锁定!
1
关注
1093
浏览

computeMatrix画图问题

为什么被折叠? 0 个回复被折叠
chatGPT机器人 机器人 用户来自于: IANA
2023-06-23 12:54
您好!根据您的问题描述,您的思路是正确的。以下是一些可能有用的代码示例,可以帮助您使用computeMatrix创建热图并对peaks进行注释。 首先,您需要一些bed文件来指定您要在热图中显示的区域。假设您已经有了三个bed文件,分别命名为unique.bed,common.bed和all.bed。 接下来,您可以使用以下命令使用computeMatrix来创建热图: ```bash computeMatrix reference-point --referencePoint center -S signal.bw -R unique.bed common.bed all.bed -a 5000 -b 5000 -o matrix.gz ``` 这将创建一个名为matrix.gz的gzip压缩文件,其中包含对unique.bed,common.bed和all.bed中的区域的信号进行计算后生成的矩阵。这里使用了参考点的中心,并且信号文件为signal.bw。您可以根据需要更改这些参数。 现在,您可以使用以下命令将该矩阵绘制为热图: ```bash plotHeatmap -m matrix.gz -out heatmap.pdf --colorList white,blue --heatmapHeight 7 --heatmapWidth 10 --refPointLabel "Center" ``` 这将创建一个名为heatmap.pdf的PDF文件,其中包含从matrix.gz生成的热图。这里使用了白色和蓝色的颜色列表,并且热图的高度为7,宽度为10。您可以根据需要更改这些参数。 最后,您可以使用以下命令将peaks注释到热图中: ```bash annotatePeaks.pl matrix.gz hg19 -size 5000 -detailedAnnotation > annotated_peaks.txt ``` 这将创建一个名为annotated_peaks.txt的文本文件,其中包含对热图中每个peak的详细注释。这里使用了hg19作为参考基因组,并且peak的大小为5000bp。您可以根据需要更改这些参数。 希望这些代码对您有所帮助!如果您有任何问题,请随时问我。
孟浩巍 超级管理员 用户来自于: 北京市
2023-06-25 16:13

G图可以,但是H,和I图明显是富集分析的结果。

应该是找到peak以后又和gene做了关联,然后做了gene的富集分析。

关于作者

lyx_2023 注册会员

这家伙很懒,还没有设置简介

问题动态

发布时间
2023-06-23 12:50
更新时间
2023-06-25 16:13
关注人数
1 人关注

相关问题

服务器linux系统上安装Rstudio和相应配置问题?
求助,运行软件时遇到perl问题
生存分析KM-plot问题
基因表达芯片的分类问题
关于hub基因的问题
染色体重叠区域问题
DoubletfFinder Pk值相关计算的问题
cox单因素分析问题
axel,wget下载数据出错的问题,没头绪
想问一个RNA-Seq的流程问题

推荐内容

请问cytoscape是怎么构建网络的?
纤维二糖功能
安装R包sf报错
【求助】怎么能画出子网络图
GWAS和WGS的选择消除分析有什么区别呀?
bowtie2比对报错((ERR): bowtie2-align exited with value 1)
不同样品中寻找特异的OTU
how to speed up following code 
染色体号是罗马数字怎么写sh循环
23年生信NCCL室间质评
All Rights Reserved Powered BY WeCenter V4.1.0 © 2024