上面是最近邮件的提问，而这个问题在Bioconductor上最有人问过, https://support.bioconductor.org/p/104676/

Hi Guangchuang,

Is there a way to control annotations in the the annotatePeak function? For instance I would like to ignore "downstream" and consider those peaks "intergenic". Or for instance combine "1st exon" with "other exon" so there is one category of "exon".

Any help would be great, thanks for doing this!

-Andrew

所以呢在发行版中，早已支持，下面就是个简单的例子：

1. > require(ChIPseeker)

2. Loading required package: ChIPseeker

3. 
   

4. ChIPseeker v1.16.0  For help: https://guangchuangyu.github.io/ChIPseeker

5. 
   

6. If you use ChIPseeker in published research, please cite:

7. Guangchuang Yu, Li-Gen Wang, Qing-Yu He. ChIPseeker: an R/Bioconductor package for ChIP peak annotation, comparison and visualization. Bioinformatics 2015, 31(14):2382-2383

8. > ff = getSampleFiles()[[4]]

9. > options(ChIPseeker.ignore_1st_exon = TRUE)

10. > options(ChIPseeker.ignore_1st_intron = TRUE)

11. > options(ChIPseeker.ignore_downstream = TRUE)

12. > options(ChIPseeker.ignore_promoter_subcategory = TRUE)

13. > x = annotatePeak(ff)

14. >> loading peak file...                         2018-01-09 11:36:33

15. >> preparing features information...            2018-01-09 11:36:33

16. >> identifying nearest features...              2018-01-09 11:36:34

17. >> calculating distance from peak to TSS...     2018-01-09 11:36:34

18. >> assigning genomic annotation...              2018-01-09 11:36:34

19. >> assigning chromosome lengths                 2018-01-09 11:36:52

20. >> done...                                      2018-01-09 11:36:52

21. Warning message:

22. In loadTxDb(TxDb) :

23.  >> TxDb is not specified, use 'TxDb.Hsapiens.UCSC.hg19.knownGene' by default...

24. > x

25. Annotated peaks generated by ChIPseeker

26. 1331/1331  peaks were annotated

27. Genomic Annotation Summary:

28.     Feature  Frequency

29. 1   Promoter 58.7528174

30. 2     5' UTR  0.3005259

31. 3     3' UTR  2.1036814

32. 4       Exon  3.0052592

33. 5     Intron  7.9639369

34. 6 Intergenic 27.8737791

我们可以通过设置 options来控制，但这只做用于总结结果，用于可视化， annotatePeak总是会给出细节性的注释。

往期精彩

• 话题投票

• 掐架的额外收获

• 使用自己的数据来给进化树上色

• 纯文本可视化FASTQ质量数据

• ggtree直播PPT第一部分

• 你们想要的：QQ跑在Linux上

• 你的专业是那一只Pokémon神兽？

• s03 - Complementing a Strand of DNA