用atac-seq数据计算的TSS enrichment score - 做专业的生信问答网站

登录

该问题已被锁定！

2: 关注

1141: 浏览

用atac-seq数据计算的TSS enrichment score

RNA_seq

我用atac-seq数据的bigwig，用computeMatrix去计算TSS score,得到的图的纵坐标不是从1开始的？按照tss score的公式，这里应该是从1开始的呀，最小值一定是1。如图所示：

好问题 0 评论收藏举报

2 回答

chatGPT机器人机器人用户来自于: IANA
2023-05-17 14:28

首先，TSS enrichment score是用来衡量某个基因的转录起始位点（TSS）在ATAC-seq数据中的富集程度的。计算TSS enrichment score时，通常会将基因的TSS定义为起点，然后计算在一定距离范围内（例如2kb）ATAC-seq信号的富集程度。

对于问题中提到的情况，纵坐标不是从1开始的，这可能是由于computeMatrix计算过程中使用的是默认参数，将所有信号值都转换为正数。因此，得到的图中信号值的最小值可能不是1。

为了确保TSS enrichment score的计算正确，可以通过以下步骤进行：

将bigwig文件中的信号值转换为正数。
将信号值归一化，使其范围在0到1之间。
将信号值按照距离TSS的远近进行排序。
计算在一定距离范围内的信号值的平均值，并将其除以全局平均信号值，得到TSS enrichment score。

通过以上步骤，可以确保得到的TSS enrichment score是正确的，同时也可以保证纵坐标从1开始。

赞同 0 0评论

孟浩巍超级管理员用户来自于: 北京市
2023-05-17 14:46

ChatGPT回答得挺好……

赞同 0 0评论

关于作者

: kkkkkk 注册会员
找一个可以记录问题的地方

2: 回答

1: 文章

9: 问题

问题动态

发布时间: 2023-05-17 14:18

更新时间: 2023-05-17 14:46

关注人数: 2 人关注

相关问题

如何计算碱基偏好性？: 664 浏览 3 关注 1 回答 0 评论

来自不同project的RNA-Seq数据可以直接合并分析吗？: 730 浏览 2 关注 1 回答 0 评论

请问芯片数据分析的一般流程及涉及的常用算法？: 1043 浏览 2 关注 3 回答 0 评论

如何针对affy表达谱芯片计算结果分析不同样本之间是否表达差异: 587 浏览 2 关注 1 回答 0 评论

Hichip得到loop数据如何注释？: 689 浏览 2 关注 2 回答 0 评论

用log2(fold_change)数据做热图，遇到想要表达的目的基因log2(fold_change)值为inf或-inf时，怎么办？: 572 浏览 2 关注 1 回答 0 评论

如何对特征数量少的空间蛋白组数据进行细胞聚类？: 1316 浏览 2 关注 1 回答 0 评论

ciriRNA表达量如何计算？（更）: 628 浏览 2 关注 1 回答 0 评论

R软件里Bray-Curtis指数计算: 426 浏览 1 关注 0 回答 0 评论

TCGA数据下载: 648 浏览 2 关注 1 回答 2 评论

推荐内容

问 seurat中如何通过坐标/细胞名删掉细胞: 1466 浏览 2 关注 1 回答 0 评论

问 RIP-seq有没有类似ChIP-seq的blacklist region合集？: 1432 浏览 2 关注 2 回答 0 评论

问转录组: 1202 浏览 2 关注 1 回答 0 评论

问请问RNA-seq采用poly A(+)策略建库，处理数据时若不去除rRNA会对后续分析有何影响？: 770 浏览 2 关注 1 回答 0 评论

问多个处理组的RNAseq分析中关于counts表格合并的问题: 1186 浏览 2 关注 1 回答 0 评论

问公司双端测序的数据R1R2处理: 961 浏览 2 关注 1 回答 0 评论

问 RNA-seq，样本中存在同一基因对应不同的FPKM值: 692 浏览 2 关注 1 回答 0 评论

问 RNA-seq分析，rlog需要生物学重复才行？有别的Normalization方法吗？: 875 浏览 1 关注 0 回答 0 评论

问 R语言作火山图问题: 764 浏览 0 关注 1 回答 0 评论

问 RNA-seq不同样本多个生物学重复不同处理条件下的如何找差异基因: 769 浏览 2 关注 2 回答 0 评论

All Rights Reserved Powered BY WeCenter V4.1.0 © 2024