[BioC] Question about ChIPpeakAnno
Zhu, Lihua (Julie)
Julie.Zhu at umassmed.edu
Wed Jun 6 21:26:48 CEST 2012
Sorry, Petros!
It is - instead of _.
http://www.bioconductor.org/help/course-materials/2011/BioC2011/LabStuff/ChI
PpeakAnno-BioC2011.pdf
For your other questions in a separate email, you need to give sessionInfo()
and tell us whether you are able to get the example work. Thanks!
Best regards,
Julie
On 6/6/12 3:22 PM, "Petros Kolovos" <p.kolovos at erasmusmc.nl> wrote:
> Dear Julie
>
> the
> http://www.bioconductor.org/help/course-materials/2011/BioC2011/BioC2011_ChIPp
> eakAnno.pdf
> doesn't work. Could you please check it?
>
> Best regards
> Petros
>
>
>> Petros,
>>
>> Totaltest specifies the total number of tests performed to obtain the list
>> of peaks.
>>
>> This is one of the frequently asked questions and many have contributed
>> wisdoms to address this question. Could you please take a look at the
>> slides
>> at
>> http://www.bioconductor.org/help/course-materials/2011/BioC2011/BioC2011_ChI
>> PpeakAnno.pdf, esp. the faq slide near the end?
>>
>> Please do cc Bioconductor list (bioconductor
>> <bioconductor at stat.math.ethz.ch>) for sharing. Thanks!
>>
>> Best regards,
>>
>> Julie
>>
>>
>> On 6/6/12 11:02 AM, "Petros Kolovos" <p.kolovos at erasmusmc.nl> wrote:
>>
>>> Dear Julie,
>>>
>>> I would like to ask you also something else.
>>>
>>> In the makeVennDiagrams function what do you mean by total tests?
>>>
>>> For example i have the following code
>>>
>>>> venn=makeVennDiagram(RangedDataList (pr, intron),NameOfPeaks =
>>>> c("peaks","intron"),maxgap = 0, totalTest = 320000)
>>>
>>> PR=my peaks has 50000 peaks
>>> intron= table aquired from UCSC with all the exons 250000
>>> peaks/coordinates
>>>
>>> So how much should I put the totalTest
>>>
>>> Is there a problem if I make a venn diagram of peaks vs introns?
>>>
>>> Thanks
>>>
>>> Best regards
>>> Petros
>>>
>>>
>>>
>>>> Dear Petros,
>>>>
>>>> You can safely ignore the warning message and proceed. I will modify
>>>> the
>>>> code for the next release to replace multiple with select internally.
>>>>
>>>> Best regards,
>>>>
>>>> Julie
>>>>
>>>>
>>>> On 6/6/12 4:04 AM, "Petros Kolovos" <p.kolovos at erasmusmc.nl> wrote:
>>>>
>>>>> Dear Dr Julie Zhu,
>>>>>
>>>>> Good morning
>>>>>
>>>>> My name is Petros Kolovos and I am a PhD student at Rotterdam,
>>>>> Netherlands.
>>>>>
>>>>> I was using your library in order to annotate some chip peaks and in
>>>>> order
>>>>> to make some venn diagrams.
>>>>>
>>>>> But I faced some problems with the venn diagrams
>>>>>
>>>>> Here is what I am doing
>>>>>
>>>>>> venn=makeVennDiagram(RangedDataList (pr, Anno_exon),NameOfPeaks =
>>>>> c("peaks","exon"),maxgap = 0, totalTest = 500000)
>>>>> Warning message:
>>>>> In findOverlappingPeaks(Peaks[[1]], Peaks[[2]], NameOfPeaks1 =
>>>>> NameOfPeaks[1], :
>>>>> Please use select instead of multiple!
>>>>>
>>>>> What should I do?
>>>>>
>>>>> Could you please help me as I am a rookie in this field
>>>>>
>>>>> Thank you in advance
>>>>>
>>>>> Yours sincerely
>>>>> Petros Kolovos
>>>>>
>>>>>
>>>>>
>>>>>
>>>>
>>>>
>>>>
>>>
>>>
>>
>>
>>
>
>
More information about the Bioconductor
mailing list