[BioC] a question about girafe

Martin Morgan mtmorgan at fhcrc.org
Sun Jul 3 14:09:36 CEST 2011


On 07/01/2011 08:01 PM, Sean Davis wrote:
> On Fri, Jul 1, 2011 at 10:55 PM, chun-jiang he<camelbbs at gmail.com>  wrote:
>> Thanks Sean. But I dont know where I can input my .sam file.
>
> Well, you'll want to convert your .sam file to a sorted .bam file and
> then index the .bam file (use samtools to do that, for instance).

or Rsamtools asBam, sortBam, indexBam which are simple wrappers around 
samtools code so have the same basic memory / speed characteristics.

Martin

> Then, you can use the readAligned() function from the ShortRead
> package.  In the girafe vignette, you can see section 2.2 for an
> example.
>
> Sean
>
>
>> On Fri, Jul 1, 2011 at 9:48 PM, Sean Davis<sdavis2 at mail.nih.gov>  wrote:
>>
>>> On Fri, Jul 1, 2011 at 6:14 PM, chun-jiang he<camelbbs at gmail.com>  wrote:
>>>> Hi,
>>>>
>>>> Does anybody use girafe to analyze .sam file extracting from bowtie.
>>>>
>>>> Some example is better.
>>>
>>> Hi, Jiang.  Have you looked at the vignette for the girafe package?
>>>
>>>
>>> http://www.bioconductor.org/packages/release/bioc/vignettes/girafe/inst/doc/girafe.pdf
>>>
>>> It does include examples.  If you have problems, send the code that
>>> you have tried, the error messages, and the output of sessionInfo().
>>>
>>> Sean
>>>
>>>
>>>> thanks so much,
>>>>
>>>> Jiang
>>>>
>>>>         [[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
>>>>
>>>
>>
>>         [[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
>>
>
> _______________________________________________
> 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


-- 
Computational Biology
Fred Hutchinson Cancer Research Center
1100 Fairview Ave. N. PO Box 19024 Seattle, WA 98109

Location: M1-B861
Telephone: 206 667-2793



More information about the Bioconductor mailing list