[BioC] Query neighboring genes
Zhu, Lihua (Julie)
Julie.Zhu at umassmed.edu
Fri May 18 01:05:55 CEST 2012
Yadav,
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[!is.na(annotatedPeak$insideFeature) &
annotatedPeak$insideFeature == "downstream" &
!is.na(annotatedPeak$shortestDistance) & annotatedPeak$shortestDistance
<=5000,]
Best regards,
Julie
On 5/17/12 4:17 PM, "Yadav Sapkota" <ysapkota at ualberta.ca> wrote:
> Hi,
>
> 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:
>
> source("http://bioconductor.org/biocLite.R")
> biocLite("ChIPpeakAnno")
> library("ChIPpeakAnno")
> data(myPeakList)
> data(TSS.human.GRCh37)
> annotatedPeak = annotatePeakInBatch (myPeakList[1:6,], AnnotationData =
> TSS.human.GRCh37)
>
> --Yadav
>
>
>
>
> On Thu, May 17, 2012 at 10:52 AM, Zhu, Lihua (Julie) <Julie.Zhu at umassmed.edu>
> wrote:
>> Yadav,
>>
>> You could try annotatePeakInBatch in ChIPpeakAnno package.
>>
>> Best regards,
>>
>> Julie
>>
>>
>> On 5/17/12 12:40 PM, "Yadav Sapkota" <ysapkota at ualberta.ca> wrote:
>>
>>> Hi,
>>>
>>> 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.
>>>
>>> Regards,
>>> Yadav Sapkota
>>> Uni of Alberta
>>>
>>> [[alternative HTML version deleted]]
>>>
>>> _______________________________________________
>>> Bioconductor mailing list
>>> Bioconductor at r-project.org
>>> https://stat.ethz.ch/mailman/listinfo/bioconductor
>>> Search the archives:
>>> http://news.gmane.org/gmane.science.biology.informatics.conductor
>>
>>
>
>
More information about the Bioconductor
mailing list