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

登录

该问题已被锁定！

2: 关注

1145: 浏览

用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 人关注

相关问题

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

Aspera数据下载报错: 1224 浏览 2 关注 2 回答 0 评论

在很多分析表达谱芯片的教程中，使用limma 包寻找差异表达基因，都是先找到差异表达探针，再把探针注释为基因id，为何不先进行注释再计算差异表达呢？这两种方法有什么区别呢？: 808 浏览 2 关注 1 回答 0 评论

分析CRISPR 高通量筛选数据: 969 浏览 1 关注 2 回答 0 评论

GWAS 数据存储网站: 454 浏览 1 关注 0 回答 0 评论

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

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

聚类数据标准化: 926 浏览 2 关注 2 回答 0 评论

人源细胞系ATAC-seq比对率较低，显示支原体污染，如何处理？: 813 浏览 2 关注 1 回答 0 评论

转录组数据样本聚类结果不理想: 696 浏览 2 关注 1 回答 0 评论

推荐内容

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

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

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

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

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

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

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

问 GATK和Samtools联合寻找变异: 1051 浏览 2 关注 2 回答 0 评论

问【坐标北京】近期准备做一些RNA-seq，想请问一下大家哪些公司比较靠谱？-w-: 753 浏览 2 关注 1 回答 0 评论

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

All Rights Reserved Powered BY WeCenter V4.1.0 © 2024