bioconductor入门第二弹--GRanges概述

Granges是genomic ranges的简称，它是bioconductor中的R包存储genomic intervals的数据结构，快速且有效，是想要深入分析数据所必须掌握的。它的重要性类似于bed注释文件，这么一比较这玩意的重要性就显而易见了

做基因组学分析时，最重要的就是intervals或者是intervals的集合。那么什么是intervals呢？它是基因组上的一些实体，我也想不出来一个合适的名词去翻译它，你可以粗鲁的称呼它为对象。

虽然翻译不出来，但是我举几个例子，一定能让你明白它究竟是个什么玩意

如下图，这是UCSC browser上随便挑选出的一段序列的的截图，我们能从上面看到DNA clusters，SNPs ，promoters， genes，这些东西都可以称之为intervals，甚至一个测序read被map到基因组上的位置也可以称之为一个intervals

许多基因组工作涉及到各种intervals的各种相互关系，如哪些promoters含有SNP位点，某个TF结合位点又与哪些promoter重合，哪些基因被测序reads覆盖。说了这么多intervals，那么它和我们要介绍的GRanges有什么关系呢？

上图是一个GRanges的结构，其中的A1，A2和A3是它的三个intervals，所以我们可以理解，intervals是GRanges对象的横向最小单位。

从纵向上看，GRanges对象由三列组成，第一列是染色体名称，第二列是intervals所在的具体位置（这个位置信息是IRanges对象所提供的），第三列则是该位置所在链的方向

除此之外，我们还能看到左下角是染色体的信息，比如seqlengths告诉了我们chr1的长度，若未定义则为NA

我们可以使用R来处理GRanges对象，但除此之外我们还可以使用bedtools处理intervals之间的相互关系。比如，我最常用的是intersectBed，它可以很快的帮助我知道我所需要的信号处于哪些位置。

IRanges是包含intervals的向量， IRanges和GRanges类似，但是GRanges具有染色体编号和正负链信息，IRanges为GRanges提供了基础功能，可以进行intervals的简单代数运算

IRanges有三列，起始，终止以及宽度。当我们给定了任意两列信息就可以由此推断剩下一列的信息。

下面使用IRanges函数构建IRanges对象

这个个Iranges对象包含三个不同的intervals，每一个intervals都称之为一个range(范围区间)，这里的ir1有三个ranges

start(), end(), width()这三个方法可以获得IRanges对象的起始列，终止列，以及宽度那一列的信息。IRanges对象的每个ranges也可以有自己的名字，如

因为是向量，因此IRanges只有一个维度，不能使用dim()方法，但是我们可以使用length()获取其长度。我们可以像操作向量一样，对IRanges对象进行操作，如

甚至可以利用c()将两个IRanges对象合并

有时候在阅读R包的资料时，你也许会看到NormalIRanges对象，它是由使用reduce函数对IRanges对象处理得到。一个NormalIRanges，他其中的ranges（就是每一 行 数据）有以下几个特点：

用文字可能比较抽象，用一张图应该能帮助大家更好的理解。上面的表示我们的其实Iranges对象，而下面的则是我们使用reduce函数处理后所得到的normal Ranges对象。

下面是reduce()函数的对象：

disjoin函数和reduce函数恰好相反，找出完全不重合的intervals，将重合的intervals进行分割。

还有一些方法可以用于批量转换IRanges对象的坐标位置，比如：

除此之外，还有一些方法对于上述三个方法做了进一步的拓展，由于篇幅原因就不在这里介绍，大家感兴趣可以自己搜搜看

reseize()函数可能是一个处理数据很有用的方法，它有一个fix参数，可以固定原始ranges中不变的列，具体用法如下:

既然IRanges是一个向量，那么肯定可以对IRanges进行集合操作，我们通常使用unnion取并集，intersect取交集：

这些函数，还有一些变种，punion(), pintersect(), psetdiff(), pgap()。它们与R基础函数中的pmax很相似，很少被用到，因此不做讲解

对测序数据进行下游分析时，经常需要寻找overlap，使用findOverlaps()函数寻找两个IRanges对象中的overlap。它的返回值是一个Hists对象，它可以描述两个Iranges之间重叠关系，是一个两列的矩阵，由IRanges的索引组成。可以由queryHist()和subjectHits()得到其索引信息。

findOverlaps的参数很多，是否有最小重叠，是否完全重叠等。我自己也记不住，所以就不在这里介绍了

countOverlaps()函数可以返回overlap的数目，该功能速度更快，占用内存更小，因为不需要记录ranges重叠的具体信息，仅仅考虑重叠的个数

nearest()函数用于寻找两个IRanges中，第一个IRanges对象中的ranges最接近第二个IRanges对象中的哪个ranges。

如果你用过Y叔的R包ChIPseeker，你可能会对这个功能比较熟悉，它将peak注释到基因的一个方法就是通过找该peak与哪个基因最接近

类似的函数还有precede，follow等，实在是太多了，而且不常用，就不一一介绍了……

相信看到这里你们对于IRanges对象已经有了一个比较深的认识，对于GRanges也有了大致的了解。下一期不出意外，我会详细写写GRanges对象的操作，希望能带你进一步认识GRanges对象