[BioC] mypeaklist in ChIPpeakAnno

Zhu, Julie Julie.Zhu at umassmed.edu
Mon Aug 16 16:34:14 CEST 2010


Dear Ramu,

Glad it worked. Thanks for letting me know!

Kind regards,

Julie


On 8/16/10 5:44 AM, "Ramu Chenna" <ramu.chenna at biotec.tu-dresden.de> wrote:

> Dear Julie
> 
> It works, thanks.  Basically what we wanted to do is, from the peak
> list, we want to know
> how many peaks are in promotor regions and further classify them into
> cpg island promotors and
> non-cpg island promoters.
> 
> Note:  The pdf documentation is not done correct there are lines that are
> truncated,  (may be a problem in setting the margin in latex?)
> 
> Best,
> Ramu
> 
> Zhu, Julie wrote:
>> Dear Ramu,
>> 
>> You can use read.table to read the data in bed file format from MACS, then
>> use BED2RangedData to convert the data to RangedData.
>> 
>> Library(ChIPpeakAnno)
>> X = read.table("mypeak.bed",sep="\t")
>> MyPeakList = BED2RangeData(X)
>> 
>> ?BED2RangedData
>> 
>> BTW, I am on vacation so may not be able to check email frequently.
>> 
>> Kind regards,
>> 
>> Julie
>> 
>> 
>> On 8/13/10 3:44 AM, "Ramu Chenna" <ramu.chenna at biotec.tu-dresden.de> wrote:
>> 
>>   
>>> Dear Julie,
>>> 
>>> The link did not work, but I found the ppt.  Thanks.
>>> 
>>> I am not a R geek, and I just need to know how to read MACS output
>>> as say 'myPeakList' as RangedData. It would be nice to know there is already
>>> a R function doing this.
>>> 
>>> Regards,
>>> Ramu
>>> 
>>> Zhu, Julie wrote:
>>>     
>>>> Dear Ramu,
>>>> 
>>>> It is a RangedData. Please take a look at the 2nd example given in the
>>>> Bioc2010 meeting in Seattle.
>>>> 
http://www.bioconductor.org/help/course-materials/2010/BioC2010/BioC2010_Ch>>>>
I
>>>> PpeakAnno.pdf
>>>> 
>>>> 
>>>>   
>>>> 
>>>>   
>>>>       
>> 
>> 
>>   
> 



More information about the Bioconductor mailing list