有小伙伴问为什么R里循环出不来图？
其实看我公众号的各位读者们，你们不应该有这个问题，如果有，那么证明我的文章没有好好看！有两篇很重要的文章《扪心自问，meme几何？》和《树变图，图变树？ ？》，特别是第二篇，和这个问题息息相关。
首先假设我们有一个图：
require(ggplot2)
p = qplot(1:10)
因为在R里，所有都是对象，而在终端里你输入一个对象回车，R会去找相应的方法去打印出来，如果是S3对象，会去寻找print方法，而如果是S4会去找show方法。
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
抛出一个话题，R是不是有个大坑，在自动化脚本中(或者在循环中)，某些脚本执行会不出结果，特别是作图脚本。ps:分步执行可以出图
其实看我公众号的各位读者们，你们不应该有这个问题，如果有，那么证明我的文章没有好好看！
有两篇很重要的文章《扪心自问，meme几何？》和《树变图，图变树？》，特别是第二篇，和这个问题息息相关。
首先假设我们有一个图：
require(ggplot2)
p = qplot(1:10)
因为在R里，所有都是对象，而在终端里你输入一个对象回车，R会去找相应的方法去打印出来，如果是S3对象，会去寻找print方法，而如果是S4会去找show方法。
所以你在终端里打p，其实等同于你打print(p)，它其实调用了ggplot2:::print.ggplot，
你打plot(p)也可以，如果你看过代码的话，画图这个动作是定义在plot.ggplot方法里的，而print.ggplot <- plot.ggplot。
首先基因下游是什么？
我们知道上游很重要啊，因为可能会调控转录，但注释的时候，没有上游这个东西，为什么呢？因为转录起始位点TSS的上下游被定义为promoter，所以啊上游被包括在promoter中，也就没有上游这个category了。
我在2018年4月份的一条推特文，https://twitter.com/guangchuangyu/status/988672880551518210，可以把UpSetR变成ggplot然后可以各种玩，比较我在推特文中所展示的，图上嵌入venn图、使用cowplot拼图等。然而过了大半年，我的pull request并没有被作者接收，也就是说这功能还不存在于官方的版本中。
如果有兴趣的话，应该安装我fork的版本。
devtools::install_github("GuangchuangYu/UpSetR")
生物狗喜欢画barplot，但毕竟只有均值和标准误差信息量太低，万一有outliers呢？万一不是正态分布呢？所以就吐槽了又吐槽，但大家还是乐此不疲，因为大家都在用，因为很多人只会画barplot，所以产生各种反barplot运动，并呼吁大家使用boxplot，boxplot有四分位数的统计量，可以反应outlier和数据的分布，比barplot高得多了。
然而事情是不断在进化的，我们的手段是要跟上潮流的，现在的潮流就是除了要有统计量，还要有原始数据，甚至于有统计推断，全集中一图中，揭开数据的红盖头，给大家呈现全方位无死角的面貌。
