[BioC] Query neighboring genes
Zhu, Lihua (Julie)
Julie.Zhu at umassmed.edu
Fri May 18 01:05:55 CEST 2012
Please set parameters output = "both", maxgap = 100000, select="all"
followed by filtering the results further using a combination of
insideFeature and shortestDistance. For detailed parameter setting, please
type help(annotatePeakInBatch) in a R session.
For example, the following will return peaks located within 5kb downstream
of a gene.
annotatedPeak$insideFeature == "downstream" &
!is.na(annotatedPeak$shortestDistance) & annotatedPeak$shortestDistance
On 5/17/12 4:17 PM, "Yadav Sapkota" <ysapkota at ualberta.ca> wrote:
> I tried this function and it works fine for me however, I could not find any
> option to define the upstream and downstream size while looking for
> neighboring genes.
> For example, I wanted to only query genes that are located 100 KB upstream and
> 50 KB downstream from a chromosomal range. Any hint would be appreciated.
> Below is the code I am using for demo:
> annotatedPeak = annotatePeakInBatch (myPeakList[1:6,], AnnotationData =
> On Thu, May 17, 2012 at 10:52 AM, Zhu, Lihua (Julie) <Julie.Zhu at umassmed.edu>
>> You could try annotatePeakInBatch in ChIPpeakAnno package.
>> Best regards,
>> On 5/17/12 12:40 PM, "Yadav Sapkota" <ysapkota at ualberta.ca> wrote:
>>> Is there any tool or package that can query neighboring genes (lets say 100
>>> KB upstream and downstream) of a specific chromosomal range (chr1:100-200)?
>>> I have thousands of these chromosomal ranges.
>>> You help will be greatly appreciated.
>>> Yadav Sapkota
>>> Uni of Alberta
>>> [[alternative HTML version deleted]]
>>> Bioconductor mailing list
>>> Bioconductor at r-project.org
>>> Search the archives:
More information about the Bioconductor