[BioC] converting transcripts to gene level

Valerie Obenchain vobencha at fhcrc.org
Thu Feb 2 17:43:01 CET 2012


You may want to check out the cummeRbund package. It doesn't do the 
conversion you want in but offers other storage and manipulation options 
of Cufflinks data that may be of interest. It's written by the authors 
of TopHat and Cufflinks.

Valerie


On 02/02/2012 04:42 AM, Michael Lawrence wrote:
> Hi,
>
> There's nothing that does exactly that. You can use
> rtracklayer::import(gtf, format = "gff2") to get the file into a RangedData
> or GRanges. And then after your manipulations write it out again with
> export(object, file, "gff2").
>
> Michael
>
> On Wed, Feb 1, 2012 at 5:15 PM, Abhishek Pratap<apratap at lbl.gov>  wrote:
>
>> Hi All
>>
>> I have cufflinks generated transcripts file in gtf format. I wondering
>> if I can use any R package to convert it to gene level gtf file. By
>> this means all the transcripts are collapsed into genes and each gene
>> contains all the exons seen in all the transcripts for that gene.
>>
>>
>> Sample Data attached.
>>
>> In the attached GTF file I have 3 made up transcripts from 2 genes.
>> What I would want to get as a output is a gtf file containing two
>> genes and all the exons associated with them.
>>
>> Thanks!
>> -Abhi
>>
>> _______________________________________________
>> 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



More information about the Bioconductor mailing list