使用GREAT对peak进行功能注释

欢迎关注”生信修炼手册”!

GREAT是一款peak区间进行基因注释的工具,除了给出peak对应的基因外,还集成了多种基因的功能分析,网址如下

http://great.stanford.edu/public/html/index.php

目前该在线工具只支持以下几个物种

  1. Human

  2. Mouse

  3. Zebrafish

在使用时,还需要注意对应的基因组版本。用法比较简单,选择对应的基因组版本,然后上传对应的BED格式的peak文件即可,示意如下

结果展示如下,给出了peak关联基因的个数和TSS距离的频数分布柱状图

除此之外,还给出以下多种基因的功能分析

  1. GO Molecular Function

  2. GO Biological Process

  3. GO Cellular Component

  4. Mouse Phenotype

  5. Human Phenotype

  6. Disease Ontology

  7. MsigDB Cancer neighborhood

  8. Placenta Disorders

  9. PANTHER Pathway

  10. BioCyc Pathway

  11. MsigDB Pathway

  12. MGI Expression

  13. MgisDB Perturbation

不同于传统的费舍尔精确检验,功能富集分析的p值是基于二项分布的计算得到的,计算过程如下所示

以GO中MF这一类别的功能注释为例,示意如下

通过GREAT可以方便的对peak关联的基因功能进行探究。

·end·

—如果喜欢,快分享给你的朋友们吧—

扫描关注微信号,更多精彩内容等着你!

  • 1
    点赞
  • 5
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
好的,这是一个很好的问题。使用Python语言可以很容易地实现读者管理功能的开发。下面是一个简单的示例代码,它可以让你开始开发你自己的读者管理系统。 ```python class Reader: def __init__(self, name, age, gender, books=[]): self.name = name self.age = age self.gender = gender self.books = books def __str__(self): return f"{self.name}, {self.age}, {self.gender}, {len(self.books)} books" class ReaderManager: def __init__(self): self.readers = [] def add_reader(self, reader): self.readers.append(reader) def remove_reader(self, reader): self.readers.remove(reader) def search_reader(self, name): for reader in self.readers: if reader.name == name: return reader return None def get_all_readers(self): return self.readers if __name__ == "__main__": reader1 = Reader("John", 25, "Male", ["Python for Dummies"]) reader2 = Reader("Jane", 30, "Female", ["The Catcher in the Rye", "To Kill a Mockingbird"]) reader3 = Reader("Bob", 35, "Male") reader_manager = ReaderManager() reader_manager.add_reader(reader1) reader_manager.add_reader(reader2) reader_manager.add_reader(reader3) print("All readers:") for reader in reader_manager.get_all_readers(): print(reader) reader_manager.remove_reader(reader2) print("All readers after removing reader2:") for reader in reader_manager.get_all_readers(): print(reader) reader4 = Reader("Alice", 20, "Female", ["The Great Gatsby"]) reader_manager.add_reader(reader4) print("Searching for reader 'Bob':") bob = reader_manager.search_reader("Bob") if bob: print(bob) else: print("Reader not found.") ``` 这个示例代码定义了两个类:`Reader`和`ReaderManager`。`Reader`类表示一个读者,其中包含读者的姓名、年龄、性别和所借阅的书籍列表。`ReaderManager`类表示一个读者管理器,它可以添加、删除、搜索和获取所有读者。在`if __name__ == "__main__":`语句中,我们创建了一些读者并将它们添加到`ReaderManager`中,然后演示了如何获取所有读者、删除一个读者和搜索一个读者。 当然,这只是一个示例代码,你可以根据自己的需求来修改和扩展它。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值